[论文解读] 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

一句话总结

本研究使用高分辨率、理想化的流体动力学模拟，基于Athena代码，探究类Perseus星系团核心中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代码进行模拟，以最小化人工耗散，保持低数值扩散。
• 使用高分辨率结构化网格，以解析星系团核心及喷流驱动外流中的小尺度结构。
• 开发新型分析技术，分别测量冷却速率、总加热速率，以及弱激波和湍流耗散的贡献。
• 模拟在类Perseus星系团环境中的动能喷流，聚焦于受控的理想化设置下的反馈机制。
• 通过比较不同流体求解器实现方式的结果，评估对数值扩散和格式特性的敏感性。
• 在网格尺度上分析热不稳定性与多相气体形成，以评估其在加热与数值收敛性中的作用。

实验结果

研究问题

• RQ1 动能喷流能否通过非对称加热在星系团核心成功抑制冷却流？
• RQ2 弱激波、湍流耗散与直接喷流加热对总加热速率的相对贡献如何？
• RQ3 数值求解器的选择如何影响模拟结果，特别是在解析热不稳定性与多相气体方面？
• RQ4 湍流在不同半径处（特别是超过50 kpc时）在多大程度上是显著的加热源？
• RQ5 当前的分辨率与物理解释是否足以实现对喷流加热星系团核心的数值收敛？

主要发现

• 冷却流的抑制主要发生在喷流轴向，该区域加热具有显著的非对称性且最有效。
• 在约50 kpc半径处，弱激波对加热有显著贡献，是该尺度下最主要的加热机制。
• 湍流仅在星系团中心区域是主要加热源，在50 kpc处贡献不足总加热的5%。
• 数值求解器的特性显著影响模拟结果，表明由于在网格尺度上未能充分解析热不稳定性，导致收敛性差。
• 由热不稳定性引发的多相气体形成在网格尺度上发展，且对所选流体求解器高度敏感，限制了数值结果的可靠性。
• 喷流反馈的定性图像与以往研究一致，但定量结果依赖于数值实现方式，凸显了模拟该物理区时仍存在未解决的挑战。