Abstract:A class of high order asymptotic preserving (AP) schemes has been developed for the BGK equation in Xiong et. al. (2015) [37], which is based on the micro-macro formulation of the equation. The nodal discontinuous Galerkin (NDG) method with Lagrangian basis functions for spatial discretization and globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) scheme as time discretization are introduced with asymptotic preserving properties. However, it is only necessary to solve the kinetic equation when the hydrodynamic description breaks down. Motivated by the recent work in Filbet and Rey (2015) [23], it is more naturally to construct a hierarchy scheme under the NDG-IMEX framework without hybridization, as the formal analysis in [37] shows that when $\epsilon$ is small, the NDG-IMEX scheme becomes a local discontinuous Galerkin (LDG) scheme for the compressible Navier-Stokes equations, and when $\epsilon=0$ it is a discontinuous Galerkin (DG) scheme for the compressible Euler equations. Moveover, we propose to combine the kinetic regime with the hydrodynamic regime including both the compressible Euler and Navier-Stokes equations. Numerical experiments demonstrate very decent performance of the new approach. In our numerics, all three regimes are clearly divided, leading to great savings in terms of the computational cost.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an efficient numerical method for solving the BGK (Bhatnagar - Gross - Krook) equation, especially being able to automatically adapt and switch between different physical models (from kinetic theory to fluid mechanics) at different scales, so as to reduce the computational cost and improve the simulation accuracy. Specifically:
1. **Treatment of multi - scale problems**: When the Knudsen number (\(\varepsilon\)) varies within a large range, how to effectively simulate the behavior of rarefied gases is a challenge. When \(\varepsilon\) approaches zero, the system should converge to a fluid - mechanics description (such as the compressible Euler or Navier - Stokes equations), while when \(\varepsilon\) is large, a kinetic model needs to be used.
2. **Requirement for efficient computation**: Directly using the kinetic model for simulation is very time - consuming, especially in the regions where the fluid - mechanics description is applicable. Therefore, researchers hope to find a method that can use the fluid - mechanics model as much as possible to reduce the computational cost while ensuring the accuracy.
3. **Adaptive switching mechanism**: To achieve this goal, the paper proposes a hierarchy scheme, which can adaptively switch solvers between different physical regions. This scheme is based on the Nodal Discontinuous Galerkin (NDG) and Implicit - Explicit (IMEX) Runge - Kutta methods, ensuring high - order accuracy and asymptotic - preserving properties (Asymptotic Preserving, AP).
4. **Specific scientific problems**: The paper mentions that for multi - scale problems in rarefied gas dynamics, traditional single models are difficult to take into account physical phenomena at different scales. For example, near shock waves or in the boundary - layer region, the fluid - mechanics description may fail, and the kinetic description is too expensive. Therefore, this research aims to develop a method that can switch seamlessly between these two descriptions, thereby simulating these complex phenomena more accurately and efficiently.
In summary, the core problem of this paper is to develop an efficient and accurate numerical method that can automatically select appropriate physical models at different scales and ensure a smooth transition from kinetic theory to fluid mechanics through asymptotic - preserving properties, thereby significantly reducing the computational cost and improving the accuracy of the simulation.