Numerical Simulations for Radiation Hydrodynamics. I. Diffusion Limit

Wenlong Dai, Paul R. Woodward

Research output: Contribution to journalArticlepeer-review

33 Scopus citations


A second order accurate finite difference scheme is proposed for multidimensional radiation hydrodynamical equations in a diffusion limit. The radiation hydrodynamical equations in the limit constitute a hyperbolic system of conservation laws plus radiative heat transfer. A Godunov scheme including linear and nonlinear Riemann solvers is proposed for the set of conservation laws. The scheme with the linear Riemann solver works well for relatively strong shocks. The nonlinear Riemann solver is specially designed for flows involving strong shocks. The radiative heat conduction is treated implicitly. The treatment possesses a number of advantages over typical implicit methods. The most notable are the second order accuracy in both space and time, quick damping of numerical errors when the size of time steps is large, iterative solver and the fast convergence, the accurate treatment for the nonlinearity, and the energy conservation. Numerical examples are given to show the features of the schemes.

Original languageEnglish (US)
Pages (from-to)182-207
Number of pages26
JournalJournal of Computational Physics
Issue number1
StatePublished - May 1 1998

Bibliographical note

Funding Information:
The work presented here has been supported by the Department of Energy through Grants DE-FG02-87ER25035 and DE-FG02-94ER25207, by the National Science Foundation through Grant ASC-9309829, by NASA through Grant USRA=5555-23=NASA, and by the University of Minnesota through its Minnesota Supercomputer Institute.


Dive into the research topics of 'Numerical Simulations for Radiation Hydrodynamics. I. Diffusion Limit'. Together they form a unique fingerprint.

Cite this