[논문 리뷰] Anisotropic thermal conduction in galaxy clusters with MHD in Gadget

Alexander Arth, Klaus Dolag외 3명|arXiv (Cornell University)|2014. 12. 19.

Galaxies: Formation, Evolution, Phenomena참고 문헌 5인용 수 14

한 줄 요약

이 논문은 GADGET 코드를 사용한 스무스드 파icle 유체역학(SPH) 시뮬레이션에 자기장에 따라 정렬된 열전도를 물리적으로 유도된 방식으로 구현한다. 이는 자기장에 따라 정렬된 열 이동을 고려한 이방성 열전도를 제안하며, 관측된 온도 변동과 일치하는 결과를 도출한다. 이는 효율 계수를 자유롭게 조정할 필요 없이 고효율 모델보다도 뛰어난 성능을 보이며, 자유 매개변수 없이 현실적인 모델링이 가능하다.

ABSTRACT

We present an implementation of thermal conduction including the anisotropic effects of magnetic fields for SPH. The anisotropic thermal conduction is mainly proceeding parallel to magnetic fields and suppressed perpendicular to the fields. We derive the SPH formalism for the anisotropic heat transport and solve the corresponding equation with an implicit conjugate gradient scheme. We discuss several issues of unphysical heat transport in the cases of extreme ansiotropies or unmagnetized regions and present possible numerical workarounds. We implement our algorithm into the GADGET code and study its behaviour in several test cases. In general, we reproduce the analytical solutions of our idealised test problems, and obtain good results in cosmological simulations of galaxy cluster formations. Within galaxy clusters, the anisotropic conduction produces a net heat transport similar to an isotropic Spitzer conduction model with an efficiency of one per cent. In contrast to isotropic conduction our new formalism allows small-scale structure in the temperature distribution to remain stable, because of their decoupling caused by magnetic field lines. Compared to observations, isotropic conduction with more than 10 per cent of the Spitzer value leads to an oversmoothed temperature distribution within clusters, while the results obtained with anisotropic thermal conduction reproduce the observed temperature fluctuations well. A proper treatment of heat transport is crucial especially in the outskirts of clusters and also in high density regions. It's connection to the local dynamical state of the cluster also might contribute to the observed bimodal distribution of cool core and non cool core clusters. Our new scheme significantly advances the modelling of thermal conduction in numerical simulations and overall gives better results compared to observations.

연구 동기 및 목표

• 전자 이동에 대한 자기장 영향을 고려하여 은하단에서 이방성 열전도를 모델링하기 위해.
• 극도로 이방성 또는 자기장이 없는 영역에서 비물리적인 열전도가 발생하는 문제를 해결하기 위해.
• 우주론적 SPH 시뮬레이션에 이방성 열전도를 구현하고 검증하기 위한 수치적 방법을 개발하기 위해.
• 이방성 열전도가 온도 구조, 자기장 및 은하단 관측량에 미치는 영향을 등방성 모델과 비교하기 위해.
• 이방성 열전도가 관측된 냉핵 및 비냉핵 은하단의 이중성 분포를 설명할 수 있는지 평가하기 위해.

제안 방법

• 자기장 선에 따라 정렬된 이방성 열전도를 위한 SPH 형식을 유도하며, 자기장에 따라 변하는 전도도를 가진 확산 방정식을 사용한다.
• 수치적 안정성을 확보하기 위해 음이 아닌 공액 기울기 방법을 사용해 유도된 열전도 방정식을 해석한다.
• 약한 또는 자기장이 없는 영역에서 비물리적인 열전도가 발생하는 것을 방지하기 위해 수치적 보완 조치를 적용한다.
• 우주론적 은하단 형성 시뮬레이션을 위해 GADGET-3 코드에 스무스드 파icle MHD(SPMHD)와 함께 알고리즘을 통합한다.
• 완전한 우주론적 시뮬레이션에 적용하기 전에 분석적 해와의 비교를 통해 3차원 이상의 이상적 시험 케이스를 사용해 알고리즘을 검증한다.
• 등방성 및 이방성 열전도 모델 간의 결과를 비교하고 효율 계수를 변화시켜 온도 프로파일과 발광도를 평가한다.

실험 결과

연구 질문

• RQ1 이방성 열전도는 우주론적 SPH 시뮬레이션에서 은하단의 온도 구조와 진화에 어떻게 영향을 미치는가?
• RQ2 이방성 열전도 모델은 등방성 모델보다 관측된 온도 변동을 더 잘 재현할 수 있는가?
• RQ3 수직 방향의 열전도는 은하간 매질의 역학적 특성과 형상에 어떤 영향을 미치는가?
• RQ4 이방성 열전도의 포함 여부는 냉핵 및 비냉핵 은하단의 형성과 특성에 어떤 영향을 미치는가?
• RQ5 이방성 열전도는 고밀도 영역에서 자기장 구조와 열에너지 분포에 얼마나 큰 영향을 미치는가?

주요 결과

• 이방성 열전도는 특히 스피츠터 열전도의 10% 이상의 효율을 가진 등방성 모델과 비교할 때 관측된 데이터와 매우 유사한 온도 변동을 생성한다.
• 이방성 열전도를 사용한 시뮬레이션은 스피츠터 값의 약 1%에 해당하는 등방성 열전도와 동일한 순열 열전도를 보이며, 미세 구조의 온도 구조를 유지한다.
• 10% 이상의 효율을 가진 등방성 열전도는 온도 프로파일을 과도하게 부드럽게 만들며, 관측 결과와 일치하지 않지만, 이방성 열전도는 현실적인 미세 구조를 유지한다.
• 특히 저자기장 또는 고밀도 영역에 있는 입자 중 상당수는 무시할 수 없는 수준의 수직 열전도를 보이며, 이 성분을 포함하는 것이 중요함을 시사한다.
• 이방성 열전도의 결과로 거대 은하단에서 반경 방향 자기장 및 온도 프로파일이 더 급격해지고, 평균 은하단 온도의 분산이 증가한다.
• 이방성 열전도 시뮬레이션의 발광도 분포는 프랭크 등(2013)의 관측 데이터와 가장 잘 일치하며, 특히 고효율 등방성 모델과 비교할 때 뚜렷한 우월성을 보인다.