[論文レビュー] Exponential Runge-Kutta schemes for inhomogeneous Boltzmann equations with high order of accuracy

Qin Li, Lorenzo Pareschi|arXiv (Cornell University)|Aug 13, 2012

Gas Dynamics and Kinetic Theory参考文献 26被引用数 9

ひとこと要約

本稿は、剛性領域における非一様ボルツマン方程式に対して、時間変動するマクスウェル平衡を用いて漸近的保存性（AP）と一様精度を達成する高次指数ルンゲ・クッタ（ExpRK）スキームを導入する。この手法は非線形解法を回避し、クヌーゼン数にかかわらず高次精度を保ち、シュウ・オシヤ REPRESENTATION を用いて正値性と強安定性を確保する。数値的検証により、3次までの精度が確認された。

ABSTRACT

We consider the development of exponential methods for the robust time discretization of space inhomogeneous Boltzmann equations in stiff regimes. Compared to the space homogeneous case, or more in general to the case of splitting based methods, studied in Dimarco Pareschi (SIAM J. Num. Anal. 2011) a major difficulty is that the local Maxwellian equilibrium state is not constant in a time step and thus needs a proper numerical treatment. We show how to derive asymptotic preserving (AP) schemes of arbitrary order and in particular using the Shu-Osher representation of Runge-Kutta methods we explore the monotonicity properties of such schemes, like strong stability preserving (SSP) and positivity preserving. Several numerical results confirm our analysis.

研究の動機と目的

• すべてのクヌーゼン数領域で安定かつ高精度を保つ、剛性のある非一様ボルツマン方程式の高次時間離散化手法の開発。
• 分割法や陰解法の限界を克服し、非線形連立方程式の解法を避ける一方で高次精度を維持する。
• 指数積分法を用いて、剛性のある運動論方程式において漸近的保存性（AP）と単調性（正値性や強安定性）を確保する。
• 領域分解や分割を用いずに、非一様ボルツマン方程式へ高次指数法を拡張する。
• ソース項を伴う運動論方程式に対して、強固で一様精度が保たれ、正値性を保つ数値フレームワークを提供する。

提案手法

• 非一様ボルツマン方程式に対して、固定平衡を用いるExpRK-Fと時間変動平衡を用いるExpRK-Vの2つのクラスの指数ルンゲ・クッタスキームを提案する。
• ルンゲ・クッタ法のシュウ・オシヤ表現を用いて、強安定性保存性（SSP）および正値性保存性を保証する。
• ボルツマン方程式を指数積分法で再定式化し、衝突項を指数時間積分を介して組み込む数値スキームを導出する。
• モーメントマッチングと漸近的解析を用いて、時間変動するマクスウェル平衡 $M$ 及びその時間微分 $\partial_t M$ を計算する。
• $\varepsilon$-正則化技術を用いて剛性のある衝突項を扱い、流体極限における一様精度を保証する。
• モーメント保存性と収縮性の議論を用いて検証し、各段階で $f^{(i)}$ と $M^{(i)}$ が同一のモーメントを持つことを証明する。

実験結果

リサーチクエスチョン

• RQ1 非一様ボルツマン方程式に対して、クヌーゼン数にかかわらず一様精度を保つ高次指数ルンゲ・クッタスキームを構築可能か？
• RQ2 剛性のある運動論方程式において、指数積分法で漸近的保存性と強安定性（SSP）をどのように達成できるか？
• RQ3 固定平衡状態と時間変動平衡状態の使用が、ExpRKスキームの精度と単調性に与える影響は何か？
• RQ4 ボルツマン方程式の高次指数積分法で、分布関数の正値性をどのように保証できるか？
• RQ5 剛性領域において、提案されたスキームは分割法やIMEX法と比較して、精度と安定性で優位性を示すか？

主な発見

• 提案されたExpRK-Vスキームは、時間変動する平衡を用いることで、すべてのクヌーゼン数領域で一様な高次精度を達成し、中間領域では固定平衡スキームを上回る性能を示す。
• シュウ・オシヤ表現により、強安定性保存性（SSP）および正値性保存性が実現され、物理的に意味のある解が保証される。
• 数値的検証により、ExpRK-FおよびExpRK-Vスキームの両方で、剛性領域においても3次収束率が確認された。
• $\varepsilon \to 0$ の極限で正しい流体極限（オイラー方程式）が得られることを確認し、漸近的保存性が裏付けられた。
• $f$ と $M$ の差が $O(\varepsilon)$ で有界であることが示され、モーメント方程式における補正項が小さく安定であることが保証された。
• カットオフ付きマクスウェル分子に対して、$d_2$-ノルムの収縮性 $\mathrm{d}_2(P(f), P(g)) \leq \mathrm{d}_2(f,g)$ が証明され、フーリエ領域における安定性が裏付けられた。