REVIEW
Introducing Enzo, an AMR Cosmology Application
Astrophysics · Submitted 1 Mar 2004 · 15 references
Pith's verdict: A tool-release paper, not a science result: Enzo is real and became a standard AMR cosmology code, but §3's comparison with Lagrangian codes is asserted without any supporting curve.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Enzo is a publicly released, parallel adaptive-mesh code that concentrates computing power only where cosmic structure forms, coupling dark matter particles to a high-order gas solver.
The reading
What carries the argument
The load-bearing structure is the block-structured adaptive mesh hierarchy of Berger and Colella, in which the root grid carries the periodic cosmological volume and child grids, refined by integer factors of two in these runs, are laid down wherever user-chosen criteria demand. Two solvers sit on this hierarchy and must agree with each other: the elliptic Poisson equation is handled by FFT on the triply periodic root grid and by multigrid relaxation on the non-periodic subgrids, with forces then interpolated to particle positions, while the gas is advanced by a cosmology-modified piecewise parabolic method that solves both the internal energy and total energy equations everywhere so that strong shocks and hypersonic bulk flows keep correct pressures and temperatures. The dual energy formalism is the specific device that lets the code survive cosmological hypersonic flows where kinetic energy dwarfs thermal energy by many orders of magnitude.
Extended reading notes
Core claim
The claim is that a single code can deliver both the large volume and the small-scale resolution that cosmological structure formation demands, without the fixed mass resolution and poor shock handling the authors attribute to smoothed particle hydrodynamics. Enzo achieves this by treating each grid patch as an object in a distributed hierarchy, with a three-zone-deep ghost-zone layer wide enough to feed the hydrodynamics stencil, and by distributing patches across processors with a simple dynamic load-balancing scheme. On top of that skeleton the authors place a full physics set: two radiative cooling models, one assuming collisional equilibrium at a fixed metallicity and one solving a 28-equation primordial reaction network for nine species, several uniform ultraviolet background models after Haardt and Madau, and a heuristic star formation and feedback prescription in which rapidly cooling, densifying gas is converted into star particles that return metals and thermal energy as winds. The quantitative demonstration is that in a $3~h^{-1}$ Mpc box with $64^3$ or $128^3$ root grids and five or six refinement levels, cells populate the entire density range rather than clustering only at high overdensity, and the mean cell mass at the end of a run is about five to ten times smaller than the root-grid mass.
Load-bearing premise
The gravity solver must produce a potential that is consistent across refinement boundaries, since the periodic root grid is handled by a Fourier transform and the subgrids by multigrid relaxation, and the paper states this in a single sentence without an error estimate, a convergence test, or a description of how subgrid boundary values are supplied from the parent grid.
Editorial extensions
If this is right
- Cosmological simulations can cover several megaparsecs while resolving kiloparsec-scale structure, which is the range galaxy formation requires.
- Low-density gas is well sampled rather than undersampled, unlike particle-based methods, so the Lyman-alpha forest and other underdense regions can be modelled with the same run that resolves collapsed objects.
- The same code base serves non-cosmological problems, since the mesh, gravity, and hydro machinery do not depend on the expansion of the universe.
- Because the code is released publicly with HDF5 output and MPI portability across IBM SP, SGI Origin, and Linux clusters, results become reproducible and extensible by groups that did not write it.
- Adding magnetohydrodynamics and flux-limited radiation diffusion, which the authors say is underway, extends the same hierarchy to magnetized and radiatively driven astrophysical flows.
Reading between the lines
- The paper never quantifies how accurately the FFT root-grid potential and the multigrid subgrid potential agree at refinement boundaries, so a reader cannot tell from this document whether particle trajectories acquire level-dependent errors; a convergence study comparing a fully refined uniform run against the adaptive run would settle it.
- Because the refinement criteria are overdensity, slope, shocks, Jeans length, and cooling time, the effective resolution depends on thresholds the user picks; the authors' own comparison of a doubled overdensity threshold shows the cell count falls while the density distribution barely changes, which suggests the mass-resolution gain is somewhat threshold-insensitive but leaves the accuracy of tha
- The star formation prescription is explicitly heuristic, converting rapidly cooling gas into particles that return metals and energy; whether such a sub-grid rule reproduces observed galaxy stellar masses and wind-driven metal enrichment is a separate question the paper does not attempt to answer here.
- The stated contrast with smoothed particle hydrodynamics rests on shock resolution and low-density sampling; a direct head-to-head on a problem with a known analytic shock or a known mass function would turn that contrast from a design argument into a measured one.
Original abstract
In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. Enzo solves dark matter N-body dynamics using the particle-mesh technique. The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. We also show results illustrating properties of the adaptive mesh portion of the code. Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ .
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. Enzo solves dark matter N-body dynamics using the particle-mesh technique. The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. We also show results illustrating properties of the adaptive mesh portion of the code. Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ .
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. | contribution 1.00 | 0.83 | 0.36 | asserted 1.00 |
| s1 | Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. | contribution 0.88 | 0.24 | 0.10 | asserted 0.93 |
| s2 | Enzo solves dark matter N-body dynamics using the particle-mesh technique. | method 0.99 | 0.57 | 0.04 | asserted 1.00 |
| s3 | The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). | method 1.00 | 0.64 | 0.03 | asserted 1.00 |
| s4 | Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. | method 0.99 | 0.67 | 0.03 | asserted 1.00 |
| s5 | Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. | method 0.66 | 0.81 | 0.04 | asserted 1.00 |
| s6 | We also show results illustrating properties of the adaptive mesh portion of the code. | result 0.78 | 0.92 | 0.06 | asserted 1.00 |
| s7 | Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ . | contribution 0.80 | 0.32 | 0.04 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:astro-ph/0403044 by the authors.
References
How later papers used it
47 later papers cite it
6 later papers by other authors are listed below: 5 used its method or tool and 1 compared its own results with it. 2 more papers by the same authors are listed: 2 used its method or tool.
Uses its method or tool 7 of 28 papers with this label · 2 by the same authors
-
Used the cosmological codes GADGET3 and Enzo v2.1.1 with initial conditions that include $v_$ self-consistently.
The impact of the supersonic baryon-dark matter velocity difference on the z~20 21cm background (2012) · 1.00 · paraphrased: the licence does not allow quoting
-
Carried out all simulations with its modified version of the multi-dimensional hydrodynamics and N-body grid-based code Enzo.
The effect of binding energy and resolution in simulations of the common envelope binary interaction (2017) · 1.00 · paraphrased: the licence does not allow quoting
-
Used these initial conditions with the cosmological codes GADGET3 and Enzo v2.1.1 with the default piecewise parabolic method to solve the equations of hydrodynamics.
The formation of the first cosmic structures and the physics of the z~20 Universe (2012) · 1.00 · paraphrased: the licence does not allow quoting
-
Used the publicly available codes, a Lagrangian TreeSPH code and a hybrid adaptive mesh refinement, grid-based code, to run the simulations in this work.
How to Zoom: Bias, Contamination, and Lagrange Volumes in Multimass Cosmological Simulations (2013) · 1.00 · paraphrased: the licence does not allow quoting
-
Created the initial conditions for the simulations using ENZO, a publicly available cosmological N-body + hydrodynamics code.
The Bispectrum of f(R) Cosmologies (2011) · 1.00 · paraphrased: the licence does not allow quoting
By the same authors
-
Used the code Enzo to simulate the formation of the first star within a periodic volume 250 comoving kpc on a side.
Accretion onto the First Stellar Mass Black Holes (2008) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Presented simulations performed with the adaptive mesh refinement code Enzo.
The Formation of Population III Binaries from Cosmological Initial Conditions (2009) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
Compares with it 1 of 2 papers with this label
-
Compared the result from the identical test, using the pressure-energy formulation without the $∇h$ terms in the EOM, and results from a high-resolution run with Enzo, a grid code.
A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems (2012) · 0.96 · paraphrased: the licence does not allow quoting
How the literature uses it
45 arXiv papers cite this in 54 sentences
Papers citing it · 16 sentences without a confident label
Uses its method or tool 26 papers
-
For the data presented in this paper, we solve the ideal MHD equations with self-gravity using the adaptive mesh refinement (AMR) code Enzo [ref] extended to MHD by [ref] .
The Two States of Star Forming Clouds · Numerical Method, Simulations, Analysis · 0.99
-
Our initial conditions are spherically-averaged halos derived from cosmological initial conditions with the Enzo adaptive mesh refinement (AMR) code [ref] that are photoionized by the star in the ZEUS-MP reactive flow radiative transfer code [ref] .
The Destruction of Cosmological Minihalos by Primordial Supernovae · Introduction · 1.00
-
In order to study the effects of a global photodissociating background, we take a single “standard” Pop III star formation simulation using the Enzo code [ref] and gradually increase the strength of the Lyman-Werner (H$_2$ photodissociating) background.
Local and Global Radiative Feedback from Population III Star Formation · Effects of a Global Photodissociating Background · 0.99
Uses its data or values 1 paper
-
< g r a p h i c s > Generating an image from pre-generated set of isodensity surfaces with the isolated galaxy test simulation of [ref] .
AstroBlend: An Astrophysical Visualization Package for Blender · Summary and Future Plans · 0.89
Compares with it 2 papers
-
We also compare the result from the identical test in , using the pressure-energy formulation but without including the $∇h$ terms in the EOM, as well as the results from a high-resolution run with Enzo [ref] , a grid code.
A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems · The “Blob” Test · 0.96
-
The physical properties of the forest obtained from the TreePM/SPH code GADGET-2 are in agreement at the percent level with those inferred from the moving-mesh code AREPO [ref] and with the Eulerian code ENZO [ref] .
IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function · Hydrodynamical Simulations · 0.82
Cites it as background 4 papers
-
The parallelization of AMR alogrithms, which should combine the cost/time savings of both methods is not straight forward however and there have been many different approaches [ref] , .
Efficient Parallelization for AMR MHD Multiphysics Calculations; Implementation in AstroBEAR · Introduction · 0.97
-
Extensive theory has since been developed for the method in different contexts, and mesh refinement algorithms have been widely adopted in the literature [ref] .
Toward a Consistent Framework for High Order Mesh Refinement Schemes in Numerical Relativity · Introduction · 0.92
-
There are many cosmological hydrodynamic simulation codes available today [ref] .
EUNHA: a new cosmological hydro simulation code · Introduction · 0.94
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.