[论文解读] Anisotropic thermal conduction in galaxy clusters with MHD in Gadget

Alexander Arth, Klaus Dolag+3 more|arXiv (Cornell University)|Dec 19, 2014

Galaxies: Formation, Evolution, Phenomena参考文献 5被引用 14

一句话总结

本文在GADGET代码的光滑粒子流体动力学（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模拟中出现的非物理解热输运问题。
• 在包含磁流体动力学（MHD）的宇宙学SPH模拟中，实现并验证一种各向异性导热的数值方案。
• 与各向同性模型相比，评估各向异性导热对温度结构、磁场和星系团可观测量的影响。
• 确定各向异性导热是否能够解释观测到的冷核与非冷核星系团的双峰分布。

提出的方法

• 推导出一种与磁场线对齐的各向异性热传导SPH形式，采用依赖磁场的电导率的扩散方程。
• 使用隐式共轭梯度法求解所得的热传导方程，以确保数值稳定性。
• 对弱磁场或零磁场区域应用数值修正，以防止出现非物理解热传递。
• 将该算法集成到GADGET-3代码中，结合光滑粒子MHD（SPMHD），用于模拟宇宙学尺度星系团的形成。
• 使用理想化的三维测试案例，在应用于完整宇宙学模拟前，将该方案与解析解进行验证。
• 对比各向同性与各向异性导热模型的结果，调整效率因子，评估温度剖面与发射率。

实验结果

研究问题

• RQ1 在宇宙学SPH模拟中，各向异性热传导如何影响星系团的温度结构和演化？
• RQ2 与各向同性导热模型相比，各向异性导热是否能更好地再现星系团中观测到的温度涨落？
• RQ3 垂直方向热传导对星系团内介质的动力学和形态有何影响？
• RQ4 在模拟中引入各向异性导热如何影响冷核与非冷核星系团的形成与特性？
• RQ5 各向异性导热在多大程度上改变了高密度区域的磁场结构与热能分布？

主要发现

• 各向异性导热产生的温度涨落与观测数据高度一致，尤其在与效率超过Spitzer导热10%的各向同性模型对比时表现更优。
• 采用各向异性导热的模拟中，净热传导量仅相当于Spitzer值约1%的各向同性导热，同时保持了小尺度温度结构。
• 效率超过10%的各向同性导热会导致温度剖面过度平滑，与观测不一致，而各向异性导热则保持了真实的亚结构。
• 大量粒子（尤其在低B场或高密度区域）表现出不可忽略的垂直方向热传导，凸显了包含该分量的重要性。
• 各向异性导热使大质量星系团的径向磁场和温度剖线更陡峭，并增加了星系团平均温度的方差。
• 与各向同性模型相比，各向异性导热模拟的发射率分布最符合Frank等人（2013年）的观测数据，尤其在高效率情况下。