Loading...

Archive

    25 January 2012, Volume 29 Issue 1 Previous Issue    Next Issue

    For Selected: Toggle Thumbnails
    Finite Volume Scheme for Diffusion Equation with Vertex Reconstruction
    ZHAO Qiang, YUAN Guangwei, DONG Zhiwei
    2012, 29(1): 1-9. 
    Abstract ( )   PDF (4678KB) ( )  
    Nine.point scheme for diffusion equation on distorted quadrilateral is discussed.Based on a classical nine。point scheme, and with continuity of diffusion flux, a formula for designed for computing Vertex unknowns. Numerical experiments show that the scheme is robust,and the accuracy is higher than that of other weighted methods for diffusion problem with continuous or discontinuous coefficient on distorted meshes.
    A Kind of Conservative Remapping Algorithms Based on RBF Interpolations
    ZHAO Xiaojie, ZHAO Ning, WANG Donghong
    2012, 29(1): 10-16. 
    Abstract ( )   PDF (4460KB) ( )  
    In order to ensure conservative precision and monotonicity of a remapping algorithm and discontinuity with high resolution, a new kind of conservative RBF remapping algorithm is developed based on radial based function interpolation. Conservative precision is tested with a conservative error remapping algorithm. Numerical results are obtained for smooth and discontinuous functions. Compared with other remapping methods,the method gives good numerical results.
    Parallel Algorithm for Brownian Configuration Fields with Finite Volume Method
    ZHANG Huifeng, OUYANG Jie, DAI Xiangyan
    2012, 29(1): 17-24. 
    Abstract ( )   PDF (398KB) ( )  
    In a FENE bead-spring chain model,a parallel algorithm for Brownian configuration fields with finite volume method is used to simulate pressure driven pipe nows, velocity-pressure driven pipe nows with same and opposite directions. Validity of the algorithm is verified with numerical examples. Velocity,stress and stretch of molecular chain of pipe nows are obtained. It shows that the proposed algorithm overcomes shortcomings of huge computation and time-consuming for Brownian configuration fields. And it has excellent expansibility. Parallel efficiency is over 85% as number of beads is greater than two. It also shows that the time of now to reach steady state is shorter with increasing number of chain segments.
    Numerical Shock Instability for 2-D Shallow Water Equations
    SHEN Zhijun, HU Lijun, YAN Wei
    2012, 29(1): 25-35. 
    Abstract ( )   PDF (5391KB) ( )  
    In calculation of multidimensional fluid mechanics problems with numerical schemes that accurately capture contact discontinuity, perturbation near shock wave may increase dramatically. This is called numerical shock instability. In this paper numerical shock instability on shallow water equations is studied. By analyzing linear stability of several numerical schemes, marginal stability of schemes are found having close relation with numerical shock instability. According to eigenvalue analysis, a hybrid method is designed to remedy nonphysical phenomenon by locally modify the original schemes. Numerical experiments show efficiency and robustness of HLLC-HLL hybrid scheme in eliminating shock instability of shallow water equations.
    Interface Treating Method for Multi-medium Flow on Triangular Meshes
    SONG Peichang, WANG Chunwu
    2012, 29(1): 36-42. 
    Abstract ( )   PDF (6903KB) ( )  
    Simulation of multi-medium compressible flows is generally on structured meshes. As triangular meshes provide a flexible way for decomposing a complex domain, we present an interface treating method for multi-medium flows on triangular meshes with a rGFM (real ghost fluid method). A level set function is used to capture interface. A Riemann problem is constructed at interface to provide fluid states. By defining interface boundary conditions with these states, calculation can be carried out as if in a single medium flow. A high-accuracy RKDG method is then employed to discretizing equations. Several examples are given to test robustness and efficiency of the algorithm. It shows that the method can capture sharp interface and shock wave accurately.
    Bubble Meshing Method for Two-parametric Surface
    ZHANG Weiwei, NIE Yufeng, WANG Lei
    2012, 29(1): 43-50. 
    Abstract ( )   PDF (6901KB) ( )  
    For mesh generation of a two-parameter surface, anisotropic and non-uniform node placement method with bubble simulation is applied to optimize node distribution in parameter area. Then the parameter area is meshed with constrained Delaunay triangulation. Finally, according to the mapping method, two-parametric surface mesh is obtained. A second order Riemann metric tensor determines distribution of nodes in the parameter area. It could be co-generated with a three-dimensional surface metric tensor and gradient of surface functions. Numerical examples show that the node placement method with bubble simulation can generate node set meeting requirements of Riemann metric in parameter area. Nodes are meshed and mapped back into the surface. A high quality surface mesh is obtained.
    A Gas-kinetic Scheme for Compressible Two-phase Flow Containing Non-Conservative Products
    SHE Bangwei, ZHAO Guiping
    2012, 29(1): 51-57. 
    Abstract ( )   PDF (308KB) ( )  
    We propose a method for BN two-phase flow containing non-conservative terms. In microscopic aspect, a GKS model which matches with the BN model is constructed. Numerical fluxes are obtained with integrating particle velocity distribution function. And non-conservative terms are explicitly introduced into evolution and construction of numerical fluxes. Non-conservative terms are well handled. This method avoids complex iterative process of Riemann solver. It has high calculation efficiency. Numerical tests show high precision.
    A Fast Communication Algorithm for Parallel Structured Mesh Applications
    LIU Xu, ZHANG Aiqing, XIAO Li, MO Zeyao
    2012, 29(1): 58-64. 
    Abstract ( )   PDF (4602KB) ( )  
    Communication algorithms need to exchange data between adjacent subdomains. The detection of adjacent subdomains is an import task for communication algorithms, which can be done by solving an intersection problem of boxes. An interval tree algorithm is proposed for intersection problems. Taking advantage of structured mesh applications, it is able to maintain 0 (Nlog N) time complexity. Numerical experiments show that the algorithm is able to achieve high computational efficiency and good scalability. It supports parallel computing of miUion-box scale.
    Numerical Study of the Effect of Viaduct Noise Barrier Height on Train Aerodynamic Characteristics
    LUO Jianbin, HU Aijun
    2012, 29(1): 65-72. 
    Abstract ( )   PDF (8966KB) ( )  
    With computational fluid dynamics method, a numerical study on effect of viaduct noise barrier height on aerodynamic characteristics of a high speed train was made. High calculation precision is obtained with suitable meshing, appropriate turbulent model and boundary conditions. It shows that as moving in leeward the side force on head train or the middle train decreases gradually with increasing of noise barrier height. A rolling moment on the head is most and it fails as noise harrier height goes up. The roiling moment on middle train running in windward is greater than that in leeward as noise barrier height is increased. Under two running conditions, the train relative position in cavity flow is different and surface pressure distribution is changed. Accordingly, aerodynamic force or moment on train body is varied.
    A 7.62 mm Rifle Bullet Impacting on a 30CrMnSiA Steel Target Plate:SPH-FEM Conversion Algorithm
    ZHANG Zhichun, QIANG Hongfu, FU Xuejin, XIA Wei
    2012, 29(1): 73-81. 
    Abstract ( )   PDF (7972KB) ( )  
    Conversion criterion and time step control are corrected for smoothed particle hydrodynamics-finite element method (SPH-FEM) conversion algorithm. Impact of a 7.62 mm rifle bullet against a special heat treated 30CrMnSiA steel target plate is simulated in full-size in 3D using the corrected SPH-FEM conversion algorithm, in which distorted finite elements of the bullet are converted into SPH particles. An elastic-plasticity model is used for bullet, and corrected Johnson-Cook and Gruneisen EOS are used for target. Two destroyed modes of disc-shape plastic deformation and plugging perforation are calculated for different initial velocities of bullet. Good agreement between numerical results and experimental observations shows that the corrected SPH-FEM conversion algorithm could use finite elements and SPH particles efficiently, which provides an effective tool for simulation of a soft core bullet impacting on a hard plate.
    Using Super-high Energy Neutrons to Detect Inertial Confinement Fusion
    LI Shu, TIAN Dongfeng, DENG Li
    2012, 29(1): 82-86. 
    Abstract ( )   PDF (278KB) ( )  
    In inertial confinement fusion (ICF) target, collision between fusion neutrons may produce super-high energy neutrons. A relationship between velocity of fusion and ratio of generated super-high energy neutron is derived. An ideal fusion model is simulated by numerical method. It indicates that the ratio of generated super-high energy neutron is increased sharply as burning volume of DT is decreased. With this information, 2D-effect and mixing-effect during compression of target are deduced with escaped super-high energy neutrons. Furthermore, a relationship between velocity of burning T and ratio of generated super-high energy neutron is concluded. It shows that using super-high energy neutron to detect ICF is effective.
    Simulation of Convective Heat Transfer in Magnetic Microencapsulated Phase Change Material Slurries in a Circular Tube
    HUANG Yong, XUAN Yimin, LI Qiang
    2012, 29(1): 87-94. 
    Abstract ( )   PDF (2571KB) ( )  
    A mathematical model is established to describe convective heat transfer in magnetic microencapsulated phase change material slurry (MMPCMS) flows in a circular tube under an external magnetic field. Influences of magnetic field strength, volume fraction of magnetic phase change materials (MMPCM), mass flow rate, and heat flux on convective heat transfer are analyzed. It shows that convective heat transfer of slurries is obviously enhanced by external magnetic field. Since distributions of MMPCM volume fraction and slurry temperature are changed clue to magnetic force on MMPCM.
    Stochastic Multiscale Model for Dynamic Yield Behavior of Magnetorheological Fluids
    LI Jie, PENG Yongbo
    2012, 29(1): 95-100. 
    Abstract ( )   PDF (369KB) ( )  
    Dynamic yield behavior of magnetorheological fluids is investigated by upscaling information of microscale interaction between particles, employing a large-scale molecular dynamical simulation technique, to macroscale bulk behavior. We conduct a stochastic muhiscale model for dynamic yield of magnetorheological fluids based on equivalence of system energy at different scales. It is revealed that the dynamic yield exhibits nonlinear and stochastic fluctuations due to heterogeneity of sequence and number of cluster-sheet reconstructions under shear fields loading, as well as Brownian motion of suspensions with initial random conditions. Meanwhile, we investigate fluctuation of microscale particle motions, relationship between stress and strain, and constitutive relationship of shear rate. It is noted that the microscale thermal fluctuation is far more than macroscalc variations since the upscaling from microscale to macroscale results in degradation of fluctuations. Besides, the macroscale variations relies on external magnetic field as in the constitutive relationship of shear rate, i.e. Bingham model, which is supposed to be considered in the design and optimization of magnetorheological devices.
    Brownian Motion of Non-spherical Particles:Fluctuating-Lattice Boltzmann Investigation
    NIE Deming, LIN Jianzhong
    2012, 29(1): 101-107. 
    Abstract ( )   PDF (421KB) ( )  
    A fluctuating-lattice Boltzmann method is used for direct numerical simulations of Brownian motion of non-spherical particles. Numerical results include mean-square velocities and velocity autocorrelation functions of elliptical and rectangular particles. It showns that equi-partition theorem applies to non-spherical particles as well and energy of each degree of freedom depends on corresponding mean-square velocity or rotational velocity. Furthermore, translational temperature of Brownian particle is identical to its rotational temperature. Finally, velocity autocorrelation function of Brownian particle decays at relative long times as~ct-1, where coefficlent c is independent of particle shape.
    Thermoelectric Analysis of Interconnect Considering Via and Fringing Effects
    WANG Ning, DONG Gang, YANG Yintang, WANG Zeng, WANG Fengjuan, DING Can
    2012, 29(1): 108-114. 
    Abstract ( )   PDF (4679KB) ( )  
    With consideration of via effect and heat fringing effect, a thermoelectric simulation method is proposed which modifies node heat flow due to temperature distribution. Based on thermoelectric duality, thermal resistance models including inner/inter-layer and vias are presented. Take advantage of feedback relationship between heat and electric, the node network heat flow model is modified with temperature distribution. Multilevel interconnects temperature distribution with polymer and silicon oxide as insulator dielectric are analyzed. Compared with results of finite element, the relative standards deviation of the proposed method can be reduced by 71.2% and 12. 9% respectively than those of available models. With consideration of via effect and heat fringing effect, we calculate peak temperature rise in different technology nodes. It shows that interconnect temperature distribution is overestimated in traditional models.
    Improved Calculation of Magnetic Field in Simulation of Guiding-Center Motion
    XIE Qing, YANG Weihong
    2012, 29(1): 115-120. 
    Abstract ( )   PDF (2426KB) ( )  
    In simulation of a charged particle's guiding-center orbit in Tokamak equilibrium configurations, an efficient way to build physical functions such as magnetic field was obtained. Discrete Fourier transform is applied to replace cubic spline interpolation and B spline interpolation methods in calculation of gulding-center's trajectory. It has same accuracy as cubic spline interpolation and facilitates calculations in tearing mode as well. The method holds naturally property of symmetry. Derivatives of 3rd order or above can be obtained easily. Moreover, effect of incorrect data is reduced to the minimal limit compared with cubic spline interpolation.
    Finite Element Method with Diffusion Synthetic Acceleration for Short Pulsed Laser in Turbid Media
    AN Wei, ZHU Tong
    2012, 29(1): 121-126. 
    Abstract ( )   PDF (371KB) ( )  
    A finite element model is developed to simulate short pulsed laser in turbid media by solving a transient radiative transfer equation. The model employs diffuse synthetic acceleration to reduce computation time and improve solution efficiency, especially for scattering dominated media, It indicates that the model should adopt analytical formulation to deal with scattering term produced by pulsed laser since it accurately simulates sharp changed wave-front without any numerical diffusion and oscillation. Computation time of diffusion synthetic acceleration is about 50% to 60% that of source iteration.
    Analysis of Terahertz Quantum Cascade Laser Beam
    WANG Jian, WU Weidong, ZHANG Xiaoli, DUAN Suqing
    2012, 29(1): 127-132. 
    Abstract ( )   PDF (2490KB) ( )  
    Threshold gain of laser mode in terahertz quantum cascade laser is calculated with finite element method. It shows that contact layer thickness and doping density have far greater impact on threshold gain than waveguide width and lasing wavelength. As thickness and density becoming smaller (greater), threshold gain of TM1 (TM0) mode becomes smaller. Far-field characteristics of the beam are analyzed using vector diffraction theory. The beam spot is basically elliptical. With increase of waveguide width or lasing wavelength, far-field divergence angles at x direction decrease or increase linearly. Though the contact layer thickness and doping density for TM0 mode are different from that for TM1, their far-field divergence angles at x direction are same. Moreover, far-field divergence angles at y direction are unaffected by waveguide width, and less affected by lasing wavelength. In terms of threshold gain and beam quality, TM1 is more desirable than TM0 mode.
    Comb Filtering Properties of Total Reflection Tunnelling in Cylindrical Photonic Crystals
    LIU Qineng
    2012, 29(1): 133-138. 
    Abstract ( )   PDF (282KB) ( )  
    To study comb filtering properties of tunnelling in total reflection of light in 1 D cylindrical photonic crystals, transition matrix method is adopted. Transmissivity of light is studied as incident angle is greater than full reflection angle. Tunnelling under condition of total reflection is found in transmitting waves. Good comb filtering characteristics is shown. Responses of tunnelling versus period, quantum number of mode, and cylinder radius are obtained. They are basis for design of high-quality comb filters.
    Structure and Potential Energy Function of PCIx (X=1,2) Molecular Ground State
    HAN Xiaoqin, XIAO Xiajie, LIU Yufang
    2012, 29(1): 139-144. 
    Abstract ( )   PDF (2433KB) ( )  
    QCISD method is used to optimize possible ground-state structures of PCl and PCl2 molecule with multiple basis sets. Dissociation energy, harmonic frequency and force constants are calculated with 6-311 + G (dr) and 6-311G (dr) basis sets, respectively. They agree with experimental observations well. Potential energy functions of PCl2 are derived from many-body expansion theory. Potential energy curves verify structure characteristics and potential depth. Cl+ PCl and P + ClCl based on molecular reaction potential energy surface are discussed. It describes molecular reaction dynamics successfully.
    Generation and Properties of Dark Solitons in Heteronuclear Two-component Bose-Einstein Condensates
    ZHOU Yanzhen, ZHANG Suying, HAN Wei
    2012, 29(1): 145-151. 
    Abstract ( )   PDF (8803KB) ( )  
    Generation and properties of dark solitons are studied in detail by means of numerical simulation of quasi one-dimensional heteronuclear two-component Bose-Einstein condensates in syntonic potential. Dark solitons are only induced in one of the condensate components by modulational instability with instantaneous conversion of repulsive interspecies interaction to attraction. Solitons pass through each other periodically in syntonie potential. In addition, the number of dark solitons is affiliated with ratio of particle mass and ratio of particles numbers.
    General Multiphase Matching for Grover Quantum Search Algorithm
    JIN Wenliang, CHEN Xiangdong
    2012, 29(1): 152-158. 
    Abstract ( )   PDF (403KB) ( )  
    Since a quantum system is inevitably influenced by some unpredictable perturbations, we thereby conclude that all the experimental realizations of Grover quantum search algorithm reported were, in fact, achieved in a three-dimensional complex subspace. We also prove that in a two-dimensional complex subspace, for any given initial superposition of basis states|γ0>=cosβ0|α>+sinβ0e|β)(β0 is a small positive real number, ζ is an arbitrary real number), there exists a set of solutions Fj={(θj,θj-1,…,θ1),(φj,φj-1,…,φ1)} such that a desired state can be found with certainty for some positive integer j≥2, where the phase rotation angles θl andθt are real numbers but not equal to 2k'π,1 ≤ 1 ≤ j,k'is an arbitrary integer. If it is only required that a desired state can be found with high success probability, then as the total number of the desired and undesired states in an unsorted database is sufficiently large the above set of solutions Fj can be written in the form  for a relatively small positive integer j.
Share: