Thermal conduction in cosmological SPH simulations

• Published in: Monthly notices of the Royal Astronomical Society Vol. 351; no. 2; pp. 423 - 435
• Main Authors: Jubelgas, Martin, Springel, Volker, Dolag, Klaus
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Science Ltd 21.06.2004; Blackwell Science
• Subjects: Astronomy; conduction; Earth, ocean, space; Exact sciences and technology; Fundamental astronomy and astrophysics. Instrumentation, techniques, and astronomical observations; galaxies: clusters: general; Galaxy clusters; Galaxy groups, clusters, and superclusters. Large-scale structure of the universe; Mathematical procedures and computer techniques; methods: numerical; Observation and data reduction techniques. Computer modeling and simulation; Stellar systems. Galactic and extragalactic objects and systems. The universe; Spherical symmetry; Plasma; galaxies: clusters: general; Local equilibrium; Gas dynamics; Hydrostatic equilibrium; Saturation; Digital simulation; Luminosity; Cooling flow; Numerical method; methods: numerical; Galaxy clusters; Intracluster matter; conduction; Hydrodynamic model; Diffusion equation; Noise temperature; Thermal conduction
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1111/j.1365-2966.2004.07801.x

Thermal conduction in the intracluster medium has been proposed as a possible heating mechanism for offsetting central cooling losses in rich clusters of galaxies. However, because of the coupled non-linear dynamics of gas subject to radiative cooling and thermal conduction, cosmological hydrodynamical simulations are required to predict reliably the effects of heat conduction on structure formation. In this study, we introduce a new formalism to model conduction in a diffuse ionized plasma using smoothed particle hydrodynamics (SPH), and we implement it in the parallel TreePM/SPH-code gadget-2. We consider only isotropic conduction and assume that magnetic suppression can be described in terms of an effective conductivity, taken as a fixed fraction of the temperature-dependent Spitzer rate. We also account for saturation effects in low-density gas. Our formulation manifestly conserves thermal energy even for individual and adaptive time-steps, and is stable in the presence of small-scale temperature noise. This allows us to evolve the thermal diffusion equation with an explicit time integration scheme along with the ordinary hydrodynamics. We use a series of simple test problems to demonstrate the robustness and accuracy of our method. We then apply our code to spherically symmetric realizations of clusters, constructed under the assumptions of hydrostatic equilibrium and a local balance between conduction and radiative cooling. While we confirm that conduction can efficiently suppress cooling flows for an extended period of time in these isolated systems, we do not find a similarly strong effect in a first set of clusters formed in self-consistent cosmological simulations. However, their temperature profiles are significantly altered by conduction, as is the X-ray luminosity.