Quick search Find article
Quick search
Find article

Comparing AMR and SPH Cosmological Simulations. I. Dark Matter and Adiabatic Simulations

FREE

Brian W. O'Shea1,2, Kentaro Nagamine1,3, Volker Springel4, Lars Hernquist3 and Michael L. Norman1

Show affiliations


We compare two cosmological hydrodynamic simulation codes in the context of hierarchical galaxy formation: the Lagrangian smoothed particle hydrodynamics (SPH) code GADGET, and the Eulerian adaptive mesh refinement (AMR) code Enzo. Both codes represent dark matter with the N-body method but use different gravity solvers and fundamentally different approaches for baryonic hydrodynamics. The SPH method in GADGET uses a recently developed "entropy conserving" formulation of SPH, while for the mesh-based Enzo two different formulations of Eulerian hydrodynamics are employed: the piecewise parabolic method (PPM) extended with a dual energy formulation for cosmology, and the artificial viscosity-based scheme used in the magnetohydrodynamics code ZEUS. In this paper we focus on a comparison of cosmological simulations that follow either only dark matter, or also a nonradiative ("adiabatic") hydrodynamic gaseous component. We perform multiple simulations using both codes with varying spatial and mass resolution with identical initial conditions.

The dark matter-only runs agree generally quite well provided Enzo is run with a comparatively fine root grid and a low overdensity threshold for mesh refinement, otherwise the abundance of low-mass halos is suppressed. This can be readily understood as a consequence of the hierarchical particle-mesh algorithm used by Enzo to compute gravitational forces, which tends to deliver lower force resolution than the tree-algorithm of GADGET at early times before any adaptive mesh refinement takes place. At comparable force resolution we find that the latter offers substantially better performance and lower memory consumption than the present gravity solver in Enzo.

In simulations that include adiabatic gasdynamics we find general agreement in the distribution functions of temperature, entropy, and density for gas of moderate to high overdensity, as found inside dark matter halos. However, there are also some significant differences in the same quantities for gas of lower overdensity. For example, at z = 3 the fraction of cosmic gas that has temperature log T > 0.5 is ~80% for both Enzo ZEUS and GADGET, while it is 40%-60% for Enzo PPM. We argue that these discrepancies are due to differences in the shock-capturing abilities of the different methods. In particular, we find that the ZEUS implementation of artificial viscosity in Enzo leads to some unphysical heating at early times in preshock regions. While this is apparently a significantly weaker effect in GADGET, its use of an artificial viscosity technique may also make it prone to some excess generation of entropy that should be absent in Enzo PPM. Overall, the hydrodynamical results for GADGET are bracketed by those for Enzo ZEUS and Enzo PPM but are closer to Enzo ZEUS.


Subject headings

cosmology: theory; galaxies: formation; methods: numerical


Dates

Issue 1 (2005 September)

Received 2003 December 29, accepted for publication 2005 May 9



  1. Comparing AMR and SPH Cosmological Simulations. I. Dark Matter and Adiabatic Simulations

    Brian W. O'Shea et al. 2005 ApJS 160 1

  2. Particle detectors in the presence of boundaries

    P C W Davies et al 1989 Class. Quantum Grav. 6 1041

  3. The TOTEM Experiment at the CERN Large Hadron Collider

    The TOTEM Collaboration et al 2008 JINST 3 S08007

  4. Resonant transmission of microwaves through a finite length subwavelength metallic slit

    J R Suckling et al 2005 New J. Phys. 7 250

  5. A bound on 6D Script N = 1 supergravities

    Vijay Kumar and Washington Taylor JHEP12(2009)050

  6. Author index with titles

    2009 J. Geophys. Eng. 6 433

  7. Autocompensating quantum cryptography

    Donald S Bethune and William P Risk 2002 New J. Phys. 4 42

  8. First-principles alloy theory in oxides

    G Ceder et al 2000 Modelling Simul. Mater. Sci. Eng. 8 311

  9. Some systematics of the coupling constant dependence of N = 4 Yang-Mills

    Anirban Basu et al JHEP09(2004)045

  10. International comparison of volume measurement standards at 50 l at the CENAM (Mexico), PTB (Germany), Measurement Canada and NIST (USA)

    J M Maldonado et al 2002 Metrologia 39 91

Users also read

What's this?
This innovative new feature generates a list of articles 'also read' by other users based on them reading the original article. Article abstracts citations and references are all considered and weighted accordingly. We hope that this will help you find relevant papers for your research.

  1. Cosmological Adaptive Mesh Refinement Magnetohydrodynamics with Enzo
  2. The Santa Barbara Cluster Comparison Project: A Comparison of Cosmological Hydrodynamics Solutions
  3. The cosmic code comparison project

View by subject




Export








Please login to access our web services, or create an account if you don't yet have one.

You must have cookies enabled in your web browser to be able to login.

Username
Password

Forgotten your password? Get a new one here.