[論文レビュー] High order and energy preserving discontinuous Galerkin methods for the Vlasov-Poisson system

Blanca Ayuso de Dios, Soheil Hajian|arXiv (Cornell University)|Sep 18, 2012

Gas Dynamics and Kinetic Theory参考文献 43被引用数 11

ひとこと要約

本稿では、全質量と全エネルギーを保存しつつ高次精度を維持する1次元Vlasov-Poisson系に対する修正された高次不連続ガレルキン(DG)法の族を提示する。この手法は位相空間における不連続な多項式近似を用い、境界条件を弱く強制するとともに、エネルギー保存を保証するための新規なフラックス定式化を採用しており、プラズマ物理学シミュレーションにおける数値ベンチマークによって検証されている。

ABSTRACT

We present a computational study for a family of discontinuous Galerkin methods for the one dimensional Vlasov-Poisson system that has been recently introduced. We introduce a slight modification of the methods to allow for feasible computations while preserving the properties of the original methods. We study numerically the verification of the theoretical and convergence analysis, discussing also the conservation properties of the schemes. The methods are validated through their application to some of the benchmarks in the simulation of plasma physics.

研究の動機と目的

• Vlasov-Poisson系に対して、特に質量およびエネルギーといった主要な物理的不変量を保存する高次不連続ガレルキンスキームの開発を目的とする。
• 粒子-フィールド法における数値ノイズの課題に対処するため、DG法に基づく決定論的で高次のEuler形式のアプローチを導入することを目的とする。
• 既存のDG定式化を修正し、計算の実行可能性を保ちつつエネルギー保存性と最適収束率を維持することを目的とする。
• 標準的なプラズマ物理学ベンチマーク上での手法の検証を通し、長時間シミュレーションにおけるロバストネスと精度を示すこと。
• ベルンシュタイン多項式を用いて電場の符号および極値を効率的に検出する手法を検討し、安定性および正値制御を支援すること。

提案手法

• 本手法は、位置および速度変数の両方において高次多項式基底関数（ラグランジュおよびベルンシュタイン）を用いた不連続ガレルキン空間離散化を採用する。
• 半離散DG定式化が用いられ、境界条件の弱い強制を通じて質量およびエネルギー保存が確保される数値フラックスが構築される。
• 電場の再構築にはベルンシュタイン多項式が用いられ、係数解析による符号変化および極値の効率的・信頼性の高い検出が可能になる。
• 計算効率の向上を図るため、既存の研究から変更を加え、特に基底変換のためのバーンスタイン行列の逆行列を事前計算する。
• 時間積分には標準的な陽的時間ステッピングが用いられ、ブロック対角質量行列構造を活用して効率を向上させる。
• 理論的分析および数値的検証により、分布関数および電場の両方において最適収束率が保証されている。

実験結果

リサーチクエスチョン

• RQ1 Vlasov-Poisson系において、エネルギー保存性を維持しつつ最適収束率を達成できる高次不連続ガレルキン法を設計可能か？
• RQ2 エネルギー保存型DGスキームの計算コストを、保存性の性質を損なわず低減可能か？
• RQ3 ベルンシュタイン多項式表現は、DGスキームにおける電場評価の効率性およびロバストネスをどの程度向上させるか？
• RQ4 修正されたDGスキームは、標準的なプラズマ物理学ベンチマークにおいて、精度および長時間安定性の観点でどの程度の性能を示すか？
• RQ5 基底関数の選択（ラグランジュ対ベルンシュタイン）が、電場再構築の条件数および正値保存性に与える影響はいかなるか？

主な発見

• 提案されたDG法は、$L^2$-ノルムにおいて分布関数および電場の両方で、$k+1$次の最適収束率を達成する。
• 局所的不連続ガレルキン定式化および境界条件の弱い強制により、全質量が正確に保存される。
• 理論的証明および数値的確認により、修正されたスキームは長時間シミュレーションでさえも系全体の全エネルギーを保存することが示された。
• ベルンシュタイン多項式の使用により、電場の符号変化および極値の効率的かつ信頼性の高い検出が可能になり、数値的ロバストネスが向上した。
• 基底変換のための逆バーンスタイン行列を事前計算することで、時間積分中の係数チェックが高速化され、計算コストが顕著に低減された。
• 数値ベンチマークにより、本手法の精度および安定性が、ランダム減衰および2ストリーム不安定性といった重要なプラズマ現象のシミュレーションにおいて確認された。