[論文レビュー] Improving Large-scale Convection Zone-to-Corona Models

W. P. Abbett, G. H. Fisher|arXiv (Cornell University)|May 4, 2010

Solar and Space Plasma Dynamics被引用数 3

ひとこと要約

本稿では、太陽の大気における大規模3次元MHDモデルを向上させる2つの新規手法を提案する。1つは光沢的厚さの放射冷却を効率的に近似する手法であり、もう1つは時系列ベクトル磁場計測データを用いたデータ同化技術である。冷却手法は、従来のパrameter化手法に比べて物理的妥当性が向上し、データ同化手法により、観測された光球磁場を用いて物理的に整合性のある活動領域シミュレーションの初期化が可能になる。

ABSTRACT

We introduce two new methods that are designed to improve the realism and utility of large, active region-scale 3D MHD models of the solar atmosphere. We apply these methods to RADMHD, a code capable of modeling the Sun's upper convection zone, photosphere, chromosphere, transition region, and corona within a single computational volume. We first present a way to approximate the physics of optically-thick radiative transfer without having to take the computationally expensive step of solving the radiative transfer equation in detail. We then briefly describe a rudimentary assimilative technique that allows a time series of vector magnetograms to be directly incorporated into the MHD system.

研究の動機と目的

• 対流圏からコロナに至る大規模な太陽MHDモデルにおいて、計算的に効率的かつ物理的に意味のある放射冷却の取り扱いを開発すること。
• 観測された時系列ベクトル磁場計測データをMHDフレームワークに組み込むことで、現実的な活動領域シミュレーションの初期化という課題に取り組むこと。
• 手動パrameter化の代わりに巨視的放射移動近似を用いることで、RADMHDシミュレーションの物理的忠実度を向上させること。
• 対流圏からコロナへの磁気エネルギーおよびヘリシティ輸送のより正確なモデリングを可能にすること。
• 過剰に指定されすぎないが、観測された磁場の時間的変化を保持する安定でデータ制約付きの初期化手法を確立すること。

提案手法

• 全放射移動方程式を解かずに、下部太陽大気における光沢的厚さの冷却をモデル化する巨視的放射移動近似を開発した。
• RADMHDコードに新しい冷却手法を組み込み、より熱力学的現実性を高めた対流圏からコロナへの遷移をシミュレートした。
• 時系列ベクトル磁場計測データを用いて、速度逆問題を介して光球境界条件を駆動するデータ同化手法を実装した。
• 光球層における観測力（磁場計測時系列から得られる力）と内部MHD力のバランスを制御するため、0から1の間の信頼度行列（ξ）を定義した。
• 磁場計測時系列から導出された誘導流れ場（u_inv）を用い、光球層の力学的制約を課した。
• 磁場に平行および垂直な成分にそれぞれ別々の処理を施すことで、光球での運動量方程式への強制を物理的に整合性を持たせた。

実験結果

リサーチクエスチョン

• RQ1 全放射移動方程式を解かずに、大規模3次元MHDモデルにおいて光沢的厚さの放射冷却をどのように効率的に近似できるか？
• RQ2 観測された時系列ベクトル磁場計測データを、3次元MHDモデルに効果的に同化することで、現実的な活動領域シミュレーションの初期化が可能か？
• RQ3 この同化手法は、観測された磁場の時間的変化をどの程度保持しつつ、数値的安定性と物理的整合性を維持できるか？
• RQ4 従来の手動パrameter化に比べ、新しい冷却近似は熱力学的現実性と計算効率の両面でどの程度向上しているか？
• RQ5 データ同化された境界条件は、コロナにおける磁気エネルギーおよびヘリシティの進化にどのような影響を与えるか？

主な発見

• 新しい放射冷却近似は、従来のパrameter化されたニュートン冷却手法に比べ、下部大気におけるエネルギー損失の物理的根拠をより強く持つ。
• 計算コストを低減しつつも、大規模領域における対流圏からコロナへの遷移をモデル化するのに十分な精度を維持している。
• データ同化手法により、NOAA AR 8210の4時間分の時系列ベクトル磁場計測データが、120 Mm³のRADMHDモデルに安定的かつ物理的に整合性のある境界力として組み込まれた。
• 信頼度行列（ξ）により、観測制約と内部MHD力学のバランスを柔軟に制御でき、過剰な指定を回避できる。
• 得られた初期状態は、発散関数を最小化するだけの単純な緩和法に比べ、MHD整合性を持つ電流から出発するため、より物理的に妥当である。
• これらの手法を用いたシミュレーションは現在進行中であり、コロナ加熱、磁場の出現、エネルギー輸送の研究における現実性の向上に期待が持てる。