Journal Article

A comparison of cosmological codes: properties of thermal gas and shock waves in large-scale structures

F. Vazza, K. Dolag, D. Ryu, G. Brunetti, C. Gheller, H. Kang and C. Pfrommer

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 418, issue 2, pages 960-985
Published in print December 2011 | ISSN: 0035-8711
Published online November 2011 | e-ISSN: 1365-2966 | DOI:
A comparison of cosmological codes: properties of thermal gas and shock waves in large-scale structures

Show Summary Details


Cosmological hydrodynamical simulations are a valuable tool for understanding the growth of large-scale structure and the observables connected with this. Yet, comparably little attention has been given to validation studies of the properties of shocks and of the resulting thermal gas between different numerical methods – something of immediate importance as gravitational shocks are responsible for generating most of the entropy of the large-scale structure in the Universe. Here, we present results for the statistics of thermal gas and the shock wave properties for a large volume simulated with three different cosmological numerical codes: the Eulerian total variations diminishing (TVD) code, the Eulerian piecewise parabolic method based code enzo and the Lagrangian smoothed particle hydrodynamics (SPH) code gadget. Starting from a shared set of initial conditions, we present convergence tests for a cosmological volume of side-length 100 Mpc h−1, studying in detail the morphological and statistical properties of the thermal gas as a function of mass and spatial resolution in all codes. By applying shock-finding methods to each code, we measure the statistics of shock waves and the related cosmic ray acceleration efficiencies, within the sample of simulations and for the results of the different approaches. We discuss the regimes of uncertainties and disagreement among codes, with a particular focus on the results at the scale of galaxy clusters. Even if the bulk of thermal and shock properties is reasonably in agreement among the three codes, yet some significant differences exist (especially between Eulerian methods and SPH). In particular, we report (a) differences of huge factors (∼10–100) in the values of average gas density, temperature, entropy, Mach number and shock thermal energy flux in the most rarefied regions of the simulations (ρ/ρcr < 1) between grid and SPH methods; (b) the hint of an entropy core inside clusters simulated in grid codes; (c) significantly different phase diagrams of shocked cells in grid codes compared to SPH and (d) sizable differences in the morphologies of accretion shocks between grid and SPH methods.

Keywords: methods: numerical; galaxies: clusters: general; intergalactic medium; large-scale structure of Universe

Journal Article.  16977 words.  Illustrated.

Subjects: Astronomy and Astrophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.