[論文レビュー] Splitting methods for Fourier spectral discretizations of the strongly magnetized Vlasov-Poisson and the Vlasov-Maxwell system

J. Ameres|arXiv (Cornell University)|Jul 11, 2019

Gas Dynamics and Kinetic Theory参考文献 40被引用数 5

ひとこと要約

本稿は、強磁場下のVlasov-PoisissonおよびVlasov-Maxwell系のフーリエスペクトル離散化に対する新規分割法を導入し、4次元および3次元で高精度なシミュレーションを可能にしている。ケルビン＝ヘルムホルツ不安定性やウェイベル流れ不安定性といった困難なテストケースにおいても安定性と幾何的保存性を示し、ストラング分割とハミルトニアン積分法を用いることで、低解像度でもエネルギー誤差が10⁻¹²未塔に抑えられる。

ABSTRACT

Fourier spectral discretizations belong to the most straightforward methods for solving the unmagnetized Vlasov--Poisson system in low dimensions. In this article, this highly accurate approach is extended two the four-dimensional magnetized Vlasov--Poisson system with new splitting methods suited for strong magnetic fields. Consequently, a comparison to the asymptotic fluid model is provided at the example of a turbulent Kelvin--Helmholtz instability. For the three dimensional electromagnetic Vlasov--Maxwell system different novel charge conserving implementations of a Hamiltonian splitting are discussed and simulation results of the Weibel streaming instability are presented.

研究の動機と目的

• 強磁場下の高次元Vlasov系におけるフーリエスペクトル離散化に対する幾何的分割法の欠如に応える。
• 強い磁場下における4次元磁場付きVlasov-Poisson系のための効率的で構造保存型時間積分法を開発する。
• ハミルトニアン分割法を3次元Vlasov-Maxwell系に拡張し、電荷およびエネルギー保存を確保する。
• ケルビン＝ヘルムホルツおよびウェイベル流れ不安定性を含むベンチマークプラズマ不安定性に対して、低解像度および高解像度の両方で手法を検証する。
• 提案手法が粗いグリッドおよび強い磁場下でも長期間にわたり安定性と幾何的精度を維持することを示す。

提案手法

• 位置空間および速度空間における移流を分離するためのハミルトニアン分割法を適用し、フーリエ変換を用いて正確に解を得る。
• 空間および速度空間の離散化にフーリエスペクトル法を用い、計算効率を高めるためにFFTWを活用する。
• 定数係数移流のフーリエ空間内での正確な積分を実現するため、指数型分割スキーム（例：スコヴェル法および指数型ボリス法）を実装する。
• Vlasov-Maxwell系における電荷保存を確保するため、元のハミルトニアン構造を保存する幾何的積分法を用いる。
• 2次精度を持つストラング分割を時間積分に適用し、エネルギーおよび運動量保存性の性質を維持する。
• ベッセル関数およびガウス求積を用いて、スコヴェル法におけるアモール方程式のジロ平均化を処理し、正確な時間積分を実現する。

実験結果

リサーチクエスチョン

• RQ1 新規分割スキームを用いることで、フーリエスペクトル法を4次元の強磁場下Vlasov-Poisisson系に効果的に拡張可能か？
• RQ2 提案された分割法は、強い磁場下における乱流不安定性（例：ケルビン＝ヘルムホルツ不安定性）のシミュレーションにおいて、どのように性能を発揮するか？
• RQ3 ハミルトニアン分割法は3次元Vlasov-Maxwell系に適応可能か？ これにより電荷およびエネルギー保存が保証されるか？
• RQ4 解像度および時間刻みが、ウェイベル型不安定性の長時間シミュレーションにおけるエネルギーおよび運動量誤差に与える影響は何か？
• RQ5 標準的な擬似スペクトル法と比較して、新規スキームは幾何的構造および安定性をどの程度維持するか？

主な発見

• 提案されたスコヴェル法の分割スキームは、漸近的極限が存在しない状況下でも、広い磁場強度範囲において高精度および安定性を維持する。
• ケルビン＝ヘルムホルツ不安定性のシナリオでは、強い磁場下におけるキネティックモデルが流体モデルとよく一致し、漸近的領域の妥当性が裏付けられる。
• ハミルトニアン分割法を用いたVlasov-Maxwellソルバーは、高解像度（128×128×128グリッドポイント）でエネルギー誤差が1.71×10⁻¹²未塔に抑えられ、高次精度GEMPIC法と同等の性能を示す。
• 低解像度（32×32×32）でも、長時間にわたり安定で、エネルギー誤差が10⁻¹²未塔に抑えられ、強靭性が示される。
• ウェイベル流れ不安定性のシミュレーションでは、キネティックエネルギーおよび磁場の成長が一貫して観測され、対称および非対称ケースでは異なるフィラメント構造が明確に再現される。
• ベッセル関数およびガウス求積の使用により、スコヴェル法におけるアモール方程式の正確な時間積分が可能となり、長期間の保存性を実現する重要な要因となる。