[論文レビュー] Simulations of Jet Heating in Galaxy Clusters: Successes and Numerical Challenges

Davidé Martizzi, Eliot Quataert+2 more|arXiv (Cornell University)|May 16, 2018

Astrophysics and Star Formation Studies参考文献 2被引用数 5

ひとこと要約

本研究では、アテナコードを用いた高解像度で理想化された流体力学的シミュレーションを用い、ペルセウス型銀河団コアにおけるAGNジェットが冷却フローをどのように制御するかを調査した。ジェットによる非等方的加熱が、冷却を効果的に抑制することを明らかにした。特に、大きな半径では弱い衝撃波が加熱を支配し、中心部近くでは乱流が寄与するが、数値ソルバーの選択が結果に強く影響する。これは、熱不安定性に起因する多相性ガスを正しく捉えるために解像度が不十分であるため、収束性が悪いことが要因である。

ABSTRACT

We study idealised hydrodynamical simulations of Perseus-like galaxy cluster cores, with the goal of studying cooling flow regulation by jets generated by active galactic nuclei (AGN). The simulations are 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 total heating rate, the heating rate from weak shocks and the heating rate from the decay of turbulent motions excited by the jet. We find that successful regulation of the cooling flow can be achieved with kinetic jets. Heating is anisotropic and is mostly distributed along the jet axis where the cooling flow is suppressed. Away from the jet axis the cooling flow is reduced by weak shock heating, but not fully suppressed. Turbulence is a significant source of heating only near the cluster centre, but this mechanism becomes inefficient at $\sim$50 kpc scales where it only represents a few % of the total heating rate and where most of the heating is provided by weak shocks. This qualitative picture is in agreement with previous numerical work. However, we find that several details of the simulations depend on the choice made for the hydro solver. This is 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.

研究の動機と目的

• 理想化された流体力学的シミュレーションを用いて、AGNジェットが銀河団コアにおける冷却フローをどのように制御するかを調査すること。
• 冷却を抑制するための加熱メカニズムの寄与度を定量化すること——弱い衝撃波、乱流の減衰、直接的なジェット加熱。
• 数値ソルバーの特性がシミュレーション結果に与える影響を評価すること、特に熱不安定性および多相性ガス形成に関連して。
• ジェット加熱銀河団コアのモデリングにおける数値収束の課題を評価すること。

提案手法

• 低数値拡산を実現するため、EulerianのGodunovに基づくコードAthenaを用いてシミュレーションを実施した。
• 小スケール構造およびジェット駆動型流れを解像するため、高解像度で構造化されたグリッドを採用した。
• 冷却率、総加熱率、弱い衝撃波および乱流の減衰からの寄与を別々に測定するための新規な解析手法を開発した。
• 制御された理想化された設定下で、ペルセウス型の銀河団環境における運動的ジェットのフィードバック機構に焦点を当てた。
• 数値拡散およびスキーム特性への感受性を評価するために、異なる流体ソルバー実装の結果を比較した。
• 熱不安定性および多相性ガス形成はグリッドスケールで分析し、それらが加熱および数値収束性に果たす役割を評価した。

実験結果

リサーチクエスチョン

• RQ1 運動的ジェットは、非等方的加熱によって銀河団コアにおける冷却フローを効果的に抑制できるか？
• RQ2 弱い衝撃波、乱流の減衰、直接的ジェット加熱が、総加熱率に占める相対的寄与度は何か？
• RQ3 数値ソルバーの選択が、特に熱不安定性および多相性ガスの解像度に与える影響は何か？
• RQ4 半径が異なる位置、特に50 kpcを超えた領域では、乱流がどの程度の加熱源として重要となるか？
• RQ5 現在の解像度および物理的実装では、ジェット加熱銀河団コアのモデリングにおいて適切な数値収束が達成可能か？

主な発見

• 冷却フローの抑制は、主にジェット軸に沿って発生し、その領域での加熱は非等方的かつ最も効果的である。
• 弱い衝撃波は、約50 kpcの半径で加熱に顕著に寄与し、これらのスケールでの加熱の大部分を占める。
• 乱流は中心部の銀河団領域でのみ主要な加熱源であり、50 kpcでは総加熱の5%未満に寄与する。
• 数値ソルバーの特性がシミュレーション結果に強く影響しており、グリッドスケールでの熱不安定性の解像度が不十分であるため、収束性が悪いことが示唆された。
• 熱不安定性に起因する多相性ガスの形成はグリッドスケールで発生し、選択された流体ソルバーに強く依存するため、数値的信頼性が制限される。
• ジェットフィードバックの定性的な画像は、先行研究と一致するが、定量的結果は数値実装に依存しており、この分野における未解決の課題が浮き彫りになった。