[論文レビュー] Iterative Splitting Methods for Coulomb Collisions in Plasma Simulations

Jürgen Geiser|arXiv (Cornell University)|Jun 21, 2017

Gas Dynamics and Kinetic Theory参考文献 16被引用数 4

ひとこと要約

本稿では、プラズマシミュレーションにおけるクーロン衝突を記述する非線形確率微分方程式を解くための反復的分割スキームを提案する。Langevin方程式を用い、減衰係数と拡散係数を含む。反復的リラクゼーションにより、より大きな時間ステップを許容することで、精度と安定性が向上し、剛性のある問題において、Euler-Maruyama や Milstein などの陽的スキームに比べ、計算効率と収束性に優れる。

ABSTRACT

In this paper, we present splitting methods that are based on iterative schemes and applied to plasma simulations. The motivation arose of solving the Coulomb collisions, which are modeled by nonlinear stochastic differential equations. We apply Langevin equations to model the characteristics of the collisions and we obtain coupled nonlinear stochastic differential equations, which are delicate to solve. We propose well-known deterministic splitting schemes that can be extended to stochastic splitting schemes, by taking into account the stochastic behavior. The benefit decomposing the different equation parts and solve such parts individual is taken into account in the analysis of the new iterative splitting schemes. Numerical analysis and application to various Coulomb collisions in plasma applications are presented.

研究の動機と目的

• プラズマシミュレーションにおけるクーロン衝突を記述する剛性があり非線形な確率微分方程式を解く課題に対処すること。
• 直接解法に比べ、数値的安定性が向上し、より大きな時間ステップを許容できる反復的分割スキームの開発。
• ウィーナ過程と非線形リラクゼーションを組み込むことで、決定論的分割法を確率的設定に拡張すること。
• Euler-Maruyama や Milstein などの陽的スキームと比較して、反復的スキームの収束性と計算効率の分析。
• 速度、角度の余弦、方位角の動的挙動を含む、実際のプラズマ衝突モデルへの本手法の有効性の提示。

提案手法

• 本手法は、Langevin方程式を、移流（減衰）項と拡散項に分割し、時間ステップごとに逐次的に解く。
• 各反復では、固定時間ステップ τ を用いて線形化された確率微分方程式を解き、前段階のサブプロブレムの解を入力として使用する。
• 反復的プロセスは、バナッハ空間における抽象的コーシー問題として定式化され、移流項と拡散項を表す作用素 A と B がそれぞれ使用され、ウィーナ過程の増分が含まれる。
• 反復的更新は、i が奇数のとき c_i(t)、i が偶数のとき c_{i+1}(t) として定義され、時間 t^n における初期条件 c^n を用いる。
• 半群理論と作用素ノルムを用いて収束性を分析し、L2ノルムにおける誤差境界を導出する。
• 数値的に実装され、Euler-Maruyama および Milstein スキームと比較され、強収束および弱収束の指標を用いて評価される。

実験結果

リサーチクエスチョン

• RQ1 反復的分割スキームは、プラズマシミュレーションにおける非線形確率微分方程式に対して、直接解法に比べ高次の精度とより優れた安定性を達成できるか？
• RQ2 クーロン衝突モデルに適用した場合、反復的分割スキームの収束性と計算コストはどのように振る舞うか？
• RQ3 反復的リラクゼーションにより、Langevinに基づくプラズマシミュレーションで、数値的精度を損なわずに、どの程度大きな時間ステップを許容できるか？
• RQ4 反復的リラクゼーションは、速度、μ、φ 成分の分散および平均二乗誤差にどのような影響を与えるか？
• RQ5 時間ステップを変化させた場合、Euler-Maruyama や Milstein などの陽的スキームと比較して、反復的スキームの計算時間はどの程度か？

主な発見

• 反復的分割スキーム（Iter1 および Iter2）は、Δt = 10⁻¹ の時間ステップで安定した解を得られるが、Euler-Maruyama や Milstein などの陽的スキームは、このような大きなステップで振動を示す。
• 反復的スキームの計算コストは陽的スキームの 2〜3 倍高いが、より大きな時間ステップが使用可能であるため、全体の計算効率が向上する。
• 強収束誤差に関して、v、μ、φ について、反復的スキームは大きな時間ステップでも陽的スキームに比べ顕著に誤差が小さく抑えられ、Δt が小さくなるに従い誤差が減少する。
• 弱収束および分散解析から、Iter1 および Iter2 は、複数回のシミュレーションにわたって、より低い平均二乗誤差と一貫性のある統計的挙動を維持している。
• 図6の三次元解のプロットは、反復的スキームが、速度、μ、φ 成分の時間的挙動を陽的スキームよりもよく捉えていることを確認している。
• 誤差の計算には、Δt = 10⁻⁵ で実行された Milstein スキームに基づく参照解が使用され、反復的スキームは、より少ない時間ステップ数で同等の精度を達成している。