paper

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

arXiv:2203.12177 · doi:10.1017/S1743921322001429

Abstract

We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used -- but less accurate -- closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities.

7 pages, 2 figures, to appear in the proceedings of IAUS362: The predictive power of computational astrophysics as discovery tool (D. Bisikalo, T. Hanawa, C. Boily, J. Stone, eds.). arXiv admin note: substantial text overlap with arXiv:2202.08778

References in corpus (4)