[논문 리뷰] Simulations of Jet Heating in Galaxy Clusters: Successes and Numerical Challenges

Davidé Martizzi, Eliot Quataert외 2명|arXiv (Cornell University)|2018. 05. 16.

Astrophysics and Star Formation Studies참고 문헌 2인용 수 5

한 줄 요약

이 연구는 페르세우스 유사 은하단 중심부의 냉각 유량을 조절하는 데 있어 AGN 제트의 역할을 고해상도, 이де알라이즈드 유체역학 시뮬레이션을 통해 탐구한다. Athena 코드를 사용한 결과, 제트가 제트 축을 따라 비등방적으로 기체를 가열함으로써 냉각을 효과적으로 억제함을 확인하였다. 특히, 더 큰 반경에서는 약한 충격파가 주요 가열 기여를 하였고, 난류는 중심부 근처에서만 중요한 기여를 하였다. 그러나 수치 해법의 선택이 결과에 큰 영향을 미치며, 열 불안정성에 의해 유도되는 다상태 기체를 포착하는 데 있어 수렴성이 열악한 편이었다.

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.

연구 동기 및 목표

• 이dealized 유체역학 시뮬레이션을 통해 AGN 제트가 은하단 중심부의 냉각 유량을 어떻게 조절하는지 탐구하기 위해.
• 냉각을 억제하는 데 기여하는 다양한 가열 메커니즘—약한 충격파, 난류 붕괴, 직접 제트 가열—의 기여도를 정량화하기 위해.
• 수치 해법의 특성이 시뮬레이션 결과에 미치는 영향, 특히 열 불안정성과 다상태 기체 형성에 대해 평가하기 위해.
• 냉각 유량과 제트 피드백의 복잡한 물리적 특성을 포착하는 데 있어 수치적 수렴 문제를 평가하기 위해.

제안 방법

• 저항성의 영향을 최소화하기 위해 낮은 수치적 확산을 가진 에ulerian Godunov 기반 코드인 Athena를 사용하여 시뮬레이션을 수행하였다.
• 은하단 중심부의 미세 구조와 제트에 의해 유도된 유출을 해상도 높은 정렬된 격자로 해상도를 확보하였다.
• 냉각 속도, 총 가열 속도, 그리고 약한 충격파와 난류 붕괴의 기여도를 별도로 측정하기 위해 새로운 분석 기법을 개발하였다.
• 통제된 이dealized 설정에서 피드백 메커니즘을 집중적으로 분석하기 위해 페르세우스 유사 은하단 환경에서의 운동적 제트를 모델링하였다.
• 수치적 확산과 해법 특성에 대한 민감도를 평가하기 위해 다양한 유체 해법 구현 방식 간의 결과를 비교하였다.
• 열 불안정성과 다상태 기체 형성의 역할을 평가하기 위해 격자 스케일에서 분석하였다.

실험 결과

연구 질문

• RQ1 운동적 제트가 비등방성 가열을 통해 은하단 중심부의 냉각 유량을 성공적으로 억제할 수 있는가?
• RQ2 약한 충격파, 난류 붕괴, 직접 제트 가열의 기여도 중 어느 것이 총 가열 속도에 가장 큰 기여를 하는가?
• RQ3 수치 해법의 선택이 시뮬레이션 결과에 미치는 영향은 어떻게 되며, 특히 열 불안정성과 다상태 기체를 해결하는 데 있어 어떤가?
• RQ4 특히 50 kpc 이상의 반경에서 난류가 중요한 가열 기여원으로 작용하는 정도는 어느 정도인가?
• RQ5 현재의 해상도와 물리적 구현 방식은 제트로 가열된 은하단 중심부를 모델링하는 데 있어 적절한 수치적 수렴성을 보장하는가?

주요 결과

• 냉각 유량 억제 효과는 주로 제트 축을 따라 발생하며, 이는 비등방성 가열이 가장 효과적으로 작용하는 영역이다.
• 반경 약 50 kpc에서는 약한 충격파가 가열에 중요한 기여를 하며, 이 스케일에서 대부분의 가열 기여를 차지한다.
• 난류는 중심 은하단 영역 내에서만 주요 가열 기여원이며, 50 kpc에서는 총 가열의 5% 미만 기여한다.
• 수치 해법의 특성이 시뮬레이션 결과에 큰 영향을 미치며, 격자 스케일에서 열 불안정성을 제대로 포착하지 못함으로써 수렴성이 열악함을 시사한다.
• 열 불안정성에 의해 유도되는 다상태 기체 형성은 격자 스케일에서 발생하며, 선택된 유체 해법에 매우 민감하여 수치적 신뢰성이 제한된다.
• 제트 피드백의 정성적 그림은 이전 연구와 일치하지만, 정량적 결과는 수치적 구현 방식에 따라 달라지며, 이 영역을 시뮬레이션하는 데 있어 해결되지 않은 과제가 여전히 존재함을 시사한다.