Abstract:We consider a linear kinetic transport equation under a diffusive scaling, that converges to a diffusion equation as the Knudsen number $\varepsilon\rightarrow0$. In [3, 21], to achieve the asymptotic preserving (AP) property and unconditional stability in the diffusive regime with $\varepsilon\ll 1$, numerical schemes are developed based on an additional reformulation of the even-odd or micro-macro decomposed version of the equation. The key of the reformulation is to add a weighted diffusive term on both sides of one equation in the decomposed system. The choice of the weight function, however, is problem-dependent and ad-hoc, and it can affect the performance of numerical simulations. To avoid issues related to the choice of the weight function and still obtain the AP property and unconditional stability in the diffusive regime, we propose in this paper a new family of AP schemes, termed as IMEX-DG-S schemes, directly solving the micro-macro decomposed system without any further reformulation. The main ingredients of the IMEX-DG-S schemes include globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) temporal discretizations with a new IMEX strategy, discontinuous Galerkin (DG) spatial discretizations, discrete ordinate methods for the velocity space, and the application of the Schur complement to the algebraic form of the schemes to control the overall computational cost. The AP property of the schemes is shown formally. With an energy type stability analysis applied to the first order scheme, and Fourier type stability analysis applied to the first to third order schemes, we confirm the uniform stability of the methods with respect to $\varepsilon$ and the unconditional stability in the diffusive regime. A series of numerical examples are presented to demonstrate the performance of the new schemes.

Energy diminishing implicit-explicit Runge–Kutta methods for gradient flows

Energy diminishing implicit-explicit Runge--Kutta methods for gradient flows

Higher-Order Energy-Decreasing Exponential Time Differencing Runge-Kutta methods for Gradient Flows

Average energy dissipation rates of additive implicit-explicit Runge-Kutta methods for gradient flow problems

Average energy dissipation rates of explicit exponential Runge-Kutta methods for gradient flow problems

High-order energy stable discrete variational derivative schemes for gradient flows

A novel high-order linearly implicit and energy-stable additive Runge-Kutta methods for gradient flow models

Original energy dissipation preserving corrections of integrating factor Runge-Kutta methods for gradient flow problems

Energy-dissipative spectral renormalization exponential integrator method for gradient flow problems

Two novel numerical methods for gradient flows: generalizations of the Invariant Energy Quadratization method

Variational Extrapolation of Implicit Schemes for General Gradient Flows

A Second-Order, Global-in-Time Energy Stable Implicit-Explicit Runge–Kutta Scheme for the Phase Field Crystal Equation

An Implicit-Explicit Relaxation Extrapolated Runge-Kutta and Energy-Preserving Finite Element Method for Klein-Gordon-Schrodinger Equations

Explicit Exponential Runge–Kutta Methods for Semilinear Time-Fractional Integro-Differential Equations

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

A Unified Design of Energy Stable Schemes with Variable Steps for Fractional Gradient Flows and Nonlinear Integro-differential Equations

High order strong stability preserving multi-derivative implicit and IMEX Runge--Kutta methods with asymptotic preserving properties

Stability analysis and error estimates of local discontinuous Galerkin methods with implicit–explicit time-marching for nonlinear convection–diffusion problems

Asymptotic preserving IMEX-DG-S schemes for linear kinetic transport equations based on Schur complement

A unified framework on the original energy laws of three effective classes of Runge-Kutta methods for phase field crystal type models