FULLY COUPLED SIMULATION OF COSMIC REIONIZATION. I. NUMERICAL METHODS AND TESTS
Abstract
We describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large ∼(100 Mpc){sup 3} cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes selfconsistently on the same mesh, as opposed to a postprocessing approach which coarsegrains the radiative transfer. We do, however, employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use to advance the system of equations. Radiation transport is done in the gray fluxlimited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties ofmore »
 Authors:
 CASS, University of California, San Diego, 9500 Gilman Drive La Jolla, CA 920930424 (United States)
 Southern Methodist University, 6425 Boaz Lane, Dallas, TX 75205 (United States)
 SDSC, University of California, San Diego, 9500 Gilman Drive La Jolla, CA 920930505 (United States)
 Center for Relativistic Astrophysics, Georgia Institute of Technology, 837 State Street, Atlanta, GA 30332 (United States)
 Publication Date:
 OSTI Identifier:
 22340107
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Astrophysical Journal, Supplement Series; Journal Volume: 216; Journal Issue: 1; Other Information: Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ACCURACY; COSMOLOGY; DIFFUSION; EQUATIONS; GEOMETRY; HYDRODYNAMICS; IONIZATION; LAWRENCE LIVERMORE NATIONAL LABORATORY; MASS; NONLUMINOUS MATTER; POINT SOURCES; RADIANT HEAT TRANSFER; RADIATION TRANSPORT; REYNOLDS NUMBER; SIMULATION; STARS; VALIDATION
Citation Formats
Norman, Michael L., So, Geoffrey C., Reynolds, Daniel R., Harkness, Robert P., and Wise, John H.. FULLY COUPLED SIMULATION OF COSMIC REIONIZATION. I. NUMERICAL METHODS AND TESTS. United States: N. p., 2015.
Web. doi:10.1088/00670049/216/1/16.
Norman, Michael L., So, Geoffrey C., Reynolds, Daniel R., Harkness, Robert P., & Wise, John H.. FULLY COUPLED SIMULATION OF COSMIC REIONIZATION. I. NUMERICAL METHODS AND TESTS. United States. doi:10.1088/00670049/216/1/16.
Norman, Michael L., So, Geoffrey C., Reynolds, Daniel R., Harkness, Robert P., and Wise, John H.. 2015.
"FULLY COUPLED SIMULATION OF COSMIC REIONIZATION. I. NUMERICAL METHODS AND TESTS". United States.
doi:10.1088/00670049/216/1/16.
@article{osti_22340107,
title = {FULLY COUPLED SIMULATION OF COSMIC REIONIZATION. I. NUMERICAL METHODS AND TESTS},
author = {Norman, Michael L. and So, Geoffrey C. and Reynolds, Daniel R. and Harkness, Robert P. and Wise, John H.},
abstractNote = {We describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large ∼(100 Mpc){sup 3} cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes selfconsistently on the same mesh, as opposed to a postprocessing approach which coarsegrains the radiative transfer. We do, however, employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use to advance the system of equations. Radiation transport is done in the gray fluxlimited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the wellknown inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. We illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with 3200{sup 3} Eulerian grid cells and dark matter particles.},
doi = {10.1088/00670049/216/1/16},
journal = {Astrophysical Journal, Supplement Series},
number = 1,
volume = 216,
place = {United States},
year = 2015,
month = 1
}

Here, we describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large similar to(100 Mpc)(3) cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes selfconsistently on the same mesh, as opposed to a postprocessing approach which coarsegrains the radiative transfer. But, we employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use tomore »Cited by 17

Fully coupled simulation of cosmic reionization. II. Recombinations, clumping factors, and the photon budget for reionization
We use a fully coupled cosmological simulation including dark matter dynamics, multispecies hydrodynamics, nonequilibrium chemical ionization, fluxlimited diffusion radiation transport, and a parameterized model of star formation and feedback (thermal and radiative) to investigate the epoch of hydrogen reionization in detail. In this paper, the first of several application papers, we investigate the mechanics of reionization from stellar sources forming in highz galaxies, the utility of various formulations for the gas clumping factor on accurately estimating the effective recombination time in the intergalactic medium (IGM), and the photon budget required to achieve reionization. We also test the accuracy of themore » 
ON IMPROVING ANALYTICAL MODELS OF COSMIC REIONIZATION FOR MATCHING NUMERICAL SIMULATION
The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emergedmore » 
A fully coupled method for massively parallel simulation of hydraulically driven fractures in 3dimensions: FULLY COUPLED PARALLEL SIMULATION OF HYDRAULIC FRACTURES IN 3D
This study describes a fully coupled finite element/finite volume approach for simulating fieldscale hydraulically driven fractures in three dimensions, using massively parallel computing platforms. The proposed method is capable of capturing realistic representations of local heterogeneities, layering and natural fracture networks in a reservoir. A detailed description of the numerical implementation is provided, along with numerical studies comparing the model with both analytical solutions and experimental results. The results demonstrate the effectiveness of the proposed method for modeling largescale problems involving hydraulically driven fractures in three dimensions.Cited by 4 
COSMIC REIONIZATION ON COMPUTERS: NUMERICAL AND PHYSICAL CONVERGENCE
In this paper I show that simulations of reionization performed under the Cosmic Reionization On Computers project do converge in space and mass, albeit rather slowly. A fully converged solution (for a given star formation and feedback model) can be determined at a level of precision of about 20%, but such a solution is useless in practice, since achieving it in productiongrade simulations would require a large set of runs at various mass and spatial resolutions, and computational resources for such an undertaking are not yet readily available. In order to make progress in the interim, I introduce a weakmore »