Geometric numerical discretization of a quasineutral hybrid model of drift-kinetic electrons and fully kinetic ions

技術要約：ドリフト・キネティック電子と完全キネティックイオンの準中性ハイブリッドモデルにおける幾何学的数値離散化

問題提起
イオンをキネティックに扱い、電子を簡略化されたモデルで近似するハイブリッド・キネティック・モデルは、プラズマ物理学において計算効率と物理的正確性の間の重要なバランスを提供する。しかし、準中中性ハイブリッド・モデルの数値解法は、重大な課題を提示している。従来の粒子インセル（PIC）法は、離散的な電場ソルバーと粒子の更新との間の不整合により、長期的な数値的不安定性、保存則の違反、および非物理的なエネルギー成長に悩まされることが多い。さらに、高周波電磁波をフィルタリングするために変位電流を無視する準中性ハイブリッド・モデルでは、電場の標準的な時間発展方程式が除去される。これは、電場成分を一貫して決定するための代替的な定式化を必要とし、しばしば大規模で疎な線形系へとつながる。GEMPIC（Geometric Electromagnetic PIC）フレームワークのような既存の構造保存手法は、完全キネティックおよびドリフト・キネティック・モデルへと拡張されてきたが、ドリフト・キネティック電子と完全キネティック・イオンを組み合わせた「準中性」ハイブリッド・モデルに特化した構造保存離散化は、完全には実現されていなかった。

手法
著者らは、準中性ハイブリッド・ブラソフ・マックスウェル方程式を解くために、幾何学的電磁PICフレームワークであるGEMPICXを拡張する。本手法は以下の柱に基づいている：

1. 変分定式化： モデルは、ドリフト・キネティック電子（ゼロ・ラーマー半径極限）と完全キネティック・イオンのラグランジアンに基づく離散作用原理から導出される。準中性を強制するため、電場エネルギー項は除去されている。
2. 構造保存離散化： 双対格子を用いたミメティック有限差分法が採用されている。この離散化はド・ラーム複体を尊重し、離散演算子（勾配、回転、発散）が正確な幾何学的恒等式を満たすことを保証する。電場変数は、ガウスの法則と磁場のソレノイダル性を正確に保存するために、原始格子および双対格子の空間（ノード、エッジ、フェイス、ボリューム）に定義される。
3. 電場の再構成： 変位電流が存在しないため、電場は $\partial_t \mathbf{E}$ を通じて進化しない。代わりに：
   • 垂直成分 ($\mathbf{E}_\perp$) は、ドリフト・キネティック電子電流とイオン電流を関連付ける、修正された準中性アンペールの式を解くことによって得られる。
   • 平行成分 ($E_\parallel$) は、ファラデーの法則とアンペールの法則の時間微分から導かれるカール・カール（curl-curl）方程式によって決定される。これは、大規模で疎な、対称正定値線形系をもたらす。
4. 時間ステップ・スキーム： モデルに固有の高周波波分枝によって導入される硬さ（スティフネス）に対処するため、著者らは3つの時間積分戦略を調査している：
   • 完全陽解法： 4次の低ストレージ・ルンゲ・クッタ（LSRK）スキーム。
   • 暗黙・陽解法（IMEX）スキーム： 「CN-RK-CN」および「RK-CN-RK」という2つの分割戦略。これらは、最も高い周波数を含む硬い項をクランク・ニコルソン法を用いて暗黙的に扱い、非硬い項をLSRKを用いて陽的に進める。

主要な貢献

• GEMPICXの拡張： 本論文は、準中性ハイブリッド・ドリフト・キネティック領域を解くためにGEMPICXフレームワークを拡張し、連続系の幾何学的構造を保存する離散ラグランジアンと半離散方程式を導出した。
• 硬さの暗黙的処理： 著者らは、ハイブリッド準中性モデルに、陽解法に対して $\Delta t \propto \Delta z^2$ という厳しい安定制約を課す、非物理的な高周波波分枝が含まれていることを特定した。彼らは、硬い項を暗黙的に扱うためにシステムを分割することで、この制約を $\Delta t \propto \Delta z$ まで緩和できることを示した。これにより、大幅な計算節約が可能となる。
• 分散関係による検証： 数値モデルは、冷たいプラズマおよび温かいプラズマの両方について、解析的な分散関係と比較することで厳密に検証されている。研究は、背景磁場に対して平行および垂直に伝播する波をカバーしている。

結果
数値シミュレーションは、GEMPICXソフトウェア・フレームワークを用いて、準1次元周期領域で行われた。

• 波スペクトル： シミュレーションは、イオンサイクロトロン波（ICW）、圧縮性アルヴェン・ウィスラー波（CAW-WHW）、および温かいプラズマにおけるベルンシュタイン波を含む主要な固有モードを正常に再現した。数値結果は、低周波数において解析的な予測と優れた一致を示した。
• 安定性と効率性： 完全陽解法は、安定性を維持するために非常に小さな時間ステップ ($\Delta t \approx 0.02$) を必要としたが、IMEXスキームは大幅に大きなステップ ($\Delta \text{t} \approx 0.75$ で 'RK-CN-RK') を許容した。これにより、IMEXスキームは完全陽解法と比較して約95%の計算節約を実現した。
• 減衰： モデルは、温かいプラズマにおけるイオンサイクロトロン波の減衰を正常にシミュレートし、初期の過渡状態の後に解析的な減衰率と一致した。
• スキーム比較： 両方のIMEXスキームが効率性を提供したが、'RK-CN-RK' スキームは、硬い項の進化における位相誤差により非物理的なフォールディング（偽の分散）を示す 'CN-RK-CN' スキームよりも、高周波におけるスペクトル特性が優れていることが示された。

意義と主張
本論文は、提案された構造保存離散化が、準中性ハイブリッドプラズマをシミュレートするための堅牢かつ効率的なフレームワークを提供すると主張している。離散作用原理からスキームを導出することで、本手法はガウスの法則やソレノイダル制約の正確な保存といった主要な幾何学的特性を継承する。主な意義は、IMEX時間ステップ法を適用して、準中性ハイブリッドモデルに典型的な厳しい時間ステップ制限を克服することにあり、これにより、完全陽解法の法外なコストをかけることなく、アルヴェン波やイオンサイクロトロン力学の現実的な研究を可能にしている。著者らは、現在の研究がデカルト格子に限定されており、偏極および磁化効果を無視しているものの、本フレームワークが、より複雑な幾何学や有限ラーマー半径効果が重要な領域への将来的な拡張のための基礎を確立するものであると述べている。