[論文レビュー] High-order Hamiltonian splitting for Vlasov-Poisson equations

Fernando Casas, Nicolas Crouseilles+2 more|arXiv (Cornell University)|Oct 7, 2015

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

ひとこと要約

本稿では、ハミルトニアンを運動エネルギー（T）と電気的エネルギー（U）の成分に分解することにより、Vlasov-Poisson系の高次ハイミルトニアン分割スキームを提案する。コンmutator解析とルンゲ＝クッタ＝ニストローム型の次数条件を用いて、1次元では6次精度の効率的スキーム、高次元では一般の高次スキームを構築し、厳密な誤差推定と収束性を証明し、古典的なストラング分割法を上回る優れた精度を示している。

ABSTRACT

We consider the Vlasov-Poisson equation in a Hamiltonian framework and derive new time splitting methods based on the decomposition of the Hamiltonian functional between the kinetic and electric energy. Assuming smoothness of the solutions, we study the order conditions of such methods. It appears that these conditions are of Runge-Kutta-Nystr{ö}m type. In the one dimensional case, the order conditions can be further simplified, and efficient methods of order 6 with a reduced number of stages can be constructed. In the general case, high-order methods can also be constructed using explicit computations of commutators. Numerical results are performed and show the benefit of using high-order splitting schemes in that context. Complete and self-contained proofs of convergence results and rigorous error estimates are also given.

研究の動機と目的

• 古典的なストラング分割法を超える精度を実現するため、Vlasov-Poisson系の高次時間積分スキームの開発。
• 系のハミルトニアン構造を活用して、カシミール不変量と位相空間構造を保存する幾何的積分法の構築。
• ルンゲ＝クッタ＝ニストローム型の形式およびリー代数を用いて、分割スキームの次数条件を導出し、簡略化する。
• 1次元において、恒等式 [[T, U], U] = 2mU を用いて段数を削減し、効率的な6次精度スキームを構築する。
• 重み付きソボレフ空間における、提案スキームの完全な収束証明と厳密な誤差推定を提供する。

提案手法

• Vlasov-Poissonのハミルトニアンを運動エネルギー（T）と電気的エネルギー（U）に分解する。
• TおよびUに対応する正確なフロー ϕτ_T および ϕτ_U の合成として時間分割スキームを構築する。
• Baker–Campbell–Hausdorff（BCH）計算およびコンミュータ展開を用いて、高次精度のための次数条件を導出する。
• 1次元では、恒等式 [[T, U], U] = 2mU を用いて次数条件を簡略化し、段数を減らした最適化された6次精度スキームを構築する。
• 高次元では、[[T, U], U] のフローを修正されたポテンシャルを用いて計算し、重み付きソボレフ空間におけるエネルギー推定を用いて収束性を導出する。
• 厳密な存在性および安定性結果（定理5.4）を適用して、適切な正則性および有界性仮定のもとで、H^r_νノルムにおける収束次数pを証明する。

実験結果

リサーチクエスチョン

• RQ1 ハミルトニアン分割とコンミュータ解析を用いて、Vlasov-Poisson系の高次分割スキームを体系的に構築可能か？
• RQ2 1次元における高次分割の簡略化された次数条件は何か？ そして、それらは段数を最小限に抑えた効率的な6次精度スキームの構築を可能にするか？
• RQ3 1次元におけるコンミュータ構造 [[T, U], U] = 2mU および高次元における類似構造は、高次スキームの設計および効率性にどのように影響を与えるか？
• RQ4 重み付きソボレフ空間におけるVlasov-Poisson系の高次分割スキームについて、厳密な収束性および誤差推定を確立可能か？
• RQ5 コンミュータフローのための修正されたポテンシャルの使用は、高次時間積分法の効率性および精度をどのように向上させるか？

主な発見

• 1次元においては、5段のみで構成される6次精度スキームが構築され、標準的な高次スキームと比較して計算コストを顕著に削減した。
• 1次元における恒等式 [[T, U], U] = 2mU により、次数条件の簡略化が可能となり、段数を減らした効率的な高次スキームの構築が可能になった。
• 高次元（d ≥ 2）においては、[[T, U], U] のフローが修正ポテンシャル K(f) を持つハミルトニアンフローに等価であることが示され、高次スキームの構築が可能になった。
• 重み付きソボレフ空間 H^r_ν において収束次数pが厳密に証明され、誤差境界 ‖(ψτ_p)^n(f) − ϕt_H(f)‖_H^r_ν ≤ Cτ^p が t = nτ ≤ T に対して成り立つ。
• ポisson積分法としての構造のおかげで、カシミール不変量（例：f の L^p ノルム）が保存され、長時間安定性が保証される。
• 数値結果により、高次分割スキームが古典的ストラング分割法に比べて精度および収束速度の面で優れていることが確認された。