A parallel, energy-stable low-rank integrator for nonlinear multi-scale thermal radiative transfer

Thermal radiative transfer models physical phenomena ranging from supernovas in astrophysics to radiation from a hohlraum striking a fusion target in plasma physics. Transport and absorption of particles in radiative transfer at different rates lead to a complex interaction between the material and particles that involves highly varying time scales. Resolving these effects can require prohibitively small step sizes, which, combined with nonlinear effects and the particle density’s high-dimensional phase space, render conventional numerical methods computationally expensive. This work presents an asymptotic–preserving, mass conservative, rank-adaptive, and parallel integrator for a macro–micro decomposition-based dynamical low-rank approximation of the thermal radiative transfer equations. The proposed integrator efficiently incorporates reflection-transmission type boundary conditions in the low-rank factors. It captures the nonlinear effects of thermal radiation and is energy stable with the step size restriction capturing both hyperbolic and parabolic CFL conditions. The efficacy of the proposed integrator is demonstrated with numerical experiments.