[論文レビュー] Transport and mixing in the radiation zones of rotating stars: I-Hydrodynamical processes

S. Mathis, Jean-Paul Zahn|ArXiv.org|Jun 18, 2004

Stellar, planetary, and galactic studies参考文献 20被引用数 104

ひとこと要約

本稿は、回転混合のモデル化を改善するもので、流体力学的取り扱いをより高次の球面調和関数に拡張し、時間微分を完全に保持することで、タコクリン循環や急速な進化段階の正確なシミュレーションを可能にする。従来のモデルの限界を克服し、特に強い剪切駆動乱流と緯度差別的回転を示す質量の大きな星において、非定常的ダイナミクスと緯度差別的回転を捉えることができる。

ABSTRACT

The purpose of this paper is to improve the modelization of the rotational mixing which occurs in stellar radiation zones, through the combined action of the thermally driven meridional circulation and of the turbulence generated by the shear of differential rotation. The turbulence is assumed to be anisotropic, due to the stratification, with stronger transport in the horizontal directions than in the vertical. The main difference with the former treatments by Zahn (1992) and Maeder & Zahn (1998) is that we expand here the departures from spherical symmetry to higher order, and include explicitly the differential rotation in latitude, to first order. This allows us to treat simultaneously the bulk of a radiation zone and its tachocline(s). Moreover, we take fully into account the non-stationarity of the problem, which will enable us to tackle the rapid phases of evolution. The system of partial differential equations, which govern the transport of angular momentum, heat and chemical elements, is written in a form which makes it ready to implement in a stellar evolution code. Here the effect of a magnetic field is deliberately ignored; it will be included in forthcoming papers.

研究の動機と目的

• 標準的な1.5次元近似を超えて、球対称性からのずれの取り扱いを拡張することで、星の放射層における回転混合のモデル化を改善すること。
• 線形近似において緯度差別的回転の効果を明示的に取り入れ、タコクリン循環の取り扱いを可能にすること。
• 動的緩和を除くすべての時間微分を保持することで、回転と組成の勾配が急峻な急速な進化段階をよりよく解像できること。
• 角運動量、熱、化学元素輸送に焦点を当て、星の進化コードへの実装に適した偏微分方程式系を提供すること。
• 今後、回転混合プロセスに磁場を組み込むための土台を築くこと。現在のモデルでは磁場は無視されている。

提案手法

• すべての流体力学的変数を、標準的な1.5次元手法を越えて高次の球面調和関数に展開し、緯度方向の変化を明示的に取り扱えるようにする。
• 温度揺らぎ（Ψₗ）に対して緩和法を用いた方程式系を構築することで、組成と回転の勾配が急峻な領域での数値解法の安定性を高める。
• 安定層上昇のため、乱流輸送を異方的（νᵥ ≪ νₕ、Dᵥ ≪ Dₕ）にモデル化し、水平方向の混合を優遇する。
• von Zeipelの研究に基づく熱的不均衡に起因する緯度循環と、角運動量輸送による回転へのフィードバックをモデルに組み込む。
• 星の表面（r = R）でポテンシャル場と連続する境界条件を適用し、Ψₗに対して径方向微分をゼロとし、Λₗに対して混合境界条件を設定する。
• 非線形項を避けるために、線形化・軸対称形式を用いて、差別的回転（Ψₗ経由）と組成不均一性（Λₗ経由）の相互作用を記述する。

実験結果

リサーチクエスチョン

• RQ1 放射層における緯度循環と剪切駆動乱流の取り扱いをどのように改善すれば、タコクリンのダイナミクスを正確に捉えることができるか？
• RQ2 高次の球面調和項を含めることで、緯度差別的回転と循環のモデル化にどのような影響を与えるか？
• RQ3 すべての時間微分を保持することで、回転と組成の勾配が急峻な急速な進化段階のシミュレーションがどのように向上するか？
• RQ4 温度揺らぎΨₗの緩和法による方程式の定式化は、星内部の急峻な勾配領域での数値的安定性をどのように向上させるか？
• RQ5 方程式系をどのように構築すれば、物理的一致性を保ちつつ、星の進化コードへの直接的な実装が可能になるか？

主な発見

• 本モデルは、対流層と放射層の間の差別的回転に起因するオクタポール型循環（タコクリン循環）を正確に再現できており、これは従来、手動で仮定されていたものである。
• 動的緩和を除くすべての時間微分を保持することで、角速度と組成の勾配が急激に変化する急速な進化段階のシミュレーションが可能になり、その精度が向上した。
• 温度揺らぎΨₗの緩和法による定式化により、組成勾配が急峻な領域での数値的感度が低下し、星の進化コードにおける安定性が向上した。
• 高次の球面調和関数の導入により、回転と循環の緯度変化のより正確な表現が可能になり、1.5次元のシェルラル近似を超えた記述が実現した。
• 角運動量、熱、化学元素輸送の自己整合的フレームワークを提供し、星の進化コードへの統合に適したものとなった。
• 本フレームワークは、磁場の将来的な組み込みを明示的に設計されており、太陽型星の平坦な回転プロファイルを説明する上で磁場が不可欠であるとされているが、現在のモデルでは中心部が急速に回転する予測となっている。