[論文レビュー] Extremum-Preserving Limiters for MUSCL and PPM

Michael Sekora, Phillip Colella|ArXiv.org|Mar 24, 2009

Computational Fluid Dynamics and Aerodynamics参考文献 15被引用数 8

ひとこと要約

本稿では、局所的極値でのみ2階微分に非線形制約を適用することにより、滑らかな極値でも高次精度を維持する、MUSCLおよびPPMスキームの極値保存型制限子を導入する。この手法は、単調性と安定性、$4^{\text{th}}$-次空間精度を保ちながら計算コストを最小限に抑え、不連続性での性能を損なわず、滑らかな領域における数値誤差を顕著に低減する。

ABSTRACT

Limiters are nonlinear hybridization techniques that are used to preserve positivity and monotonicity when numerically solving hyperbolic conservation laws. Unfortunately, the original methods suffer from the truncation-error being first-order accurate at all extrema despite the accuracy of the higher-order method. To remedy this problem, higher-order extensions were proposed that relied on elaborate analytic and geometric constructions. Since extremum-preserving limiters are applied only at extrema, additional computational cost is negligible. Therefore, extremum-preserving limiters ensure higher-order spatial accuracy while maintaining simplicity. This report presents higher-order limiting for (i) computing van Leer slopes and (ii) adjusting parabolic profiles. This limiting preserves monotonicity and accuracy at smooth extrema, maintains stability in the presence of discontinuities and under-resolved gradients, and is based on constraining the interpolated values at extrema (and only at extrema) by using nonlinear combinations of second derivatives. The van Leer limiting can be done separately and implemented in MUSCL (Monotone Upstream-centered Schemes for Conservation Laws) or done in concert with the parabolic profile limiting and implemented in PPM (Piecewise Parabolic Method). The extremum-preserving limiters elegantly fit into any algorithm which uses conventional limiting techniques. Limiters are outlined for scalar advection and nonlinear systems of conservation laws. This report also discusses the fourth-order correction to the point-valued, cell-centered initial conditions that is necessary for implementing higher-order limiting.

研究の動機と目的

• 伝統的なMUSCLおよびPPMスキームの制限子が極値で$4^{\text{th}}$-次精度を失う問題に対処すること。
• 不連続性での性能を損なわず、滑らかな極値でのみ単調性と精度を保証する高次制限戦略の開発。
• 制限プロセスが滑らかな領域で$4^{\text{th}}$-次空間精度を維持しつつ、計算的に効率的であることを保証すること。
• 特性に基づく制限または成分ごとの適用を用いて、非線形双曲型保存則へのこのアプローチの拡張。

提案手法

• 極値における補間値に、2階微分の組み合わせを用いた非線形制約を適用し、極値保存を実現する。
• MUSCLではvan Leer制限子を勾配再構築に用い、$4^{\text{th}}$-次精度を達成するため極値保存型補正を組み合わせる。
• PPMでは$4^{\text{th}}$-次および$6^{\text{th}}$-次パラボリック再構築を実装し、極値保存型制限子は極値でのみ適用する。
• セル中心初期条件に対する$4^{\text{th}}$-次補正を構築し、MUSCLおよびPPMの両方で高次制限を可能にする。
• 非線形系では、単一成分の変数に成分ごとまたは特性場に制限子を適用することで、単調性と精度を維持する。
• 極値保存型制限子を既存のMUSCLおよびPPMフレームワークに最小限のアルゴリズム的変更で統合し、シンプルさと効率性を保つ。

実験結果

リサーチクエスチョン

• RQ1 極値保存型制限子は、MUSCLおよびPPMスキームで滑らかな極値における$4^{\text{th}}$-次精度を回復させつつ、単調性や安定性を損なわずに行えるか。
• RQ2 滑らかでない解および不連続解において、極値保存型制限子は従来の制限子と比較して誤差低減にどの程度優れているか。
• RQ3 極値保存型制限子は制限定数$C_{VL}$および$C_{PPM}$の選択に対してどれほど感度を示すか。
• RQ4 非線形双曲型保存則に、高次精度と単調性を維持したまま、極値保存型制限子を効果的に拡張できるか。

主な発見

• ガウス波（G）では、$N=128$セルにおける$L_1$誤差が約2.5倍低減され、収束次数$R_1 \approx 2.5$となり、極値における$4^{\text{th}}$-次精度の向上が示された。
• 半円波（SC）では、$N=128$における$L_1$誤差が約2.4倍低減され、収束次数$R_1 \approx 1.2$となり、滑らかな領域での中程度の改善が確認された。
• 方形波（S）では、滑らかな極値が存在しないため$L_1$誤差の低減は顕著ではなく、$R_1 \approx 0.8$となり、$2^{\text{nd}}$-次挙動と整合的であった。
• 制限定数$C_{VL}$および$C_{PPM}$に対して極値保存型制限子は感度が低く、$C$を1.25から5に変更しても$L_\infty$誤差はわずか4%しか変化しなかった。
• 非線形系では、Shu-Osher問題において極値保存型制限子はWENO3およびWENO5と同等の解像度を示し、高精度な性能を発揮した。
• アダプティブメッシュリファインメントを用いたChomboにおけるWoodward-Colellaのラムプ問題の数値結果から、極値保存型制限子の多変数衝撃捕捉におけるロバスト性と安定性が確認された。