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.

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.

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

High Order Hierarchical Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes For The BGK Equation

A Discontinuous Galerkin Method Based On A Bgk Scheme For The Navier-Stokes Equations On Arbitrary Grids

A BGK-Based Discontinuous Galerkin Method for the Navier-Stokes Equations on Arbitrary Grids

High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling

Asymptotic Analysis of IMEX-RK Methods for ES-BGK Model at Navier-Stokes level

Time-Implicit High-Order Accurate Positivity-Preserving Discretizations for the Navier–Stokes and Navier–Stokes–Korteweg Equations

A high-order gas-kinetic Navier–Stokes flow solver

An Asymptotic Preserving Scheme for the ES-BGK model

A New Gas-Kinetic Scheme Based on Analytical Solutions of the BGK Equation.

A robust implicit high-order discontinuous Galerkin method for solving compressible Navier-Stokes equations on arbitrary grids

A High-order Gas-Kinetic Navier-Stokes Solver I: One-dimensional Flux Evaluation

Comparison of the high-order Runge-Kutta discontinuous Galerkin method and gas-kinetic scheme for inviscid compressible flow simulations

A DGBGK scheme based on WENO limiters for viscous and inviscid flows

Stability and Consistency of Kinetic Upwinding for Advection–diffusion Equations

On high-order/low-order and micro-macro methods for implicit time-stepping of the BGK model

A New High-Order Multidimensional Scheme

A Gas-Kinetic BGK Scheme for the Navier–Stokes Equations and Its Connection with Artificial Dissipation and Godunov Method

Discontinuous Galerkin BGK Method for Viscous Flow Equations: One-Dimensional Systems

High-order discontinuous Galerkin schemes with subcell finite volume limiter and adaptive mesh refinement for a monolithic first-order BSSNOK formulation of the Einstein-Euler equations

A comparative study of an asymptotic preserving scheme and unified gas-kinetic scheme in continuum flow limit.