[論文レビュー] Toward detailed prominence seismology - I. Computing accurate 2.5D magnetohydrodynamic equilibria

J. W. S. Blokland, Rony Keppens|arXiv (Cornell University)|Jun 24, 2011

Ionosphere and magnetosphere dynamics参考文献 24被引用数 17

ひとこと要約

本稿では、太陽のフレア帯においてローレンツ力、重力、圧力勾配、ねじれた磁場をバランスさせる2.5次元磁気力学(2.5D MHD)平衡を高精度で計算する数値的手法を提示する。有限要素法コード（FINESSE）を用いて、4次収束を達成し、高温のキャビティ内に冷却凝縮を有する現実的な平衡状態と多層構造を生成する。これにより、将来的なMHDモードの精密なセイズモロジー解析が可能となる。

ABSTRACT

Context. Prominence seismology exploits our knowledge of the linear eigenoscillations for representative magnetohydro- dynamic models of filaments. To date, highly idealized models for prominences have been used, especially with respect to the overall magnetic configurations. Aims. We initiate a more systematic survey of filament wave modes, where we consider full multi-dimensional models with twisted magnetic fields representative of the surrounding magnetic flux rope. This requires the ability to compute accurate 2.5 dimensional magnetohydrodynamic equilibria that balance Lorentz forces, gravity, and pressure gradients, while containing density enhancements (static or in motion). Methods. The governing extended Grad-Shafranov equation is discussed, along with an analytic prediction for circular flux ropes for the Shafranov shift of the central magnetic axis due to gravity. Numerical equilibria are computed with a finite element-based code, demonstrating fourth order accuracy on an explicitly known, non-trivial test case. Results. The code is then used to construct more realistic prominence equilibria, for all three possible choices of a free flux-function. We quantify the influence of gravity, and generate cool condensations in hot cavities, as well as multi- layered prominences. Conclusions. The internal flux rope equilibria computed here have the prerequisite numerical accuracy to allow a yet more advanced analysis of the complete spectrum of linear magnetohydrodynamic perturbations, as will be demonstrated in the companion paper.

研究の動機と目的

• フレア帯構造においてローレンツ力、重力、圧力勾配をバランスさせる2.5次元MHD平衡を数値的に高精度に計算する手法の開発。
• 拡張されたグラド・シャフランォフ方程式を解くことにより、高温キャビティ内に冷却凝縮を有する現実的なフレア帯構造および多層構造の密度増強をモデル化すること。
• 異なるフラックス関数の仮定（温度、密度、またはエントロピー依存性）のもとで、重力が平衡構造、特にシャフランォフシフトに与える影響を定量化すること。
• 高精度な平衡モデルを提供することで、将来の線形MHD安定性および固有モード解析の基盤を築くこと。
• 既知の解析解との比較を通じて数値手法の妥当性を検証し、非自明なテストケースにおいて4次収束を示すこと。

提案手法

• トランスレーション不変性を有する2.5次元MHD平衡を計算する目的で、拡張されたグラド・シャフランォフ方程式を導出し、有限要素法に基づくコード（FINESSE）を用いて数値的に解く。
• 重力、圧力勾配、ねじれた磁場を考慮し、フラックス関数の選択として、温度、密度、またはエントロピー依存性の3つの異なる選択肢を採用する。
• 重力下における円形フラックスロープのシャフランォフシフトについて、解析的予測を導出し、数値結果の妥当性を検証する。
• 等温で重力平衡の状態におけるDungey（1953）の解析解と比較することで、数値解の4次収束を実証する。
• 太陽の条件を反映するように変化する重力パラメータを用いて、単層および二重層構造のフレア帯に対して平衡を生成する。
• 磁場線および密度等高面の3次元可視化を用いて、フラックスロープのヘリカルかつねじれた構造と凝縮の層状構造を示す。

実験結果

リサーチクエスチョン

• RQ1 有限要素法コードは、フレア帯構造においてローレンツ力、重力、圧力勾配をバランスさせる2.5次元MHD平衡をどの程度高精度に計算できるか？
• RQ2 重力は磁場軸（シャフランォフシフト）をどの程度シフトさせるか。また、フラックス関数の依存性に関する異なる仮定がこのシフトに与える影響はいかほどか？
• RQ3 高温キャビティ内に冷却凝縮を有する現実的で、多層構造を有するフレア帯構造は、自己自己一貫的に平衡状態として生成可能か？
• RQ4 強力な重力下において、フラックス関数が温度、密度、またはエントロピーに依存する場合、平衡の性質（圧力、密度、ベータ）にどのような差が生じるか？
• RQ5 FINESSEコードは、複雑な非線形MHD平衡をどの程度の数値的精度で解明できるか。また、解析的ベンチマークと比較してその性能はいかがなものか？

主な発見

• FINESSEコードは、解析的Dungey（1953）解との比較において4次収束を達成しており、中程度の解像度でも高い数値的精度を確認している。
• 弱い重力下では、数値的に得られたシャフランォフシフトが解析的予測とよく一致する。一方、強い重力下では、径方向力のバランスに基づく数値的手法のみが正確にシフトを捉えている。
• 冷却フレア帯では、重力を増加させると、最大圧力および最大密度の位置が下方にずれる。この効果は、温度がフラックス関数に依存する場合に最も顕著に現れる。
• フラックス関数の3つの選択肢（T, n, S）は、強い重力下で圧力および密度分布に顕著な差をもたらし、観測的制約により正しい物理的仮定を区別できる可能性を示唆する。
• 重力ポテンシャルが十分に強い場合、明確な密度および圧力層を有する二重層構造のフレア帯が、平衡状態として成功裏に生成された。
• コードは、ヘリカルな磁場線および層状の凝縮を示す現実的な3次元可視化を生成しており、将来的な安定性およびセイズモロジー研究へのモデルの妥当性を確認している。