Realizability-preserving monolithic convex limiting in continuous Galerkin discretizations of the M1 model of radiative transfer

We discretize the M1 model of radiative transfer using continuous finite elements and propose a tailor-made monolithic convex limiting (MCL) procedure for enforcing physical realizability. The M1 system of nonlinear balance laws for the zeroth and first moments of a probability distribution function is derived from the linear Boltzmann equation and equipped with an entropy-based closure for the second moment. To ensure hyperbolicity and physical admissibility, evolving moments must stay in an invariant domain representing a convex set of realizable states. We first construct a low-order method that is provably invariant domain preserving (IDP). Introducing intermediate states that represent spatially averaged exact solutions of homogeneous Riemann problems, we prove that these so-called bar states are realizable in any number of space dimensions. This key auxiliary result enables us to show the IDP property of a fully discrete scheme with a diagonally implicit treatment of reactive terms. To achieve high resolution, we add nonlinear correction terms that are constrained using a two-step MCL algorithm. In the first limiting step, local bounds are imposed on each conserved variable to avoid spurious oscillations and maintain positivity of the scalar-valued zeroth moment (particle density). The second limiting step constrains the magnitude of the vector-valued first moment to be realizable. The flux-corrected finite element scheme is provably IDP. Its ability to prevent nonphysical behavior while attaining high-order accuracy in smooth regions is verified in a series of numerical tests. The developed methodology provides a robust simulation tool for dose calculation in radiotherapy.