English

Numerical analysis of a spherical harmonic discontinuous Galerkin method for scaled radiative transfer equations with isotropic scattering

Numerical Analysis 2023-09-19 v1 Numerical Analysis

Abstract

In highly diffusion regimes when the mean free path ε\varepsilon tends to zero, the radiative transfer equation has an asymptotic behavior which is governed by a diffusion equation and the corresponding boundary condition. Generally, a numerical scheme for solving this problem has the truncation error containing an ε1\varepsilon^{-1} contribution, that leads to a nonuniform convergence for small ε\varepsilon. Such phenomenons require high resolutions of discretizations, which degrades the performance of the numerical scheme in the diffusion limit. In this paper, we first provide a--priori estimates for the scaled spherical harmonic (PNP_N) radiative transfer equation. Then we present an error analysis for the spherical harmonic discontinuous Galerkin (DG) method of the scaled radiative transfer equation showing that, under some mild assumptions, its solutions converge uniformly in ε\varepsilon to the solution of the scaled radiative transfer equation. We further present an optimal convergence result for the DG method with the upwind flux on Cartesian grids. Error estimates of (1+O(ε))hk+1\left(1+\mathcal{O}(\varepsilon)\right)h^{k+1} (where hh is the maximum element length) are obtained when tensor product polynomials of degree at most kk are used.

Keywords

Cite

@article{arxiv.2309.09394,
  title  = {Numerical analysis of a spherical harmonic discontinuous Galerkin method for scaled radiative transfer equations with isotropic scattering},
  author = {Qiwei Sheng and Cory Hauck and Yulong Xing},
  journal= {arXiv preprint arXiv:2309.09394},
  year   = {2023}
}