A higher-order Godunov method for radiation hydrodynamics: Radiation subsystem

Michael David Sekora, James McLellan Stone

Research output: Contribution to journalArticlepeer-review

11 Scopus citations


A higher-order Godunov method for the radiation subsystem of radiation hydrodynamics is presented. A key ingredient of the method is the direct coupling of stiff source term effects to the hyperbolic structure of the system of conservation laws; it is composed of a predictor step that is based on Duhamel's principle and a corrector step that is based on Picard iteration. The method is second-order accurate in both time and space, unsplit, asymptotically preserving, and uniformly well behaved from the photon free streaming (hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and to the strong equilibrium diffusion (hyperbolic) limit. Numerical tests demonstrate second-order convergence across various parameter regimes.

Original languageEnglish (US)
Pages (from-to)135-152
Number of pages18
JournalCommunications in Applied Mathematics and Computational Science
Issue number1
StatePublished - 2009

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Computational Theory and Mathematics
  • Applied Mathematics


  • Asymptotic preserving methods
  • Godunov methods
  • Hyperbolic conservation laws
  • Radiation hydrodynamics
  • Stiff relaxation
  • Stiff source terms


Dive into the research topics of 'A higher-order Godunov method for radiation hydrodynamics: Radiation subsystem'. Together they form a unique fingerprint.

Cite this