[論文レビュー] 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

ひとこと要約

本研究では、AGB星のヘリウムシェルフラッシュ対流境界上部における水素の混合を3次元、完全圧縮流体として$4\pi$幾何構造で初めてシミュレートした。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星におけるヘリウムシェルフラッシュ対流境界上部での陽子を多く含む燃料の混合を調査すること。この現象は、サクラリの星のような後AGB星におけるHインジェストフラッシュに関連している。
• 標準の混合長モデルが不適切となる、深く硬い対流境界環境における定量的混合率を特定すること。
• 星内部の対流境界混合における3次元流体ダイナミクスシミュレーションの数値的信頼性を検証すること。特に分解能収束性と安定性に焦点を当てる。
• PPM+PPB多流体移流スキームが、圧縮性で乱流的な流れにおける鋭い流体界面を高分解能で正確かつ低拡散で追跡できることを確立すること。

提案手法

• ヘリウムシェルフラッシュ対流と上部境界での水素を多く含む燃料の混合をモデル化するため、$4\pi$球対称幾何構造における3次元、完全圧縮流体ダイナミクスシミュレーションを実施する。
• 各セル内で圧力と温度が平衡に達していると仮定し、各流体の密度を別々に保存しないことで計算コストを削減する。水素とヘリウムを追跡するための混合分率$f$を用いた多流体アプローチを採用する。
• 混合分率$f$のセル界面を越えた高次元で安定な移流を達成するため、PPM+PPB（区分的放物線法＋放物線的プロファイル制限）再構成技術を適用する。
• 個々の流体密度を明示的に保存しないで、混合密度と$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$から正確な密度再構成が可能となり、計算コストを著しく削減しながらも精度を維持した。