Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

DG-IMEX method for a two-moment model for radiation transport in the $$\mathscr{O}$$($$v$$/$$c$) limit

Journal Article · · Journal of Computational Physics
 [1];  [2];  [1];  [3];  [3]
  1. Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
  2. Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States); University of Tennessee Knoxville, TN (United States)
  3. University of Tennessee Knoxville, TN (United States)

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $$\mathscr{O}$$($$v$$/$$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $$\mathscr{O}$$($$v$$/$$c$) Eulerian-frame number density and is consistent, to $$\mathscr{O}$$($$v$$/$$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

Research Organization:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC); USDOE National Nuclear Security Administration (NNSA); National Science Foundation (NSF)
Grant/Contract Number:
AC05-00OR22725
OSTI ID:
2538479
Journal Information:
Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 520; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English

References (53)

Moment closure hierarchies for kinetic theories journal June 1996
Study of conservation and recurrence of Runge–Kutta discontinuous Galerkin schemes for Vlasov–Poisson systems journal January 2013
Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae journal November 2020
Neutrino transport in general relativistic neutron star merger simulations journal February 2023
Relativistic transport theory journal May 1966
On the Solution of the Time-Dependent inertial-Frame Equation of radiative Transfer in Moving Media to journal August 1982
Efficient implementation of essentially non-oscillatory shock-capturing schemes journal August 1988
Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes II journal July 1989
Maximum entropy Eddington factors journal December 1978
Radiation hydrodynamics in the fluid frame journal September 1979
Radiation transfer in the fluid frame journal November 1983
Relating Eddington factors to flux limiters journal February 1984
Radiation transfer in the fluid frame: A covariant formulation journal April 1986
Radiation transfer in the fluid frame: A covariant formulation journal April 1986
Comoving-frame radiation transport for nonrelativistic fluid velocities journal May 2001
Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations journal November 1997
On maximum-principle-satisfying high order schemes for scalar conservation laws journal May 2010
On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes journal November 2010
A realizability-preserving discontinuous Galerkin method for the M1 model of radiative transfer journal July 2012
Low-storage implicit/explicit Runge–Kutta schemes for the simulation of stiff high-dimensional ODE systems journal April 2015
Realizability-preserving DG-IMEX method for the two-moment model of fermion transport journal July 2019
A numerical model for multigroup radiation hydrodynamics journal May 2011
Neutrino opacities in nuclear matter journal October 2006
Towards dense linear algebra for hybrid GPU accelerated manycore systems journal June 2010
Theory of core-collapse supernovae journal April 2007
Flash-X: A multiphysics simulation software instrument journal July 2022
Runge-Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems journal September 2001
Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations journal November 2002
Radiative Transfer in Spherically Symmetric Flows journal December 1972
A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Spacetime
  • Liebendorfer, Matthias; Messer, O. E. Bronson; Mezzacappa, Anthony
  • The Astrophysical Journal Supplement Series, Vol. 150, Issue 1 https://doi.org/10.1086/380191
journal January 2004
Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics journal April 2006
A new Multi-Dimensional General Relativistic Neutrino Hydrodynamic code for Core-Collapse Supernovae. i. Method and code Tests in Spherical Symmetry journal June 2010
Three-Dimensional Boltzmann Hydro code for core Collapse in Massive Stars. i. Special Relativistic Treatments journal September 2014
Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers journal February 2015
A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics journal September 2015
Conservative formulations of general relativistic kinetic theory journal July 2003
Conservative special relativistic radiative transfer for multidimensional astrophysical simulations: Motivation and elaboration journal August 2005
Conservative 3 + 1 general relativistic variable Eddington tensor radiation transport equations journal May 2013
Conservative 3 + 1 general relativistic Boltzmann equation journal July 2013
Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables journal October 1988
Asymptotic Analysis of Upwind Discontinuous Galerkin Approximation of the Radiative Transport Equation in the Diffusive Limit journal January 2010
Anderson Acceleration for Fixed-Point Iterations journal January 2011
A Regularized Entropy-Based Moment Method for Kinetic Equations journal January 2019
Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity journal June 2011
Iterative Procedures for Nonlinear Integral Equations journal October 1965
Explosion Mechanisms of Core-Collapse Supernovae journal November 2012
Discontinuous Finite Element Transport Solutions in Thick Diffusive Problems journal March 2001
Fornax: A Flexible Code for Multiphysics Astrophysical Simulations journal February 2019
Chimera: A Massively Parallel Code for Core-collapse Supernova Simulations journal April 2020
A DG-IMEX Method for Two-moment Neutrino Transport: Nonlinear Solvers for Neutrino–Matter Coupling* journal April 2021
Multigroup Radiation Magnetohydrodynamics Based on Discrete Ordinates including Compton Scattering journal October 2022
Discontinuous Galerkin methods for the one-dimensional Vlasov-Poisson system journal November 2011
High-order entropy-based closures for linear transport in slab geometry journal January 2011

Similar Records

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity
Journal Article · Sun Feb 02 23:00:00 EST 2025 · Physical Review. D. · OSTI ID:2538323

Realizability-preserving DG-IMEX method for the two-moment model of fermion transport
Journal Article · Wed Mar 27 00:00:00 EDT 2019 · Journal of Computational Physics · OSTI ID:1513400

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model
Journal Article · Fri Apr 15 00:00:00 EDT 2022 · Journal of Computational Physics · OSTI ID:2538513