The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate me...The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation.展开更多
The homotopy analysis method (HAM), as a new mathematical tool, has been employed to solve many nonlinear problems. As a fundamental equation in non-equilibrium statistical mechanics, the Boltzmann integro-differentia...The homotopy analysis method (HAM), as a new mathematical tool, has been employed to solve many nonlinear problems. As a fundamental equation in non-equilibrium statistical mechanics, the Boltzmann integro-differential equation (BE) describing the movement of particles is of strong nonlinearity. In this work, HAM is preliminarily applied to dilute granular flow which is relatively simple. By choosing the Maxwell velocity distribution function as the initial solution, the concrete expression of the first-order approximate solution to BE with collision term being the BGK model is given. Furthermore it is consistent with the solution using Chapman-Enskog method but does not rely on little parameters.展开更多
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.展开更多
基金Supported by the National Natural Science Foundation of China (Grant Nos. 90205009 and 10321002)the National Parallel Computing Center
文摘The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation.
基金Supported by the National Basic Research Program of China (Grant No. 2007CB714101)Research Fund of the State Key Laboratory for Hydroscience and Engineering in Tsinghua University (Grant No. 2008-ZY-6)
文摘The homotopy analysis method (HAM), as a new mathematical tool, has been employed to solve many nonlinear problems. As a fundamental equation in non-equilibrium statistical mechanics, the Boltzmann integro-differential equation (BE) describing the movement of particles is of strong nonlinearity. In this work, HAM is preliminarily applied to dilute granular flow which is relatively simple. By choosing the Maxwell velocity distribution function as the initial solution, the concrete expression of the first-order approximate solution to BE with collision term being the BGK model is given. Furthermore it is consistent with the solution using Chapman-Enskog method but does not rely on little parameters.
基金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.