[論文レビュー] Conservative discontinuous Galerkin/Hermite Spectral Method for the Vlasov-Poisson System

Francis Filbet, Tao Xiong|arXiv (Cornell University)|Apr 6, 2020

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

ひとこと要約

本稿では、速度空間でエルミート多項式、位相空間で不連続ガラーキン離散化を用いる、Vlasov-Poisson系の保存的不連続ガラーキン/エルミートスペクトル法を提案する。高次精度を達成するとともに、質量および全エネルギーの解析的保存を保証する。数値結果は、ベンチマークとしてのプラズマ不安定性における高次収束性と機械精度に近い保存量の維持を確認した。

ABSTRACT

We propose a class of conservative discontinuous Galerkin methods for the Vlasov-Poisson system written as a hyperbolic system using Hermite polynomials in the velocity variable. These schemes are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in general are hard to achieve within other numerical method frameworks for simulating the Vlasov-Poisson system. The proposed scheme employs discontinuous Galerkin discretization for both the Vlasov and the Poisson equations, resulting in a consistent description of the distribution function and electric field. Numerical simulations are performed to verify the order of accuracy and conservation properties.

研究の動機と目的

• 質量および全エネルギーを厳密に保存する高次数値スキームの開発を目的とする。
• 粒子-場（PIC）法の欠点（ノイズが $1/\sqrt{N}$ のスケーリングを示す）と、スペクトル法の欠点（周期的でない条件下でギブス振動を発生する）を克服することを目的とする。
• 不連続ガラーキン法の利点（局所的保存性と高次精度）と、エルミートスペクトル法の利点（効率的な速度空間表現）を統合することを目的とする。
• 振動を抑えるフィルタリング技術を導入し、保存性を損なわずに強い非線形領域における数値的安定性を確保することを目的とする。
• ランダム減衰、2ストリーム、およびバンプ・オン・テイル不安定性といった複雑なキネティック現象の正確なシミュレーションを可能とすることを目的とする。

提案手法

• 速度変数においてエルミート多項式を基底関数として用いることで、Vlasov-Poisson系を双曲型系に再定式化し、速度空間におけるスペクトル的精度を実現する。
• Vlasov方程式およびPoisson方程式の両方に対して不連続ガラーキン（DG）法を適用し、物理空間および速度空間における局所的保存性と高次精度を確保する。
• 半離散DGスキームは、弱形式および数値フラックスを用いて、離散的に質量および全エネルギーを保存するように構築され、保存性が証明されている。
• 不適切な振動を低減するために、特にHou-Liフィルタを適用し、保存性に影響を与えないようにする。
• 時間積分には、空間離散化の保存性と精度を維持する強安定性保存性ルンゲ・クッタ法を用いる。
• 本手法は1次元位相空間に実装され、$N_x \times N_H$ 個のグリッドポイントを用い、$N_H$ はエルミートモード数を表す。半ラグランジュ有限差分法（SLFDM）と比較される。

実験結果

リサーチクエスチョン

• RQ1 不連続ガラーキン/エルミートスペクトル法は、Vlasov-Poisson系において高次精度を達成するとともに、質量および全エネルギーを保存できるか？
• RQ2 提案スキームは、バンプ・オン・テイル不安定性や2ストリーム不安定性といった強い非線形プラズマ不安定性をシミュレートする際に、どのように性能を発揮するか？
• RQ3 フィルタリング技術は、保存性を劣化させることなく、不適切な振動をどの程度抑制できるか？
• RQ4 同等の解像度下で、本スキームの精度および保存性は、標準的な半ラグランジュ有限差分法（SLFDM）と比較してどの程度優れているか？
• RQ5 長時間にわたるシミュレーションにおいて、質量、運動量、エネルギーを機械精度で維持できるか？

主な発見

• 提案されたDG/エルミートスペクトル法は、すべてのベンチマークで理論的次数と一致する収束率を達成し、高次精度を実現した。
• 質量、運動量、全エネルギーは、すべてのシミュレーション中、機械精度で保存された。特に、$\alpha = 0.04$ のバンプ・オン・テイル不安定性のような強い非線形領域でも同様であった。
• 電場の$L^2$および$L^\infty$ノルムは、SLFDMの基準解と一貫した振るまいを示し、定性的および定量的に類似していた。
• 全電場エネルギーは時間経過に伴い良好に保存され、物理的に不適切な増大は観察されず、スキームの長期的安定性を裏付けた。
• 分布関数の表面プロットから、DGFSMは、$T=200$ の遅い時刻においても位相空間構造を正確に捉えており、数値的爆発や過剰な拡散の兆候は見られなかった。
• SLFDMの基準（$N_v = 320$）よりも少ない速度モード数（$N_H = 512$）で、同等またはより高い位相空間構造の解像度を達成した。これにより、本手法の効率性が示された。