[论文解读] A calibration of the mixing-length for solar-type stars based on hydrodynamical simulations - I. Methodical aspects and results for solar metallicity

H.‐G. Ludwig, B. Freytag+1 more|arXiv (Cornell University)|Nov 11, 1998

Stellar, planetary, and galactic studies参考文献 1被引用 35

一句话总结

本论文通过太阳型恒星的二维辐射流体动力学（RHD）模拟，标定对流表层的混合长度参数（α_MLT）。其得出的α_MLT值在F型矮星处约为1.3，在K型亚巨星处约为1.75，太阳附近呈平台状，且通过日震学数据验证了该方法，内部精度达到±0.05，并为恒星模型应用提供了拟合公式。

ABSTRACT

Based on detailed 2D numerical radiation hydrodynamics (RHD) calculations of time-dependent compressible convection, we have studied the dynamics and thermal structure of the convective surface layers of solar-type stars. The RHD models provide information about the convective efficiency in the superadiabatic region at the top of convective envelopes and predict the asymptotic value of the entropy of the deep, adiabatically stratified layers. This information is translated into an effective mixing-length parameter alpha suitable to construct standard stellar structure models. We validate the approach by a detailed comparison to helioseismic data. The grid of RHD models for solar metallicity comprises 58 simulation runs with a helium abundance of Y=0.28 in the range of effective temperatures 4300K < Teff < 7100K and gravities 2.54 < logg < 4.74. We find a moderate, nevertheless significant variation of alpha between about 1.3 for F-dwarfs and 1.75 for K-subgiants with a dominant dependence on Teff. In the close neighbourhood of the Sun we find a plateau where alpha remains almost constant. The internal accuracy of the calibration of alpha is estimated to be +/- 0.05 with a possible systematic bias towards lower values. An analogous calibration of the convection theory of Canuto & Mazzitelli gives a different temperature dependence but a similar variation of the free parameter. For the first time, values for the gravity-darkening exponent beta are derived independently of mixing-length theory: beta=0.07...0.10. We show that our findings are consistent with constraints from stellar stability considerations and provide compact fitting formulae for the calibrations.

研究动机与目标

• 通过直接的流体动力学模拟标定太阳型恒星的混合长度参数（α_MLT），摆脱启发式假设。
• 量化α_MLT在赫罗图上的变化，特别是在太阳附近的区域，以改进恒星结构建模。
• 通过与日震学数据对比，验证基于模拟的标定方法，确保与观测到的太阳结构一致。
• 独立于混合长度理论，推导重力变暗指数β，提供新的物理约束。
• 为标准恒星演化模型提供紧凑且可直接使用的α_MLT和α_CMT（Canuto-Mazzitelli理论）的拟合公式。

提出的方法

• 对具有太阳金属丰度（Y=0.28）的58个时间依赖、可压缩的二维辐射流体动力学（RHD）模拟进行太阳型恒星研究，覆盖4300 K ≤ T_eff ≤ 7100 K 和 2.54 ≤ log g ≤ 4.74。
• 从RHD模拟中提取深层绝热层的熵（s_env），利用状态方程关联热力学变量。
• 通过将模拟的对流效率与标准混合长度理论相匹配，将s_env转换为有效混合长度参数α_MLT。
• 对Canuto & Mazzitelli（1991, 1992）对流理论（CMT）应用相同方法，推导用于比较的α_CMT。
• 基于归一化的温度（T̃ = (T_eff - 5770)/1000）和重力（g̃ = log(g/27500)）构建辅助拟合函数，生成紧凑且数值可用的公式。
• 通过将预测的太阳结构（s_env, α_MLT）与日震学约束对比，验证标定结果，确保在小误差范围内一致。

实验结果

研究问题

• RQ1 对于具有太阳金属丰度的太阳型恒星，混合长度参数α_MLT在赫罗图上的变化规律如何？
• RQ2 RHD模拟能否准确再现日震学约束下的太阳结构，从而验证模拟方法的可靠性？
• RQ3 α_MLT对有效温度和表面重力的依赖关系如何？与标准混合长度理论的假设相比有何差异？
• RQ4 能否从RHD模拟中独立于混合长度理论推导出重力变暗指数β？
• RQ5 如何将RHD模拟结果紧凑地参数化，以用于标准恒星结构模型？

主要发现

• 混合长度参数α_MLT从F型矮星的约1.3变化至K型亚巨星的1.75，主要依赖于有效温度。
• 在太阳附近（T_eff ≈ 5770 K），α_MLT保持近乎恒定的平台状，值约为1.587，与日震学约束一致。
• α_MLT标定的内部精度估计为±0.05，可能存在略微偏低的系统性偏差。
• 独立推导出的重力变暗指数β为0.07至0.10，与恒星稳定性考虑一致。
• 对Canuto & Mazzitelli（CMT）对流理论的标定显示，其温度依赖关系不同，但自由参数的变化趋势与MLT相似。
• 通过归一化温度和重力变量，提供了α_MLT和α_CMT的紧凑拟合公式，可直接集成至恒星演化模型中。