Simulations of jet heating in galaxy clusters: Successes and challenges

Davide Martizzi*, Eliot Quataert, Claude André Faucher-Giguère, Drummond Fielding

*Corresponding author for this work

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealized, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated with multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at ∼50 kpc scales where it only represents a few per cent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multiphase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.

Original languageEnglish (US)
Pages (from-to)2465-2486
Number of pages22
JournalMonthly Notices of the Royal Astronomical Society
Volume483
Issue number2
DOIs
StatePublished - Jan 1 2019

Fingerprint

galaxies
heating
cooling
simulation
Godunov method
shock heating
kinetics
thermal instability
gases
gas
active galactic nuclei
physics
flow velocity
compression
grids
rate
expansion
high resolution
configurations
method

Keywords

  • Galaxies: active
  • Galaxies: clusters: general
  • Galaxies: jets
  • Methods: numerical

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

Martizzi, Davide ; Quataert, Eliot ; Faucher-Giguère, Claude André ; Fielding, Drummond. / Simulations of jet heating in galaxy clusters : Successes and challenges. In: Monthly Notices of the Royal Astronomical Society. 2019 ; Vol. 483, No. 2. pp. 2465-2486.
@article{9c804f5d693e4aa2bf609f75c58b1c34,
title = "Simulations of jet heating in galaxy clusters: Successes and challenges",
abstract = "We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealized, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated with multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at ∼50 kpc scales where it only represents a few per cent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multiphase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.",
keywords = "Galaxies: active, Galaxies: clusters: general, Galaxies: jets, Methods: numerical",
author = "Davide Martizzi and Eliot Quataert and Faucher-Gigu{\`e}re, {Claude Andr{\'e}} and Drummond Fielding",
year = "2019",
month = "1",
day = "1",
doi = "10.1093/mnras/sty3273",
language = "English (US)",
volume = "483",
pages = "2465--2486",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "2",

}

Simulations of jet heating in galaxy clusters : Successes and challenges. / Martizzi, Davide; Quataert, Eliot; Faucher-Giguère, Claude André; Fielding, Drummond.

In: Monthly Notices of the Royal Astronomical Society, Vol. 483, No. 2, 01.01.2019, p. 2465-2486.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Simulations of jet heating in galaxy clusters

T2 - Successes and challenges

AU - Martizzi, Davide

AU - Quataert, Eliot

AU - Faucher-Giguère, Claude André

AU - Fielding, Drummond

PY - 2019/1/1

Y1 - 2019/1/1

N2 - We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealized, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated with multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at ∼50 kpc scales where it only represents a few per cent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multiphase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.

AB - We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealized, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated with multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at ∼50 kpc scales where it only represents a few per cent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multiphase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.

KW - Galaxies: active

KW - Galaxies: clusters: general

KW - Galaxies: jets

KW - Methods: numerical

UR - http://www.scopus.com/inward/record.url?scp=85060552345&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85060552345&partnerID=8YFLogxK

U2 - 10.1093/mnras/sty3273

DO - 10.1093/mnras/sty3273

M3 - Article

AN - SCOPUS:85060552345

VL - 483

SP - 2465

EP - 2486

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 2

ER -