[论文解读] On the structure of the turbulent interstellar atomic hydrogen. I- Physical characteristics

P. Hennebelle, E. Audit|CERN Bulletin|Dec 29, 2006

Astrophysics and Star Formation Studies参考文献 1被引用 91

一句话总结

本文通过高分辨率二维数值模拟研究了湍流热双稳态原子氢（HI），结果表明，暖中性介质（WNM）中的亚音速湍流可自然形成冷中性介质（CNM）团块，其质量函数谱符合𝒩(M) ∝ M⁻¹.⁷，与观测到的CO团块谱非常接近。结果表明，CNM结构源于亚音速、近等温流中的密度涨落，其速度功率谱由于涡旋分量占主导且团块受压强约束，符合Kolmogorov标度，与超音速等温湍流有本质不同。

ABSTRACT

{We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}

研究动机与目标

• 研究湍流热双稳态原子氢（HI）在空间尺度3–4个数量级范围内的统计特性。
• 理解在HI中观测到的小尺度CNM团块的起源与结构，特别是其质量分布和动力学特性。
• 确定亚音速、近等温WNM中的湍流如何导致CNM碎片的形成与稳定。
• 将模拟结果与观测特征（如微小原子结构TSAS和低柱密度云）进行比较。
• 评估数值分辨率和热传导在塑造流动统计特性中的作用。

提出的方法

• 对热双稳态HI流进行高分辨率2D数值模拟，空间尺度分辨率比以往研究提高10倍。
• 采用亚音速、近等温WNM模拟湍流，CNM通过热不稳定性与碎片化形成。
• 应用Press & Schechter形式化方法推导亚音速气体中密度涨落的CNM质量函数理论预测。
• 计算密度与速度场的功率谱，以分析湍流特性及能量在不同尺度上的分布。
• 系统性地改变数值分辨率与热导率，以评估数值收敛性与物理结果的鲁棒性。
• 对柱密度、质量与速度弥散分布进行统计分析，以与观测数据对比。

实验结果

研究问题

• RQ1 观测到的HI中CNM团块质量谱的起源是什么？其与CO团块质量谱相比有何异同？
• RQ2 在亚音速、热双稳态HI湍流中，速度谱与密度谱如何区别于超音速等温湍流？
• RQ3 为何该模型中CNM碎片在高湍流环境下仍能长期存在？压强平衡在此过程中起何作用？
• RQ4 动能在WNM与CNM之间的分布如何随空间尺度变化？
• RQ5 数值分辨率与热传导在多大程度上影响模拟HI结构的统计特性？

主要发现

• CNM质量谱符合𝒩(M) ∝ M⁻¹.⁷，与观测到的CO团块质量谱高度一致，表明二者可能起源于稀薄原子气体中的共同机制。
• 基于Press & Schechter形式化的理论模型预测：在2D中𝒩(M) ∝ M⁻⁵⁄³，在3D中𝒩(M) ∝ M⁻¹⁶⁄⁹，为观测到的质量函数提供了物理解释。
• 速度功率谱符合Kolmogorov预测（k⁻⁵⁄³），且主要由涡旋分量（不可压缩分量）主导，与超音速等温流中的行为截然不同。
• 在大于~1 pc的尺度上，动能主要由WNM贡献；但在小于~1 pc的尺度上，动能则主要来自CNM碎片的平动运动。
• 密度功率谱相对平坦，与超音速等温湍流一致，但结构受压强约束且寿命较长，与超音速流中由激波约束的结构有本质区别。
• 热传导对小尺度结构有轻微影响，但对大尺度流动的统计特性无显著影响。