Physical Modeling and Numerical Studies of Three-dimensional Non-equilibrium Multi-temperature Flows

TL;DR

This paper develops a high-order three-dimensional multi-temperature gas-kinetic scheme for accurately and efficiently simulating non-equilibrium flows, validated through various complex flow scenarios and comparisons with experimental and other numerical methods.

Contribution

It introduces a novel high-order 3D multi-temperature GKS within a two-stage fourth-order framework, improving accuracy and efficiency over existing methods.

Findings

High-order accuracy and robustness confirmed through validation.

Significant efficiency improvements over UGKS and DSMC.

Effective simulation of complex non-equilibrium flow phenomena.

Abstract

For increasingly rarefied flowfields, the Navier-Stokes (NS) equations lose accuracy partially due to the single temperature approximation. To overcome this barrier, a continuum multi-temperature model based on the Bhatnagar-Gross-Krook (BGK) equation coupled with the Landau-Teller-Jeans relaxation model has been proposed for two-dimensional hypersonic non-equilibrium multi-temperature flow computation. In recent study, a two-stage fourth-order gas-kinetic scheme (GKS) has been developed for equilibrium flows, which achieves a fourth-order accuracy in space and time as well as high efficiency and robustness. In this paper, targeting for accurate and efficient simulation of multi-temperature non-equilibrium flows, a high-order three-dimensional multi-temperature GKS is implemented under the two-stage fourth-order framework, with the fourth-order Simpson interpolation rule for the newly emerged source term. Simulations on decaying homogeneous isotropic turbulence, low-density nozzle flow, rarefied hypersonic flow over a flat plate, and type IV shock-shock interaction are used to validate the multi-temperature model through the comparison with experimental measurements. The unified gas kinetic scheme (UGKS) results, and the Direct simulation Monte Carlo (DSMC) solutions will be used as well in some cases for validation. Computational results not only confirm the high-order accuracy and quite robustness of this scheme, but also show the significant improvement on computational efficiency compared with UGKS and DSMC, especially in the near continuum flow regime.