Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

技術的概要：非圧縮性ホール MHD 方程式に対する発散保存型混合有限要素法の誤差解析

問題提起
非圧縮性ホール・磁気流体力学（ホール MHD）系は、電導性流体における流体運動と磁場の相互作用をモデル化し、イオンスキン深度に匹敵するスケールで顕著となるホール効果を組み込んでいる。この効果は、誘導方程式に剛性が高く、高度に非線形な項（$\eta \text{curl}(\text{curl} B \times B)$）を導入し、分散挙動、急峻な勾配、高周波振動（ホイッスラー波）を引き起こす。これらの特徴は、分析的および計算的な課題を大幅に引き起こす。具体的には、3 次元における強解の全球的存在は未解決の問題であり、標準的な数値手法は安定性と、磁場が発散フリー（$\text{div } B = 0$）であるという厳密な要件の両方においてしばしば困難に直面する。この発散フリー条件の違反は、非物理的な磁気単極子を誘発し、深刻な数値的不安定性を引き起こす可能性がある。さらに、スカラー補助変数（SAV）法やスペクトル法に基づく既存の数値的手法は、しばしば真の物理的エネルギーではなく修正されたエネルギー汎関数を散逸させ、離散レベルで発散フリー条件を厳密に強制できず、あるいは厳密な誤差解析や幾何学的柔軟性を欠いている。

手法
これらの課題に対処するため、著者らは Voigt 正則化されたホール MHD 系に対する構造保存型、線形、完全離散混合有限要素法（FEM）を提案する。Voigt 正則化は、運動量方程式における時間微分を空間フィルタリングされたバージョン $(I - \alpha_1 \Delta)\partial_t u$ に置き換え、一般化オームの法則に電子慣性項 $\alpha_2 \partial_t J$ を保持する。この修正により、元の系の本質的な物理的構造とエネルギー法則を保持しつつ、適切性と数値的安定性が向上する。

数値スキームは、以下の構成要素を用いて構築される：

1. 混合定式化: 電場 $E$ と電流密度 $J$ を主要変数として用いて問題を再定式化し、運動量、誘導、オームの法則、アンペールの法則の方程式からなる結合系を導出する。
2. 空間離散化: 発散フリー条件を厳密に満たすことを保証するために、有限要素外微分幾何（FEEC）を採用する。
   • 速度 $u$ と圧力 $p$ は、MINI 有限要素ペア（最低次数の inf-sup 安定）を用いて近似される。
   • 磁場 $B$ と電流密度 $J$ は、それぞれ最低次数のラビエト＝トーマス（$RT_0$）空間とネデレク（$X_0$）空間を用いて近似される。
   • これらの空間は離散的なド・ラーム複体を形成し、$\text{div } B_h = 0$ が離散レベルで厳密に保存されることを保証する。
3. 時間離散化: 歪対称な線形陰的時間進行スキームが用いられる。対流項とローレンツ項は、線形性と無条件のエネルギー安定性を維持するために、前時間ステップの値を用いて半陰的に処理される。
4. 誤差解析: 著者らは、誤差を近似誤差（ストークス射影と直交射影を用いる）と離散誤差に分解することで、最適収束率を確立する。非線形結合項を処理するために、離散的ガフニー不等式とエネルギー評価を利用する。

主な貢献

• 構造保存スキーム: 本論文は、磁場に対するソレノイダル拘束を厳密に強制する線形・完全離散 FEM を提示する。これにより、精度や安定性を低下させる可能性のある射影ステップやペナルティ法の必要性を回避する。
• 無条件エネルギー安定性: 提案されたスキームは、剛性のホール項が存在する場合であっても、連続的な物理的エネルギー法則を反映し、無条件に離散的エネルギー散逸法則を満たす。
• 厳密な誤差評価: 著者らは、Voigt 正則化された問題に対する完全な誤差解析を提供し、適切な正則性仮定の下で $O(h^s + \tau)$ の順序の最適収束率を確立する。
• 正則化されていない系への拡張: 決定的なことに、この解析は正則化されていないホール MHD 系（$\alpha_1 = \alpha_2 = 0$）に拡張される。著者らは、抵抗率 $\sigma > 0$ である限り、正則化されていない系に対しても収束が証明可能であることを示す。非抵抗領域（$\sigma = 0$）では、Voigt 正則化が誤差評価において決定的な役割を果たすことが示される。
• 2.5 次元および 3 次元実装: この枠組みは、2.5 次元（並進対称性）および完全な 3 次元設定の両方で検証され、複雑な幾何学と物理現象に対する手法の適用可能性を実証する。

結果
2.5 次元および 3 次元の数値シミュレーションは、理論的知見を確認する：

• 収束: スキームは、速度、磁場、電流密度に対して予測された空間および時間収束率を達成する。速度と電流密度の収束率は、ホール項に起因する解の正則性の限界により、いくつかのケースでわずかに低下するが、磁場は最適収束を維持する。
• 発散保存: 最大発散誤差（$\max | \text{div } B_h |$）はシミュレーション全体を通じてソルバの許容誤差のレベルに留まり、ソレノイダル拘束の厳密な保存を確認する。
• エネルギー安定性: 系の全エネルギーは、理論的なエネルギー散逸法則と一致して、時間とともに単調に減少する。
• 物理的忠実度: アーノルド・ベルトラーミ・チャイルドレス（ABC）流、オルサグ・タン・ボルテックス、およびハリスシートにおけるホール駆動再結合のシミュレーションは、磁気再結合、磁気島やフラックスロープの形成、ホール MHD 力学に特徴的な四極子電流パターンなど、複雑な物理現象を成功裡に捉える。このスキームは、比較的粗いメッシュであってもこれらの特徴を正確に解像する。

意義
本論文は、ホール MHD 系のシミュレーションに対して数学的に厳密かつ計算的に堅牢な枠組みを提供すると主張している。Voigt 正則化と発散保存型混合有限要素定式化を組み合わせることで、著者らはホール項の非線形性と厳密な発散フリー要件という二重の課題を克服する。この手法は、安定性、収束性、複雑な幾何学への処理能力を備えた物理的に整合した近似を提供し、厳密な誤差解析や厳密な発散保存を欠いていた既存の手法が残したギャップを埋める。その結果は、ホール効果が支配的な天体物理および実験室プラズマ現象の研究において、このアプローチが実用的なツールとなり得ることを示唆している。