气体化学反应流动的DSMC／EPSM混合算法研究

发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的：EPSM耦合模型，并通过混合参数进行流区的自动识别，将：EPSM方法与蒙特卡罗直接模拟方法(OSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的二维高超音速竖板绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

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

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

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

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

高速流场气动光学RANS/DSMC混合数值模拟算法研究

提出了适用于高速流场气动光学数值模拟的RANS/DSMC混合算法.通过RANS对全局时均流场进行数值模拟,再对其中局部流场的脉动量采用DSMC进行数值模拟,以Maxwellian速度分布实现RANS宏观量信息向DSMC微观量信息的传递.采用超声速环境下尖劈模型对该混合算法进行校验,对比试验结果论证了算法的正确性.在计算...     

绕圆柱体跨声速流场

本文将保角曲线坐标方法应用于绕物体外部的跨声速流动.文中讨论了圆柱体跨声速绕流,计算了来流马赫数M_∞。为亚临界,超临界和M_∞为l的流场.给出了不同来流马赫数下柱面马赫数和压力分布,柱面前端中心流线上马赫数分布.给出了超临界绕流时不同来流马赫数下的声速线和来流马赫数M_∞从0.2,0.3到1的等马赫数线分布.本文部分结果与已有的理论结果进行了比较.本方法计算简单,精度好,是计算厚体跨声速绕流的有效方法,适用于不同形状的柱体和机翼的跨声速绕流.     

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

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

梯形断面明渠丁坝绕流水力特性三维数值模拟

采用两相流混合模型,并选取RNG k-ε湍流模型封闭两相流时均方程,对梯型断面明渠非淹没式丁坝绕流水力特性进行了三维数值模拟。采用有限体积法离散计算区域,求解速度与压力耦合方程组时使用半隐式SIMPLE(Semi-Implicit Method for Pressure-Linked Equations)算法,模拟自由水面时采用了VOF(Volume of Fluid)法。对丁坝后不同的回流长度进行了分析比较,并将模拟结果与实测资料进行了对比验证,结果表明两者吻合较好,相对误差小于8%,说明该模型能够很好地模拟明渠丁坝绕流的水力特性分布规律。     

大尺度网格下DSMC仿真碰撞概率的修正方法

从DSMC碰撞概率的实现途径出发,通过引入最小二乘法及亚松弛的思想给出一种大尺度网格下分子仿真碰撞概率的修正方法,将该方法应用于非结构四面体网格下的DSMC模拟中,计算分析了文献中钝锥模型的绕流结果,数值结果初步验证了本文方法的合理性和有效性。     

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模拟值等结果对比分析,验证统一算法模拟自由分子流到连续流再入过程高超声速绕流问题的可靠性与精度.     

超音速混合层稳定性分析及增强混合的研究

利用流动稳定性提高超音速混合层的混合效率,对于提高超音速流的高效混合是一个有效途径。研究结果表明,有展向曲率的三维混合层中,三维扰动的增长率很大,且法向的掺混能力也较强,可以有效地增强混合。对于高马赫数来流的超音速混合层,这一特性依然存在,这将有利于提高高超音速混合层的混合能力。     

A hybrid DSMC/Navier–Stokes frame to solve mixed rarefied/nonrarefied hypersonic flows over nano‐plate and micro‐cylinder

We extend a hybrid DSMC/Navier–Stokes (NS) approach to unify the DSMC and the NS simulators in one framework capable of solving the mixed non‐equilibrium and near‐equilibrium flow regions efficiently. Furthermore, we use a one‐way state‐based coupling (Dirichlet–Dirichlet boundary‐condition coupling) to transfer the required information from the continuum region to the rarefied one. The current hybrid DSMC–NS frame is applied to the hypersonic flows over nanoflat plate and microcylinder cases. The achieved solutions are compared with the pure DSMC and NS solutions. The results show that the current hybrid approach predicts the surface heat transfer rate and shear stress magnitudes very accurately. Some important conclusions can be drawn from this study. For example, although the shock wave region would be a non‐equilibrium region, it is not necessary to use a pure DSMC simulator to solve it entirely. This is important when the researchers wish to predict the surface properties such as velocity slip, temperature jump, wall heat flux rate, and friction drag magnitudes accurately. Our investigation showed that our hybrid solution time would be at least 40% (for the flat plate) and 35% (for the cylinder) of the time that must be spent by a pure DSMC solver to attain the same accuracy.Copyright © 2013 John Wiley & Sons, Ltd.     

A novel efficient hybrid DSMC–dynamic collision limiter algorithm for multiscale transitional flows

A new 2D parallel multispecies polyatomic particle–based hybrid flow solver is developed by coupling the Direct Simulation Monte Carlo (DSMC) method with a novel Dynamic Collision Limiter (DCL) approach to solve multiscale transitional flows. The hybrid DSMC‐DCL solver can solve nonequilibrium multiscale flows with length scales ranging from continuum to rarefied. The DCL method, developed in this work, dynamically assigns different number of collisions in cells, which is based on the local value of K‐S parameter such that the number of collisions per time step is limited in near‐equilibrium flow regions. Present hybrid solver uses the Kolmogorov‐Smirnov statistical test as the continuum breakdown parameter, based on which, the solution domain is decomposed into near‐equilibrium and nonequilibrium flow regions. Direct Simulation Monte Carlo is used where nonequilibrium flow regions are encountered, while the DCL method is used where flow regions are found to be in near‐equilibrium state. In this work, we have studied hypersonic flow of nitrogen over a blunt body with an aerospike and supersonic flow of argon through a micronozzle. The results obtained by the hybrid DSMC‐DCL solver are compared and shown to agree well with the experimental data and with those obtained from DSMC, with significant savings in the computational cost.     

二维自适应非结构网格DSMC并行算法研究

研究了二维自适应非结构网格DSMC并行算法实现的过程.首先提出了一类非结构网格自适应策略,有效降低了网格尺度对计算结果的影响,提高了流场的分辨率;然后基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,利用分区并行思想,设计了非结构网格DSMC并行算法,节约了计算时间.利用For-tran90的动态分配内存技术编制了通用计算程序;最后对过渡流域高超声绕流进行了数值模拟,计算结果初步验证了算法的可行性与有效性.     

微喷管流动粒子仿真与Navier-Stokes数值模拟比较研究

运用DSMC（Direct Simulation Monte—Carlo）方法从分子运动论层次对大膨胀比、喉部转角为尖角的微喷管流动现象进行模拟,考察来流总压对喷管性能的影响,并与Navier—Stokes方程运算结果、实验结果进行比较。研究表明：在模拟微型喷管的流动现象时,DSMC方法比N—S方程更加适用。     

Extension of the DSMC method to high pressure flows

A collision-limiter method, designated as equilibrium direct simulation Monte Carlo (eDSMC), is proposed to extend the DSMC technique to high pressure flows. The method is similar to collision-limiter schemes considered in the past with the important distinction that for inviscid flows, equilibrium is enforced in the entire flow by providing a sufficient number of collisions, based on pre-simulation testing. To test the method with standard DSMC and Navier–Stokes (NS) methods, axi-symmetric nozzle and embedded-channel flows are simulated and compared with experimental temperature data and pre-existing calculations, respectively. The method is shown to agree with third-order Eulerian nozzle flows and first-order channel flows. Chapman–Enskog theory is utilized to predict the range of initial conditions where eDSMC is potentially useful for modeling flows that contain viscous boundary layer regions. Comparison with supersonic nozzle data suggests that the eDSMC method is not adequate for capturing the large variation in flow length scales occurring in supersonic expansions into a vacuum. However, when eDSMC is used in combination with the baseline-DSMC method a near-exact solution is obtained with a considerable computational savings compared to the exact DSMC solution. Viscous flow channel calculations are found to agree well with an exact Navier–Stokes (NS) calculation for a small Knudsen number case as predicted by Chapman–Enskog theory.     

不同形状微管道中的气体流动

不同形状微尺度管道(圆形、六边形、半圆形、不同宽高比的矩形)中的气体流动特性是微机电系统设计最为关心的问题之一.文中利用信息保存(IP)方法和直接模拟Monte Carlo(DSMC)方法进行研究,给出两种方法的计算结果相互符合,并与其它研究者的BGK模型方程计算结果进行了比较.对于微尺度管道中关心的低Mach数流动, IP方法的统计收敛效率明显优于DSMC方法.通过拟合IP和DSMC结果,给出了圆形、六边形、半圆形、不同宽高比的矩形截面情况下无量纲质量流率与等效Knudsen数的关系.     

Investigation on a coupled Navier–Stokes–Direct Simulation Monte Carlo method for the simulation of plume flowfield of a conical nozzle

To assess the plume effects of space thrusters, the accurate plume flowfield is indispensable. The plume flow of thrusters involves both continuum and rarefied flow regimes. Coupled Navier–Stokes–Direct Simulation Monte Carlo (NS–DSMC) method is a major approach to the simulation of continuum‐rarefied flows. An axisymmetric coupled NS–DSMC solver, possessing adaptive‐interface and two‐way coupling features, is investigated in this paper for the simulation of the nozzle and plume flows of thrusters. The state‐based coupling scheme is adopted, and the gradient local Knudsen number is used to indicate the breakdown of continuum solver. The nitrogen flows in a conical nozzle and its plume are chosen as the reference case to test the coupled solver. The threshold value of the continuum breakdown parameter is studied based on both theoretical kinetic velocity sampling and coupled numerical tests. Succeeding comparisons between coupled and full DSMC results demonstrate their conformities, meanwhile, the former saves 58.8% computational time. The pitot pressure evaluated from the coupled simulation result is compared with the experimental data proposed in literatures, revealing that the coupled method makes precise predictions on the experimental pitot pressure. Copyright © 2014 John Wiley & Sons, Ltd.     

Computation of head–disk interface gap micro flowfields using DSMC and continuum–atomistic hybrid methods

This paper discusses computational modeling of micro flow in the head–disk interface (HDI) gap using the direct simulation Monte Carlo (DSMC) method. Modeling considerations are discussed in detail both for a stand‐alone DSMC computation and for the case of a hybrid continuum–atomistic simulation that couples the Navier–Stokes (NS) equation to a DSMC solver. The impact of the number of particles and number of cells on the accuracy of a DSMC simulation of the HDI gap is investigated both for two‐ and three‐dimensional configurations. An appropriate implicit boundary treatment method for modeling inflow and outflow boundaries is used in this work for a three‐dimensional DSMC micro flow simulation. As the flow outside the slider is in the continuum regime, a hybrid continuum–atomistic method based on the Schwarz alternating method is used to couple the DSMC model in the slider bearing region to the flow outside the slider modeled by NS equation. Schwarz coupling is done in two dimensions by taking overlap regions along two directions and the Chapman–Enskog distribution is employed for imposing the boundary condition from the continuum region to the DSMC region. Converged hybrid flow solutions are obtained in about five iterations and the hybrid DSMC–NS solutions show good agreement with the exact solutions in the entire domain considered. An investigation on the impact of the size of the overlap region on the convergence behavior of the Schwarz method indicates that the hybrid coupling by the Schwarz method is weakly dependent on the size of the overlap region. However, the use of a finite overlap region will facilitate the exchange of boundary conditions as the hybrid solution has been found to diverge in the absence of an overlap region for coupling the two models. Copyright © 2009 John Wiley & Sons, Ltd.