Loading...

Archive

    25 June 1994, Volume 11 Issue 2 Previous Issue    Next Issue

    For Selected: Toggle Thumbnails
    INTERACTION OF SPECIAL SOLITONS IN TODA LATTICE
    Duan Wenshan, Lu Kepu, Zhao Jinbao
    1994, 11(2): 129-133. 
    Abstract ( )   PDF (248KB) ( )  
    The amplitude of particle velocity of one soliton in Toda lattice is given. The paper also provides the dependence of the maximum velocity in the interacting process of both solitions on the initial amplitude prior to collision.
    COMPUTATION OF IMPACT FORCE AND ANALYSIS OF DYNAMIC RESPONSE ON IMPACTED ELASTIC PLATES WINKLER FOUNDATION-SUPPORTED
    Huang Wei, Zou Yida
    1994, 11(2): 134-140. 
    Abstract ( )   PDF (402KB) ( )  
    Discussion covess dynamic response problem of an elastic plate with arbitrary boundary shapes on a half-space linear viscoelastic Winkler foundation impacted by a moving body.Undes the condition of a known initial velocity of the body, a non-linear Volterra integral equation on impact force F(t) is deduced, and the dynamic transverse deflection W(x,y,t) of the plate is put forward. To solve the integral equation, an effective numerical computation method is used. At the end, results of two practical examples are shown and discussed.
    THE MODIFIED TEMAM SCHEME ON NONUMIFORM STAGGERED MESH IN NUMERICAL SIMULATION OF DRIVEN CAVITY FLOW
    Wu Yadan, Huang Lanjie
    1994, 11(2): 141-148. 
    Abstract ( )   PDF (474KB) ( )  
    An efficient method for multiple-scale problems is presented, which extends the authors' modified Temam scheme on the staggered mesh [5] to nonuniform staggered mesh. The nonuniform mesh is generated by a continuous coordinate transformation without singlarity according to boundary layer thickness so as to ensure accuracy. The constraint condition, i.e. continuity equation, is realized by the pressure correction projection method; the related Poisson equation on nonuniform mesh is solved by the generalized cyclic reduction program from FISHPACK. Numerical simulation of driven flow in a square cavity was done for steady-state solutions for Re=100, 400, 1000, 5000. These solutions accurately match those given in literature; for high Reynold numbers, tertiary vortices near the corners in the cavity are resolved. These results shows that with suitable nonuniform mesh, the above method remarkably improves computational efficiency.
    FURTHER DISCUSSION ON “AN OPTIMAL ALGORITHM FOR MATRIX MULTIPLICATION”
    Jiang Changjun, Wu Zhehui
    1994, 11(2): 149-153. 
    Abstract ( )   PDF (284KB) ( )  
    Three improved ways are given for the algorithm of matrix multiplication in the paper [1] and the new algorithm has the advantages of usefulness and parallilism.
    A NUMERICAL TREATMENT FOR TURBULENT NONLINEARITY
    Hu Shuanghui, Huang Lin, Qiu Xiaoming
    1994, 11(2): 154-160. 
    Abstract ( )   PDF (402KB) ( )  
    The paper dercribes a numerical method for the resistivity gradient driven turbulence in the Tokamak edge. The mode selection for the mode-mode coupling or nonlinear interaction is discussed. The corresponding computer code for three-dimension magnetic hydrodynamics equations has been constructed, which can save the computation time effectively.
    THE SPECTRAL METHOD FOR A PERTURBED SINE-GORDON EQUATION
    Guo Boling, Wu Xianghui
    1994, 11(2): 161-166. 
    Abstract ( )   PDF (302KB) ( )  
    A pertured sine-Gordon equation utt-uxx+αsinu-βuxxt=g(u),t>0,-∞< x< ∞,with the periodic initial and boundary value problem is studied. A spectral method is proposed,a priori estimates are established and the existnce and uniqueness are proved for the generalized solution.The error estimates are given and the convergence is proved for the approximate solution of the spectral method.
    THE PROBLEM ABOUT HIGH ORDER SINGULAR INTEGRAL OF BOUNDARY ELEMENT METHOD FOR ELASTODYMMICS
    Chen Guoroug
    1994, 11(2): 167-171. 
    Abstract ( )   PDF (279KB) ( )  
    The problem on singular integral of the Boundary Element Method for elastodynamics is studied, and a calculating method of the high order singular integral is presented.The method can indirectly find the singular integral by making use of the special solution of elastodynamic eguation. As a result, the difficulties for singular integral will be avoided. Finally a example is given, in the paper.
    COMPUTER SIMULATION OF RADIATION DAMAGE IN THE FIRST WALL CAUSED BY FUSION NEUTRON
    Li Rongwu, Pan Zhengying, Shao Qigun, Gao Xinghua
    1994, 11(2): 172-178. 
    Abstract ( )   PDF (491KB) ( )  
    A neutron transport code and a charged-ion radiation damage code have been linked to simulate neutron induced defects in the static, multi-component amorphous target. In simulating collision cascades, replacement collisions have been taken into account, especially. It corrected the value of DPA about 10%~20%. Improper replacements were recorded also. The radiation damage in two kinds of materials candidated for the first wall in fusion reactor is calculated. Calculation gives the production rates of PKA and DPA, gas production rates and their depth distributions. The results show that damage distributions are almost uniform body distributions, the DPA value of PE-16 is a little larger than that of SS316, the hydrogen production rates is nearly two times as much as that of SS316, the helium production rates is 1.45 times as much as that of SS316. It means that the radiation damage of PE-16 is more serious than that of SS316. The cross section data from MC(87) n library is used in the calculation.
    NUMERICAL SIMULATION OF TURBULENT FLOW AND MELT TOPOGRAPHY IN ALUMINUM REDUCTION CELLS
    Huang Zhaolin, Yang Zhifeng, Wu Jianghang
    1994, 11(2): 179-184. 
    Abstract ( )   PDF (355KB) ( )  
    A mathematical model is proposed for numerical simulation of the flow of the two superposed liquids (the electrolyte and molten aluminum) and of the electrolyte-metal interface wave. By employing a modified scheme based on the SIMPLE algorithm, velocities within both liquids and melt topography are calculated for a 280kA cell. Satisfactory results of numerical prediction are given.
    THE DIFFRACTION OF THE SOLITARY WAVE BY A HUMP AT THE SEA BOTTOM AND PRACTICAL COMPARISONS OF THE DIFFERENT OBC
    Zou Guangyuan, Wu Junxiao
    1994, 11(2): 185-194. 
    Abstract ( )   PDF (638KB) ( )  
    Open boundary condition (OBC) is one of the most challenge subjects in the field of computational method for infinite domain. There are already many different kinds of OBCs in wave problems. But, they still lack practical comparisons of effectiveness. This article devotes to such comparisons of several important OBCs at the artificial boundaries for the nonlinear wave equation Boussinesq equation. A test example is given on the diffraction of solitary wave by a hump at the sea bottom. 3-D scattering wave diagrams are also given. The effectiveness is analyzed with instructive conclusion.
    APPlication of A TVD Technique to Oblique Shock wave Reflection in Weakly Ionizing, Chemically Relaxing Gas Flow
    Hu Guangchu, Le Jialin, Cao Wenxiang
    1994, 11(2): 195-202. 
    Abstract ( )   PDF (445KB) ( )  
    In this paper several planar oblique shock reflection in weakly ionizing,chemically relaxing gas flow have been successfully simulated using a second-order upwind Total Variation Diminishing(TVD) scheme. In present nonequilibrium gas flow calculations the gas-dynamics and species-continuity equations are solved by fully coupling manner,and the fully coupling procedure is described in detail.Numerical solutions are presented with detailed discussion for three examples of oblique shock wave reflection.The numerical results obtained here agree well with the experiments in Refe.[4] by I.I.Glass and calculation results in Refe.[1] H.M.Glaz.
    THE SURFACE STRUCTURE CALCULATED OF CO ADSORBED ON THE TRANSITION METALS SURFACES WITH LOW ENERGY ELECTRON DIFFRACTION SPECTRA
    Lan Tian, Xu Feng
    1994, 11(2): 203-211. 
    Abstract ( )   PDF (543KB) ( )  
    This paper studied the sarface structure of CO molecules adsorbed on the transition metals surfaces by Low-energy electron diffraction spectra. It is shown that for CO/Ni {100}, the atomic layer distances dNi-C=0.1810±0.0001nm, dC-O=0.1090±0.0001nm; the CO molecule bond length lC-O=0.1150nm, and the tilted angle β=18.59°±0.15°, and that the structure Parameters are close to the figure in Ni(CO)4, It is found that tbr CO/Rh{111}, the atomic layer distanes dRh-C=0.1950±0.0001nm, dC-O=0.105±0.002nm, the CO molecule bond length lC-O=0.1170nm, the tilted angle β=26.18°±2.2ånd that the structure parameters are close to the figure in Rh2(CO)8. It is also found that for CO/Cu{100}, the atomic layer distances dCu-C=0.192±0.002nm, dC-O=0.1063±0.0001nm; the CO molecule bond length lC-O=0.1129nm,the tilted angle β=19.69°±0.15°; for CO/Pd{100}, the atomic layer distances dPd-C=0.157±0.003nm, dC-O=0.1126±0.0004nm; the CO molecule bond length lC-O=0.1129nm, the tilted angle β=4.18°±2.2°. And for CO/Ti{0001}, the atomic layer distances dTi-C=0.195±0.003nm,dC-O=0.1191±0.0001nm; the CO molecule bond length lC-O=0.1129nm, and the tilted angle β=7.63°±0.15°.
    A IMPROVED EXPLICIT SYMPLECTIC ALGORITHM
    Liao Xinhao, Liu lin
    1994, 11(2): 212-218. 
    Abstract ( )   PDF (411KB) ( )  
    Based on the property that the great majority of celestial motions in solar system can be described by the model of two-body problem with perturbations, we use the composition of algorithm to construct out a improved explicit symplectic algorithm with higher accuracy. Numerical results have shown that the algorithm have advantages in the computations of celestial orbits compared with usual symplectic algorithms.
    NUMERICAL SOLUTION FOR CRITICAL PARAMETER OF THERMAL EXPLOSION
    Duan Qing sheng, Qin Cheng sen
    1994, 11(2): 219-224. 
    Abstract ( )   PDF (361KB) ( )  
    A variational method is presented for the numerical solution of the critical parameter of thermal explosion. The method is based on the variational principle [1] for the critical parameter of thermal explosion. The calcuated results show that the method is simple,convenient and more accurate.
    TWO-SOURCE-COUPLED-INPUTED GUASI SECOND ORDER SOLITON AND ITS COMMUNICATION APPLICATION
    Chen Lujun, Liang Changhong, Meng Shengwei, Guo Yao
    1994, 11(2): 225-229. 
    Abstract ( )   PDF (316KB) ( )  
    Based on nonlinear Schrodinger equation, the properties of two-source-coupled-inputed soliton and the interaction between such soli tons and the first order solitons are numerically investigated by the split-step Fourier transformation method. A tentative idea of triply optical soliton communication system is proposed and disscussed, which made the relay distance further extended with the information increased by 1.6 factor. Besides, the high power input of the light source is not necessary by using such solitons.
    CONVERGENCE AND CHOICE OF OPTIMUM PARAMETERS OF BAORJ ALGORITHM
    Hu Jiagan
    1994, 11(2): 230-236. 
    Abstract ( )   PDF (346KB) ( )  
    The convergence of BAORJ algorithm is considered and the optimum paramenters and the corresponding rate of convergence for the model problem are obtained. These results will benefit the use of BAORJ algorithm.
    HIGHER ORDER TWO PARAMETER PARALLEL JACOBI-TYPE METHOD AND ITS CONVERGENCE
    Hu Jiagan, Liu xingping
    1994, 11(2): 237-243. 
    Abstract ( )   PDF (393KB) ( )  
    The higher ordar parallel Jacobi-Type method for solving system of linear algebraic equations is proposed, the convergence of the method is analysed and the optimum parameters and the corresponding spectral radius of the iterative matrix for the model proplem and the like are given. In the end there are some numerical examples to illustrate the effectiveness of our method.
    A GROUP THEORY METHOD ON CALCULATION OF NUCLEUS-NUCLEUS COLLISION AT HIGH ENERGY
    Chen Kezhong, Guo Bingying
    1994, 11(2): 244-247. 
    Abstract ( )   PDF (238KB) ( )  
    In view of the serious complexity of the calculation of nucleus-nucleus colli sion, a group theory method is proposed, by which the integrals with same value in the amplitude equation can be identified, and the calculation would become simple and easy to proceed on computer. The α-αelastic scattering at 5.07 GeV/c is calculated concretely and the result is satisfactory.
    THE MCHF PACKAGE AND COMPUTATIONAL CODES OF EIGENCHANNEL QUANTUM DEFECT PARAMETERS
    Qian Xingzhong, Wang Zhiwen
    1994, 11(2): 248-251. 
    Abstract ( )   PDF (248KB) ( )  
    The program for computing EQDT parameters is achieved by using the program source ready in original MCHF package on basis of the eigen-channel quantum defect theory with electrostatic and spin-orbit interactions.
    MICROCOMPUTER PROGRAM FOR RADIAL DISTRIBUTION FUNCTION OF NONCRYTALLINE MATERIALS
    Wang Junqiao, Chen Jizhou
    1994, 11(2): 252-255. 
    Abstract ( )   PDF (242KB) ( )  
    A new RDF program with wide adaptability and good reliability is suggested. It can be used both to prepare experimental radial distribution qunction (RDF) from wide angle X-ray scattering (WAXS) and to calculate theoretical RDF from model structure. In addition, it can also determine structure parameters from RDF.
Share: