[論文レビュー] On the velocity space discretization for the Vlasov-Poisson system: comparison between Hermite spectral and Particle-in-Cell methods. Part 2: fully-implicit scheme

Enrico Camporeale, Gian Luca Delzanno+2 more|arXiv (Cornell University)|Dec 17, 2013

Magnetic confinement fusion research参考文献 68被引用数 9

ひとこと要約

本稿では、速度空間を直交ヘルミート基底を用いて離散化し、配置空間をフーリエ分解によって離散化する、完全陰解法によるヘルミートスペクトル法を提案する。この手法により、電荷、運動量、エネルギーの正確な保存が可能となる。数値実験では、計算コストの僅か数分の1で、PIC法に比べて最大10⁷倍の精度を達成しており、標準的なプラズマ問題において、高解像度PICシミュレーションでさえも大きく上回っている。

ABSTRACT

We describe a spectral method for the numerical solution of the Vlasov-Poisson system where the velocity space is decomposed by means of an Hermite basis, and the configuration space is discretized via a Fourier decomposition. The novelty of our approach is an implicit time discretization that allows exact conservation of charge, momentum and energy. The computational efficiency and the cost-effectiveness of this method are compared to the fully-implicit PIC method recently introduced by Markidis and Lapenta (2011) and Chen et al. (2011). The following examples are discussed: Langmuir wave, Landau damping, ion-acoustic wave, two-stream instability. The Fourier-Hermite spectral method can achieve solutions that are several orders of magnitude more accurate at a fraction of the cost with respect to PIC. This paper concludes the study presented in Camporeale et al. (2013) where the same method has been described for a semi-implicit time discretization, and was compared against an explicit PIC.

研究の動機と目的

• Vlasov-Poissonシミュレーションにおいて、電荷、運動量、エネルギーの正確な保存を保証する、フーリエ-ヘルミートスペクトル法に対する完全陰解法時間離散化の開発。
• 最近開発された完全陰解法PIC法と比較して、完全陰解法ヘリクテスペクトル法の計算効率と精度を評価すること。
• スペクトル的速度空間離散化が、エコー や高調波生成といった弱非線形な運動論的現象を、PIC法が捉えきれない場合に、それらを解像できるかを評価すること。
• 標準的なプラズマベンチマーク（ラムバル波、ランダム減衰、イオン音波、2ストリーム不安定性）において、ヘルミートスペクトル法とPIC法のコストパフォーマンスを比較して評価すること。

提案手法

• 速度空間は直交ヘルミート基底を用いて離散化され、分布関数がヘルミート多項式の級数展開として表現可能となる。
• 配置空間はフーリエスペクトル分解により離散化され、空間的表現の高精度が実現される。
• Crank-Nicolson時間離散化が採用され、2次精度と無条件安定性が保証される。
• 連立方程式を同時に解く完全陰解法の定式化により、電荷、運動量、エネルギーの正確な保存が実現される。
• 非線形性を有する時間積分に起因する非線形系を解くために、ジャコビアンフリー・ニュートン・クリロフ（JFNK）ソルバが用いられ、高い収束精度が達成される。
• 粒子サンプリングではなくスペクトル表現を用いるため、PICに内在する粒子ノイズを回避する。

実験結果

リサーチクエスチョン

• RQ1 完全陰解法ヘルミートスペクトル法は、Vlasov-Poissonシミュレーションにおいて、電荷、運動量、エネルギーの正確な保存を達成できるか？
• RQ2 標準的なプラズマベンチマークにおいて、完全陰解法ヘルミートスペクトル法の精度と計算コストは、完全陰解法PIC法と比べてどうなるか？
• RQ3 ヘリクテスペクトル法は、エコー や高調波生成といった弱非線形運動論的効果を、高粒子数でも捉えきれない標準PICシミュレーションが見逃す現象を解像できるか？
• RQ4 衝突なしプラズマの長期的シミュレーションにおいて、誤差/単位CPU時間という指標で、ヘルミートスペクトル法とPIC法の相対的有効性は何か？

主な発見

• イオン音波ベンチマークにおいて、完全陰解法ヘルミートスペクトル法は、主要モードの高次高調波を正確に捉えているが、PICシミュレーションは51,200粒子/セルでさえもそれらを解像できない。
• 約550秒のCPU時間で、ヘルミートスペクトル法の解は、最も優れたPIC解よりも約10⁷倍も高精度である。これは、コストパフォーマンスにおける顕著な優位性を示している。
• ヘルミートスペクトル法は、電荷、運動量、エネルギーの正確な保存を達成しているが、PIC法でさえも運動量を保存せず、継続的な監視が必要となる。
• 2ストリーム不安定性の線形成長率を正確に解像し、非線形遷移を高精度で捉えている。精度と計算効率の両面でPIC法を上回っている。
• ラムバル波、ランダム減衰、イオン音波、2ストリーム不安定性の全テストケースにおいて、ヘルミートスペクトル法の性能優位性は一貫しており、その頑健さと一般性が確認された。
• 結果から、ヘルミートスペクトル法は、PIC法ですら高粒子数を用いても、ノイズ制限により到達できない現象（例：エコー、波-波結合）を解像可能であることが確認された。