Self-Consistent Solution of Cosmological Radiation-Hydrodynamics and Chemical Ionization

Abstract

We consider a PDE system comprising compressible hydrodynamics, flux-limited diffusion radiation transport and chemical ionization kinetics in a cosmologically-expanding universe. Under an operator-split framework, the cosmological hydrodynamics equations are solved through the Piecewise Parabolic Method, as implemented in the Enzo community hydrodynamics code. The remainder of the model, including radiation transport, chemical ionization kinetics, and gas energy feedback, form a stiff coupled PDE system, which we solve using a fully-implicit inexact Newton approach, and which forms the crux of this paper. The inner linear Newton systems are solved using a Schur complement formulation, and employ a multigrid-preconditioned conjugate gradient solver for the inner Schur systems. We describe this approach and provide results on a suite of test problems, demonstrating its accuracy, robustness, and scalability to very large problems.

36 pages, 13 figures, submitted to Journal of Computational Physics; added NSF grant acknowledgement; final accepted version including revisions

References in corpus (4)

• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• The Physics and Early History of the Intergalactic Medium
• Lyman-alpha radiative transfer during the Epoch of Reionization: contribution to 21-cm signal fluctuations
• Late Reheating of the IGM by Quasars: A Radiation Hydrodynamical Simulation of Helium II Reionization

Cited by in corpus (15)

• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• RAMSES-RT: Radiation hydrodynamics in the cosmological context
• Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests
• A Radiation Transfer Solver for Athena using Short Characteristics
• Direct Collapse to Supermassive Black Hole Seeds with Radiative Transfer: Isolated Halos
• Optimized Multi-Frequency Spectra for Applications in Radiative Feedback and Cosmological Reionization
• A Hybrid Godunov Method for Radiation Hydrodynamics
• Direct Collapse to Supermassive Black Hole Seeds with Radiation Transfer: Cosmological Halos
• A novel approach for accurate radiative transfer in cosmological hydrodynamic simulations
• Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
• Smoothed Particle Radiation Hydrodynamics: Two-Moment method with Local Eddington Tensor Closure
• Fully Coupled Simulation of Cosmic Reionization. III. Stochastic Early Reionization by the Smallest Galaxies
• Three-Dimensional Modelling of Ionized Gas. I. Did very massive stars of different metallicities drive the second cosmic reionization?
• High-Performance Astrophysical Simulations and Analysis with Python
• Star Formation with Adaptive Mesh Refinement Radiation Hydrodynamics