[论文解读] Hydrodynamic simulations of H entrainment at the top of He-shell flash convection

Paul R. Woodward, Falk Herwig+1 more|arXiv (Cornell University)|Jul 15, 2013

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

一句话总结

本研究首次在 $4\pi$ 球对称几何下，对渐近巨星支（AGB）恒星中氦壳层闪光对流区顶部的氢包覆进行三维全压缩流体动力学模拟，采用先进的多流体迎风格式结合PPM+PPB重构方法，实现了高分辨率、稳定的模拟。关键结果是在 $1536^3$ 分辨率下获得收敛的包覆速率为 $4.38 \pm 1.48 \times 10^{-13}\,M_{\odot}/\text{s}$，证明了该方法在恒星内部强对流边界模拟中的数值保真度。

ABSTRACT

We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4π$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 \pm 1.48 imes 10^{-13}M_\odot \mathrm{/s}$.

研究动机与目标

• 研究渐近巨星支（AGB）恒星中氦壳层闪光对流区顶部质子富集燃料的包覆过程，该过程与萨库赖星（Sakurai’s Object）等后-AGB恒星中的氢吸入闪光密切相关。
• 在标准混合长度模型失效的深层、强刚性对流边界环境中，定量测定包覆速率。
• 验证三维流体动力学模拟在恒星内部对流边界混合过程中的数值保真度，重点关注分辨率收敛性与稳定性。
• 证实PPM+PPB多流体迎风格式能够实现对可压缩、湍流流动中锐利流体界面的精确、低扩散追踪。

提出的方法

• 在 $4\pi$ 球对称几何下开展三维全压缩流体动力学模拟，以模拟氦壳层闪光对流区在顶部边界处氢富集燃料的包覆过程。
• 采用多流体方法，引入混合分数 $f$ 来追踪氢与氦，假设每个网格单元内压力与温度保持平衡，从而避免为每种流体单独存储密度。
• 应用PPM+PPB（分段抛物线法结合抛物线轮廓限制）重构技术，实现对单元界面间混合分数 $f$ 的高阶、稳定迎风格式。
• 通过利用平均分子量比值，将流体密度从混合密度和 $f$ 插值重构，严格保证质量守恒，避免显式存储各流体的独立密度。
• 在多个分辨率（$1024^3$ 至 $1536^3$）下进行模拟，以评估包覆速率的收敛性，并确定数值保真度的分辨率阈值。
• 以高分辨率 $1536^3$ 模拟作为主要结果，基于网格细化和包覆速率稳定性进行收敛性分析。

实验结果

研究问题

• RQ1 氢富集物质进入氦壳层闪光对流区顶部边界的包覆速率是多少？其数值结果如何依赖于模拟分辨率？
• RQ2 大尺度对流运动与剪切流如何影响强刚性、深层恒星对流区中对流边界的稳定性和动力学行为？
• RQ3 PPM+PPB多流体迎风格式是否能够在高分辨率下实现对可压缩、湍流恒星对流中锐利流体界面的稳定、低扩散模拟？
• RQ4 在厚壳层对流中，全局或近全局尺度的对流结构在多大程度上主导了包覆过程？
• RQ5 随着分辨率提高，数值模拟是否收敛于一个稳定的包覆速率？达到定量可靠性所需的分辨率是多少？

主要发现

• 基于 $1536^3$ 分辨率模拟，包覆速率被定量确定为 $4.38 \pm 1.48 \times 10^{-13}\,M_{\odot}/\text{s}$，该结果显示出最清晰的收敛迹象。
• 分辨率为 $1024^3$ 及以下的模拟无法达到足够的数值保真度，仅 $1536^3$ 分辨率能提供可靠的包覆速率估计。
• 相干的、近全局尺度的对流运动主导了流场，且在无边界层分离的情况下，上边界剪切流对开尔文-赫姆霍兹不稳定性具有稳定性。
• 大尺度对流中反向流的边界层分离会诱发开尔文-赫姆霍兹不稳定性，从而将全局对流动力学与包覆过程联系起来。
• PPM+PPB方法可在仅2–3个网格单元厚度的界面处实现平滑、高分辨率的混合分数 $f$ 跟踪，有效避免数值振荡与扩散。
• 局部压力与温度平衡的假设使得仅通过混合密度和 $f$ 即可准确重构密度，显著降低计算成本，同时保持高精度。