A hybrid Godunov method for radiation hydrodynamics
Abstract
From a mathematical perspective, radiation hydrodynamics can be thought of as a system of hyperbolic balance laws with dual multiscale behavior (multiscale behavior associated with the hyperbolic wave speeds as well as multiscale behavior associated with source term relaxation). With this outlook in mind, this paper presents a hybrid Godunov method for onedimensional radiation hydrodynamics that is 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. Moreover, one finds that the technique preserves certain asymptotic limits. The method incorporates a backward Euler upwinding scheme for the radiation energy density E{sub r} and flux F{sub r} as well as a modified Godunov scheme for the material density {rho}, momentum density m, and energy density E. The backward Euler upwinding scheme is firstorder accurate and uses an implicit HLLE flux function to temporally advance the radiation components according to the material flow scale. The modified Godunov scheme is secondorder accurate and directly couples stiff source term effects to the hyperbolic structure of the system of balance laws. This Godunov technique is composed of a predictor step that is based on Duhamel's principle and a corrector stepmore »
 Authors:
 Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ 08544 (United States)
 (United States)
 Publication Date:
 OSTI Identifier:
 21417245
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 229; Journal Issue: 19; Other Information: DOI: 10.1016/j.jcp.2010.05.024; PII: S00219991(10)002858; Copyright (c) 2010 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.
 Country of Publication:
 United States
 Language:
 English
 Subject:
 97 MATHEMATICAL METHODS AND COMPUTING; APPROXIMATIONS; ASYMPTOTIC SOLUTIONS; BALANCES; DIFFUSION; ENERGY DENSITY; EQUILIBRIUM; FUNCTIONS; HYDRODYNAMICS; MATHEMATICAL MODELS; ONEDIMENSIONAL CALCULATIONS; CALCULATION METHODS; FLUID MECHANICS; MATHEMATICAL SOLUTIONS; MEASURING INSTRUMENTS; MECHANICS; WEIGHT INDICATORS
Citation Formats
Sekora, Michael D., Email: sekora@math.princeton.ed, Stone, James M., and Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544. A hybrid Godunov method for radiation hydrodynamics. United States: N. p., 2010.
Web. doi:10.1016/j.jcp.2010.05.024.
Sekora, Michael D., Email: sekora@math.princeton.ed, Stone, James M., & Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544. A hybrid Godunov method for radiation hydrodynamics. United States. doi:10.1016/j.jcp.2010.05.024.
Sekora, Michael D., Email: sekora@math.princeton.ed, Stone, James M., and Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544. 2010.
"A hybrid Godunov method for radiation hydrodynamics". United States.
doi:10.1016/j.jcp.2010.05.024.
@article{osti_21417245,
title = {A hybrid Godunov method for radiation hydrodynamics},
author = {Sekora, Michael D., Email: sekora@math.princeton.ed and Stone, James M. and Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544},
abstractNote = {From a mathematical perspective, radiation hydrodynamics can be thought of as a system of hyperbolic balance laws with dual multiscale behavior (multiscale behavior associated with the hyperbolic wave speeds as well as multiscale behavior associated with source term relaxation). With this outlook in mind, this paper presents a hybrid Godunov method for onedimensional radiation hydrodynamics that is 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. Moreover, one finds that the technique preserves certain asymptotic limits. The method incorporates a backward Euler upwinding scheme for the radiation energy density E{sub r} and flux F{sub r} as well as a modified Godunov scheme for the material density {rho}, momentum density m, and energy density E. The backward Euler upwinding scheme is firstorder accurate and uses an implicit HLLE flux function to temporally advance the radiation components according to the material flow scale. The modified Godunov scheme is secondorder accurate and directly couples stiff source term effects to the hyperbolic structure of the system of balance laws. This Godunov technique is composed of a predictor step that is based on Duhamel's principle and a corrector step that is based on Picard iteration. The Godunov scheme is explicit on the material flow scale but is unsplit and fully couples matter and radiation without invoking a diffusiontype approximation for radiation hydrodynamics. This technique derives from earlier work by Miniati and Colella (2007) . Numerical tests demonstrate that the method is stable, robust, and accurate across various parameter regimes.},
doi = {10.1016/j.jcp.2010.05.024},
journal = {Journal of Computational Physics},
number = 19,
volume = 229,
place = {United States},
year = 2010,
month = 9
}

We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas and radiation variables including the nonstiff source terms, and a local implicit method to integrate the stiff source terms. We adopt the M{sub 1} closure relation and include all leading source terms to O({beta}{tau}). We employ the reduced speed of light approximation (RSLA) with subcycling of the radiation variables in order to reduce computational costs. Our code is dimensionally unsplitmore »

A GODUNOV METHOD FOR MULTIDIMENSIONAL RADIATION MAGNETOHYDRODYNAMICS BASED ON A VARIABLE EDDINGTON TENSOR
We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusionlike approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiationmodified linear acoustic and magnetosonic waves, the structure of radiationmodified shocks, and twodimensional tests of photon bubble instability andmore » 
Implicitexplicit hybrid method for Lagrangian hydrodynamics
We describe a new implicitexplicit hybrid method for solving the equations of hydrodynamics. The scheme is an extension of the explicit secondorder piecewiseparabolic method (PPM) which is unconditionally stable. The scheme is thus of the Godunov type. It is conservative, accurate to second order in both space and time, and makes use of a nonlinear Riemann solver to obtain fluxes of the conserved quantities. The hybrid character of the method provides increased accuracy and computational efficiency. Switching between implicit and explicit formulations occurs smoothly and in a natural way and is performed separately for each characteristic family of waves. Themore »