[論文レビュー] On Runge-Kutta discontinuous Galerkin schemes for Vlasov-Poisson systems

Yingda Cheng, Irene M. Gamba+1 more|arXiv (Cornell University)|Sep 28, 2012

Gas Dynamics and Kinetic Theory参考文献 49被引用数 3

ひとこと要約

本稿では、衝突なしプラズマをモデル化する1次元Vlasov-Poisson系を解く高次Runge-Kutta不連続ガラーキン（RKDG）スキームを提示する。この手法は保存性を保証し、制限子を用いて分布関数の正値性を維持し、フーリエ解析と多項式近似を通じて、Landau減衰や2ストリーム不安定性といったベンチマーク問題において高精度な解を得る。

ABSTRACT

In this paper we consider Runge-Kutta discontinuous Galerkin (RKDG) schemes for Vlasov-Poisson systems that model collisionless plasmas. One-dimensional systems are emphasized. The RKDG method, originally devised to solve conservation laws, is seen to have excellent conservation properties, be readily designed for arbitrary order of accuracy, and capable of being used with a positivity-preserving limiter that guarantees positivity of the distribution functions. The RKDG solver for the Vlasov equation is the main focus, while the electric field is obtained through the classical representation by Green's function for the Poisson equation. A rigorous study of recurrence of the DG methods is presented by Fourier analysis, and the impact of different polynomial spaces and the positivity-preserving limiters on the quality of the solutions is ascertained. Several benchmark test problems, such as Landau damping, two-stream instability and the KEEN (Kinetic Electrostatic Electron Nonlinear) wave, are given.

研究の動機と目的

• 1次元衝突なしプラズマ力学におけるVlasov-Poisson系を解く高次・保存的数値スキームの開発を目的とする。
• 正値性保持制限子を用いて、運動論的シミュレーションにおける分布関数の正値性を保証することを目的とする。
• フーリエ解析を用いてDG法の再発現および安定性特性を分析することを目的とする。
• Landau減衰や2ストリーム不安定性といった標準的ベンチマーク問題に対するスキームの妥当性を検証することを目的とする。
• 多項式次数および制限子の選択が解の品質および精度に与える影響を調査することを目的とする。

提案手法

• Runge-Kutta不連続ガラーキン（RKDG）法をVlasov方程式に適用し、高次精度と局所的保存性を実現する。
• ポアソン方程式のグリーン関数表現を用いて、電場を高精度にポテンシャル再構築する。
• 分布関数の負値化を防ぐために、RKDGフレームワークに正値性保持制限子を統合する。
• DG離散化の再発現およびスペクトル特性を調査するために、フーリエ解析を用いる。
• 多項式空間の次数を変化させることで、運動論方程式の解における収束性と精度を評価する。
• 数値実験を通じて、安定性、精度、保存性特性をベンチマーク問題に対して評価する。

実験結果

リサーチクエスチョン

• RQ1 異なる多項式空間が、Vlasov-Poisson系に対するRKDGスキームの精度および安定性にどのように影響を与えるか？
• RQ2 正値性保持制限子は、長時間シミュレーションにおいて物理的整合性をどの程度維持するか？
• RQ3 DG離散化における再発現の役割は何か？ また、解の品質にどのように影響するか？
• RQ4 RKDGスキームは、Landau減衰や2ストリーム不安定性といった主要なプラズマ現象をどの程度正確に捉えられるか？
• RQ5 高次空間的および時間的離散化は、解の精度および保存性特性にどのような影響を与えるか？

主な発見

• RKDGスキームは空間および時間において高次精度を示し、滑らかな解に対して理論的期待と一致する収束率を達成する。
• 正値性保持制限子は、すべてのテストされたベンチマークで非負の分布関数を的確に維持した。
• フーリエ解析により、RKDG法が良好な再発現特性を示し、不自然な振動を最小限に抑えることが明らかになった。
• 電場の期待される指数的減衰を保持しながら、最小限の数値拡散でLandau減衰を正確に再現した。
• 2ストリーム不安定性に関しては、不安定性の増大および非線形飽和を高解像度で解像した。
• KEEN波解は正確にシミュレートされ、本スキームが複雑な非線形運動論的現象を扱える能力を確認した。