[論文レビュー] A Simple and Effective High-Order Shock-Capturing Limiter for Discontinuous Galerkin Methods

Scott Moe, James A. Rossmanith+1 more|arXiv (Cornell University)|Jul 10, 2015

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

ひとこと要約

本稿では、バート・ジスペルセンおよびチャン・シュウの最大原理保存フレームワークのアイデアを組み合わせた、不連続ガラーキン（DG）法における単純で高次の衝撃捕捉型リミッターを提案する。局所的に隣接要素のデータを用いて解変数を制限し、高次補正を最小限に減衰させることで、衝撃近傍での振動を効果的に抑止しながら、滑らかな領域では高次精度を維持する。この手法は、構造化および非構造化メッシュ上での1次元および2次元ベンチマークで検証された。

ABSTRACT

The discontinuous Galerkin (DG) finite element method when applied to hyperbolic conservation laws requires the use of shock-capturing limiters in order to suppress unphysical oscillations near large solution gradients. In this work we develop a novel shock-capturing limiter that combines key ideas from the limiter of Barth and Jespersen [AIAA-89-0366 (1989)] and the maximum principle preserving (MPP) framework of Zhang and Shu [Proc. R. Soc. A, 467 (2011), pp. 2752--2776]. The limiting strategy is based on traversing the mesh element-by-element in order to (1) find local upper and lower bounds on user-defined variables by sampling these variables on neighboring elements, and (2) to then enforce these local bounds by minimally damping the high-order corrections. The main advantages of this limiting strategy is that it is simple to implement, effective at shock capturing, and retains high-order accuracy of the solution in smooth regimes. The resulting numerical scheme is applied to several standard numerical tests in both one and two-dimensions and on both Cartesian and unstructured grids. These tests are used as benchmarks to verify and assess the accuracy and robustness of the method.

研究の動機と目的

• 不連続性の近傍で物理的に意味のない振動を回避する、強固で高次の衝撃捕捉型リミッターをDG法に開発すること。
• 衝撃近傍での散乱を最小限に抑えながら、滑らかな領域での高次精度を維持すること。
• リミッティングプロセスに検出を埋め込むことで、明示的な衝撃検出ステップの必要性を排除すること。
• 実装が単純で、構造化および非構造化メッシュの両方で効果的であること。
• MPPフレームワークの適用範囲を正の保存にとどまらず、強い衝撃に対応できるように拡張すること。

提案手法

• リミッターは要素ごとに作用し、隣接要素からの点値を用いて、ユーザー定義変数の局所的な上限および下限を計算する。
• 解の品質を評価する滑らかさ指標を用い、勾配がパラメータαに比例するh^1.5の閾値を超えた場合にのみリミッティングを発動する。
• 高次多項式補正の減衰を最小限に抑え、セル平均を維持することで高次精度を保持する。
• エントロピー層などの複雑な流れの特徴を示す領域では、原始変数の代わりに保存変数のセル平均を用いて境界推定を行う。
• 感度の高い領域（例：前方に突き出したステップの上部）では、修正されたリミッター係数α = 20h^1.5を用いることで、誤ったエントロピー層を低減する。
• SSPRK4時間積分法とLLF数値フラックスを用い、直交座標および非構造メッシュ上にP3要素で実装する。

実験結果

リサーチクエスチョン

• RQ1 明示的な衝撃検出を避けつつ、滑らかな領域での精度を維持する単純で高次のDGリミッターを設計可能か？
• RQ2 バート・ジスペルセンおよびMPP原理の組み合わせは、精度の順序を低下させることなく、衝撃近傍での振動を効果的に抑えるか？
• RQ3 構造化および非構造化メッシュの両方で、顕著な修正なしに安定して適用可能か？
• RQ4 エントロピー層を伴う前方に突き出したステップ問題のような複雑な衝撃構造を、このリミッターはいかに正しく捉えるか？
• RQ5 高勾配領域における境界推定に、原始変数の代わりに保存変数を用いる影響は何か？

主な発見

• 1次元および2次元のエーラー方程式において、衝撃近傍での物理的に意味のない振動が効果的に抑制された。特に挑戦的な前方に突き出したステップ問題においても同様の効果が得られた。
• 前方に突き出したステップ問題において、特にステップ上部の領域でα = 20h^1.5の高分解能リミッター係数を用いることで、誤ったエントロピー層が顕著に低減された。
• 直交座標メッシュのシミュレーションでは顕著なエントロピー層が観察されたが、先行研究で一般的に観測される範囲内に収まり、リミッターの強健性を示した。
• 非構造メッシュの結果は直交座標メッシュの結果と良好に一致し、本手法が非構造メッシュ上でも有効であることを確認した。
• 接触不連続性の鋭い解像度と滑らかな領域の滑らかなプロファイルから、滑らかな領域での高次精度が保持されていることが裏付けられた。
• リミッティングの意思決定が局所的境界推定プロセスに内蔵されているため、追加の衝撃検出ステップは不要である。