General synthetic iterative scheme for nonlinear gas kinetic simulation of multi-scale rarefied gas flows

Zhu, Lianhua and Pi, Xingcai and Su, Wei and Li, Zhi-Hui and Zhang, Yonghao and Wu, Lei (2021) General synthetic iterative scheme for nonlinear gas kinetic simulation of multi-scale rarefied gas flows. Journal of Computational Physics, 430. 110091. ISSN 0021-9991

[thumbnail of Zhu-etal-JCP-2020-General-synthetic-iterative-scheme-for-nonlinear-gas-kinetic-simulation] Text (Zhu-etal-JCP-2020-General-synthetic-iterative-scheme-for-nonlinear-gas-kinetic-simulation)
Zhu_etal_JCP_2020_General_synthetic_iterative_scheme_for_nonlinear_gas_kinetic_simulation.pdf
Accepted Author Manuscript
Restricted to Repository staff only until 29 December 2021.
License: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 logo

Download (4MB) | Request a copy from the Strathclyde author

    Abstract

    The general synthetic iterative scheme (GSIS) is extended to find the steady-state solution of the nonlinear gas kinetic equation, resolving the long-standing problems of slow convergence and requirement of ultra-fine grids in near-continuum flows. The key ingredient of GSIS is the tight coupling of gas kinetic and macroscopic synthetic equations, where the constitutive relations explicitly contain Newton's law of shear stress and Fourier's law of heat conduction. The higher-order constitutive relations describing rarefaction effects are calculated from the velocity distribution function; however, their constructions are simpler than our previous work (Su et al., 2020 [28]) for linearized gas kinetic equations. On the other hand, solutions of macroscopic synthetic equations are used to accelerate the evolution of gas kinetic equation at the next iteration step. A rigorous linear Fourier stability analysis of the present schemes in periodic system shows that the error decay rate of GSIS can be smaller than 0.5, which means that the deviation to steady-state solution can be reduced by 3 orders of magnitude in 10 iterations. Other important advantages of the GSIS are: (i) it does not rely on the specific form of Boltzmann collision operator, and (ii) it can be solved by sophisticated techniques in computational fluid dynamics, making it amenable to large scale engineering applications. In this paper, the efficiency and accuracy of GSIS are demonstrated by a number of canonical test cases in rarefied gas dynamics, covering different flow regimes.

    ORCID iDs

    Zhu, Lianhua ORCID logoORCID: https://orcid.org/0000-0003-1615-7371, Pi, Xingcai, Su, Wei, Li, Zhi-Hui, Zhang, Yonghao and Wu, Lei;