[論文レビュー] Optimized Runge-Kutta Methods with Automatic Step Size Control for Compressible Computational Fluid Dynamics

Hendrik Ranocha, Lisandro Dalcín+2 more|arXiv (Cornell University)|Apr 14, 2021

Numerical methods for differential equations参考文献 84被引用数 51

ひとこと要約

本稿では、固定CFL数の代わりに誤差ベースの制御を用いる、圧縮性計算流体力学（CFD）のための最適化された明示的ルンゲ・クッタ法と適応的ステップサイズ制御を提案する。この手法は、安定性制限領域では近似的に最適な効率性を達成し、精度制限領域では正確な時間誤差制御を提供する。工業的CFDベンチマークおよび困難な粘性衝撃波シミュレーションにおいて、従来のCFLベースの手法を上回る性能を発揮する。

ABSTRACT

We develop error-control based time integration algorithms for compressible fluid dynamics (CFD) applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime. Focusing on discontinuous spectral element semidiscretizations, we design new controllers for existing methods and for some new embedded Runge-Kutta pairs. We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice. We compare a wide range of error-control-based methods, along with the common approach in which step size control is based on the Courant-Friedrichs-Lewy (CFL) number. The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances, while additionally providing control of the temporal error at tighter tolerances. The numerical examples include challenging industrial CFD applications.

研究の動機と目的

• 圧縮性CFDシミュレーションにおけるCFLベースの時間刻み制御の非効率性と手動によるCFLチューニングの負担を解消すること。
• 安定性制限領域および精度制限領域の両方で高い効率性を維持する誤差制御時間積分アルゴリズムの開発。
• 不連続スペクトル有限要素法に特化した最適化された埋め込みルンゲ・クッタ対およびPIDベースのステップサイズ制御器の設計。
• 誤差ベース制御が、緩い許容誤差では自然に最大安定CFL数に近づき、きつい許容誤差では時間誤差を制御できることを示すこと。
• 実世界のCFDアプリケーションにおいて安定性と効率性を保証するための、制御器パラメータ選定の実用的ガイドラインの提供。

提案手法

• h-p適応スペクトル有限要素法に適した、最適化された安定性領域と埋め込み誤差推定器を備えた新しい埋め込みルンゲ・クッタ対の設計。
• 局所切捨て誤差推定に基づいてステップサイズを調整するPIDベースのステップサイズ制御器の実装。これにより、安定性と正確性の両方が確保される。
• ボーサー表記法を用いて、低記憶、最初のステップと最後のステップが同じ（FSAL）のルンゲ・クッタ対と埋め込み誤差推定を定義。
• 圧縮性エーラー方程式およびナビエ・ストークス方程式の不連続ガラーキン半離散化と時間積分器を統合。
• PETScに基づくSSDCソルバーフレームワーク内に制御器を統合し、スケーラブルで非構造的、曲線グリッドを用いたシミュレーションを可能に。
• 数値実験を用いて制御器パラメータ（例：比例、積分、微分ゲイン）をキャリブレーションし、ロバスト性と安定性を確保。

実験結果

リサーチクエスチョン

• RQ1 誤差ベースのステップサイズ制御は、手動によるCFLチューニングを必要とせずに、安定性制限領域におけるCFLベース制御と同等の効率性を達成できるか？
• RQ2 最適化された埋め込みルンゲ・クッタ対に適応的制御器を組み合わせた手法は、異なる許容誤差において、標準的手法と比較して正確性と計算コストの点でどのように差をつけるか？
• RQ3 高次元、不連続スペクトル有限要素離散化において、ステップサイズ制御の安定性とロバスト性を保証する制御器パラメータ設定は何か？
• RQ4 誤差ベース制御は、安定性限界付近で近似的に最適なステップサイズを維持しつつ、きつい許容誤差でも時間誤差を制御できるか、その程度はどの程度か？
• RQ5 新規手法は、粘性衝撃波や乱流といった困難な産業的CFD問題において、どのように性能を発揮するか？

主な発見

• 最適化された誤差制御手法は、安定性制限領域において近似的に最適な効率性を達成し、緩い許容誤差では最大安定CFL数に近いステップサイズを自動的に採用する。
• きつい許容誤差では、信頼性の高い時間誤差制御が可能であり、誤差は有界に保たれ、期待される次数で収束することが、粘性衝撃波テストにおけるL2誤差ノルムで示された。
• p=4の多項式において、RK5(4)10F[3S*+]スキームは許容誤差10−6でL2誤差5.42×10−5を達成し、関数評価回数はたったの135,081回で、高い効率性を示した。
• RK4(3)9F[3S*+]スキームは許容誤差10−8で相対L2誤差1.00×10−7を達成し、関数評価回数108,879回で、優れた正確性と安定性を示した。
• 提案された制御器は、特に精度制限領域においてCFLベース制御を上回る効率性とロバスト性を示した。CFL制御は誤差を効果的に低減できない場合がある。
• 制御器パラメータ（例：β = (0.45, -0.13, 0.00)）は安定性にとって極めて重要であり、不適切な設定はステップ拒否や発散を引き起こした。これにより、体系的なキャリブレーションの必要性が強調された。