An Adjoint Projection Formulation for Enforcing the divergence-free Constraint in Smoothed Particle Magnetohydrodynamics

技術要約：SPMHDにおける $\nabla \cdot \mathbf{B} = 0$ 制約を強制するための随伴射影定式化

問題提起
磁気流体力学（MHD）において、ソルレノイダル制約 $\nabla \cdot \mathbf{B} = 0$ は基本的な物理的要求である。格子ベースのコードでは、制約付き輸送（CT）を用いてこの制約を離散的に保持できるが、平滑化粒子磁気流体力学（SPMHD）のようなメッシュレス手法には、CTを直接適用するための自然な面–エッジ・トポロジーが欠如している。その結果、SPMHにおける主流のアプローチは、補助的なスカラー場を進化させてダイバージェンス誤差を伝播・減衰させる「ダイバージェンス・クリーニング」となっている。しかし、クリーニング手法は物理パラメータ（伝播速度や減衰タイムスケール）を導入するため、特に大きなダイナミックレンジを持つ系においては、それらがシミュレーション結果に影響を与える可能性がある。PriceとMonaghanによる過去の射影の試みは、補正のために重力のような総和を用いているが、誤差を測定するために使用されるダイバージェンス演算子と、補正を計算するために使用される勾配演算子との間に不整合が生じており、短波長のノイズに対して性能が低いという課題があった。

手法
本論文では、各MHD更新の後に楕円問題を解くことで、離散的なダイバージェンスフリー制約を強制する射影法を提案する。核心となる革新は、誤差を測定するために用いられる離散ダイバージェンス演算子の厳密な**随伴（アジョイント）**となる補正演算子を、体積重み付き計量内で構築することである。

1. 定式化: 本手法は、離散的な磁気エネルギーノルム（粒子体積 $V_i = m_i/\rho_i$ によって定義される）の変化を最小化しつつ、制約 $D\mathbf{B}^{n+1} = 0$ を満たす補正された磁場 $\mathbf{B}^{n+1}$ を求める。ここで $D$ は離散SPHダイバージェンス演算子である。
2. 随伴勾配: ラグランジュ未定乗数場 $\pi$ を導入することで、補正は $\mathbf{B}^{n+1} = \mathbf{B}^* - \mathcal{G}\pi$ として導かれる。演算子 $\mathcal{G}$ は、$D$ の体積計量による随伴 $\mathcal{G} = M_V^{-1} D^T$ として明示的に構築される。これにより、補正空間が体積内積に関してダイバージェンスフリー部分空間に対して直交することが保証される。
3. 離散化: 本論文では、$D$ とその随伴勾配 $\mathcal{G}$ の両方の明示的なSPH公式を導出している。決定的なのは、随伴勾配が、近傍リストの非対称性（変化する平滑化長や密度に起因するもの）を考慮し、特定の近傍リストの対称性と体積重み付けを使用している点である。
4. 数値解法: 結果として得られる射影方程式 $L\pi = r$ （ここで $L = D\mathcal{G}$）は、対称正定値線形システムを形成する。これにより、グローバル行列を明示的に組み立てることなく、前処理付き共役勾配法（PCG）を用いて効率的にシステムを解くことができる。実装においては、活性な粒子（$\pi$ が解かれる対象）と固定境界粒子を区別している。

主な貢献

• 一貫した随伴定式化: 本論文は、勾配補正が体積重み付き計量の下で、選択された離散ダイバージェンス演算子の厳密な随伴として構築された、最初のSPMHD射影定式化を提供している。これにより、従来の射影の試みで見られた離散化レベルでの不整合が解消された。
• エネルギー最小化特性: 体積計量は磁気エネルギーにとって自然な計量であることが示されている。したがって、この射影は離散的な磁気エネルギーを増加させないエネルギー最小化補正をもたらす。これは物理的な忠実度にとって極めて重要な特性である。
• 効率的なソルバー: 定式化は、PCGによって解ける対称正定値システムを導き、従来の総和法による $O(N^2)$ のコストを回避している。
• 実用的な停止基準: 著者らは、毎ステップでマシン精度までの収束を要求するのではなく、射影ステップ間の蓄積されたダイバージェンス誤差に基づく実用的な停止基準を導入しており、精度と計算コストのバランスをとっている。

結果
本手法は、2次元のDedner型ダイバージェンステストおよび3次元の磁化重力崩壊シミュレーションを用いてテストされた。

• ダイバージェンステスト: 摂動を加えた格子およびランダムな粒子分布を用いた2Dテストにおいて、本手法は十分なPCG反復によって、離散ダイバージェンス残差を浮動小数点丸め誤差レベル（$\sim 10^{-15}$）まで減少させることに成功した。特筆すべきは、従来のメソッドが苦戦していた平滑化長（$h$）と同程度の波長における誤差も補正できた点である。
• 崩壊シミュレーション: 3D崩壊ランにおいて、射影法は標準的な双曲型ダイバージェンス・クリーニングを大幅に上回る性能を示した。
  • 誤差減少: 実用的な停止基準（蓄積された誤差をMHD更新誤差の端数まで減少させる）を用いることで、射影法はダイバージェンス・クリーニングによるものよりも十分に低い正規化ダイバージェンス誤差を抑制した。
  • 堅牢性: 射影間隔（$N_{interval} = 1, 10, 100$）を変化させても、密度およびプラズマ-$\beta$ の構造は一貫していた。対照的に、ダイバージェンス・クリーニングを用いた実行では、流出構造やプラズマ $\beta$ に定量的な差異が見られた。
  • コスト: 実用的な間隔（$N_{interval} \approx 10$）において、射影の計算コストは、SPMHD更新コストの約1%から10%という控えめな範囲であった。極端なストレステストにおいても、本スキームは誤差を丸め誤差レベルまで駆動することが可能であった。

意義と主張
本論文は、この随伴射影定式化が、SPMHDおよび関連する粒子またはメッシュレスMHDスキームにおける $\nabla \cdot \mathbf{B}$ 誤差を制御するための、堅牢で魅力的な代替案であることを主張している。

強調されている主な利点は以下の通りである：

1. パラメータ独立性: 物理パラメータ（波速、減能スケール）に依存するダイバージェンス・クリーニングとは異なり、射影法のパラメータ（停止基準、間隔）は無次元であり、テストされた範囲内では解に最小限の影響を与えることが示されている。
2. 物理的一貫性: 補正が人工的な磁気エネルギーを注入しないことを保証する。
3. スケーラビリティ: この手法はSPMHDに限定されない。著者らは、適切な修正を加えることで、演算子ベースの射影のアイデアが他のメッシュレスまたは移動メッシュスキーム（GIZMO, AREPOなど）にも適用可能であることを示唆している。

著者らは、終了と頻度のための自由パラメータを選択する必要はあるものの、ダイバージェンス・クリーニングの感度の問題なしに一貫した高忠実度の結果を生み出す能力により、本手法が粒子ベースのMHDにおける数値的ダイバージェンスを制御するための優れた選択肢であると結論付けている。