[論文レビュー] Buoyant Bubbles in a Cooling Intracluster Medium I. Hydrodynamic Bubbles

A. Gardini|ArXiv.org|Nov 14, 2006

Galaxies: Formation, Evolution, Phenomena参考文献 101被引用数 11

ひとこと要約

本研究では、高解像度3次元流体力学的シミュレーションを用いて、相対論的プラズマバブル（ラジオローブとしてモデル化）がクラスターコアにおける冷却流れに与える影響を調査した。乱流運動を引き起こすが、冷却率の低下は最大で25%にとどまり、純粋な流体力学的メカニズムだけでは、クールコアクラスタにおける観測された冷却流れの抑制を説明できないことが示された。

ABSTRACT

Over the past several years, numerous examples of X-ray cavities coincident with radio sources have been observed in so-called "cool core" clusters of galaxies. Motivated by these observations, we explore the evolution and the effect of cavities on a cooling intracluster medium (ICM) numerically, adding relevant physics step by step. In this paper we present a first set of hydrodynamical, high resolution (1024^3 effective grid elements), three-dimensional simulations, together with two-dimensional test cases. The simulations follow the evolution of radio cavities, modeled as bubbles filled by relativistic plasma, in the cluster atmosphere while the ICM is subject to cooling. We find that the bubble rise retards the development of a cooling flow by inducing motions in the ICM which repeatedly displace the material in the core. Even bubbles initially set significantly far from the cluster center affect the cooling flow, although much later than the beginning of the simulation. The effect is, however, modest: the cooling time is increased by at most only 25%. As expected, the overall evolution of pure hydrodynamic bubbles is at odds with observations, showing that some additional physics has to be considered in order to match the data.

研究の動機と目的

• 冷却ガス中における冷却流れの発展に、浮力を持つ相対論的プラズマバブルが与える影響を調査すること。
• X線のキャビティとして観測されたように、純粋な流体力学的バブルがクールコアクラスタにおける冷却流れを抑制できるかどうかを特定すること。
• バブルのダイナミクス（上昇運動やヴォーテックスリングの形成など）が、ICM内でのエネルギーおよび運動量の再分配に果たす役割を評価すること。
• 解像度や座標系の違いに対して、シミュレーション結果の頑健性を評価し、純粋な流体力学モデルの限界を同定すること。
• 磁場、粘性、熱伝導、AGN活動の段階といった追加の物理的要因を組み込む今後の研究の基盤を確立すること。

提案手法

• 冷却ICM内でのバブル進化をモデル化するため、$1024^3$ 個の有効グリッド要素を用いた高解像度3次元流体力学的シミュレーションを実施。
• 圧縮性Euler方程式に冷却項と相対論的プラズマ仮定を組み込んだFLASHコードを用いて計算。
• ヴォーテックスリングの形成やコア物質の混合を焦点に、2次元のテストケースを実施し、3次元の結果を検証および拡張。
• 初期位置（例：中心からの距離が10、20、30 kpc）を変化させた複数のバブル配置をシミュレーションし、効果の空間的依存性を評価。
• グリッドに起因する数値不安定性を軽減し、結果の頑健性を向上させるために、初期条件にランダムノイズを導入。
• 解像度、座標系、初期摂動の違いに対して結果を比較し、感度を評価。

実験結果

リサーチクエスチョン

• RQ1 浮力を持つ相対論的プラズマバブルは、クラスターコアにおける冷却流れの形成をどの程度抑制できるか？
• RQ2 バブルの初期位置がクラスターコアからの距離に依存する場合、冷却流れの発生を遅らせたり、破壊したりする能力にどのような影響を与えるか？
• RQ3 ヴォーテックスリングや混合された物質は、ICM内でのエネルギーおよび運動量の再分配にどのような役割を果たすか？
• RQ4 グリッドに起因する数値アーティファクト（例：数値不安定性）は、バブル進化シミュレーションの信頼性にどのように影響するか？
• RQ5 なぜ純粋な流体力学的シミュレーションでは、観測された冷却流れの抑制を完全に再現できないのか？さらにどのような追加の物理的メカニズムが必要か？

主な発見

• 浮力を持つバブルはICM内に乱流運動を引き起こし、冷却の急激な進行を遅らせるが、元の状態と比較して最大で25%の遅延にとどまる。
• 中心からの距離が遠い（例：30 kpc）初期位置に置かれたバブルでもコアに影響を与えるが、長時間の遅れを経て、近いバブルに比べて影響は弱い。
• SB1011ケース（初期位置10 kpc）では、冷却時間が約25%遅延するが、最終的には依然として急激な冷却が発生する。
• SB1030ケース（初期位置30 kpc）では、降下する物質がコアを圧縮し、冷却を加速させ、$t=380$ Myrには急激な冷却が発生する。
• バブルのダイナミクスはヴォーテックスリングの形成に類似しており、バブルから上昇する細い物質の塊（混合物）が形成され、コアに質量を堆積させるが、コアの構造を破壊しない。
• 数値的結果は解像度や座標系に対して頑健であるが、2次元シミュレーションは丸め誤差に対してより敏感であり、3次元結果の信頼性を裏付ける。