[论文解读] The density distribution and physical origins of intermittency in supersonic, highly magnetised turbulence with diverse modes of driving

James R. Beattie, Philip Mocz+2 more|arXiv (Cornell University)|Sep 22, 2021

Solar and Space Plasma Dynamics被引用 5

一句话总结

本研究通过36组不同驱动模式的MHD模拟，研究了高磁化超音速湍流中的间歇性。结果表明，非高斯密度涨落——在流体湍流中典型存在——在高马赫数（ℳ ≳ 4）时消失，转而呈现高斯统计特性，与流体动力学预期相反。间歇性的主要驱动机制是与大尺度磁场B₀对齐的单个强激波，而磁压通过自相似性和时序均质化抑制了高ℳ下的极端密度涨落。

ABSTRACT

The probability density function (PDF) of the logarithmic density contrast, $s=\ln (ρ/ρ_0)$, with gas density $ρ$ and mean density $ρ_0$, for hydrodynamical supersonic turbulence is well-known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, $\mathcal{M}$. By studying the mass- and volume-weighted $s$-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of $s$. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetised plasma for $\mathcal{M} \lesssim 4$. However, for $\mathcal{M} \gtrsim 4$, the non-Gaussian signatures disappear, leaving approximately Gaussian $s$-statistics -- exactly the opposite of hydrodynamical turbulence in the high-$\mathcal{M}$ limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the $\mathcal{M} \lesssim 4$ non-Gaussian features we use one-dimensional (1D) pencil beams to explore the dynamics along and across the large-scale magnetic field, $\mathbf{B}_0$. We discuss kinetic, density and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along $\mathbf{B}_0$. We discuss the Gaussianisation of the $\mathcal{M} \gtrsim 4$ $s$-fields through the lens of two phenomenologies: the self-similarity of the $s$-field and homogenisation of the dynamical timescales between the over- and under-dense regions in the compressible gas.

研究动机与目标

• 理解与星际介质相关的高磁化、超音速湍流中非高斯密度涨落（间歇性）的物理起源。
• 研究不同驱动模式（压缩型与涡旋型）对对数密度对比度s = ln(ρ/ρ₀)统计特性的影响。
• 确定为何在高湍流马赫数（ℳ ≳ 4）时s-PDF近似呈高斯分布，与纯流体动力学湍流行为相反。
• 识别主导非高斯密度PDF特征的物理机制，如激波与波模式。
• 探讨磁场通过自相似性和动力学时序均质化抑制极端密度涨落的作用。

提出的方法

• 对超音速湍流开展36组三维、等温、平均场MHD模拟，覆盖从亚阿夫达马赫数到跨阿夫达马赫数的湍流马赫数（ℳ）和阿夫达马赫数（ℳ_A）范围。
• 通过驱动参数b控制，采用不同驱动模式（压缩型或涡旋型）注入压缩或涡旋强迫，系统研究其对s-PDF的影响。
• 计算质量加权与体积加权的s-PDF，并应用四种独立的非高斯性度量（如峰度、偏度及高阶累积量）量化间歇性。
• 采用与大尺度磁场B₀平行和垂直的一维铅笔状光束，分离并分析局部动力学，包括激波、速度与磁场涨落。
• 开展一维数值激波实验，分离单个激波与多激波相互作用对非高斯s-PDF组分的贡献。
• 评估两种现象学解释在高ℳ下实现高斯化的合理性：s场的自相似性，以及稠密与稀疏区域之间动力学时序的均质化。

实验结果

研究问题

• RQ1 湍流马赫数ℳ如何影响高磁化、超音速湍流中s-PDF的非高斯性？
• RQ2 哪些物理机制——如激波、波或磁场结构——主导了间歇性非高斯密度涨落的生成？
• RQ3 为何在高ℳ（ℳ ≳ 4）时s-PDF近似呈高斯分布，与纯流体动力学湍流行为相反？
• RQ4 不同驱动模式（压缩型与涡旋型）如何影响s-PDF中非高斯特征的振幅与结构？
• RQ5 磁场涨落与波模式（如快磁声波、剪切阿夫达波）在间歇性或极端密度涨落抑制中起到多大程度的贡献？

主要发现

• 当ℳ ≲ 4时，s-PDF表现出显著的非高斯特征，间歇性主要由与大尺度磁场B₀对齐的单个强激波驱动。
• 当ℳ ≳ 4时，s-PDF近似呈高斯分布，逆转了流体动力学湍流中的预期趋势，表明磁场抑制了极端密度涨落。
• s-PDF的非高斯组分随压缩型驱动模式增强（b值升高）而单调增加，证实压缩强迫增强了间歇性。
• 一维铅笔状光束分析表明，沿B₀方向的激波产生体积稀少、稀薄的s-结构（1.5–2σ速度事件），是导致非高斯性的主要贡献者。
• 一维激波实验确认，单个激波而非多激波相互作用，是模拟中非高斯s-PDF组分的主要来源。
• 高ℳ下实现高斯化的最佳解释是s场的自相似性，以及稠密与稀薄区域之间动力学时序的均质化，其中磁压在抑制极端涨落中起关键作用。