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

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

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

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

MEMS稀薄气体内部流动模拟中的信息保存法

首先综述了处理低速稀薄气体流动的一些方法: 线化Boltzmann方程方法、Lattice Boltzmann方法(LBM)、加滑移边界的Navier-Stokes方程、以及DSMC方法, 并讨论它们在模拟MEMS中过渡领域低速流动特别是内部流动所遇到的困难, 其中表明了LBM现有方案不适合模拟过渡领域中的MEMS流动问题. 信息保存(IP)法通过保存一个模拟分子所代表的大量分子的平均信息,克服了流速低使得信息噪声比小而引起统计模拟的困难. 本文给出了方法的一些理论证实. MEMS中内部流动的特点, 即流速低和大的长宽比的特点, 引起椭圆性问题, 即出入口边界条件相互影响需要协调的问题. 通过对(长约几千微米的)微槽道流动应用IP方法的算例,演示了采用守恒形式的质量守恒方程和超松弛法可成功地解决这一问题. 借助同样的方法,用IP方法求解了真实长度(1\,000\,$\mu$m)硬盘驱动器读写头在过渡领域的薄膜支撑问题, 压力分布与具有严格气体动理论基础的概括化Reynolds方程完全相符, 而在此之前, DSMC方法只对短的读写头(5\,$\mu$m)与Reynolds方程做了校验. 作者建议将原来用于求解读写头润滑问题的Reynolds方程退化来求解过渡领域中的微槽道流动问题, 从而提供了一个有严格气体动理论品性的检验方法来验证求解MEMS内部流动的各种方法.      

气体动理学统一算法中相容性条件不满足引起的数值误差及其影响研究

求解玻尔兹曼(Boltzmann) 模型方程的气体动理学统一算法(unified gas kinetic scheme,UGKS) 是为模拟存在显著稀薄气体效应流动而建立的. 在该方法中,如果速度空间离散采用传统的离散速度坐标法(discreteordinate method,DOM),将会导致相容性条件得不到严格满足,从而引入数值误差. 本文从理论分析及数值试验两方面说明了该数值误差,正比于来流马赫数,反比于来流努森数. 引入了守恒型的离散速度坐标法(conservativediscrete ordinate method,CDOM),在离散层面上确保了相容性条件得到严格满足. 圆柱绕流计算结果表明,来流马赫数较高、努森数较小时,相容性条件满足与否对计算结果影响较大,采用CDOM 可以在较稀的速度空间网格上得到网格无关解,缩减计算量最大可达2/3.      

稀薄流非线性模型方程离散速度坐标法有限差分解

从一般非线性Bo ltzm ann方程出发,发展并实现了一套适于大范围K nudsen数稀薄流问题数值模拟的统一算法。采用BGK模型和Shakov模型近似碰撞项,进而引入两个二速度无量纲简化分布函数,通过关于分子速度第三分量取矩积分,将三速度单一模型方程变换为二速度微分方程组。基于G auss-H erm ite积分公式和正交多项式G auss积分公式,借助离散速度坐标法消除简化模型方程对分子速度空间的连续依赖性,从相空间到物理空间得到一组带源项双曲守恒离散方程,并给出其显式和隐式二阶迎风TVD有限差分解。以二维圆柱A r气体超声速绕流算例,验证了数值算法的有效性,比较分析了漫反射和镜面反射两种气体分子壁面反射模型的计算结果。     

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

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

基于Shan-Chen模型的格子Boltzmann方法在微流动模拟研究中的应用

对格子Boltzmann方法的本质及Shan-Chen模型的核心机制进行了全面阐述, 并从应用实例角度对基于Shan-Chen模型的格子Boltzmann方法在微流动模拟方面的有效性、适应性进行了详细分析. 结果表明, Shan-Chen模型易于耦合微观条件下占主导作用的微观力, 拓宽了格子Boltzmann方法在微流动模拟方面的应用. 同时, Shan-Chen模型在润湿性边界条件表征方面的优势, 使得这种方法在微结构表面的滑移效应模拟方面具有很好的应用前景.     

基于多块对接网格的隐式气体运动论统一算法应用研究

基于玻尔兹曼模型方程的气体运动论统一算法(gas kinetic unified algorithm,GKUA) 给出了一种能模拟从连续流到自由分子流跨流域空气动力学问题的途径. 该算法采用传统计算流体力学技术将分子运动和碰撞解耦处理,若采用显式格式将受格式稳定条件限制,在模拟超声速流动尤其是近连续流和连续流区的流动时计算效率较低. 为了提高计算效率,扩展其工程实用性,采用上下对称高斯-赛德尔(LU-SGS) 方法和有限体积法构造了求解玻尔兹曼模型方程的隐式方法,同时在物理空间采用能处理任意连接关系的多块对接网格技术. 通过模拟近连续过渡区并排圆柱绕流问题,计算结果与直接模拟蒙特卡洛方法模拟值吻合较好,验证了该方法用于跨流域空气动力计算的可靠性与可行性.      

参考焓值法在高速槽道湍流中的修正

利用高来流马赫数为3, 5, 6, 7, 10的槽道湍流直接数值模拟(direct numerical simulations, DNS)数据, 评估和修正经典的参考焓值法. 研究表明在高来流马赫数槽道湍流中, 经典参考焓值法预测的壁面热流与DNS结果相差很大, 需要作适当的修正.修正参考焓值法Ⅰ和Ⅱ的预测结果明显优于经典参考焓值法;并且修正参考焓值法Ⅱ更加适用于高马赫数流动, 其壁面热流与DNS结果的相对误差在10%以内. 同时, 修正参考焓值法Ⅱ的普适性在超声速燃烧室隔离段热环境试验中得到了验证.     

Gas-kinetic numerical schemes for one- and two-dimensional inner flows

Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.     

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.     

Gas-kinetic unified algorithm for plane external force-driven flows covering all flow regimes by modeling of Boltzmann equation

The nonequilibrium steady gas flows under the external forces are essentially associated with some extremely complicated nonlinear dynamics, due to the acceleration or deceleration effects of the external forces on the gas molecules by the velocity distribution function. In this article, the gas-kinetic unified algorithm (GKUA) for rarefied transition to continuum flows under external forces is developed by solving the unified Boltzmann model equation. The computable modeling of the Boltzmann equation with the external force terms is presented at the first time by introducing the gas molecular collision relaxing parameter and the local equilibrium distribution function integrated in the unified expression with the flow state controlling parameter, including the macroscopic flow variables, the gas viscosity transport coefficient, the thermodynamic effect, the molecular power law, and molecular models, covering a full spectrum of flow regimes. The conservative discrete velocity ordinate (DVO) method is utilized to transform the governing equation into the hyperbolic conservation forms at each of the DVO points. The corresponding numerical schemes are constructed, especially the forward-backward MacCormack predictor-corrector method for the convection term in the molecular velocity space, which is unlike the original type. Some typical numerical examples are conducted to test the present new algorithm. The results obtained by the relevant direct simulation Monte Carlo method, Euler/Navier-Stokes solver, unified gas-kinetic scheme, and moment methods are compared with the numerical analysis solutions of the present GKUA, which are in good agreement, demonstrating the high accuracy of the present algorithm. Besides, some anomalous features in these flows are observed and analyzed in detail. The numerical experience indicates that the present GKUA can provide potential applications for the simulations of the nonequilibrium external-force driven flows, such as the gravity, the electric force, and the Lorentz force fields covering all flow regimes.     

A coordinate transformation method for numerical solutions of the electric double layer and electroosmotic flows in a microchannel

The electric double layer (EDL) and electroosmotic flows (EOFs) constitute the theoretical foundations of microfluidics. Numerical solution is one of the effective means of analysis in microfluidics. In general, it is difficult to obtain an accurate numerical solution of complex EOFs because of multiphysical interactions and locally high gradients. In this paper, a new coordinate transformation method is proposed to numerically solve the Poisson–Boltzmann, Navier–Stokes and Nernst–Planck equations to study the EDL and complex EOFs in a microchannel. A series of numerical examples is presented including cases of a homogeneous, discontinous wall electric potential and a locally high wall potential. A systematic comparison of numerical solutions with and without the coordinate transformation is carried out. The numerical results indicate that the coordinate transformation effectively decreases the gradient of the electric potential, ion concentration and electroosmotic velocity in the vicinity of the solid wall, and greatly improves the stability and convergency of the solution. In a transformed coordinate system with a coarse grid, the numerical solutions can be as accurate as those in the original coordinate system with a refined grid. Copyright © 2011 John Wiley & Sons, Ltd.     

A high‐order compact finite‐difference lattice Boltzmann method for simulation of steady and unsteady incompressible flows

A high‐order compact finite‐difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth‐order compact FD scheme, and the temporal term is discretized with the fourth‐order Runge–Kutta scheme to provide an accurate and efficient incompressible flow solver. A high‐order spectral‐type low‐pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also conducted to evaluate the effects of grid size, filtering, and procedure of boundary conditions implementation on accuracy and convergence rate of the solution. The accuracy and efficiency of the proposed solution procedure based on the CFDLBM method are also examined by comparison with the classical LBM for different flow conditions. Two test cases considered herein for validating the results of the incompressible steady flows are a two‐dimensional (2‐D) backward‐facing step and a 2‐D cavity at different Reynolds numbers. Results of these steady solutions computed by the CFDLBM are thoroughly compared with those of a compact FD Navier–Stokes flow solver. Three other test cases, namely, a 2‐D Couette flow, the Taylor's vortex problem, and the doubly periodic shear layers, are simulated to investigate the accuracy of the proposed scheme in solving unsteady incompressible flows. Results obtained for these test cases are in good agreement with the analytical solutions and also with the available numerical and experimental results. The study shows that the present solution methodology is robust, efficient, and accurate for solving steady and unsteady incompressible flow problems even at high Reynolds numbers. Copyright © 2014 John Wiley & Sons, Ltd.     

Efficient simulation of multidimensional continuum and non-continuum flows by a parallelised unified gas kinetic scheme solver

We develop an efficient, parallel, gas-kinetic solver for computing both continuum and non-continuum flows over non-Cartesian geometries by utilising the unified gas kinetic scheme (UGKS). UGKS, however, requires the computationally expensive update of a six-dimensional phase space at each time step restricting its application to canonical, laminar problems and simple geometries. In this paper, we demonstrate that the applications of UGKS can be increased by parallelising it and combining it with a recently developed, Cartesian grid method (UGKS-CGM). We demonstrate that our Cartesian grid methodology as well as UGKS parallelization perform and scale well on a range of numerical test cases even for a very large number of cores. Finally, we demonstrate that the solver accurately computes canonical turbulence at low Knudsen numbers. These results demonstrate that the parallelised UGKS code can be utilised to effectively study the non-equilibrium effects of rarefaction on laminar and turbulent non-continuum flows.     

Difference scheme and numerical simulation based on mixed finite element method for natural convection problem

ntroductionLetΩ R2 beaboundeddomain .Weconsiderthefollowingnon_stationarynaturalconvectionproblem :Problem (Ⅰ )　Findu =(u1,u2 ) ,p ,andTsuchthat,foranyt1>0 ,ut- μΔu +(u· )u + p=λjT　　 ((x ,y ,t) ∈Ω× (0 ,t1) ) ,divu =0　　　　　　　　　 ((x ,y,t) ∈Ω× (0 ,t1) ) ,Tt-ΔT +λu· T =0　　 ((x,y,t) ∈Ω× (0 ,t1) ) ,u =0 ,T =0　　　　　　 ((x,y,t)∈ Ω× (0 ,t1) ) ,u(x ,y ,0 ) =0 ,　T(x,y,0 ) =f(x,y)　　 ((x,y) ∈Ω) ,whereuisthefluidvelocityvectorfield ,pthepressurefield ,Tthet…