Realizability-Preserving Discontinuous Galerkin Method for Spectral Two-Moment Radiation Transport in Special Relativity

TL;DR

This paper introduces a novel numerical method for spectral two-moment radiation transport in special relativity that preserves the physical realizability of moments, ensuring nonnegative phase-space densities in simulations.

Contribution

The paper develops a realizability-preserving discontinuous Galerkin method with specialized fluxes, solvers, and limiters for relativistic radiation transport models.

Findings

Method accurately simulates relativistic radiation transport.

The scheme maintains moment realizability under various conditions.

Numerical tests confirm robustness and high accuracy.

Abstract

We present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin (DG) phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this DG-IMEX method.