[論文レビュー] Numerical simulations of stellar winds: polytropic models

Rony Keppens, JP Goedbloed|arXiv (Cornell University)|Jan 27, 1999

Solar and Space Plasma Dynamics参考文献 3被引用数 29

ひとこと要約

本稿では、球対称および軸対称の幾何構造において、多項指数磁流体力学（MHD）方程式を解くためにVersatile Advection Code（VAC）を用いて、定常状態の超音速星風の数値的検討を提示する。1次元のWeber-Davis解を2次元に拡張し、星の表面で質量流量を指定し、流れと磁力線の整合性を強制することで、風領域とデッドゾーンの両方を物理的に整合的かつ事前に密度や磁場を指定せずにモデル化することに成功し、境界条件に依存しない自己自己一貫性のある解を得た。

ABSTRACT

We discuss steady-state transonic outflows obtained by direct numerical solution of the hydrodynamic and magnetohydrodynamic equations. We make use of the Versatile Advection Code, a software package for solving systems of (hyperbolic) partial differential equations. We proceed stepwise from a spherically symmetric, isothermal, unmagnetized, non-rotating Parker wind to arrive at axisymmetric, polytropic, magnetized, rotating models. These represent 2D generalisations of the analytical 1D Weber-Davis wind solution, which we obtain in the process. Axisymmetric wind solutions containing both a `wind' and a `dead' zone are presented. Since we are solving for steady-state solutions, we efficiently exploit fully implicit time stepping. The method allows us to model thermally and/or magneto-centrifugally driven stellar outflows. We particularly emphasize the boundary conditions imposed at the stellar surface. For these axisymmetric, steady-state solutions, we can use the knowledge of the flux functions to verify the physical correctness of the numerical solutions.

研究の動機と目的

• 2次元軸対称における定常状態・超音速星風をシミュレートするための堅牢な数値フレームワークの開発と検証。
• 1次元の解析的Weber-Davis風解を、軸対称的・多項指数的・磁化・回転を伴う流れに一般化すること。
• 自己一貫性のある数値解を用いて、星風の噴出において「風」領域と「デッド」領域の形成を調査すること。
• 密度や磁場を事前に指定せずに、質量流量と流れ・磁力線の整合性に基づく物理的に整合性のある星の境界条件の規定を確立すること。
• 今後の時間依存的・3次元のコロナ加熱および風加速のシミュレーションに向け、定常状態の多項指数モデルの妥当性を検証することによる基盤を構築すること。

提案手法

• 時間依存の多項指数MHD方程式をVersatile Advection Code（VAC）を用いて数値的に解き、時間微分をゼロに設定することで定常状態の解を得る。
• エネルギー方程式を解かずに系を閉じるため、多項指数状態方程式 $ p \sim \rho^\gamma $ を用いる。
• MHD方程式の混合型特性を考慮して、時間積分を完全に陰的に処理することで、定常状態への収束を効率的に達成する。
• 星の表面における境界条件を、指定された質量流量と、極面内での速度と磁力線の整合性を強制することで実装する。
• 数値解の物理的妥当性を検証するために、フラックス関数とストリームライン保存則を用いる。
• 段階的な解の構築：1次元の等温および多項指数風から出発し、徐々に複雑性を増す2次元軸対称モデル（磁化・回転を含む）に進む。

実験結果

リサーチクエスチョン

• RQ1 1次元のWeber-Davis風解を、軸対称的2次元配置に一般化するには、どのように物理的一致性を保つことができるか？
• RQ2 磁化・回転を伴う星風において、風領域とデッドゾーンの両方が自己一貫的に形成されるための最適な星表面境界条件は何か？
• RQ3 人工的な制約なしに、多項指数MHD風の流れで音速から超音速への遷移を数値的に正確に捉えることは可能か？
• RQ4 境界で密度や磁場ではなく質量流量を指定する選択が、解の物理的妥当性および収束性に与える影響は何か？
• RQ5 MHD方程式の混合型性質が、超音速星風シミュレーションにおける数値的解法戦略に与える影響は何か？

主な発見

• 数値解は、軸対称極限において解析的1次元Weber-Davis風解を正確に再現し、コードと手法の妥当性を検証した。
• 全テストされた多項指数モデルにおいて、衝撃を伴わず、滑らかで連続的な加速が、亜音速から超音速への流れに達している。
• 星の表面で密度や磁場を固定しないことで、開放的（風）および閉じた（デッド）磁力線領域の両方を自己一貫してモデル化できた。
• 質量流量と速度・磁力線の整合性に基づく境界条件の規定により、正しい回転結合と物理的一致性が確保された。
• 解は、質量流量および磁束が極面ストリームラインに沿って保存されることを確認し、数値的手法の正しさを裏付けた。
• このフレームワークは、時間依存的・3次元シミュレーションへの拡張が可能であり、今後のコロナ加熱および衝撃形成の研究の基盤を提供する。