An asymptotic-preserving dynamical low-rank method for the multi-scale multi-dimensional linear transport equation

被引:36
作者
Einkemmer, Lukas [1 ]
Hu, Jingwei [2 ]
Wang, Yubo [2 ]
机构
[1] Univ Innsbruck, Dept Math, A-6020 Innsbruck, Austria
[2] Purdue Univ, Dept Math, W Lafayette, IN 47907 USA
关键词
Dynamical low-rank integrator; Linear transport equation; Macro-micro decomposition; Diffusion limit; Asymptotic preserving; Implicit-explicit Runge-Kutta scheme (IMEX); PROJECTOR-SPLITTING INTEGRATOR; KINETIC-EQUATIONS; TIME INTEGRATION; AP SCHEMES; ALGORITHM;
D O I
10.1016/j.jcp.2021.110353
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We propose a dynamical low-rank method to reduce the computational complexity for solving the multi-scale multi-dimensional linear transport equation. The method is based on a macro-micro decomposition of the equation; the low-rank approximation is only used for the micro part of the solution. The time and spatial discretizations are done properly so that the overall scheme is second-order accurate (in both the fully kinetic and the limit regime) and asymptotic-preserving (AP). That is, in the diffusive regime, the scheme becomes a macroscopic solver for the limiting diffusion equation that automatically captures the low-rank structure of the solution. Moreover, the method can be implemented in a fully explicit way and is thus significantly more efficient compared to the previous state of the art. We demonstrate the accuracy and efficiency of the proposed low-rank method by a number of four-dimensional (two dimensions in physical space and two dimensions in velocity space) simulations. (C) 2021 The Author(s). Published by Elsevier Inc.
引用
收藏
页数:21
相关论文
共 27 条
[1]   IMPLICIT-EXPLICIT RUNGE-KUTTA SCHEMES FOR HYPERBOLIC SYSTEMS AND KINETIC EQUATIONS IN THE DIFFUSION LIMIT [J].
Boscarino, S. ;
Pareschi, L. ;
Russo, G. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2013, 35 (01) :A22-A51
[2]  
Chandrasekhar S., 1960, Radiative Transfer
[3]  
DAVISON B, 1973, NEUTRON TRANSPORT TH
[4]   Asymptotic-Preserving methods and multiscale models for plasma physics [J].
Degond, Pierre ;
Deluzet, Fabrice .
JOURNAL OF COMPUTATIONAL PHYSICS, 2017, 336 :429-457
[5]  
Ding Z., 2019, ARXIV190704247
[6]   A low-rank projector-splitting integrator for the Vlasov-Maxwell equations with divergence correction [J].
Einkemmer, Lukas ;
Ostermann, Alexander ;
Piazzola, Chiara .
JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 403
[7]   A LOW-RANK ALGORITHM FOR WEAKLY COMPRESSIBLE FLOW [J].
Einkemmer, Lukas .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2019, 41 (05) :A2795-A2814
[8]   A QUASI-CONSERVATIVE DYNAMICAL LOW-RANK ALGORITHM FOR THE VLASOV EQUATION [J].
Einkemmer, Lukas ;
Lubich, Christian .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2019, 41 (05) :B1061-B1081
[9]   A LOW-RANK PROJECTOR-SPLITTING INTEGRATOR FOR THE VLASOV-POISSON EQUATION [J].
Einkemmer, Lukas ;
Lubich, Christian .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2018, 40 (05) :B1330-B1360
[10]   An almost symmetric Strang splitting scheme for nonlinear evolution equations [J].
Einkemmer, Lukas ;
Ostermann, Alexander .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2014, 67 (12) :2144-2157