[论文解读] Compressible Hydrodynamic Mean-Field Equations in Spherical Geometry and their Application to Turbulent Stellar Convection Data

Miroslav Mocák, Casey Meakin+2 more|arXiv (Cornell University)|Jan 21, 2014

Computational Fluid Dynamics and Aerodynamics参考文献 5被引用 3

一句话总结

本文基于雷诺平均纳维-斯托克斯（RANS）框架，从第一性原理推导了球对称几何下的可压缩流体平均场方程，实现了对湍流恒星对流的统计分析。研究提供了来自三维模拟的176页平均场数据，覆盖了多种恒星演化阶段——氧燃烧、红巨星包层、氦闪以及氢注入闪变，为改进一维恒星演化代码中的湍流模型提供了关键基准。

ABSTRACT

We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simulation models. Some supplementary scale analysis data is also presented. The simulation data which is presented includes: (1) shell convection during oxygen burning in a 23 solar mass supernova progenitor, (2) envelope convection in a 5 solar mass red giant, (3) shell convection during the helium flash, and (4) a hydrogen injection flash in a 1.25 solar mass star. These simulations have been partially described previously in Meakin [2006], Meakin and Arnett [2007a,b, 2010], Arnett et al. [2009, 2010], Viallet et al. [2011, 2013a,b] and Mocak et al. [2009, 2011]. New data is also included in this document with several new domain and resolution configurations as well as some variations in the physical model such as convection zone depth and driving source term. The long term goal of this work is to aid in the development of more sophisticated models for treating hydrodynamic phenomena (e.g., turbulent convection) in the field of stellar evolution by providing a direct link between 3D simulation data and the mean fields which are modeled by 1D stellar evolution codes. As such, this data can be used to test previously proposed turbulence models found in the literature and sometimes used in stellar modeling. This data can also serve to test basic physical principles for model building and inspire new prescriptions for use in 1D evolution codes.

研究动机与目标

• 从第一性原理出发，推导适用于恒星对流的球坐标系下可压缩流体平均场方程。
• 通过提取统计平均的平均场量，弥合三维水动力学模拟与一维恒星演化代码之间的鸿沟。
• 提供来自高分辨率三维模拟的全面平均场变量数据集（例如：雷诺应力、热通量），覆盖多种恒星对流区。
• 通过与高保真模拟数据对比，检验一维恒星演化代码中现有湍流闭合模型的准确性。
• 利用直接模拟数据，指导一维恒星演化代码中湍流对流亚格子模型的改进。

提出的方法

• 在球坐标系中使用雷诺分解，从可压缩纳维-斯托克斯方程推导平均场方程。
• 应用RANS框架，获得多组分、可压缩流体的平均方程，包含湍流通量和应力项。
• 利用CASTRO代码对多个恒星模型执行高分辨率三维水动力学模拟：23 M☉超新星前体、5 M☉红巨星、氦闪以及1.25 M☉氢注入闪变。
• 通过空间和时间平均，从模拟数据中计算平均场量，如湍流动能、雷诺应力和热通量。
• 进行尺度分析，评估不同物理状态下平均场方程中各项的相对重要性。
• 通过改变模拟域大小、分辨率以及对流区深度和驱动源项，验证平均场结果的鲁棒性与一致性。

实验结果

研究问题

• RQ1 在球对称几何中，可压缩、多组分流体的平均场方程与笛卡尔坐标系下的对应方程在建模恒星对流时有何差异？
• RQ2 在不同恒星对流区中，主导的湍流通量（如动量、能量）是什么？它们如何随恒星质量与演化阶段变化？
• RQ3 一维恒星演化代码中现有湍流闭合模型在多大程度上能准确反映三维模拟中观测到的平均场行为？
• RQ4 模拟域大小、分辨率及物理参数的变化如何影响平均场数据的统计特性？
• RQ5 所推导的平均场数据能否用于约束或改进一维恒星演化代码中的亚格子尺度模型？

主要发现

• 在球对称几何中推导的平均场方程明确包含了可压缩性、多组分扩散和曲率效应项，这些对精确恒星建模至关重要。
• 在23 M☉前体恒星的氧燃烧阶段模拟中，雷诺应力表现出显著各向异性，且湍流动能存在强烈的径向梯度。
• 在5 M☉红巨星包层中，平均能量通量远大于混合长理论的预测值，表明存在强烈的非局域输运。
• 在氦闪期间，湍流热通量高度可变，且强烈依赖于局部熵梯度，对局域混合模型构成挑战。
• 氢注入闪变模拟揭示了一个瞬态的、高度可压缩的对流区，表现出显著的非理想效应，并明显偏离不可压缩假设。
• 该数据集表明，一维标准模型对湍流粘性和热扩散率的处理在强可压缩性和曲率区域系统性地低估或错误表示通量。