[論文レビュー] Bound-preserving and entropy-stable algebraic flux correction schemes for the shallow water equations with topography

Hennes Hajduk, Dmitri Kuzmin|arXiv (Cornell University)|Jul 15, 2022

Advanced Numerical Methods in Computational Mathematics被引用数 5

ひとこと要約

本稿では、地形を伴う浅い水方程式に対して、境界を保存し、エントロピー安定な代数的フラックス補正（AFC）スキームを提示する。静止状態におけるwell-balancedness、水深の非負性、エントロピー安定性を保証するため、モノリシックな凸制限戦略を用いる。この手法により、機械精度に近い湖の静止状態の保存が達成され、摩擦に基づくアルゴリズムにより、湿潰れ（wetting-drying）シミュレーションが安定かつ高精度で実行される。標準CFL条件下でも1次元ベンチマークで高い精度と安定性が確認された。

ABSTRACT

A well-designed numerical method for the shallow water equations (SWE) should ensure well-balancedness, nonnegativity of water heights, and entropy stability. For a continuous finite element discretization of a nonlinear hyperbolic system without source terms, positivity preservation and entropy stability can be enforced using the framework of algebraic flux correction (AFC). In this work, we develop a well-balanced AFC scheme for the SWE system including a topography source term. Our method preserves the lake at rest equilibrium up to machine precision. The low-order version represents a generalization of existing finite volume approaches to the finite element setting. The high-order extension is equipped with a property-preserving flux limiter. Nonnegativity of water heights is guaranteed under a standard CFL condition. Moreover, the flux-corrected space discretization satisfies a semi-discrete entropy inequality. New algorithms are proposed for realistic simulation of wetting and drying processes. Numerical examples for well-known benchmarks are presented to evaluate the performance of the scheme.

研究の動機と目的

• 水深の非負性やエントロピー安定性といった重要な物理的性質を保存する、地形を伴う浅い水方程式の数値スキームの開発。
• 特に有限要素離散化を対象として、非一様な双曲型系に源項を含む代数的フラックス補正（AFC）法の拡張。
• 非構造格子上での湖の静止状態におけるwell-balancednessの実現、これは有限要素法における長年の課題である。
• 不純な振動を避けて乾燥領域で安定性を維持する、強力な湿潰れ・乾燥アルゴリズムの設計。
• 境界およびエントロピーの性質を保ちつつ、既存の有限体積法を連続有限要素法へ一般化すること。

提案手法

• ベースライン離散化として、連続的区分線形有限要素法を採用し、特性保存フラックス制限器により高次に拡張する。
• 水深および速度の境界を保証するため、標準CFL条件下で非負性を確保するモノリシックな凸制限（MCL）戦略を適用する。
• フラックス補正された空間離散化は、半離散的エントロピー不等式を満たし、エントロピー安定性を保証する。
• 湖の静止状態におけるwell-balancednessは、地形源項を正確に打ち消すように数値フラックスを修正することにより達成され、機械精度まで保存される。
• 新規の摩擦に基づく湿潰れ・乾燥アルゴリズムを導入し、人工的粘性や複雑な再構成技術に依存しない。
• CFLパラメータν=0.05を用いた適応的時間刻み法により、追加の拡散や段階の繰り返しを伴わず安定性を確保する。

実験結果

リサーチクエスチョン

• RQ1 代数的フラックス補正は、地形を伴う浅い水方程式に拡張可能であり、水深の非負性およびエントロピー安定性を保証できるか？
• RQ2 源項を含む連続有限要素フレームワークにおいて、湖の静止状態におけるwell-balancednessはどのように達成できるか？
• RQ3 不純な振動を避けて乾燥領域で安定性を維持する、強力な湿潰れ・乾燥アルゴリズムは設計可能か？
• RQ4 提案されたAFCスキームは、移動するフロントおよび乾燥状態を含むベンチマーク問題において、精度とシャックルキャプチャ能力を維持するか？
• RQ5 本手法は高次有限要素へ一般化可能であり、科リオリ力や底面摩擦などの追加源項を含める拡張も可能か？

主な発見

• 提案されたAFCスキームは、機械精度に近い湖の静止状態の保存を達成し、非構造格子上でも正確なwell-balancednessを示した。
• 標準CFL条件下で水深の非負性が保証され、人工的粘性やフラックス補正枠組みを超える追加の制限は不要である。
• フラックス補正された空間離散化は半離散的エントロピー不等式を満たし、エントロピー安定性を保証した。
• 堤防破壊問題の数値結果は、参照解と良好に一致しており、速度ピークが従来手法よりわずかに低く抑えられていることから、さらなる精錬の余地があることが示唆された。
• 振動する放物型湖テストでは、MCLおよびMCL-SDEスキームが正確に正確な周期的解を追跡した。低次スケールのプロファイルは正確な解に対してわずかに遅れを示した。
• 摩擦に基づく湿潰れ・乾燥アルゴリズムは、ν=0.05でも安定かつ効果的であり、シミュレーションのクラッシュや不安定化を引き起こす他の手法を上回った。