基于Boltzmann模型方程的气体运动论HPF并行算法

从修正的BGK-Blotzmann模型方程出发，引入离散速度坐标技术对气体分子速度分量进行离散降维，基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，发展直接求解气体分子速度分布函数的气体运动论有限差分数值格式，提出了一套能有效模拟各流域三维绕流问题的气体运动论统一算法，在分析研究气体运动论数值算法内在并行度的基础上，开展各流域三维绕流问题统一算法的HPF（高性能FORTRAN）并行化程度设计，建立一套能有效模拟各流域复杂外形体绕流的HPF并行算法软件，并进行了不同Knudsen(克努森）数下三维球体绕流及类“神舟号”返回舱外形体绕流的初步数值实验，将计算结果与过渡区有关实验数据及各流域气体绕流现象进行比较分析，证实了发展的气体运动论HPF并行算法在求解稀薄流到连续流不同流域复杂绕流问题方面的可行性。     

离散统一气体动理学格式稳态辐射输运应用

工程应用中的介质热辐射问题是典型的多尺度问题. 基于Boltzmann输运方程建立的各类气体动理学格式, 在多尺度瞬态问题中得到了广泛应用. 为了克服显式求解方案中CFL条件等的限制, 文章通过气体动理学格式实现稳态辐射输运方程的直接求解. SDUGKS格式由离散统一气体动理学格式(discrete unified gas kinetic scheme, DUGKS)的核心思想发展而来, 应用于稳态问题计算. 将SDUGKS格式进一步拓展到多尺度的稳态热辐射输运计算. SDUGKS格式继承了DUGKS格式沿特征线离散实现的界面重构, 并通过隐式增量格式的单元更新实现对辐射强度的较正, 采用逐次迭代法将辐射强度渐近收敛到稳定值. 选用多组一维和二维不同尺度的辐射传热算例, 通过与特定的解析解以及其他数值方法结果对比, 检验了SDUGKS的计算精度和计算效率, 并论证了它在多尺度问题中的渐进保持性质.      

稀薄气体动力学矩方法研究综述

临近空间位于航天器入轨与返回的必经区域, 也是临近空间高超声速飞行器长航时飞行空域, 空间环境的特殊性决定了飞行器在穿越时必须考虑稀薄大气环境对飞行器气动力防隔热通讯及控制的影响.Boltzmann方程作为描述气体分子速度分布函数演化规律的微分-积分形式, 在一定条件下能够描述从自由分子流到连续流全流域流动现象.作为Boltzmann方程的宏观表达形式, 矩方程这一经典流体力学方程形式涵盖了Euler方程N-S方程Burnett方程super-Burnett方程及近年来发展的广义流体力学方程——非线性本构关系模型等.由于成熟的CFD数值计算理论及有限矩方程较高的计算效率, 滑移过渡流矩方法相比粒子仿真与Boltzmann模型方程方法具有十分显著的优势和巨大的工程应用潜力.因此, 对近年来传统及新型矩方法研究所取得的进展进行归纳总结, 并针对关键科学问题开展理论与数值计算方法研究, 具有十分重要的理论与工程应用价值.      

表征体元尺度渗流的离散统一动理学格式

结合表征体元尺度的通用渗流模型,提出离散统一动理学格式（DUGKS）渗流方法,分别用均匀网格和非均匀网格计算二维Poiseuille、Couette、方腔流等经典渗流问题,检验DUGKS渗流方法的有效性和非均匀网格应用的优势,将DUGKS渗流方法应用到裂缝系统中.     

跨流域高超声速绕流环境Boltzmann模型方程统一算法研究

如何准确可靠地模拟从外层空间高稀薄流到近地面连续流的航天器高超声速绕流环境与复杂流动变化机理是流体物理的前沿基础科学问题. 基于对Boltzmann方程碰撞积分的物理分析与可计算建模, 确立了可描述自由分子流到连续流区各流域不同马赫数复杂流动输运现象统一的Boltzmann模型速度分布函数方程, 发展了适于高、低不同马赫数绕流问题的离散速度坐标法和直接求解分子速度分布函数演化更新的气体动理论数值格式, 建立了模拟复杂飞行器跨流域高超声速飞行热环境绕流问题的气体动理论统一算法. 对稀薄流到连续流不同Knudsen数0.002 ≤Kn_∞ ≤1.618、不同马赫数下可重复使用卫星体再入过程(110–70 km)中高超声速绕流问题进行算法验证分析, 计算结果与典型文献的Monte Carlo直接模拟值及相关理论分析符合得较好. 研究揭示了飞行器跨流域不同高度高超声速复杂流动机理、绕流现象与气动力/热变化规律, 提出了一个通过数值求解介观Boltzmann模型方程, 可靠模拟高稀薄自由分子流到连续流跨流域高超声速气动力/热绕流特性统一算法.     

基于气体动理论的二维Karman涡街数值模拟

基于从稀薄流到连续流的跨流域气体动理论统一算法(gas-kinetic unified algorithm,GKUA),通过数值求解考虑转动自由度激发的Boltzmann-Rykov模型方程,得到了一种跨流域非定常流动数值模拟的方法.该求解方法以Boltzmann模型方程为控制方程,在常温状态下如果考虑转动能激发的情况...     

圆柱界面不稳定性的多组分气体动理学数值计算

在假定单元内各组分同温同速的条件下,采用气体动理学格式(Gas-Kinetic Scheme,GKS)对空气/He和空气/R22圆柱界面不稳定性进行了数值计算,得到不同时刻的密度分布以及界面上特征位置的位移历史和平均速度。当激波穿过界面后,界面上特征位置的位移随时间逐渐增大,特征位置的平均速度与前人的实验结果和数值模拟结果吻合很好。对比结果表明,从微观气体运动角度出发的GKS方法对于界面不稳定性问题具有良好的模拟能力。     

考虑转动能的一维/二维Boltzmann-Rykov模型方程数值算法

研究考虑转动能的Boltzmann-Rykov模型方程,基于转动自由度对气体分子速度分布函数矩积分,引入约化速度分布函数,应用离散速度坐标法与数值积分技术,将气体运动论模型方程化为在离散速度坐标点处关于三个约化速度分布函数的联立方程组.应用拓展计算流体力学有限差分方法,数值计算考虑转动自由度的双原子气体一维、二维Boltzmann模型方程,得到高、低Knudsen数一维激波管内流动和二维竖直平板绕流问题的流场,分析验证考虑转动能的Boltzmann-Rykov模型方程全流域统一算法求解一维/二维气体流动问题的可靠性.结果表明,气体稀薄程度与分子内自由度对流场具有较大影响,且Knudsen数较高的稀薄气体流动呈现严重的非平衡流动特点.     

频率啁啾对强场下真空正负电子对产生的增强效应

在强背景场下真空产生正反粒子对的研究中,频率啁啾对增强粒子对的产生起着关键作用.本文介绍了狄拉克-海森伯-维格纳(Dirac-Heisenberg-Wigner)、求解量子弗拉索夫方程(quantum Vlasov equation)和计算量子场论等方法,并详细综述了它们如何应用到空间非均匀场、均匀含时场以及外部势场中正负电子对产生的研究.通过研究各种不同的场得到了不同参数(如场强和基准频率)下产生的粒子动量谱和粒子对产额,发现当频率啁啾形式或/和啁啾强度改变时结果受到显著影响.在低频场下啁啾增强的数密度可提高2—3个数量级,这主要是因为啁啾增加了场的高频成分,从而低频强场和高频弱场相结合的动力学辅助机制起到了很大的作用.一般来说在高频情况下数密度只有几倍的提高,说明动力学辅助作用被大大地抑制了.在有空间变化的场情形下,对于小空间尺度变化的场,无啁啾时本身的数密度不高,但啁啾可以对数密度有数量级的提高;对于大空间尺度变化的场,数密度逐渐趋于空间均匀的结果,啁啾也能对数密度有几倍的提高.通过Wentzel-Kramer-Brillouin近似和转变点结构的物理分析和讨论可以对相关的数值结...     

反应流体的移动网格动理学格式

在移动网格上构造一种反应流的动理学格式.首先利用BGK模型推导含化学反应的流体力学方程组,并利用其积分形式构造移动网格上离散格式,再利用自适应移动网格方法得到网格速度,最后利用时间精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.一维与二维的数值实验表明这种格式同时具有高精度、高分辨率的特点.     

Entropy considerations in numerical simulations of non-equilibrium rarefied flows

Non-equilibrium rarefied flows are encountered frequently in supersonic flight at high altitudes, vacuum technology and in microscale devices. Prediction of the onset of non-equilibrium is important for accurate numerical simulation of such flows. We formulate and apply the discrete version of Boltzmann’s H-theorem for analysis of non-equilibrium onset and accuracy of numerical modeling of rarefied gas flows. The numerical modeling approach is based on the deterministic solution of kinetic model equations. The numerical solution approach comprises the discrete velocity method in the velocity space and the finite volume method in the physical space with different numerical flux schemes: the first-order, the second-order minmod flux limiter and a third-order WENO schemes. The use of entropy considerations in rarefied flow simulations is illustrated for the normal shock, the Riemann and the two-dimensional shock tube problems. The entropy generation rate based on kinetic theory is shown to be a powerful indicator of the onset of non-equilibrium, accuracy of numerical solution as well as the compatibility of boundary conditions for both steady and unsteady problems.     

气-固界面作用物理模型的确定性实现方法

针对最新发展的气-固界面作用物理模型,提出一种确定性计算方法.算法的核心是:针对离散速度空间分布函数,采用确定性算法计算出单轮气-固碰撞作用的散射核函数矩阵,并通过考虑吸引势阱作用和迭代累加多次气-固碰撞作用,物理地反映气体分子入射/反射速度分布函数在气-固界面上的变化.与现有基于Monte Carlo随机采样的实现方...     

A unified gas-kinetic scheme for continuum and rarefied flows

With discretized particle velocity space, a multiscale unified gas-kinetic scheme for entire Knudsen number flows is constructed based on the BGK model. The current scheme couples closely the update of macroscopic conservative variables with the update of microscopic gas distribution function within a time step. In comparison with many existing kinetic schemes for the Boltzmann equation, the current method has no difficulty to get accurate Navier–Stokes (NS) solutions in the continuum flow regime with a time step being much larger than the particle collision time. At the same time, the rarefied flow solution, even in the free molecule limit, can be captured accurately. The unified scheme is an extension of the gas-kinetic BGK-NS scheme from the continuum flow to the rarefied regime with the discretization of particle velocity space. The success of the method is due to the un-splitting treatment of the particle transport and collision in the evaluation of local solution of the gas distribution function. For these methods which use operator splitting technique to solve the transport and collision separately, it is usually required that the time step is less than the particle collision time. This constraint basically makes these methods useless in the continuum flow regime, especially in the high Reynolds number flow simulations. Theoretically, once the physical process of particle transport and collision is modeled statistically by the kinetic Boltzmann equation, the transport and collision become continuous operators in space and time, and their numerical discretization should be done consistently. Due to its multiscale nature of the unified scheme, in the update of macroscopic flow variables, the corresponding heat flux can be modified according to any realistic Prandtl number. Subsequently, this modification effects the equilibrium state in the next time level and the update of microscopic distribution function. Therefore, instead of modifying the collision term of the BGK model, such as ES-BGK and BGK–Shakhov, the unified scheme can achieve the same goal on the numerical level directly. Many numerical tests will be used to validate the unified method.     

尖化前缘的稀薄气体化学非平衡流动和气动加热相似律研究

以近空间尖前缘高超声速巡航飞行器的研制为背景, 作者在前一阶段采用模型理论分析方法, 陆续研究了沿微钝前缘驻点线的化学非平衡流动和气动加热相似律, 文章是上述研究的综合回顾和深化讨论.稀薄条件下, 驻点附近流动和传热出现一系列与连续流动模型不同的新特征, 超出了经典气动热预测理论的适用范围.作者建立了一个沿驻点线能量传递和转化的广义模型, 并分别推导了具有实际物理意义的边界层外离解非平衡流动判据和边界层内复合非平衡流动判据.基于这些判据构建了预测非平衡流动驻点气动加热的桥函数, 并讨论了稀薄非平衡真实气体流动和气动加热的相似律, 发现新型近空间尖前缘飞行器遭遇的气动热环境不同于传统大钝头航天器再入问题, 传统的天地换算相似准则将会失效.这些理论分析结果可为稀薄非平衡化学反应流及气动加热的实验和计算提供一个标模检验的手段.      

A hybrid particle approach for continuum and rarefied flow simulation

A hybrid particle scheme is presented for the simulation of compressible gas flows involving both continuum regions and rarefied regions with strong translational nonequilibrium. The direct simulation Monte Carlo (DSMC) method is applied in rarefied regions, while remaining portions of the flowfield are simulated using a DSMC-based low diffusion particle method for inviscid flow simulation. The hybrid scheme is suitable for either steady state or unsteady flow problems, and can simulate gas mixtures comprising an arbitrary number of species. Numerical procedures are described for strongly coupled two-way information transfer between continuum and rarefied regions, and additional procedures are outlined for the determination of continuum breakdown. The hybrid scheme is evaluated through a comparison with DSMC simulation results for a Mach 6 flow of N₂ over a cylinder, and good overall agreement is observed. Large potential efficiency gains (over three orders of magnitude) are estimated for the hybrid algorithm relative to DSMC in a simple example involving a rarefied expansion flow through a small nozzle into a vacuum chamber.     

A unified gas kinetic scheme with moving mesh and velocity space adaptation

There is great difficulty for direct Boltzmann solvers to simulate high Knudsen number flow due to the severe steep slope and high concentration of the gas distribution function in a local particle velocity space. Local mesh adaptation becomes necessary in order to make the Boltzmann solver to be a practical tool in aerospace applications. The present research improves the unified gas-kinetic scheme (UGKS) in the following two aspects. First, the UGKS is extended in a physical space with moving mesh. This technique is important to study a freely flying object in a rarefied environment. Second, the adaptive quadtree method in the particle velocity space is implemented in the UGKS. Due to the new improvements in the discretization of a gas distribution function in the six dimensional phase space, the adaptive unified gas kinetic scheme (AUGKS) is able to deal with a wide range of flow problems under extreme flying conditions, such as the whole unsteady flying process of an object from a highly rarefied to a continuum flow regime. After validating the scheme, the capability of AUGKS is demonstrated in the following two challenge test cases. The first case is about the free movement of an ellipse flying at initial Mach number 5 in a rarefied flow at different Knudsen numbers. The force on the ellipse and the unsteady trajectory of the ellipse movement are fully captured. The gas distribution function around the ellipse is analyzed. The second case is about the study of unsteady flight of a nozzle under a bursting process of the compressed gas expanding into a rarefied environment. Due to the strong expansion wave and the huge density difference between interior and exterior regions around the nozzle, the particle distribution function changes dramatically in the particle velocity space. The use of an adaptive velocity space in the AUGKS becomes necessary to simulate such a flow and to control the computational cost to a tolerable level. The second test is a challenge problem for any existing rarefied flow solver.     

Accuracy analysis of high-order lattice Boltzmann models for rarefied gas flows

In this work, we have theoretically analyzed and numerically evaluated the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows. In the incompressible limit, the LB equation is shown to be able to reduce to the linearized Bhatnagar–Gross–Krook (BGK) equation. Therefore, when the same Gauss–Hermite quadrature is used, LB method closely resembles the discrete velocity method (DVM). In addition, the order of Hermite expansion for the equilibrium distribution function is found not to be directly correlated with the approximation order in terms of the Knudsen number to the BGK equation for incompressible flows. Meanwhile, we have numerically evaluated the LB models for a standing-shear-wave problem, which is designed specifically for assessing model accuracy by excluding the influence of gas molecule/surface interactions at wall boundaries. The numerical simulation results confirm that the high-order terms in the discrete equilibrium distribution function play a negligible role in capturing non-equilibrium effect for low-speed flows. By contrast, appropriate Gauss–Hermite quadrature has the most significant effect on whether LB models can describe the essential flow physics of rarefied gas accurately. Our simulation results, where the effect of wall/gas interactions is excluded, can lead to conclusion on the LB modeling capability that the models with higher-order quadratures provide more accurate results. For the same order Gauss–Hermite quadrature, the exact abscissae will also modestly influence numerical accuracy. Using the same Gauss–Hermite quadrature, the numerical results of both LB and DVM methods are in excellent agreement for flows across a broad range of the Knudsen numbers, which confirms that the LB simulation is similar to the DVM process. Therefore, LB method can offer flexible models suitable for simulating continuum flows at the Navier–Stokes level and rarefied gas flows at the linearized Boltzmann model equation level.     

微管道气体流动的蒙特卡洛直接模拟

采用蒙特卡洛直接模拟(DSMC)方法,数值模拟了微管道中压力驱动的气体流动,结果表明固壁边界存在速度滑移,稀薄气体效应明显;整个流场温度变化很小,流动马赫数很小,密度、压力流向变化非常大而横向几乎不变;可压缩性导致压力随流向的非线性分布,但这种效应随Knudsen数增大而减弱.