Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity
Journal Article
·
· Physical Review. D.
- The Ohio State Univ., Columbus, OH (United States)
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States); Univ. of Tennessee, Knoxville, TN (United States)
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Here 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 phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit 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 discontinuous-Galerkin-implicit-explicit method.
- Research Organization:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Organization:
- National Science Foundation (NSF); USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC)
- Grant/Contract Number:
- AC05-00OR22725
- OSTI ID:
- 2538323
- Journal Information:
- Physical Review. D., Journal Name: Physical Review. D. Journal Issue: 4 Vol. 111; ISSN 2470-0010
- Publisher:
- American Physical Society (APS)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit
A Realizability-Preserving Discontinuous Galerkin Method for the $M_1$ Model of Radiative Transfer
Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator
Journal Article
·
Wed Oct 02 20:00:00 EDT 2024
· Journal of Computational Physics
·
OSTI ID:2538479
A Realizability-Preserving Discontinuous Galerkin Method for the $M_1$ Model of Radiative Transfer
Journal Article
·
Sat Dec 31 23:00:00 EST 2011
· Journal of Computational Physics
·
OSTI ID:1048144
Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator
Journal Article
·
Sun Aug 09 20:00:00 EDT 2020
· Nuclear Fusion
·
OSTI ID:1659475