High Order Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes for the BGK Equation

Tao Xiong,Juhi Jang,Fengyan Li,Jing-Mei Qiu
DOI: https://doi.org/10.1016/j.jcp.2014.12.021
2014-06-18
Abstract:In this paper, we develop high-order asymptotic preserving (AP) schemes for the BGK equation in a hyperbolic scaling, which leads to the macroscopic models such as the Euler and compressible Navier-Stokes equations in the asymptotic limit. Our approaches are based on the so-called micro-macro formulation of the kinetic equation which involves a natural decomposition of the problem to the equilibrium and the non-equilibrium parts. The proposed methods are formulated for the BGK equation with constant or spatially variant Knudsen number. The new ingredients for the proposed methods to achieve high order accuracy are the following: we introduce discontinuous Galerkin (DG) discretization of arbitrary order of accuracy with nodal Lagrangian basis functions in space; we employ a high order globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) scheme as time discretization. Two versions of the schemes are proposed: Scheme I is a direct formulation based on the micro-macro decomposition of the BGK equation, while Scheme II, motivated by the asymptotic analysis for the continuous problem, utilizes certain properties of the projection operator. Compared with Scheme I, Scheme II not only has better computational efficiency (the computational cost is reduced by half roughly), but also allows the establishment of a formal asymptotic analysis. Specifically, it is demonstrated that when $0<\epsilon \ll 1$, Scheme II, up to $\mathcal{O}(\epsilon^2)$, becomes a local DG discretization with an explicit RK method for the macroscopic compressible Navier-Stokes equations, a method in a similar spirit to the ones in [Bassi \& Rabey 1997, Cockburn \& Shu 1998]. Numerical results are presented for a wide range of Knudsen number to illustrate the effectiveness and high order accuracy of the methods.
Numerical Analysis
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to develop high - order Asymptotic Preserving (AP) numerical methods for solving the BGK equation. The BGK equation is a simplified model of the Boltzmann equation that describes gas dynamics and can be reduced to macroscopic models, such as the Euler equation and the compressible Navier - Stokes equation, under different Knudsen numbers ($\varepsilon$). ### Specific Problem Description 1. **Model Conversion under Different Knudsen Numbers**: - When the Knudsen number $\varepsilon \to 0$, the BGK equation should be reduced to the Euler equation that describes mass, momentum, and energy conservation. - When the Knudsen number is small but not zero, the BGK equation should be reduced to the compressible Navier - Stokes equation that contains viscosity and heat conduction correction terms. 2. **Efficiency and Accuracy of Numerical Methods**: - The developed numerical method needs to maintain high - order accuracy in a wide range of Knudsen numbers and be able to correctly capture the limiting behavior of the Euler equation when $\varepsilon \to 0$. - At the same time, when the Knudsen number is small but not zero, the numerical method should be able to effectively approximate the compressible Navier - Stokes equation. ### Main Contributions of the Paper - **Micro - Macro Decomposition Framework**: By introducing micro - macro decomposition, the BGK equation is decomposed into an equilibrium part and a non - equilibrium part, thus dealing with problems at different scales more effectively. - **High - Order Discontinuous Galerkin (DG) Discretization**: The nodal - type discontinuous Galerkin method with arbitrary - order accuracy is used for spatial discretization, ensuring high - order accuracy in space. - **Implicit - Explicit Runge - Kutta (IMEX - RK) Time Discretization**: The globally stiffly - accurate IMEX - RK method is adopted for time discretization, ensuring high - order accuracy in time. - **Proposal of Two Schemes**: - **Scheme I**: Based on the direct micro - macro decomposition of the BGK equation. - **Scheme II**: Utilize the properties of projection operators to further optimize computational efficiency and establish a formal asymptotic analysis. ### Advantages of Scheme II - **Higher Computational Efficiency**: Compared with Scheme I, the computational cost of Scheme II is approximately reduced by half. - **More Complete Asymptotic Analysis**: Scheme II not only is reduced to the local DG discretization of the Euler equation when $\varepsilon \to 0$, but also can approximate the compressible Navier - Stokes equation with an accuracy of $O(\varepsilon^2)$. ### Numerical Verification Through a series of numerical experiments, the effectiveness and high - order accuracy of the proposed high - order AP scheme under different Knudsen numbers are verified, especially its superior performance in smooth and non - smooth test cases. In summary, this paper aims to develop a high - order Asymptotic Preserving numerical method applicable to a wide range of Knudsen numbers to effectively solve the BGK equation and ensure that the corresponding macroscopic models can be correctly captured under different limit conditions.