Unified gas-kinetic wave-particle methods V: Diatomic molecular flow

JOURNAL OF COMPUTATIONAL PHYSICS

卷:442

Article

Xu, Xiaocong¹  Chen, Yipei¹  Liu, Chang¹  Li, Zhihui²  Xu, Kun^1,3 

[1] Hong Kong Univ Sci & Technol, Dept Math, Hong Kong, Peoples R China

[2] Beijing Univ Aeronaut & Astronaut, Natl Lab Computat Fluid Dynam, Beijing 100191, Peoples R China

[3] Hong Kong Univ Sci & Technol, Shenzhen Res Inst, Shenzhen, Peoples R China

关键词: Unified gas-kinetic wave-particle method;    Multiscale transport;    Diatomic gas;    Hypersonic non-equilibrium flow;   

DOI  :  10.1016/j.jcp.2021.110496

来源: Elsevier

【 摘 要 】

In this paper, the unified gas-kinetic wave-particle (UGKWP) method is further developed for diatomic gas with the energy exchange between translational and rotational modes for flow study in all regimes. The multiscale transport mechanism in UGKWP is coming from the direct modellingin a discretized space, where the cell's Knudsen number, defined by the ratio of particle mean free path over the numerical cell size, determines the flow physics simulated by the wave particle method. The non-equilibrium distribution function in UGKWP is tracked by the discrete particle and analytical wave. The weights of distributed particle and wave in different regimes depend on cell's Knudsen number, where distinguishable macroscopic flow variables of particle and wave are updated inside each control volume. The UGKWP becomes a particle method in the highly rarefied flow regime and converges to the gas-kinetic scheme (GKS) for the Navier-Stokes solution in the continuum flow regime without particles. In comparison with discrete velocity method (DVM)-based unified gas-kinetic scheme (UGKS), the computational cost and memory requirement in UGKWP have been reduced by several orders of magnitude for the high speed and high temperature flow simulation, where the translational and rotational non-equilibrium can be captured accurately in the transition and rarefied regime. As a result, for the hypersonic flow around a flying vehicle, the computation can be conducted using a personal computer for the studies in all regimes. The UGKWP method for diatomic gas will be validated in various cases from one dimensional shock structure to three dimensional flow over a sphere, and the numerical solutions will be compared with the reference DSMC results and experimental measurements. (C) 2021 Elsevier Inc. All rights reserved.

【 授权许可】

Free