[논문 리뷰] Buoyant Bubbles in a Cooling Intracluster Medium I. Hydrodynamic Bubbles

A. Gardini|ArXiv.org|2006. 11. 14.

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.

연구 동기 및 목표

• 냉각 유동의 형성에 냉각된 내부군집 매질(ICM)에서 부력 있는 상대론적 플라즈마 기포의 영향을 조사하기 위해.
• 관측된 X선 공극으로부터 제안된 바와 같이, 순수 유체역학적 기포만으로도 냉각핵 군집의 냉각 유동을 억제할 수 있는지 판단하기 위해.
• 기포의 역학, 예를 들어 상승 및 소용돌이 고리 형성과 같은 요소들이 ICM 내 에너지와 운동량 재분포에 미치는 역할를 평가하기 위해.
• 다양한 해상도와 좌표계에 따른 시뮬레이션 결과의 탄력성과 순수 유체역학 모델의 한계를 평가하기 위해.
• 자기장, 점성, 열전도 및 AGN 활동 단계와 같은 추가 물리적 요소를 포함하는 향후 연구를 위한 기초를 마련하기 위해.

제안 방법

• 냉각 ICM 내에서 기포의 진화를 모델링하기 위해 $1024^3$ 개의 효과적 그리드 요소를 사용하는 고해상도 3차원 유체역학 시뮬레이션을 수행하기 위해.
• 압축성 오일러 방정식을 냉각 항과 상대론적 플라즈마 가정을 포함하여 해결하기 위해 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차원 결과의 신뢰성을 뒷받침한다.