This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity dis...This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.展开更多
In this paper,a gas-kinetic unified algorithm(GKUA)is developed to investigate the non-equilibrium polyatomic gas flows covering various regimes.Based on the ellipsoidal statistical model with rotational energy excita...In this paper,a gas-kinetic unified algorithm(GKUA)is developed to investigate the non-equilibrium polyatomic gas flows covering various regimes.Based on the ellipsoidal statistical model with rotational energy excitation,the computable modelling equation is presented by unifying expressions on the molecular collision relaxing parameter and the local equilibrium distribution function.By constructing the corresponding conservative discrete velocity ordinate method for this model,the conservative properties during the collision procedure are preserved at the discrete level by the numerical method,decreasing the computational storage and time.Explicit and implicit lower-upper symmetric Gauss-Seidel schemes are constructed to solve the discrete hyperbolic conservation equations directly.Applying the new GKUA,some numerical examples are simulated,including the Sod Riemann problem,homogeneous flow rotational relaxation,normal shock structure,Fourier and Couette flows,supersonic flows past a circular cylinder,and hypersonic flow around a plate placed normally.The results obtained by the analytic,experimental,direct simulation Monte Carlo method,and other measurements in references are compared with the GKUA results,which are in good agreement,demonstrating the high accuracy of the present algorithm.Especially,some polyatomic gas non-equilibrium phenomena are observed and analysed by solving the Boltzmann-type velocity distribution function equation covering various flow regimes.展开更多
In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys....In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys.285(2015),265-279] on uniform meshes,in this paper,in order to obtain the boundary fluxes based on the framework of unified gas kinetic scheme(UGKS),we use the real multi-dimensional reconstruction for the initial data and the macro-terms in the equation of the gray transfer equations.We can prove that the scheme is asymptotic preserving,and especially for the distorted quadrilateral meshes,a nine-point scheme [SIAM J.SCI.COMPUT.30(2008),1341-1361] for the diffusion limit equations is obtained,which is naturally reduced to standard five-point scheme for the orthogonal meshes.The numerical examples on distorted meshes are included to validate the current approach.展开更多
To directly incorporate the intermolecular interaction effects into the discrete unified gas-kinetic scheme(DUGKS)for simulations of multiphase fluid flow,we developed a pseudopotential-based DUGKS by coupling the pse...To directly incorporate the intermolecular interaction effects into the discrete unified gas-kinetic scheme(DUGKS)for simulations of multiphase fluid flow,we developed a pseudopotential-based DUGKS by coupling the pseudopotential model that mimics the intermolecular interaction into DUGKS.Due to the flux reconstruction procedure,additional terms that break the isotropic requirements of the pseudopotential model will be introduced.To eliminate the influences of nonisotropic terms,the expression of equilibrium distribution functions is reformulated in a moment-based form.With the isotropy-preserving parameter appropriately tuned,the nonisotropic effects can be properly canceled out.The fundamental capabilities are validated by the flat interface test and the quiescent droplet test.It has been proved that the proposed pseudopotential-based DUGKS managed to produce and maintain isotropic interfaces.The isotropy-preserving property of pseudopotential-based DUGKS in transient conditions is further confirmed by the spinodal decomposition.Stability superiority of the pseudopotential-based DUGKS over the lattice Boltzmann method is also demonstrated by predicting the coexistence densities complying with the van der Waals equation of state.By directly incorporating the intermolecular interactions,the pseudopotential-based DUGKS offers a mesoscopic perspective of understanding multiphase behaviors,which could help gain fresh insights into multiphase fluid flow.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 91016027 and 91130018)
文摘This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.
基金supported by the Project of manned space engineering technology(2018-14)“Large-scale parallel computation of aerodynamic problems of irregular spacecraft reentry covering various flow regimes”the National Natural Science Foundation of China(91530319).
文摘In this paper,a gas-kinetic unified algorithm(GKUA)is developed to investigate the non-equilibrium polyatomic gas flows covering various regimes.Based on the ellipsoidal statistical model with rotational energy excitation,the computable modelling equation is presented by unifying expressions on the molecular collision relaxing parameter and the local equilibrium distribution function.By constructing the corresponding conservative discrete velocity ordinate method for this model,the conservative properties during the collision procedure are preserved at the discrete level by the numerical method,decreasing the computational storage and time.Explicit and implicit lower-upper symmetric Gauss-Seidel schemes are constructed to solve the discrete hyperbolic conservation equations directly.Applying the new GKUA,some numerical examples are simulated,including the Sod Riemann problem,homogeneous flow rotational relaxation,normal shock structure,Fourier and Couette flows,supersonic flows past a circular cylinder,and hypersonic flow around a plate placed normally.The results obtained by the analytic,experimental,direct simulation Monte Carlo method,and other measurements in references are compared with the GKUA results,which are in good agreement,demonstrating the high accuracy of the present algorithm.Especially,some polyatomic gas non-equilibrium phenomena are observed and analysed by solving the Boltzmann-type velocity distribution function equation covering various flow regimes.
基金supported by the Science and Technology Development foundation of China Academy of Engineering Physics(Grant Nos.2015B0202041,2015B0202040)the Science and Technology Development foundation of China Academy of Engineering Physics(Grant 2015B0202040)+2 种基金the Science and Technology Development foundation of China Academy of Engineering Physics(Grant No.2015B0202033)for LiNSFC(Grant No.11371068)for SunNSFC(Grant No.11371068)for Zeng
文摘In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys.285(2015),265-279] on uniform meshes,in this paper,in order to obtain the boundary fluxes based on the framework of unified gas kinetic scheme(UGKS),we use the real multi-dimensional reconstruction for the initial data and the macro-terms in the equation of the gray transfer equations.We can prove that the scheme is asymptotic preserving,and especially for the distorted quadrilateral meshes,a nine-point scheme [SIAM J.SCI.COMPUT.30(2008),1341-1361] for the diffusion limit equations is obtained,which is naturally reduced to standard five-point scheme for the orthogonal meshes.The numerical examples on distorted meshes are included to validate the current approach.
基金National Numerical Wind Tunnel Project,the National Natural Science Foundation of China(No.11902266,11902264,12072283)111 Project of China(B17037).
文摘To directly incorporate the intermolecular interaction effects into the discrete unified gas-kinetic scheme(DUGKS)for simulations of multiphase fluid flow,we developed a pseudopotential-based DUGKS by coupling the pseudopotential model that mimics the intermolecular interaction into DUGKS.Due to the flux reconstruction procedure,additional terms that break the isotropic requirements of the pseudopotential model will be introduced.To eliminate the influences of nonisotropic terms,the expression of equilibrium distribution functions is reformulated in a moment-based form.With the isotropy-preserving parameter appropriately tuned,the nonisotropic effects can be properly canceled out.The fundamental capabilities are validated by the flat interface test and the quiescent droplet test.It has been proved that the proposed pseudopotential-based DUGKS managed to produce and maintain isotropic interfaces.The isotropy-preserving property of pseudopotential-based DUGKS in transient conditions is further confirmed by the spinodal decomposition.Stability superiority of the pseudopotential-based DUGKS over the lattice Boltzmann method is also demonstrated by predicting the coexistence densities complying with the van der Waals equation of state.By directly incorporating the intermolecular interactions,the pseudopotential-based DUGKS offers a mesoscopic perspective of understanding multiphase behaviors,which could help gain fresh insights into multiphase fluid flow.