Loading...

Archive

    31 December 1992, Volume 9 Issue S1 Previous Issue    Next Issue

    For Selected: Toggle Thumbnails
    THE RESEARCH ON THE REMAINDER EFFECTS OF FINITE DIFFERENCE SCHEMES
    Liu Ruxun
    1992, 9(S1): 479-484. 
    Abstract ( )   PDF (437KB) ( )  
    In the paper, based on the Modified PDE of finite difference schemes, a general and systematical research and analysis on the remainder effects of finite difference schemes are advanced. The dissipasion relationship and the dispersion relationship are defined. By the use of the two relationships, the stability, the dissipation, the dispersion and the group velocity of finite difference schemes are all disccused.
    A KIND OF GRIDDING FOR SCATTERED POTENTIAL DATA
    Cai Zongxi, Cai Yuehong
    1992, 9(S1): 485-486. 
    Abstract ( )   PDF (160KB) ( )  
    In this paper, the solution of magnetization intensity J(α,β)on imaginary magnetic dipole layer was obtained by the represantation of dipole layer field density component which was regarded as the first kind Fredholm intergral equation and the potential data on the rectangle grid point was computed by the solution J(α,β)which was taken into the represantation. It realizes gridding the scattered potential data which keep the physical character of potential data.
    THE STUDIES OF GROUNDWATER QUALITY MODELS AND ITS NUMERICAL METHOD IN KARST AREA OF JINAN
    Yiang Tianxing, Wang Yunguo, Yao Leihua
    1992, 9(S1): 487-491. 
    Abstract ( )   PDF (320KB) ( )  
    In this paper, double media of groundwater quality models of SO42- and total hardness in karst area of Jian have been built and the combination of Galerkin finite element method with generalized upwind-balance scheme is presented. This method has overcomed the numerical dispersion and the vibration phenomenon. Goodness results are obtained.
    MARKED ZEROS ALGORITHM WITH PROGRAM SOFTWARE
    Zhao Yupeng, Ma Sui, WorksXu Yumin, Deng Feiqi
    1992, 9(S1): 492-494. 
    Abstract ( )   PDF (220KB) ( )  
    In this paper, we point out the essential which has been discussed for several years in some defective papers on solution of the assignment problem and present a new independent algorithm, i.e. the "marked zeros algorithm". Besides the theoretical analyses on the algorithm, we give the program software. We also illustrate the efficiency of the algorithm by examples.
    CONVERGENCE OF BPSD METHOD FOR T(q.r)MATRIX
    Hu Jia gan
    1992, 9(S1): 495-497. 
    Abstract ( )   PDF (224KB) ( )  
    In this paper, we obtain a necessary and sufficient condition, when the coefficient matrix A of the equation Ax=f considered is an T(1, 1) matrix,a sufficient condition, when A is an T(1, 2) or T(2, 1) matrix for the convergence of BPSD method.we also obtain the optimum parameters and the optimum rate of convergence of BPSD method, when A is an T(1, 1) matrix and a necessary and sufficient condition, when A is postive definite and we point out that the necessary and sufficient condition in [1] and [9] is only sufficient. The derivation is simple and interesting. Our method can also be used to prove the convergence of some other iterative methods.
    A NEW AND UNITE METHOD OF DEALING SINGULAR INTEGRAL ABOUT BOUDNARY ELEMENT ANALYSIS OF AXISYMMETRIC HEAT CONDUCT PROBLEM
    Yao ShouGuang, Zhu Deshu
    1992, 9(S1): 498-500. 
    Abstract ( )   PDF (240KB) ( )  
    In this paper, an new and unite method about singular integral of boundary element ana-lysis for steady axisymmetric potential problem is presented. On this basis, an new simple method to solve above problem is developed. The results of numenrial analysis and application show that the method developed by author is successful.
    NUMERICAL SIMULATION OF PHYSICAL PROBLEMS WITH HIGH ORDER ACCURATE SCHEMES
    Fu Dexun, Ma Yanwen
    1992, 9(S1): 501-505. 
    Abstract ( )   PDF (341KB) ( )  
    Inorder to simulate multiscale physical problems constructing high order accurate difference schemes is needed. The fourth order accurate compact scheme is reconstructed using the method of diffusion analogy developed by authors in order to make the scheme dissipative. The modified compact scheme with fourth order accuracy in space and third order accuracy in time is used to compute the 2-D shock reflection problems. The computed results show the scheme developod to correctly simulate multiscaler physical problems also can be used to capture the shock well.
    NUMERICAL CALCULATION OF A CLASS OF HIGHLY OSCILLATORY INTEGRALS WITH THE MATHIEU FUNCTION
    Long Yongxing
    1992, 9(S1): 506-510. 
    Abstract ( )   PDF (305KB) ( )  
    This paper describes a method for computing highly oscillatory integrals with the Mathieu function. The practice proves that not only the results are highly satisfactory, but also the method is time-saving.
    CALCULATION OF SYMMETRY COEFFICIENTS IN COMPLEX CLUSTER
    Gong Xingao, Zheng Qingqi, Pan Wei
    1992, 9(S1): 512-514. 
    Abstract ( )   PDF (232KB) ( )  
    A method to calculate the symmetry coefficents in quantum chemistry calculation is in-troduced, which is very useful for calculations of complex cluster. Some results of C60 molecule are aslo represented.
    A SURVEY OF AGE METHOD COMPARING WITH SEVERAL FRACTIONAL STEP METHODS FOR SOLVING 2-D PARABOLIC EQUATIONS
    Wu Sheng Chang, Guo Ge
    1992, 9(S1): 515-519. 
    Abstract ( )   PDF (326KB) ( )  
    In this paper, we discuss the method of Alternating Group Explicit (AGE) for solving 2-D parabolic equation, and give out the results of AGE method comparing with several fractional step methods.
    NEW NUMERICAL METHODS AND THEIR THEORETICAL ANALYSES FOR THE TWO-PHASE DISPLACEMENT PROBLEMS
    Liang Dong
    1992, 9(S1): 521-525. 
    Abstract ( )   PDF (339KB) ( )  
    In this paper, for the two-phase immiscible and miscible displacement problems we propose two kinds-of new numerical methods: the viscosity splitting method and the upwind difference method over arbitraty triangulation pratitions, and discuss their theoretical analyses. The numerical results are very efficient and reasonable for model problem.
    STUDY OF THE NON-STATIONARY DESIGNING THEORY FOR SOLAR GREEN HOUSE
    Zhao Bingjun, Hu Guang zhi, Zhang Jian, Hu Guangyao, Lui Wukui
    1992, 9(S1): 527-529. 
    Abstract ( )   PDF (228KB) ( )  
    This artical uses the non-stationary heat-transfer theory to carry out dynamic simulation for solar green house, establishes the mathematical model of solar green house, uses computer to design and calculate solar house system, analyses the influence to solar green house by various parameters, and taking the area of Henan Province as a example, makes out the optimum designing datum which possess extensive practical value.
    THE SIMULATION OF HIGHER OVERVOLTAGE INDUCED BY NONSYNCHRONOUS CLOSING
    Lai Dingwen
    1992, 9(S1): 531-532. 
    Abstract ( )   PDF (153KB) ( )  
    As a continuation of our previous research this paper puts the stress on numerical simulation of tigher overvoltage induced by nonsynchronous closing. Main result as follow: Because of ordinary differential equation which is the mathematical model of preceding closing phases has two large harmonic solutions, so it implys possibility transformed from little harmonic solution to large one. Lagging currents i1(t1), i2(t1) play an important role to help intensify the billows and Waves in this transformation process.
    FINITE ELEMENT METHOD IN TEMPERATURE FIELD ANALYSIS IN SOLIDIFICATIONS OF METALS
    Chen Yongjin, Wang Honggang, Cheng Heming
    1992, 9(S1): 537-538. 
    Abstract ( )   PDF (144KB) ( )  
    By using Linearization method in little time interval, from the heat conduition equation with the 1st-order phasetransformations, the foundation equations of Finite Elememts in solidifications of Metals is established in the paper. The numerical analog on symetric problem in solidification of a grayiron in sandform is obtained by means of triangular element, and the results tallies the corresponding test well.
    AN INTERVAL METHOD FOR SEEKING CONSTRAINED NONLINEAR PROGRAMMING OF MANY VARIABLES
    Wang Haijing, Zhang Nailiang, Liu Yunhua
    1992, 9(S1): 539-541. 
    Abstract ( )   PDF (218KB) ( )  
    An interval method is given for solving constrained nonlinear programming. Continuous interval extension and the function test are used to delete all the unnecessary elements. The optimum solution and optimum value are obtained at the same time during the process of iteration. Numerical example is given and the result show that the method works well.
    CORRECTIVE METHODS OF GENERALIZED GALERKIN METHOD
    Zhang Dakai
    1992, 9(S1): 542-544. 
    Abstract ( )   PDF (186KB) ( )  
    A kind of the methods with high accurate of the generalized Galerkin method-inteprotation corrective method and error corrctive method are proposed in this paper.
    MONTE CARLO METHOD FOR FINDING ROOTS OF A TRANSCENDENTAL EQUATION
    Gong Ye
    1992, 9(S1): 547-548. 
    Abstract ( )   PDF (160KB) ( )  
    In the paper two kinds of Monte Carlo method are presented for finding roots of a transcendental equation. All real roots may by found without having to choose an initial value. Therefore, Monte Carlo method often is efficient for the problems that are inefficient using the other methods.
    THE USE OF RAY-TRACING TECHNIQUE FOR MONTE CARLO METHOD
    Zhang Liwu, Zhang Yuqin
    1992, 9(S1): 549-550. 
    Abstract ( )   PDF (170KB) ( )  
    The Ray-tracing technique is useful to decrease computing time.
    MONTE CARLO SIMULATION OF DEPTH DISTRIBUTION OF DISPLACED ATOMS IN AMORPHOUS TARGETS
    Shao Qiyun, Pan Zhengying, Huo Yukun
    1992, 9(S1): 551-554. 
    Abstract ( )   PDF (297KB) ( )  
    A Monte Carlo code is described which simulates the displacement process and the transmissions of energetic particles in solids by ion bombardment, including the operating principle, the physics model and the computing examples. Based on the binary collision approximation, the depth distributions of displasced atoms in amorphous targets induced by α particles are investigated, and the relations between the depth distributions of displaced atoms and energy deposition in solids are discussed.
    PROSPECT ON GENENAL SOFTWARE OF MONTE CARLO MATHOD
    Pei Lucheng
    1992, 9(S1): 560-562. 
    Abstract ( )   PDF (229KB) ( )  
    This is a short paper on the prospect of Monte Carlo genenal software. The content consists of cluster sampling method, zerovariance technique, self-improved method, and vectorized Monte Carlo method.
    SOME PROBLEMS ON MONTE CARLO METHOD DEVELOPMENT
    Pei Lucheng
    1992, 9(S1): 563-566. 
    Abstract ( )   PDF (332KB) ( )  
    This is a short paper on some problems of Monte Carlo method development. The content consists of deep-penetration problems, unbounded estimate problems, limitation of Mdtropolis' method, dependency problem in Metropolis'method, random error interference problems and random equations, Intellectualization and vectorization problems of general software.
    THE MONTE CARLO COMPUTATION OF DAMAGE PROBABILITY FOR THE TARGET BOMBED BY SHRAPNEL
    Huang Qijin, Liu Xingping
    1992, 9(S1): 567-568. 
    Abstract ( )   PDF (152KB) ( )  
    In this paper the Monte Carlo method in 2D and 3D is used to calculate the damage probability resulted from shrapnel-attack on the target, in addition, an exact analytical formula for 2D uniform distribution is also proposed.
    METHODS IN FINDING THE Q2-DEPENDENT NUCLEON PARTON DISTRIBUTION FUNCTIONS IN THE CONSTITUENT QUARK MODEL
    Qian Lu
    1992, 9(S1): 570-572. 
    Abstract ( )   PDF (221KB) ( )  
    The deconvolution problem of the Constituent Quark Model (CQM) in finding the partondistribution function in valon at Q02 is reduced to an optimum problem and solved successfully with Newton-Raphson method. Moreover, by discretizing the QCD evolution eqs into ordinary simultaneous differential eqs, we obtain the parton distribution functions in valon at arbitary Q2 > Q02, and thus, through convolution, we finally find the Q2-dependent nucleon parton distribution functions in CQM picture.
    HE REDUCE LANGUAGE IS APPLIED TO CALCULATE ANHARMONIC OSCILLATOR
    Jiang Saozhou
    1992, 9(S1): 573-577. 
    Abstract ( )   PDF (274KB) ( )  
    The Reduce language is a useful tool for deriving analytical formula on the computer. We calculate the effect of the high-power anharmonic terms in the anharmonic oscillator.
    SNENM/RTSP:A THREE-DIMENSIONAL NODAL AND SUBCHANNEL ANALYSIS CODE PACKAGE FOR THE CORE TRANSIENT IN PWPS
    Cheng Pingdong, Zhao Huimin, Shao Tianwei, Zhou Qianfu, Shen Sen
    1992, 9(S1): 579-582. 
    Abstract ( )   PDF (316KB) ( )  
    The paper introduces a three-dimensional nodal and subchannel analysis code package for the core transient in PWPs-SNENM/RTSP, including the model of neutronics and thermohydraulics used in the package. The preliminary results of applying this code package to the analysis of Steam-Line Break accident for Qinshan NPP are given and also evaluated briefly.
    A EFFECTIVE METHOD OF BURNUP CALCULATION
    Zhong Wenfa, Luo Jingyu, Zhang Ruiyin
    1992, 9(S1): 583-586. 
    Abstract ( )   PDF (304KB) ( )  
    A method of burnup calulation is introduced in calculation of reactor physic in this paper. Burnup calculation is completed by interpolation table of macroscopic cross section vs burnup, interpolated method and solving diffusion equation. The calculated result compared with shearon Harri s's result is satisfied. The method has characteristic of use friendly, less calculated time and occuracy demand. Now it is used in the practical calculation of reactor physic. A effective tool is provided for engineering design.
    THE NUMERICAL SIMULATION ON LOW-LEVEL RADIOACTIVE WASTE WATER, LOW-TEMPERATURE COOLING WATER DRAINED EFFECT OF IMPINGEMENT FROM THE DAYAWAN NUCLEAR POWER STATION
    Zhang Chunlin, Huang Zuzhan, Kong Lingfeng, Peng Jianwei
    1992, 9(S1): 587-592. 
    Abstract ( )   PDF (408KB) ( )  
    In this paper, we calculated the radioactive Concentration distribution of radioactive waste water, the temperature distribution of cooling water which drained and the effect of impingement from the Dayawan Nuclear Power Station on nearby waters range, discussed and analysed some problems of Computational results and computatin with Alternating Direction Implicit Method (ADI). The contents of the article indued: the establish of two-dimension tidal current equation, radioactive waste water pollutant dispersion equation and cooling water heat convection diffusion equation, the numerical difference calculation model of tidal current field, concentration field as well as temperature field, effect impingement with ADI method, numerical calcuation results.The result of research showed that: when the Dayawan Nuclear Power Station is on normal operation and after the low Level radioactive waste water and Low temperature cooling water have been drainde off into the sea, it is found that the radioactive concentration is near and enen Lower than the natural background radiation just outsite the area of about 4km2 round the water outlet, the temperature difference is not obrious just outside the area of about 6 km2 round the water outlet, the effect of impingement of abont 4.3%.
    THE STATUS AND PROSPECT OF REACTOR ANALYSIS METHODS FOR LWR
    Zhang Zongyao
    1992, 9(S1): 593-596. 
    Abstract ( )   PDF (386KB) ( )  
    This paper gives a review of status and prospect of reactor physics analysis methods and codes on lattice physics and core physics for LWR
    BUCKLING FEEDBACK OF THE SPECTRAL CALCULATIONS
    Jing Xingqing, Shan Wenzhi, Luo Jingyu
    1992, 9(S1): 597-601. 
    Abstract ( )   PDF (344KB) ( )  
    This paper studies the problems about buckling feedback of spectral calculations in physical calculations of the reactor and presents a useful method by which the buckling feedback of spectral calculations is implemented. The effect of the buckling feedback in spectra and the broad group cross scetion, comvergence of buckling feedback iteration and the effect of the spcetral zones dividing are discussed in the calculations. This method has been used for the physical design of HTR-10MW Test Module.
    STEADY-AND TRANIENT-STATE TEMPERATURE FIELD CALCULATIONS OF A HE-COOLED FIRST WALL FOR FUSION BREEDER
    Wang Xueren, Lu Xiaolan, Huang Jinhua
    1992, 9(S1): 602-606. 
    Abstract ( )   PDF (384KB) ( )  
    In this paper, a finite differential heat transfer computer code HEATING-5 has been used to model the thermal hydraulic problems of a He-cooled first wall for fusion breeder. The numerical techniquws of solving steady and transient state heat conduction equations are presented. The numerical calculations have been done. The temperature distributions of the coolant in different flow locations have been calcultaed by analytical method as the iuput parameters of the HEATING-5 code. in order go conform the fluid terms in the energy balance equation which are not included in the original code. Meanwhile, another finiteelement heat transfer computer code AYER, which includes forced convection flow terms in the heat transfer equation, is also used to solve the same problem. It shoes that the results by using HEATING-5 code to model coolant flow heat transfer are well consistent with that by using AYER code.
    AN ITERATIVE METHOD FOR FITTING MULTI-COMPONENT DECAY CURVES WITH COUNT RATE CORRECTION
    Li Weisheng, Li Wenxin, Yin Xinmin
    1992, 9(S1): 607-610. 
    Abstract ( )   PDF (338KB) ( )  
    A conversion lormula of the count rate in the measuring interval for the radioactive decay and an iterative program for fitting multi-component decay curves wity count rate correction are reported, the proposed program gas been used successfully to analyse artificial decay curves and experimetal decay curves. The comparison with conventional medhod shows that the present medhod gives satisfactory results.
    A NUMERICAL METHOD FOR SOLVING THREE-TEMPERATURE EQUATIONS OF RADIATION VIA DOMAIN DECOMPOSITION
    Feng Tinggui
    1992, 9(S1): 611-614. 
    Abstract ( )   PDF (253KB) ( )  
    A numerical method for solving three-temperature equations of radiation viadomain decomposition is presented. The method is based on dividing the domain into two or several domains and solving either three -temperature equations or single-temperature euqation in each domain. Also we present a scheme coupling interface flux between one domain and the next and analyse the convergence of the iteration method for solving the difference equations.
    A STUDY OF HIGH PARTIAL WAVE COLLISION STRENGTH
    Fang Quanyu, Cai Wei, Li Ping, Shen Zhijun, Zou Yu
    1992, 9(S1): 615-618. 
    Abstract ( )   PDF (293KB) ( )  
    The problem of the slow convergence has been resolved in partial wave expansion for the excitations of electron-ion collisions, with three approximations including plane wave, coulomb Bethe and equimultiple decrease series. As an example, we have computed some important transitions of Ne-like Ge ion and made comparisons between the results obtained by three approximations. The results slow-that all the three methods are acceptable.
    CROSS SECTION FOR THE EXCITATION OF Al10+ ION BY LOW-ENERGY-ELECTRON IMPACT
    Zou Yu, Fang Quanyu
    1992, 9(S1): 619-620. 
    Abstract ( )   PDF (182KB) ( )  
    Cross section for the excitation of Al10+ ion by electron impact from ground state to 2p,3s and 3p states has been investigated with 2-state does-coupling (2cc) and 5-state close -coupling(5cc) approximations. The indirect coupling effect has been analyzed, comparing with the results of distorted wave approximation, it is found that the coupling effect is less important for 2s-2p transition, but for 2s-3s and 2s-3p transition it should be taken into account.
    THE NUMERICAL METHOD OF LASER INDUCED VAPORIZATION OF METAL
    Tao Yingxue, Guan Jili, Shen Longjiun, Lew Chenghai
    1992, 9(S1): 621-623. 
    Abstract ( )   PDF (220KB) ( )  
    For numerical solution of laser induced vaporization of metal, we have to vesearch of calcolation method of fluid dynamics including phase transition. Because the vaporization zone is movement as time, so that, this problem is different from general fluid dynamics with a moving zone.This paper considers the energy deposition as zone deposition. According to the conditions of first class phase transition, outside the vaporization zone, we adppt the von Neumann method, and inert the vaporization zone, the another method is used. We have established the program named LHAP -1DVG, in which this method is applied.The results of a lot of calcolations indicate that this method is successful.
    NUMERICAL MODELLING CALCULATION AND ANALYSIS OF GAIN COEFFICIENT FOR CW CHEMICAL OXYGEN IODINE LASERS
    Chen Liyin, Chu Zexiang
    1992, 9(S1): 624-628. 
    Abstract ( )   PDF (377KB) ( )  
    The relation between components, pressure, flow velocity and gain coefficient of CW OICL is theoretically analyzed by using a simple model throug digital computation. The operating region concerned is [O2*]/[I2]=10~200,[O2*]/[O2]=2~9,P=1~100 torr,U=5-100 m/s, The necessity of flow conservation equation for a CW OICL is detially discussed, It is shown that the available gain region will by enlarged proportionally to flow velocity while the peak gain value almost remains inveriable when flow velocity is increased.
    NUMERICAL METHOD FOR THE KINETIC BALLOONING MODE EIGENEQUATION WITH THREE COMPONENTS
    ui Guofang, Zhou Hailin
    1992, 9(S1): 629-631. 
    Abstract ( )   PDF (218KB) ( )  
    A numerical method is described for solving the kinetic bollooning mode eigenequation with three components. The effect of high velocity particles on ballooning mode instability is investigated. The obtained results are in good agreement with experiments.
    A NEW ALGORITHM FOR 1-D INVERSE PROBLEM IN SEISMIC EXPLORATION
    Li Zhibin
    1992, 9(S1): 639-640. 
    Abstract ( )   PDF (142KB) ( )  
    In this paper a new algorithm for 1-D inverse problem of wave equations in seismic exploration is given. The inverse problem is changed to be a system of nonlinear equations by the method of difference. The system is linearized by Newton's method, then it is solved iteratively by LMS's method.
    THE COMPUTION OF RADON TRANSFORM AMD ITS APPLICATION ON GEOPHYSICS
    Li Xinxiang
    1992, 9(S1): 643-646. 
    Abstract ( )   PDF (324KB) ( )  
    According to the seismic data processing, we give a algorithm for computing Radon transform and its inverse transform, then we state the principle of appling the transform for processing the seismic data and illustrate the result of post-stack section noise-surpressing by this method.
    METHOD OF SIMULATING THE STELLAR MERIDIAN OBSERVATION AND ITS APPLICATIONS
    Zhu Zi, Xu Jin
    1992, 9(S1): 647-648. 
    Abstract ( )   PDF (158KB) ( )  
    In the paper, the mathematical equations for the output signals of the stellar meridian observation are presented via analysing the processes of the observation. The output signals of the stellar meridian observation using a photon counter are obtained by means of combining the computed signals wiry recording noises. Two practical problems in meridian observation are dicussed and solved with the simulating method.
    EFFECTS OF GENERAL RELATIVITY ON SUPERNOVA CORE COLLAPSE
    Wang Yiren, Yao Jin, Wang Weizhong
    1992, 9(S1): 649-653. 
    Abstract ( )   PDF (391KB) ( )  
    As development of equation of state supernova core collapse will reach high density above nuclear density, the general relativitic effects are most important. The density of core collapse is higher in the general relativitic (GR) core than in the Newtonian (NR) case. The higher the density reached the deeper the shock wave digs into the gravitational well and hance the shock wave will be launched with a larger energy. The formation mass point for shock wave is further from core boundary in GR case than in NR case. The disintegration of iron on the way of shock wave to propagate out will be increased. Once the core centre density reaches it maximum, the imner parts of the cross mass point for two velocity distribution in GR case and in NR case (which slight greater than the sonic mass point in NR case). The velocity value in GR case is greater than one in NR case, but the outer parts of the velocity cross point the velocity value is smaller in GR case than one in NR case, when the shock wave reaches the outer parts where the density and velocity is lower in GR case than in NR, which makes it easier for the shock wave to propagate.
Share: