[论文解读] The Stagger-grid: A grid of 3D stellar atmosphere models - III. The relation to mixing-length convection theory

Z. Magic, Achim Weiß+1 more|Mar 5, 2014

Stellar, planetary, and galactic studies参考文献 30被引用 3

一句话总结

本文利用 Stagger-grid 的 3D 辐射流体动力学（RHD）模拟，校准了 1D 恒星大气模型中的混合长度参数（αMLT），表明 αMLT 随恒星参数系统性变化——有效温度越高、表面重力越低、金属丰度越高，αMLT 越小。关键结果是提出一个基于物理的、可变的 αMLT（1.7–2.4），替代原有的恒星标定常数值，从而提升恒星模型的准确性。

ABSTRACT

We investigate the relation between 1D atmosphere models that rely on the mixing length theory and models based on full 3D radiative hydrodynamic (RHD) calculations to describe convection in the envelopes of late-type stars. The adiabatic entropy value of the deep convection zone, s_bot, and the entropy jump, Δs, determined from the 3D RHD models, are matched with the mixing length parameter, α_MLT, from 1D hydrostatic atmosphere models with identical microphysics (opacities and equation-of-state). We also derive the mass mixing length, α_m, and the vertical correlation length of the vertical velocity, C[v_z,v_z], directly from the 3D hydrodynamical simulations of stellar subsurface convection. The calibrated mixing length parameter for the Sun is α_MLT (s_bot) = 1.98. For different stellar parameters, α_MLT varies systematically in the range of 1.7 - 2.4. In particular, α_MLT decreases towards higher effective temperature, lower surface gravity and higher metallicity. We find equivalent results for α_MLT (Δs). Also, we find a tight correlation between the mixing length parameter and the inverse entropy jump. We derive an analytical expression from the hydrodynamic mean field equations that motivates the relation to the mass mixing length, α_m, and find that it exhibits qualitatively a similar variation with stellar parameter (between 1.6 and 2.4) with a solar value of α_m = 1.83. The vertical correlation length scaled with the pressure scale height yields for the Sun 1.71, but displays only a small systematic variation with stellar parameters, the correlation length slightly increasing with Teff. We derive mixing length parameters for various stellar parameters that can be used to replace a constant value. Within any convective envelope, α_m and related quantities vary a lot.

研究动机与目标

• 将 1D 恒星大气模型中使用的恒定混合长度参数 αMLT 替换为基于 3D RHD 模拟导出的、具有物理解释的可变值。
• 研究 αMLT 如何依赖于有效温度、表面重力和金属丰度等恒星参数。
• 建立混合长度参数与流体动力学量（如质量混合长度参数 αm 和垂直速度关联长度）之间的物理解释联系。
• 通过使用相同微观物理输入匹配 3D RHD 结果（s_bot 和 Δs）与 1D MLT 模型，验证校准结果的一致性。
• 从平均场流体动力学方程推导出解析表达式，解释 αMLT 随恒星参数变化的物理机制。

提出的方法

• 利用 Stagger-grid 的 3D 辐射流体动力学（RHD）模拟，计算恒星对流区的等熵熵（s_bot）和熵跃迁（Δs）。
• 通过调整 αMLT 使 1D 静力大气模型与 3D RHD 模拟在相同辐射截面和状态方程下匹配 s_bot 或 Δs，从而实现 αMLT 的校准。
• 直接从 3D RHD 模拟中提取质量混合长度参数 αm 和垂直速度关联长度 C[vz,vz]，以关联到物理对流特性。
• 基于平均场流体动力学方程推导出解析表达式，将 αMLT 与 αm 和熵跃迁的倒数关联，提供理论依据。
• 在对流区采用深度依赖的截断处理，以避免在底部边界附近计算平均关联长度时引入噪声。
• 系统性地改变恒星参数（T_eff, log g, [Fe/H]），在赫罗图上绘制 αMLT 及相关量的分布。

实验结果

研究问题

• RQ1 混合长度参数 αMLT 如何随有效温度、表面重力和金属丰度等恒星参数变化？
• RQ2 混合长度参数 αMLT 与流体动力学量（如质量混合长度参数 αm 和垂直速度关联长度）之间存在何种物理解释联系？
• RQ3 能否使用来自 3D RHD 模拟的可变 αMLT 替代 1D 模型中恒定的太阳标定值，同时保持一致性？
• RQ4 αMLT 与熵跃迁 Δs 的倒数之间是否存在系统性相关性？该相关性能否由平均场理论解释？
• RQ5 垂直速度的垂直关联长度与 αMLT 和 αm 相比如何？其随恒星参数的变化是否具有系统性？

主要发现

• 太阳的校准混合长度参数为 αMLT^☉(s_bot) = 1.98，αMLT 在不同恒星参数下变化范围为 1.7–2.4。
• αMLT 随有效温度升高而减小，随表面重力降低而减小，随金属丰度升高而减小，表现出对恒星参数的系统性依赖。
• 发现 αMLT 与熵跃迁的对数倒数之间存在紧密相关性，满足 αMLT ∼ −lnΔs，表明存在强烈的物理解释关联。
• 质量混合长度参数 αm 的取值范围为 1.6 至 2.4，太阳值为 αm^☉ = 1.83，其变化趋势与 αMLT 定性相似。
• 以压力标高归一化的垂直速度关联长度为 1.71（太阳），随有效温度略有增加，与 αMLT 的变化趋势相反。
• 其他变量（熵、温度、压力）的关联长度较低（约 1.3），而密度的关联长度接近于 1，表明不同热力学变量的标度行为存在显著差异。