论文标题

在所有流程度中,具有热力学非平衡流量的简化多尺度数值通量的统一气体运动方案

Unified gas-kinetic scheme with simplified multi-scale numerical flux for thermodynamic non-equilibrium flow in all flow regimes

论文作者

Zhang, Rui, Liu, Sha, Zhong, Chengwen, Zhuo, Congshan

论文摘要

在本文中,提出了一种具有简化多尺度数值通量的统一气动方案(UGK),用于用于热力学非平衡流量模拟,涉及所有流动方案中分子振动自由度的激发。目前的UGK将宏观流量变量的基本保护定律和离散空间中的微观气体分布函数保留。为了提高UGK的效率,简化的多尺度数值通量直接从动力学模型方程的特征差解构建。此外,提出了一种新的BGK型动力学模型来描述高温热力学非平衡效应,这是一种现象学弛豫模型,具有旋转和振动能的连续分布模式。在当前模型中,通过在麦克斯韦分布周围使用多维的遗产扩展来构建平衡分布函数,以实现正确的prandtl数量和适当的热通量放松率。此外,非结构化离散速度空间(DVS)和简单集成误差校正的应用大大降低了速度网格的数量,并使当前的方法是所有流程度中流量模拟的有效工具。在一系列情况下,检查了新方案,例如SOD的冲击管,高度平衡的冲击结构,围绕圆形圆柱体的高超音速流动,带有Knudsen(KN)数字kN = 0.01,以及在平板上的稀有超声流动,具有锋利的前沿。目前的UGK结果与DSMC的基准数据和其他经过验证的方法非常吻合。

In this paper, a unified gas-kinetic scheme (UGKS) with simplified multi-scale numerical flux is proposed for the thermodynamic non-equilibrium flow simulation involving the excitation of molecular vibrational degrees of freedom in all flow regimes. The present UGKS keep the basic conservation laws of the macroscopic flow variables and the microscopic gas distribution function in a discretized space. In order to improve the efficiency of the UGKS, a simplify multi-scale numerical flux is directly constructed from the characteristic difference solution of the kinetic model equation. In addition, a new BGK-type kinetic model for diatomic gases is proposed to describe the high-temperature thermodynamic non-equilibrium effect, which is a phenomenological relaxation model with the continuous distribution modes of rotational and vibrational energies. In present model, the equilibrium distribution functions is constructed by using a multi-dimensional Hermitian expansion around the Maxwellian distribution to achieve the correct Prandtl number and proper relaxation rate of heat fluxes. Furthermore, the application of the unstructured discrete velocity space (DVS) and a simple integration error correction reduce the number of velocity mesh significantly and make the present method be a efficient tool for simulations of flows in all flow regimes. The new scheme are examined in a series of cases, such as Sod's shock tube, high non-equilibrium shock structure, hypersonic flow around a circular cylinder with Knudsen (Kn) number Kn = 0.01, and the rarefied hypersonic flow over a flat plate with a sharp leading edge. The present UGKS results agree well with the benchmark data of DSMC and the other validated methods.

扫码加入交流群

加入微信交流群

微信交流群二维码

扫码加入学术交流群,获取更多资源