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

by   Zhichao Peng, et al.

We consider a linear kinetic transport equation under a diffusive scaling, that converges to a diffusion equation as the Knudsen number ε→0. In [3, 21], to achieve the asymptotic preserving (AP) property and unconditional stability in the diffusive regime with ε≪ 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 ε and the unconditional stability in the diffusive regime. A series of numerical examples are presented to demonstrate the performance of the new schemes.



There are no comments yet.


page 14


Stability-enhanced AP IMEX1-LDG method: energy-based stability and rigorous AP property

In our recent work [22], a family of high order asymptotic preserving (A...

High order asymptotic preserving discontinuous Galerkin methods for gray radiative transfer equations

In this paper, we will develop a class of high order asymptotic preservi...

Conservative DG Method for the Micro-Macro Decomposition of the Vlasov-Poisson-Lenard-Bernstein Model

The micro-macro (mM) decomposition approach is considered for the numeri...

Accuracy and stability analysis of the Semi-Lagrangian method for stiff hyperbolic relaxation systems and kinetic BGK model

In this paper, we develop a family of third order asymptotic-preserving ...

Uniformly accurate numerical schemes for a class of dissipative systems

We consider a class of relaxation problems mixing slow and fast variatio...

Asymptotic preserving schemes for the FitzHugh-Nagumo transport equation with strong local interactions

This paper is devoted to the numerical approximation of the spatially-ex...

Block preconditioning methods for asymptotic preserving scheme arising in anisotropic elliptic problems

Efficient and robust iterative solvers for strong anisotropic elliptic e...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.