[論文レビュー] High order semi-Lagrangian methods for the BGK equation

Maria Groppi, Giovanni Russo+1 more|arXiv (Cornell University)|Nov 28, 2014

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

ひとこと要約

本稿では、特性に沿ってBDFおよびDIRK時間積分法を用いた高次半ラグランジュスキームをBGK方程式に適用し、希薄流および流体領域の両方で漸近保存性と高い精度を達成した。数値実験において、特に補間を伴わないBDF3法が、計算効率と頑健性に優れ、DIRKスキームよりもコストパフォーマンスの面で優れていた。

ABSTRACT

A new class of high-order accuracy numerical methods for the BGK model of the Boltzmann equation is presented. The schemes are based on a semi-lagrangian formulation of the BGK equation; time integration is dealt with DIRK (Diagonally Implicit Runge Kutta) and BDF methods; the latter turn out to be accurate and computationally less expensive than the former. Numerical results and examples show that the schemes are reliable and efficient for the investigation of both rarefied and fluid regimes in gasdynamics.

研究の動機と目的

• 希薄流および流体領域の両方で安定かつ高精度を維持するBGK方程式の高次精度数値スキームの開発。
• 半ラグランジュ定式化におけるBDFおよびDIRK時間積分法の性能比較。
• 特に速度空間における空間補間の影響が精度および計算コストに与える影響の調査。
• Knudsen層の解像度を必要とせず、ε → 0 の流体極限においても精度を保持する漸近保存スキームの設計。
• 実用的な気体力学的応用を想定し、空間1次元、速度空間3次元に、Chu縮約を用いて拡張すること。

提案手法

• 半ラグランジュ的手法を用いて、時間発展を空間グリッドの制約から分離するため、BGK方程式を特性に沿って再定式化した。
• 時間積分にはBDF（後退差分公式）およびDIRK（対角的陰解法ルンゲ＝クッタ法）スキームを用い、BDF3法がより優れた計算効率を示した。
• 不均一な分布関数に対しては、WENO再構成（WENO23およびWENO35）を用いて高次空間精度を達成し、不連続付近での振動を最小限に抑えた。
• 位相空間グリッド上で直接計算を行うスキームの変種を提案し、速度空間における計算オーバーヘッドを削減した。
• Chu縮約を用いて、空間1次元および速度空間3次元に拡張し、次元を低減しながらも精度を維持した。
• 質量、運動量、エネルギーの保存およびエントロピー散逸を満たすことで、物理的一致性を保証した。

実験結果

リサーチクエスチョン

• RQ1 高次BDFおよびDIRKスキームを半ラグランジュ法と効果的に組み合わせることで、BGK方程式を高精度に解くことが可能か？
• RQ2 時間積分法の選択（BDF対DIRK）が、希薄流および流体領域における計算コストと精度にどのように影響するか？
• RQ3 速度空間における空間補間を回避することで、計算効率がどの程度向上し、精度に悪影響を与えないか？
• RQ4 提案スキームは流体極限（ε → 0）においても漸近保存性を維持するか、特にKnudsen数が小さい場合に有効か？
• RQ5 WENO再構成技術は、分布関数の不連続付近での安定性および解像度をどのように向上させるか？

主な発見

• 補間を伴わないBDF3法が、すべてのテストケースにおいて精度と計算コストの観点で最良の総合的性能を示した。
• BDFに基づくスキームは、精度は同等ながら、DIRKスキームよりも計算コストが低く抑えられた。
• WENO再構成を用いたスキームは、衝撃波や不連続付近で高次精度を達成し、不適切な振動を抑制した。
• スキームは漸近保存性を維持し、εが小さい場合でも正しく流体極限（圧縮性ナビエ＝ストークス挙動）を捉えた。
• 速度空間における補間の不在は、特に粗い速度グリッドを用いた問題において、計算コストを顕著に低減した。
• 数値結果により、希薄流および近平衡流体力学の高解像度シミュレーションにおいて、本スキームの頑健性が確認された。