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

Source: [Self-Consistent Solution of Cosmological Radiation-Hydrodynamics and Chemical Ionization](<https://www.sota2.com/research/paper/68fa2eef1e209131ce867811>)

Collection index: [llms.txt](<https://www.sota2.com/research/llms.txt>)

Data updated: 2025-10-28T13:53:31.692Z

Authors: Daniel R. Reynolds, John C. Hayes, Pascal Paschos, Michael L. Norman

Published: 2009-01-08T19:26:19.000Z

## 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.

## Recorded benchmark results

A research result applies to its stated dataset, split, metric and evaluation configuration. Compare like-for-like settings, preserve the metric's unit and direction, and distinguish paper-reported results from independent reproduction. The best result recorded here is not a claim of exhaustive, current global coverage.

## Original sources and implementation

[Original paper](<https://arxiv.org/abs/0901.1110>)

[PDF](<https://arxiv.org/pdf/0901.1110.pdf>)