Loading...

Archive

    25 January 2018, Volume 35 Issue 1 Previous Issue    Next Issue

    For Selected: Toggle Thumbnails
    A Robust Riemann Solver Without Artificial Intervention
    REN Jian, SHEN Zhijun, YAN Wei, YUAN Guangwei
    2018, 35(1): 1-12.  DOI: 10.19596/j.cnki.1001-246x.7630
    Abstract ( )   HTML ( )   PDF (20342KB) ( )  
    In order to improve robustness of multi-fluid-channel scheme on average of volume (MFCAV) and overcome its artificial intervention in actual application, a new HLLCM scheme was designed on a moving mesh, which restrains non-physical mesh tangling without artificial intervention. Numerical results show that the HLLCM scheme maintains one-dimensional spherical symmetry on equal-angle-zoned grids which has better numerical effects in keeping mesh quality and energy conservation in complex applications than MFCAV scheme.
    An Improved Third-order WENO-Z+3 Scheme and Its Application
    XU Weizheng, KONG Xiangshao, WU Weiguo
    2018, 35(1): 13-21.  DOI: 10.19596/j.cnki.1001-246x.7575
    Abstract ( )   HTML ( )   PDF (12882KB) ( )  
    To improve dissipation characteristic and resolution for flow features of classical third-order WENO scheme (WENO-JS3). An improved third-order WENO-Z+ scheme (NWENO-Z+3) was constructed with a different global smoothness indicator. Several classical one-dimensional Riemann problems and double Mach reflection case were simulated to investigate computing performance of NWENO-Z+3 scheme. It indicated that NWENO-Z+3 scheme have better characteristics than schemes such as WENO-JS3, WENO-Z and WENO-Z+3 scheme with low numerical dissipation and high resolution to flow features. NWENO-Z+3 scheme was also used to investigate propagation and wave evolution of blast waves generated by cylindrical high pressure gas in single cabin and connected cabins. It indicated that NWENO-Z+3 scheme was suitable for simulating evolution of blast waves, which contains high pressure ratio and high density ratio.
    Second-order Local-Bound-Preserving Conservative Remapping on Unstructured Polyhedral Meshes
    XU Xihua, LIU Na, CHEN Yibing
    2018, 35(1): 22-28.  DOI: 10.19596/j.cnki.1001-246x.7573
    Abstract ( )   HTML ( )   PDF (4845KB) ( )  
    We present conservatively remapping cell-centered variables from one mesh to another with second-order accuracy and boundary-preservation. It is generally applicable to any polyhedral source or target mesh. The algorithm consists of four parts:A least square based polynomial reconstruction of physical gradient; an octree-based fast donor-cell searing algorithm; a convex hull algorithm for intersection of polyhedrons and a modifying procedure for local bound preservation. The remapping scheme is scalable, second-order accurate and enjoys bound preservation property. Various benchmark problems demonstrate these properties. Numerical results show that it takes hundreds seconds to remap physical variables on tessellation with hundreds thousands to millions polyhedrons.
    Euler Numerical Methods for Reactive Flow with General Equation of States in Two Dimensions
    LIU Jingjing, ZENG Xianyang, NI Guoxi
    2018, 35(1): 29-38.  DOI: 10.19596/j.cnki.1001-246x.7594
    Abstract ( )   HTML ( )   PDF (7010KB) ( )  
    We present an efficient method to simulate reactive flow for general equation of states in two dimensions. Two kinds of nonideal equation of states for compressible fluid model coupling with reactive rate equation are concerned. The important aspect is to deal with mixture of different phase in one cell, which inevitably happens in Euler method for reactive flows. Physical variables such as pressure,velocity and speed of sound in each cell are reconstructed which results in nonlinear algebra equations. They are used to obtain flux by HLLC Riemann solver. Numerical examples of detonation in two dimensions with different equation of states demonstrate accuracy and robustness of the method.
    MRT-LBM Analysis of Acoustic Streaming in Standing Waves Between Two-dimensional Flat Plates
    ZHOU Tian, LI Xuemin, LIU Feng
    2018, 35(1): 39-46.  DOI: 10.19596/j.cnki.1001-246x.7558
    Abstract ( )   HTML ( )   PDF (6188KB) ( )  
    Acoustic streaming generated by standing waves between two-dimensional flat plates is analyzed based on multiple relaxation time lattice Boltzmann method (MRT-LBM). Simulation results are in good agreement with approximate analytical solutions of Rayleigh streaming. Effects of viscosity and width between two plates on acoustic streaming in standing waves are studied. Dimensionless velocity distributions of x component with different viscosity at x=L/4 and dimensionless velocity distributions of y component with different viscosity at x=L/2,and effects of widths on boundary layer vortex thickness and formation progress of acoustic streaming in standing waves are revealed. It demonstrates availability of MRT-LBM model in simulating acoustic streaming in standing waves.
    Numerical Simulation of Free-Rising of Bubbles in a Tube with Fins Using Front Tracking Method
    LU Min, DONG Boheng, ZHANG Ying, ZENG Liang, LIN Xiangquan, DU Peng
    2018, 35(1): 47-54.  DOI: 10.19596/j.cnki.1001-246x.7588
    Abstract ( )   HTML ( )   PDF (2121KB) ( )  
    Direct numerical simulation of FTM (Front Tracking Method) was used to study free rising motion of single bubble in a tube with fins in gravity. Ratio of fin height and pitch to tube length of a side are chosen as geometric characteristics. Morton number is used as characteristic parameter of fluid. With analysis on free rising motion of bubbles at different fin height, pitch and Morton number effects of geometric characteristics of tube and parameters of fluid on motion trajectory of bubbles are studied. It shows that under effect of fins bubble shape is asymmetry and bubbles have lateral displacement like a "snake".This phenomenon is associated with geometric structure of fin and fluid around bubble. The higher fin height is, the more obvious phenomenon becomes. The smaller the viscosity of fluid around bubble is, the greater the horizontal displacement is.
    Nonlinear Percolation Law in Low Permeability Fissure Cave Reservoir with Fractal Dimension
    LIU Huapu, LIU Huiqing, WANG Jing
    2018, 35(1): 55-63.  DOI: 10.19596/j.cnki.1001-246x.7589
    Abstract ( )   HTML ( )   PDF (1795KB) ( )  
    Fractal approach is introduced. Mathematical model of fractal triple medium fissure cave reservoir with low velocity non-Darcy effect is established. By means of Laplace transform and Stehfest numerical algorithm, bottom hole pressure solution is obtained. Seepage characteristics analysis and nonlinear parameter sensitivity analysis are conducted. Finally, combined with actual well data fractal model is verified. It shows that seepage process of fractal triple medium fissure cave reservoir can be divided into 6 seepage stages. They are presence of wellbore storage, transition flow, channeling of caves to fractures, pseudo radial flow, channeling of matrix to fractures and caves, total radial flow, respectively. Furthermore, effect of velocity non-Darcy effect on percolation law is gradually increasing with time. The higher starting pressure gradient, the higher pressure dynamic curves in total radial flow stage bend upward. Fractal coefficient affects whole percolation process. With increase of fractal coefficient, fracture tortuosity and flow resistance increase and whole pressure dynamic curves move up.
    Temperature Distribution Modeling in SAGD Circulation Phase
    QI Peng, LIU Huiqing, PANG Zhanxi, LIU Huapu, CHEN Yu
    2018, 35(1): 64-70.  DOI: 10.19596/j.cnki.1001-246x.7583
    Abstract ( )   HTML ( )   PDF (1694KB) ( )  
    Temperature pattern in circulation phase of SAGD process is studied. Two in-situ points, three-dimensional heat conduct and partial penetration of horizontal well into reservoirs, are considered. The model is processed by dimensionless treatment, Laplace transformation, separation of variables and Stehfest numerical inversion. The results are in accordance with simulations by CMG software which verifies accuracy of the model. It shows relationship between temperature and heating time and concludes that the smaller the radial distance, the less the time is needed. Calculation plate can be helpful to design heating time and amount of heat.
    Mechanical Properties of Graphene/Hydroxyapatite Composite Materials: Numerical Study
    TANG Qiuming, ZHEN Tianyi, LI Danyun, GAO Qiang
    2018, 35(1): 71-76.  DOI: 10.19596/j.cnki.1001-246x.7595
    Abstract ( )   HTML ( )   PDF (1884KB) ( )  
    We study effect of addition of graphene on mechanical properties (elastic modulus and Poisson's ratio) of hydroxyapatite composites. Firstly, generating algorithm of random distribution geometric model of graphene/hydroxyapatite composites was created and model's automatic generation program was developed. Then, finite element model of graphene/hydroxyapatite composite was established. Finite element analysis of uniaxial tensile tests of composite materials was carried out. Finally, mechanical properties of different mass fraction of graphene of composite were calculated. Validity of the algorithm is verified by experimental data. It shows that elastic modulus of composite increased 12% to 50% with addition of 0.25% to 1.25% (wt.) graphene. Mechanical properties of hydroxyapatite can be effectively improved by addition of graphene.
    Mixed Displacement-Pressure Collocation Method for Plane Elastic Problems
    WANG Zhaoqing, XU Zikang
    2018, 35(1): 77-86.  DOI: 10.19596/j.cnki.1001-246x.7585
    Abstract ( )   HTML ( )   PDF (777KB) ( )  
    Introducing a pressure variable, governing equations of elasticity are expressed as displacements and pressure coupled system of partial differential equations. Barycentric interpolation is applied to approximate unknown functions. Matrix-vector forms of discrete expressions of governing equations for plane elastic problems are obtained by using barycentric interpolation differentiation matrices. Discrete boundary conditions of displacements and pressure are obtained by using barycentric interpolation. Boundary conditions are imposed by additional method to form an over-constrained linear algebra equation system of plane elastic problem. Numerical solutions of displacement for plane elastic problem are solved with least-square method. Numerical examples illuminate efficiency and computing precision of the method.
    Calculation of Kinetic Parameters with Continuous Energy Adjoint Weighted Monte Carlo Method
    PENG Gang
    2018, 35(1): 87-94.  DOI: 10.19596/j.cnki.1001-246x.7570
    Abstract ( )   HTML ( )   PDF (677KB) ( )  
    A study of continuous energy adjoint weighted Monte Carlo method is carried out in several kinetic parameters, including effective delayed neutron fraction, prompt neutron generation time and prompt neutron decay constant, which is based on Iterated Fission Probability (IFP) method and selection of adjoint flux is extented. Difference among tracklength, collision, absorption and tally estimator are compared. Covariance and variance weight are also adopted based on these estimators. Meanwhile, different distribution domain of kinetic parameters are provided. IFP method is studied in more detail. Next event estimate and IFP method with different latent generation number are compared and influence on result are obtained. Fitness number of latent generation number is obtained. Particle kind and distribution information are automatically listed in output file.
    Monte Carlo Simulation of Interaction of Electrons with Anode in Microstructure X-Ray Tubes
    XU Xiaodong, ZHOU Bin, GUO Jinchuan, YI Minghao, XIAO Feihu
    2018, 35(1): 95-102.  DOI: 10.19596/j.cnki.1001-246x.7586
    Abstract ( )   HTML ( )   PDF (1540KB) ( )  
    With multiple scattering theory of electrons we studied energy deposition of electrons within target and energy distribution on spot of X-ray production by Monte Carlo method as energetic electrons impinge into target vertically. It shows that spread of electrons within target is relevant to incident electron energy and target material. And 99% electron energy is deposited in a cylindrical zone. Energy deposition of electrons along incident direction increased with depth first, and then to a certain depth it decreased rapidly, which is in accordance with electron backscattering theory. These results illustrate that X-ray production zone is within a certain depth under target surface. In addition, simulation results of diamond indicate that diamond film is not recommended as an electron absorption grating, but as an excellent heat sink of target. It can be used as reference for design of microstructure X-ray tubes.
    First-principles Calculation of Mn-doped LiMgN Diluted Magnetic Semiconductor
    LI Peiyuan, WU Zhimin, YE Qian, CHEN Bo, WANG Chaoqiang, XU Jian, DU Chengxu
    2018, 35(1): 103-111.  DOI: 10.19596/j.cnki.1001-246x.7579
    Abstract ( )   HTML ( )   PDF (11507KB) ( )  
    With first-principles density functional theory, geometric structures of pure LiMgN, Mn-doped LiMgN, and Mn-doped LiMgN with excess or deficient of Li are geometrically optimized. Electronic structures, magnetic properties, and optical properties were calculated. It shows that Mn doping produces spin polarized impurity bands, which makes materials exhibit half metallic properties. Their properties are affected by stoichiometry of Li. Deficient of Li makes width of impurity band, net magnetic moments and Curie temperature decrease, while half metallic increase. Excess of Li improves width of impurity band, conductivity and Curie temperature, which makes half metallic and band gap decrease. Mn-doped systems have a new dielectric peak in low-energy region, which enhances absorption of low-frequency electromagnetic waves, and shows red-shift. There is an obvious change in complex refractive index function. Energy loss decreases and blue-shift appears only in Li-deficient compound.
    Structure and magnetic properties of Ni-doped ZnO clusters
    CHEN Hongxia, DU Sijie, ZHUANG Guoce
    2018, 35(1): 112-118.  DOI: 10.19596/j.cnki.1001-246x.7572
    Abstract ( )   HTML ( )   PDF (4034KB) ( )  
    Structural and magnetic properties of (ZnO)12 clusters doped with one (monodoped) and two (bidoped) Ni atoms were studied with a first-principles method. Substitutional, exohedral, and endohedral dopings are considered. Exohedral isomers are found the most favorable for both monodoped and bidoped clusters. Magnetic coupling between Ni atoms at the nearest neighbor position is mainly governed by competition between direct Ni-Ni antiferromagnetic interaction and ferromagnetic interaction between two Ni atoms via O atom due to strong p-d hybridization. Most importantly, exohedral and endohedral bidoped clusters favor ferromagnetic state, which has potential applications in nanoscale quantum devices.
    Key Nodes Identification of Power Grid Considering Local and Global Characteristics
    WANG Yi, ZOU Yanli, HUANG Li, LI Ke
    2018, 35(1): 119-126.  DOI: 10.19596/j.cnki.1001-246x.7580
    Abstract ( )   HTML ( )   PDF (4032KB) ( )  
    In order to identify key nodes in network effectively, a comprehensive index of node importance recognition which considers local and global characteristics of network is proposed. Nodes in weighted standard test system IEEE39 and IEEE118 are sorted according to importance based on this index. Ranking results are compared with those based on betweenness centrality method and point weight method. Meanwhile these important methods of node recognition are compared using structure-based network efficiency analysis and dynamics-based loss-of-synchronization diffusion time and synchronization capability. It shows that ranking of node importance based on comprehensive index is more reasonable, which is superior to node recognition method based on betweenness centrality and point weight.
Share: