High order asymptotic preserving nodal discontinuous Galerkin IMEX schemes for the BGK equation

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 < e ? 1 , Scheme II, up to O ( e 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 and Rebay (1997) 3, Cockburn and Shu (1998) 16. Numerical results are presented for a wide range of Knudsen number to illustrate the effectiveness and high order accuracy of the methods.

[1]  Lorenzo Pareschi,et al.  Implicit-explicit runge-kutta schemes and applications to hyperbolic systems with relaxation , 2010, 1009.2757.

[2]  Carlo Cercignani,et al.  Mathematical Methods in Kinetic Theory , 1970 .

[3]  Gabriella Puppo,et al.  Microscopically implicit-macroscopically explicit schemes for the BGK equation , 2012, J. Comput. Phys..

[4]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[5]  T. G. Cowling,et al.  The mathematical theory of non-uniform gases : an account of the kinetic theory of viscosity, thermal conduction, and diffusion in gases , 1954 .

[6]  Gabriella Puppo,et al.  Implicit–Explicit Schemes for BGK Kinetic Equations , 2007, J. Sci. Comput..

[7]  Tao Xiong,et al.  Analysis of Asymptotic Preserving DG-IMEX Schemes for Linear Kinetic Transport Equations in a Diffusive Scaling , 2013, SIAM J. Numer. Anal..

[8]  Shi Jin,et al.  A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources , 2009, J. Comput. Phys..

[9]  George Em Karniadakis,et al.  A discontinuous Galerkin method for the Navier-Stokes equations , 1999 .

[10]  J. Tinsley Oden,et al.  A discontinuous hp finite element method for the Euler and Navier–Stokes equations , 1999 .

[11]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[12]  S. Rebay,et al.  Numerical evaluation of two discontinuous Galerkin methods for the compressible Navier–Stokes equations , 2002 .

[13]  Shi Jin,et al.  An Asymptotic Preserving Scheme for the ES-BGK Model of the Boltzmann Equation , 2010, J. Sci. Comput..

[14]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[15]  Tai-Ping Liu,et al.  Boltzmann Equation: Micro-Macro Decompositions and Positivity of Shock Profiles , 2004 .

[16]  C. Cercignani.  Rarefied Gas Dynamics: From Basic Concepts to Actual Calculations , 2000 .

[17]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[18]  SEBASTIANO BOSCARINO.  Error Analysis of IMEX Runge-Kutta Methods Derived from Differential-Algebraic Systems , 2007, SIAM J. Numer. Anal..

[19]  Chi-Wang Shu,et al.  High order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD , 2001 .

[20]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[21]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[22]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[23]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[24]  Shi Jin.  ASYMPTOTIC PRESERVING (AP) SCHEMES FOR MULTISCALE KINETIC AND HYPERBOLIC EQUATIONS: A REVIEW , 2010 .

[25]  Kun Xu,et al.  A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method , 2001 .

[26]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[27]  Tao Xiong,et al.  High order asymptotic preserving DG-IMEX schemes for discrete-velocity kinetic equations in a diffusive scaling , 2013, J. Comput. Phys..

[28]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[29]  F. Golse,et al.  Fluid dynamic limits of kinetic equations. I. Formal derivations , 1991 .

[30]  Serge Marguet,et al.  The Boltzmann Equation , 2017 .

[31]  L. Mieussens.  Discrete-Velocity Models and Numerical Schemes for the Boltzmann-BGK Equation in Plane and Axisymmetric Geometries , 2000 .

[32]  Carlo Cercignani.  The Boltzmann Equation , 1988 .

[33]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[34]  Luc Mieussens,et al.  Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier-Stokes asymptotics , 2008, J. Comput. Phys..