Loading...

Archive

    25 July 2003, Volume 20 Issue 4 Previous Issue    Next Issue

    For Selected: Toggle Thumbnails
    Some New Progress in the Fields of Computational Petroleum Geology and Others
    YUAN Yi-rang
    2003, 20(4): 283-290. 
    Abstract ( )   PDF (396KB) ( )  
    Numerical simulation in oil-field exploration and exploitation with computational mathematics and dynamics of fluids in porous media is discussed. The numerical method of nuclear waste disposal contamination, prediction and prevention of sea-water intrusion and the numerical simulation of transient behavior of semi-conductor devices are also studied. The mathematical model is the nonlinear convection-dominated diffusion coupled systems with initial-boundary value problem. Characteristic finite difference method, characteristic finite element method, fractional step numerical method and its theoretical analysis are mainly dealed with.
    Discrete Scheme Based on Variation Principle for Three Dimensional Diffusion Equation
    CHEN Guang-nan, LI De-yuan, WAN Zheng-su
    2003, 20(4): 291-297. 
    Abstract ( )   PDF (232KB) ( )  
    Numerical simulation for three dimensional diffusion equation is studied.Using variation principle,the difference discrete scheme with diffusion flux form on the irregular hexahedral grid is constructed.The flux of diffusion as an unknown function is obtained by minimizing the functional and solved with the equation of temperature function simultaneously.Computational formulas are deduced and numerical test is given.
    Adaptive Generalized Difference/Finite Volume Computations for Two Dimensional Nonlinear Parabolic Equation
    WU Hai-jun, LI Yong-hai, LI Rong-hua
    2003, 20(4): 298-306. 
    Abstract ( )   PDF (712KB) ( )  
    A posteriori error estimates of residual type for generalized difference method (finite volume method) are derived for 2-D nonlinear parabolic equation on triangular meshes and an adaptive scheme suitable to problems with sharp moving fronts is designed.A data structure (not a tree structure) and a flexible local coarsening method for triangular meshes suitable to adaptive computations for evolution problems are also presented.
    Differential Quadrature Domain Decomposition Method for a Kind of Free Boundary Problems of Parabolic Equations
    WU Xiong-hua, WU Yun
    2003, 20(4): 307-310. 
    Abstract ( )   PDF (174KB) ( )  
    A new numerical computational Method,Differential Quadrature Domain Decomposition Method,is presented based on the advantages of Differential Quadrature Method and Domain Decomposition Method.Numerical examples find the method is more convenient and efficient for a kind of free boundary problems of parabolic equation with weak singularity in the initial value.
    Comparison and Application of Two Kinds of Finite Difference Scheme in Frontogenesis Simulation
    JI Zhong-zhen, YANG Hong-wei, WANG Bin
    2003, 20(4): 311-314. 
    Abstract ( )   PDF (146KB) ( )  
    Based on the model of ideal frontogenesis in two dimensions, two kinds of finite difference scheme is tested. One scheme is a total energy conservation scheme that is used widely in atmosphere simulation, and the other one is an artificial compression method (ACM) that is used to capture the contact discontinuity in high-speed flows. Both of them have great abilities to simulate frontogenesis. But the effect of the ACM is better than the other to some extent.
    Unified Zonal Fortified Solutions for the Unsteady Flow Fields Induced by a Body Overtaking Shock Waves
    JIANG Juan, ZHOU Wei-jiang, MA Han-dong
    2003, 20(4): 315-320. 
    Abstract ( )   PDF (248KB) ( )  
    Numerical simulations of the flow field around a double-wedge shaped projectile overtaking shocks, which move in two different directions, are carried out by solving unsteady compressible Navier-Stokes equations. In order to have enough grid resolution, the unified zonal fortified solution algorithm is used. Two test cases, the pursuing overtaking and the colliding overtaking are considered and discussed. The complex structure of wave systems and the aerodynamic forces are discussed. By comparing with the results acquired with the original methods, it is shown that the developed unified zonal method can be an effective and simple tool for handling complex bodies and enhancing the accuracy of numerical simulation.
    Symplectic Algorithm and Simulation of Solitons for Two-dimensional Non-stationary Sine-Gordon Equation
    JIANG Chang-jin
    2003, 20(4): 321-325. 
    Abstract ( )   PDF (162KB) ( )  
    A 2×7992-order nonlinear Hamiltonian system of two-dimensional non-stationary Sine-Gordon equation is introduced when the five point difference scheme is used to discretize the differential operator L=(ə2)/(əx2)+(ə2)/(əy2) in the rectangle [-a,a]×[-a,a]. An iterative method is designed to solve the nonlinear system, which is formed by using the centered Euler scheme for the Hamiltonian system. The condition and the velocity of convergence for this method are given. Numerical examples for evaluating one-soliton and two-soliton of the Sine-Gordon equation show that the symplectic method is an efficient algorithm.
    Absorption Boundary Condition of the Electromagnetic Wave in the Finite-difference Time-domain Method
    XU Han, CHANG Wen-wei, YIN Yan
    2003, 20(4): 326-330. 
    Abstract ( )   PDF (210KB) ( )  
    The Lindman's absorption boundary condition, which has been widely used in the FDTD method, is generalized to a special form suitable to the electromagnetic field calculation from the starting point of the absorption principle. Our 2D numerical simulations indicate that the P-polarized and S-polarized electromagnetic waves with a wide varity of incident angles can be well absorbed by left and right boundaries. The reflectivity of the electromagnetic energy at either boundary is less than 0.5%.
    Velocity Acceptance Performance of Superconducting Cavities in Proton Linacs
    YU Qing-chang
    2003, 20(4): 331-334. 
    Abstract ( )   PDF (190KB) ( )  
    The proton energy gain in multicell π-mode accelerating cavities and transit time factors of the cavities are given. On the basis the velocity acceptance of the superconducting cavities is studied by means of analytical method and numerical method, and the problem about the section-division of superconducting proton linacs are discussed.
    Applications of Turbulence Models in Simulation of Complex Flows
    XIAO Zhi-xiang, LI Feng-wei, E Qin
    2003, 20(4): 335-340. 
    Abstract ( )   PDF (437KB) ( )  
    Four turbulence models:algebraic Baldwin-Lomax model with Degani-Schiff modification, two versions of Johnson-King model (J-K90A and J-K92) and two-equation k-g model, are described and evaluated for missile supersonic flow, NASA TN D-712 standard model and the civil airplane wing-body configuration (two-block C-O mesh) transonic flows. The 3-D compressible Reynolds averaged Navier-Stokes (RANS) equations are integrated numerically by central-difference with artificial viscosity scheme, finite volume formulation and explicit multi-step Runge-Kutta algorithm. The turbulence equations of k-g model are explicitly solved in the same way as the RANS. Results show that all models can perform well for attached and mildly separated flows. For large angle-of-attack separated flow over the missile, k-g and J-K92 models match the experimental results much better than those with B-L and J-K90A models. According to the B-L model with some modification, the shock location logs behind those with the other three models for the civil airplane geometry. Model k-g, because of its advantages such as without using normal-to-wall distance, simple source terms and straightforward boundary conditions, performs best in the four models for complex configuration with multi-block mesh system or multi-wall interference. B-L and J-K models, despite of cheapness and robustness, both base on empirical Prandtl mixing-length hypothesis and need to calculate the distance normal to the wall, undoubtedly limiting the application scope of these models for multi-block mesh system and complex geometry.
    Hybrid Krylov Subspace Method and Its Application
    LIU Xing-ping, HANG Xu-deng, FU Shang-wu
    2003, 20(4): 341-344. 
    Abstract ( )   PDF (211KB) ( )  
    A new type of hybrid iterative method is presented, which is competent in solving the large scale sparse linear systems derived from the 2-dimensional 3-temperature radiation dynamic energy equations. Numerical results show that the new method is as 4 times fast as the old ones.Especially on the cases the old one doesn't converge, the new method can easily get the solution to the precision required. It can successfully complete the simulation and the final physical parameters of the simulation coincide with the theory.
    Method for Solving Transient Heat and Moisture Transfer in Multilayered Building Elements
    SU Xiang-hui, ANG Hai-song, XU Feng
    2003, 20(4): 345-350. 
    Abstract ( )   PDF (229KB) ( )  
    Finite-difference method is proposed for the first time for solving the transient heat and moisture transfer problems in 1-D multilayered building elements.Transient boundary conditions are specified for both heat and moisture transfer at external surfaces, while constant values are assumed for material properties of each layer of the structure and the discontinuity of material properties for each layer is considered.The calculations for a typical three-layer wall are conducted to obtain the dynamic distributions of the relative humidity at the interfaces of every layer during 24 hours by use of the mentioned method,and the results agree with those of the transfer function method which is empolyed in reference[6].The method can be extended to 2-D and 3-D coupled heat and moisture transfer problems.
    Numerical Simulation of Current and Temperature Fields of Aluminum Reduction Cells Based on ANSYS
    LI Jie, CHENG Ying-jun, LAI Yan-qing, ZHOU Nai-jun
    2003, 20(4): 351-355. 
    Abstract ( )   PDF (321KB) ( )  
    Finite element equations to calculate current and temperature fields of aluminum reduction cells are deduced using Galerkin method.The finite element model of anode and molten electrolyte is built according to multiple elements and multiple properties of ANSYS software.With reasonable assumption of boundary conditions,the current and temperature fields of 160?KA prebaked reduction cells are computed and the temperature,voltage and electric current distributions of the cells are analyzed.The simulation results of the model well coincide with the design data,and therefore provide foundations for optimizing current aluminum electrolysis cells and developing new type cells.
    Lattice Boltzmann Model for Nonlinear Chemical Waves in the Excitable Media
    YAN Guang-wu
    2003, 20(4): 356-358. 
    Abstract ( )   PDF (213KB) ( )  
    Lattice Boltzmann model for the nonlinear chemical waves is proposed. Using Chapman-Enskog expansion and multiscale technique, higher order moments of an equilibrium distribution function is obtained. As an example, the Selkov reaction is simulated with one scratching mark, obtaining a classical chemical patterns, the spiral wave. Numerical test shows that the method can be used to simulate the nonlinear chemical waves.
    Alternating Direction Implicit Method for Solving the Two-dimensional Time-dependent Photon Diffusion Equation
    ZHANG Zhi, LUO Qing-ming, ZENG Shao-qun, ZHANG Xin-yu, HUANG De-xiu
    2003, 20(4): 359-362. 
    Abstract ( )   PDF (175KB) ( )  
    An alternating direction implicit (ADI) method for solving the two-dimensional time-dependent photon diffusion equation is studied. Numerical comparisons with analytical solution and the results of Monte Carlo simulation are achieved. Numerical results show that the method has high accuracy and good stability.
    Lattice Dynamics of Body-centered Cubic Transition Metals with Analytic EAM Interatomic Potentials
    KONG Yi, HUANG Yang-cheng, SHU Xiao-lin, WANG Ling-ling, HU Wang-yu
    2003, 20(4): 363-368. 
    Abstract ( )   PDF (220KB) ( )  
    Using the Modified Analytic Embedded Atom Method(MAEAM),the phonon spectra and lattice specific heat of bcc metals(Cr,Fe,W,Mo,Ta,V,Nb) are calculated.The obtained results are in good agreement with the experimental data.It is improved greatly as compared with the calculated results by the Johnson model.And the influence on EAM model of the characteristics of bonding orientation of these transition bcc metals is discussed.
    Theoretical Calculation of CoⅠ-like Rh18+ Ions 3p63d9-3p53d10,3p63d9-3p63d84p Transition
    CHEN Ming-lun
    2003, 20(4): 369-371. 
    Abstract ( )   PDF (113KB) ( )  
    Transitions and oscillator strength for highly ionized CoⅠ-like Rh18+ ions 3p63d9-3p53d10,3p63d9-3p63d84p are calculated by using Cowan's multiconfiguration self-consistent field method together with the fitting formula proposed here.The calculations are in agreement with observations.
    Brief Introduction of MEDICI Software and Its Application in Ionization Radiation Effects
    GUO Hong-xia, CHEN Yu-sheng, ZHOU Hui, ZHANG Yi-men, GONG Ren-xi, LÜ Hong-liang
    2003, 20(4): 372-376. 
    Abstract ( )   PDF (251KB) ( )  
    The MEDICI software of two dimensional simulation of semiconductor device is briefly introduced.Its characteristics and use are described.With the help of MEDICI,effects of total dose for MOSFET and dose rate for pn junction are simulated.Their physical models are set up.A theory approach is provided for the research of ionization radiation effects.
Share: