TY - JOUR
T1 - An asymptotic-preserving scheme for the kinetic equation with anisotropic scattering
T2 - Heavy tail equilibrium and degenerate collision frequency
AU - Wang, Li
AU - Yan, Bokai
N1 - Publisher Copyright:
©2019 Society for Industrial and Applied Mathematics.
PY - 2019
Y1 - 2019
N2 - We design an asymptotic preserving (AP) scheme for the linear kinetic equation with anisotropic scattering that leads to a fractional diffusion limit. This limit may be attributed to two reasons: A heavy tail equilibrium or a degenerate collision frequency, both of which are considered in this paper. Our scheme builds on the ideas developed in [L. Wang and B. Yan, J. Comput. Phys., 312 (2016), pp. 157-174] but with two major variations. One is a new splitting of the system that accounts for the anisotropy in the scattering cross section by introducing two extra terms. We then showed, via detailed calculation, that the scheme enjoys a relaxed AP property as opposed to the one step AP for the isotropic scattering. Another contribution is for the degenerate collision frequency case, which brings in additional stiffness. We propose to integrate a "body" term, which appears to be the main component in the diffusion limit. This term is precomputed once with a prescribed accuracy, via a change of variable that alleviates the stiffness. Numerical examples are presented to validate its efficiency in both kinetic and fractional diffusion regimes.
AB - We design an asymptotic preserving (AP) scheme for the linear kinetic equation with anisotropic scattering that leads to a fractional diffusion limit. This limit may be attributed to two reasons: A heavy tail equilibrium or a degenerate collision frequency, both of which are considered in this paper. Our scheme builds on the ideas developed in [L. Wang and B. Yan, J. Comput. Phys., 312 (2016), pp. 157-174] but with two major variations. One is a new splitting of the system that accounts for the anisotropy in the scattering cross section by introducing two extra terms. We then showed, via detailed calculation, that the scheme enjoys a relaxed AP property as opposed to the one step AP for the isotropic scattering. Another contribution is for the degenerate collision frequency case, which brings in additional stiffness. We propose to integrate a "body" term, which appears to be the main component in the diffusion limit. This term is precomputed once with a prescribed accuracy, via a change of variable that alleviates the stiffness. Numerical examples are presented to validate its efficiency in both kinetic and fractional diffusion regimes.
KW - Anisotropic scattering
KW - Asymptotic-preserving scheme
KW - Degenerate collision frequency
KW - Fractional diffusion
KW - Heavy tail equilibrium
UR - http://www.scopus.com/inward/record.url?scp=85062976951&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85062976951&partnerID=8YFLogxK
U2 - 10.1137/17M1138029
DO - 10.1137/17M1138029
M3 - Article
AN - SCOPUS:85062976951
SN - 1064-8275
VL - 41
SP - A422-A451
JO - SIAM Journal on Scientific Computing
JF - SIAM Journal on Scientific Computing
IS - 1
ER -