基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

Gas-kinetic numerical method for solving mesoscopic velocity distribution function equation

A gas-kinetic numerical method for directly solving the mesoscopic velocity distribution function equation is presented and applied to the study of three-dimensional complex flows and micro-channel flows covering various flow regimes. The unified velocity distribution function equation describing gas transport phenomena from rarefied transition to continuum flow regimes can be presented on the basis of the kinetic Boltzmann–Shakhov model equation. The gas-kinetic finite-difference schemes for the velocity distribution function are constructed by developing a discrete velocity ordinate method of gas kinetic theory and an unsteady time-splitting technique from computational fluid dynamics. Gas-kinetic boundary conditions and numerical modeling can be established by directly manipulating on the mesoscopic velocity distribution function. A new Gauss-type discrete velocity numerical integration method can be developed and adopted to attack complex flows with different Mach numbers. HPF parallel strategy suitable for the gas-kinetic numerical method is investigated and adopted to solve three-dimensional complex problems. High Mach number flows around three-dimensional bodies are computed preliminarily with massive scale parallel. It is noteworthy and of practical importance that the HPF parallel algorithm for solving three-dimensional complex problems can be effectively developed to cover various flow regimes. On the other hand, the gas-kinetic numerical method is extended and used to study micro-channel gas flows including the classical Couette flow, the Poiseuille- channel flow and pressure-driven gas flows in two-dimensional short micro-channels. The numerical experience shows that the gas-kinetic algorithm may be a powerful tool in the numerical simulation of micro-scale gas flows occuring in the Micro-Electro-Mechanical System (MEMS). The project supported by the National Natural Science Foundation of China (90205009 and 10321002), and the National Parallel Computing Center in Beijing. The English text was polished by Yunming Chen.     

GAS-KINETIC UNIFIED ALGORITHM FOR BOLTZMANN MODEL EQUATION IN ROTATIONAL NONEQUILIBRIUM AND ITS APPLICATION TO THE WHOLE RANGE FLOW REGIMES

基于过去开展稀薄自由分子流到连续流气体运动论统一算法框架,采用转动惯量描述气体分子自旋运动,确立含转动非平衡效应各流域统一玻尔兹曼模型方程.基于转动能量对分布函数守恒积分,得到计及转动非平衡效应气体分子速度分布函数方程组,使用离散速度坐标法对分布函数方程所依赖速度空间离散降维;应用拓展计算流体力学有限差分方法,构造直接求解分子速度分布函数的气体动理论数值格式;基于物面质量流量通量守恒与能量平衡关系,发展计及转动非平衡气体动理论边界条件数学模型及数值处理方法,提出模拟各流域转动非平衡效应玻尔兹曼模型方程统一算法.通过高、低不同马赫数1：5～25氮气激波结构与自由分子流到连续流全飞行流域不同克努森数（9×10^-4～10）Ramp制动器、圆球、尖双锥飞行器、飞船返回舱外形体再入跨流域绕流模拟研究,将计算结果与有关实验数据、稀薄流DSMC模拟值等结果对比分析,验证统一算法模拟自由分子流到连续流再入过程高超声速绕流问题的可靠性与精度.     

转动非平衡玻尔兹曼模型方程统一算法与全流域绕流计算应用

基于过去开展稀薄自由分子流到连续流气体运动论统一算法框架,采用转动惯量描述气体分子自旋运动,确立含转动非平衡效应各流域统一玻尔兹曼模型方程.基于转动能量对分布函数守恒积分,得到计及转动非平衡效应气体分子速度分布函数方程组,使用离散速度坐标法对分布函数方程所依赖速度空间离散降维;应用拓展计算流体力学有限差分方法,构造直接求解分子速度分布函数的气体动理论数值格式;基于物面质量流量通量守恒与能量平衡关系,发展计及转动非平衡气体动理论边界条件数学模型及数值处理方法,提出模拟各流域转动非平衡效应玻尔兹曼模型方程统一算法.通过高、低不同马赫数1:5～25氮气激波结构与自由分子流到连续流全飞行流域不同克努森数（9×10-4～10）Ramp制动器、圆球、尖双锥飞行器、飞船返回舱外形体再入跨流域绕流模拟研究,将计算结果与有关实验数据、稀薄流DSMC模拟值等结果对比分析,验证统一算法模拟自由分子流到连续流再入过程高超声速绕流问题的可靠性与精度.      

Investigation of the rarefied gas flow around a circular cylinder in stationary and oscillatory regimes

Numerically, on the basis of the Krook kinetic equation, the rarefied gas flow around a circular cylinder is investigated in stationary and oscillatory regimes. The flows around a rotating cylinder and a cylinder with a nonuniformly heated surface are considered. The Knudsen numbers at which the lift force acting on the rotating cylinder changes sign are calculated. It is shown that at low Knudsen numbers a lift force acts on the nonuniformly heated cylinder.     

气体运动论数值算法在微槽道流中的应用研究

简要介绍基于Boltzmann模型方程的气体运动论数值算法基本思想及其对二维微槽道流动问题数值计算的推广，并阐述适用于微尺度流动问题的气体运动论边界条件数值处理方法。通过对压力驱动的二维微槽道流动问题进行数值模拟，将不同Knudsen数下的微槽道流计算结果分别与有关DSMC模拟值和经滑移流理论修正的N—S方程解进行比较分析，表明基于Boltzmann模型方程的气体运动论数值算法对微槽道气体流动问题具有很好的模拟能力。     

Kinetic analysis of longitudinal Couette flow between coaxial cylinders

On the basis of a model kinetic equation, the rarefied gas flow between coaxial circular cylinders, of which the outer one is at rest while the inner one travels along its symmetry axis at a constant velocity, is studied. The problem is solved numerically in both the linear and nonlinear formulations by an implicit conservative method of second-order accuracy. The effect of the rarefaction, the cylinder radius ratio, and the inner cylinder velocity on the flow parameters is investigated. The limits of applicability of the linearized kinetic equation are established.     

Numerical investigation from rarefied flow to continuum by solving the Boltzmann model equation

Based on the Bhatnagar–Gross–Krook (BGK) Boltzmann model equation, the unified simplified velocity distribution function equation adapted to various flow regimes can be presented. The reduced velocity distribution functions and the discrete velocity ordinate method are developed and applied to remove the velocity space dependency of the distribution function, and then the distribution function equations will be cast into hyperbolic conservation laws form with non‐linear source terms. Based on the unsteady time‐splitting technique and the non‐oscillatory, containing no free parameters, and dissipative (NND) finite‐difference method, the gas kinetic finite‐difference second‐order scheme is constructed for the computation of the discrete velocity distribution functions. The discrete velocity numerical quadrature methods are developed to evaluate the macroscopic flow parameters at each point in the physical space. As a result, a unified simplified gas kinetic algorithm for the gas dynamical problems from various flow regimes is developed. To test the reliability of the present numerical method, the one‐dimensional shock‐tube problems and the flows past two‐dimensional circular cylinder with various Knudsen numbers are simulated. The computations of the related flows indicate that both high resolution of the flow fields and good qualitative agreement with the theoretical, DSMC and experimental results can be obtained. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical Analysis of the Influence of Thermal Stresses on the Nonlinear Thermomolecular Pressure Difference Effect

On the basis of a numerical analysis of the non-Navier-Stokes gas-dynamic equations for slow non-isothermal gas flows, the nonlinear thermomolecular pressure difference effect due to a large temperature gradient along the lateral surface of a capillary is investigated. It is shown that the magnitude of the effect is substantially different from the values calculated using the Navier-Stokes equations. For two models of molecular interaction (Maxwell molecules and hard spheres), the possibility of a quasi-one-dimensional interpretation of the effect for experimental estimation purposes is demonstrated. The solutions of the relaxation kinetic equation for flow in a plane capillary at small Knudsen numbers and the gas-dynamic equations for slow non-isothermal flows are compared and the range of their applicability is estimated.     

A novel slightly compressible model for low Mach number perfect gas flow calculation

By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows. The project supported by the Basic Research on Frontier Problems in Fluid and Aerodynamics in China and the National Natural Science Foundation of China (19772069)     

A hybrid Eulerian–Lagrangian method to simulate the dispersed phase in turbulent gas-particle flows

This paper presents a methodology to combine stochastic Lagrangian approach and continuum model to simulate the dispersed phase in gas-particle turbulent flows using that both approaches are based on the same Boltzmann-like kinetic equation governing the joint fluid-particle probability density function (pdf). The proposed hybrid method is based on the separate application of each approach in two adjacent domains and their coupling at the interface via flux boundary conditions. Validation of the method is carried out for non-colliding solid particles suspended in homogeneous turbulent shear flow without two-way coupling.     

Unified gas kinetic scheme combined with Cartesian grid method for intermediate Mach numbers

We develop a method to seamlessly simulate flows over a wide range of Knudsen numbers past arbitrarily shaped immersed boundaries. To achieve seamless computation, ie, not use any zone division to distinguish between continuum and non‐continuum regions, we use the unified gas kinetic scheme (UGKS), which is based on the Bhatnagar‐Groos‐Krook (BGK) approximation of the Boltzmann equation. We combine UGKS with an appropriately designed Cartesian grid method (CGM) to allow us to compute flows past arbitrary boundaries. The CGM we use here satisfies boundary conditions at the wall by using a constrained least square interpolation procedure. However, it differs from the usual, continuum CGMs in 2 ways. Firstly, to allow us capture non‐continuum effects at the boundaries, the CGM used herein interpolates the microscopic velocity distribution function in addition to the macroscopic variables. Secondly, even for the macroscopic variables, we use a gas kinetic method–based density interpolation procedure at the boundaries that allows the CGM to interface well with the UGKS method. We demonstrate the robustness and efficacy of the method by testing it on stationary immersed boundaries at various Knudsen numbers ranging from continuum to transition regimes.     

Nonlinear nonequilibrium kinetic model of the Boltzmann equation for a gas with power-law interaction

A model kinetic equation approximating the Boltzmann equation on a wide range of the intensities of nonequilibrium states of gases is derived to describe rarefied gas flows. The kinetic model is based on a distribution function dependent on the absolute velocity of gas particles. Themodel kinetic equation possesses a high computational efficiency and the problem of shock wave structure is solved on its basis. The calculated and experimental data for argon are compared.     

Multiphase flow model developed for simulating gas hydrate transport in horizontal pipe

The hydrate formation or dissociation in deep subsea flow lines is a challenging problem in oil and gas transport systems. The study of multiphase flows is complex while necessary due to the phase changes (i.e., liquid, solid, and gas) that occur with increasing the temperature and decreasing the pressure. A one-dimensional multiphase flow model coupled with a transient hydrate kinetic model is developed to study the characteristics of the multiphase flows for the hydrates formed by the phase changes in the pipes. The multiphase flow model is derived from a multi-fluid model, while has been widely used in modelling multiphase flows. The heat convection between the fluid and the ambient through the pipe wall is considered in the energy balance equation. The developed multiphase flow model is used to simulate the procedure of the hydrate transport. The results show that the formation of the hydrates can cause hold-up oscillations of water and gas.     

Weakly rarefied gas flows between nonuniformly heated parallel plates

The aim of this research is to establish the validity of the predictions of the theory of slow nonisothermal flows, to study the limits of applicability (with respect to the Knudsen number) of the conclusions reached and to determine the effect of the Knudsen layers on these flows on the basis of a numerical investigation of slow nonisothermal weakly rarefied gas flow in a plane infinite channel with weakly nonequilibrium heating of the walls and a finite wall temperature difference. The gas flow is described by a relaxation transport equation. The results obtained show how quickly, as the Knudsen number decreases, the solutions of the transport equation outside the Knudsen layers tend to the solution of the equations of gas dynamics of slow nonisothermal flows (and not to the solution of the Navier-Stokes equations).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 115–121, January–February, 1988.     

Numerical investigation of rarefied diatomic gas flows through a plane channel into a vacuum

Rarefied flows through a plane microchannel into a vacuum are numerically investigated on the basis of the model kinetic equation for a diatonic gas (nitrogen). The dependence of the gas flow rate through the channel on the Knudsen number, the wall temperature, and the channel length is determined. The energy flux transferred to the cold diatomic gas from the hot channel walls is calculated. The results for diatomic and monatomic gases are compared.