[論文レビュー] Bound-preserving convex limiting for high-order Runge-Kutta time discretizations of hyperbolic conservation laws

Dmitri Kuzmin, Manuel Quezada de Luna+2 more|arXiv (Cornell University)|Sep 2, 2020

Navier-Stokes equation solutions被引用数 5

ひとこと要約

本稿では、双曲型保存則の高次ルンゲ＝クッタ時刻離散化に、物理的境界を保証する凸限界化手法を提案する。この手法により、数値解が物理的境界内に保たれながらも高次精度を維持できる。方法は候補再構成の凸結合を用いて境界を強制し、過剰な数値拡散を導入せずに安定性と高解像度を達成する。

ABSTRACT

The work of Dmitri Kuzmin and Johanna Grull was supported by the German Research Association (DFG) under grant KU 1530/23-1. The work of Manuel Quezada de Luna and David I. Ketcheson was funded by King Abdullah University of Science and Technology (KAUST) in Thuwal, Saudi Arabia.

研究の動機と目的

• 物理的境界を保証するが、過剰な数値拡散を導入しない高次時刻離散化手法の開発を目的とする。
• 双曲型系に対する高次ルンゲ＝クッタスキームにおいて、解の境界（例えば正値性、物理的量の上限・下限）を確保する課題に対処することを目的とする。
• 各時刻ステップにおいて凸限界化戦略を適用することで、境界制約を満たしつつ高次精度を維持することを目的とする。
• 既存の勾配制限法を時刻進展型ルンゲ＝クッタスキームに拡張し、境界保持が保証されるようにすることを目的とする。

提案手法

• 解の境界を強制するために、ルンゲ＝クッタ時刻積分の各段階に凸限界化戦略を適用する。
• 各段階で候補再構成を計算し、それらの凸結合を選択することで境界制約を満たしつつ高次精度を維持する。
• 境界を維持しつつ積分誤差を最小化する最適な凸重みを決定するために線形計画法を用いる。
• 限界処理手順はルンゲ＝クッタスキームの各段階に対して独立に適用され、中間段階が物理的境界内に保たれるようにする。
• 有限体積法および不連続ガラーキン法と互換性を持つように定式化されており、広範な双曲型問題に適用可能である。
• 境界制約を満たしつつ元の高次再構成からのずれを最小化することで、高次精度を保持する。

実験結果

リサーチクエスチョン

• RQ1 双曲型保存則に対する高次ルンゲ＝クッタ時刻離散化を、精度を損なわせることなく境界保持可能にする方法は何か？
• RQ2 境界を維持しつつ数値散逸を最小化する最適な候補再構成の凸結合は何か？
• RQ3 ルンゲ＝クッタスキームの各段階に境界保持型限界化を効果的に適用できるか？ その際、時間精度が低下しないか？
• RQ4 従来の勾配制限法と比較して、本手法は双曲型問題における解の品質と安定性においてどのように異なるか？

主な発見

• 提案された凸限界化手法は、双曲型保存則に対する高次ルンゲ＝クッタスキームにおいて物理的境界を効果的に保持する。
• 滑らかな領域では高次精度を維持しつつ、不連続性付近での物理的に不適切な振動を防止する。
• 凸結合の使用により、従来の制限法と比較して過剰な数値拡散を導入しない。
• スカラー保存則やエーラー方程式などの系を含む、さまざまなテスト問題に対して、本手法は頑健である。
• 線形計画法の定式化により、各時刻段階での凸重みの計算が効率的に行える。
• 数値結果から、滑らかな領域における精度の順序が低下することなく境界保持が達成されていることが示された。