[論文レビュー] Projective and telescopic projective integration for the nonlinear BGK and Boltzmann equations

Ward Melis, Thomas Rey+1 more|arXiv (Cornell University)|Dec 18, 2017

Gas Dynamics and Kinetic Theory参考文献 64被引用数 10

ひとこと要約

本稿では、非線形BGKおよびボルツマン方程式に対して、剛性成分を減衰させるために内部時間ステップを活用し、外部のルンゲ＝クッタ段階で解を進める高次完全明示的テレスコピックプロジェクト型積分（TPI）スキームを提示する。この手法は剛性にほとんど依存しない計算コストを達成し、2次元においてBGKでは100を超える高速化、ボルツマンでは5.9の高速化を達成した。これは剛性に応じてプロジェクト型レベルが対数的に増加するためである。

ABSTRACT

We present high-order, fully explicit projective integration schemes for nonlinear collisional kinetic equations such as the BGK and Boltzmann equation. The methods first take a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution. Then, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. The procedure can be recursively repeated on a hierarchy of projective levels to construct telescopic projective integration methods. Based on the spectrum of the linearized collision operator, we deduce that the computational cost of the method is essentially independent of the stiffness of the problem: with an appropriate choice of inner step size, the time step restriction on the outer time step, as well as the number of inner time steps, is independent of the stiffness of the (collisional) source term. In some cases, the number of levels in the telescopic hierarchy depends logarithmically on the stiffness. We illustrate the method with numerical results in one and two spatial dimensions.

研究の動機と目的

• 平均自由行程が小さいことによって生じる非線形運動論的方程式（BGKおよびボルツマン）における剛性の課題に対処すること。
• 陰解法を避ける完全明示的時間積分スキームを、高次精度を維持したまま開発すること。
• 多次元位相空間における剛性衝突演算子の計算ボトルネックを克服すること。
• 高剛性（εが小さい）の流体領域でも安定性と効率性を保証すること。
• 階層的プロジェクト型時間ステッピングにより、スケーラブルかつ適応的時間積分を可能とすること。

提案手法

• 内部ステップに明示的エイラー法を用いて剛性成分を減衰させ、外部ステップに高次ルンゲ＝クッタ法を用いて解を進める、プロジェクト型レベルの階層を採用する。
• 内部ステップからの時間微分の推定を用いて、外部時間積分を明示的に構築することで、高次精度を達成する。
• テレスコピックプロジェクト型積分（TPI）を、複数のレベルにわたってプロジェクト型ステップを再帰的に適用することで、外部ステップ制限を低減する。
• 線形化衝突演算子のスペクトル理論を活用し、安定性を保証するとともにステップサイズ選択を支援する。
• 衝突演算子の評価には高速スペクトル法を、急勾配を解像するための精度を確保するためWENO2空間離散化を用いる。
• 内部ステップサイズを h₀ = ε に設定することで、外部ステップと剛性を分離し、εに依存しないコストを実現する。

実験結果

リサーチクエスチョン

• RQ1 陰解法を用いずに、剛性のある非線形BGKおよびボルツマン方程式に対して完全明示的かつ高次精度の時間積分が可能か？
• RQ2 テレスコピックプロジェクト型積分において、プロジェクト型レベルの数は剛性（ε）に対してどのように変化するか？
• RQ3 TPIの計算コストは標準的明示スキームと比較してどの程度で、剛性に依存しないか？
• RQ4 衝撃管のような強い勾配とマルチスケールダイナミクスを有する問題において、TPIスキームは精度を維持できるか？
• RQ5 剛性領域における安定性と効率性に、内部ステップサイズと外部ステップサイズの選択がどのように影響するか？

主な発見

• テレスコピックプロジェクト型積分法の計算コストは、剛性εにほとんど依存しない。これは最適な内部ステップサイズ h₀ = ε によるものである。
• BGK方程式では100を超える高速化が達成された。2次元ボルツマン方程式では、擬似マクスウェル粒子を用いて5.9の高速化を達成した。
• 必要となるプロジェクト型レベルの数は剛性に応じて対数的に増加するため、非常に小さなεに対しても効率的である。
• 外側スキームにPRK4、内側にFEを用いることで高次精度が維持され、特に渦を解像するためのWENO2空間離散化が不可欠であった。
• 標準的明示スキームがCFL制限により失敗するような流体領域（ε = 5·10⁻⁵）においても、本手法は安定かつ正確に動作した。
• このフレームワークは、将来の効率向上のため、適応的かつ空間依存のプロジェクト型レベル階層へ自然に拡張可能である。