[論文レビュー] A well-balanced reconstruction for wetting/drying fronts

Andreas Bollermann, Guoxian Chen+2 more|arXiv (Cornell University)|Dec 11, 2014

Computational Fluid Dynamics and Aerodynamics参考文献 31被引用数 3

ひとこと要約

本稿では、連続的で区分線形の床地形再構成を用いた、浅水域方程式に対するバランスの取れた正性を保証する有限体積スキームを提示する。主な革新点は、湿り/乾きセルに対して新規の再構成戦略を採用し、『静止湖』および『乾き湖』の定常状態を正確に保持することであり、局所的排水時間におけるフラックス制限により正性を維持することで、時間ステップ制限なしに湿り/乾き前線の安定シミュレーションが可能になる。

ABSTRACT

In this paper, we construct a well-balanced, positivity preserving finite volume scheme for the shallow water equations based on a continuous, piecewise linear discretization of the bottom topography. The main new technique is a special reconstruction of the flow variables in wet-dry cells, which is presented in this paper for the one dimensional case. We realize the new reconstruction in the framework of the second-order semi-discrete central-upwind scheme from (A. Kurganov and G. Petrova, Commun. Math. Sci., 2007). The positivity of the computed water height is ensured following (A. Bollermann, S. Noelle and M. Lukáčová-Medviďová, Commun. Comput. Phys., 2010): The outgoing fluxes are limited in case of draining cells.

研究の動機と目的

• 浅水域シミュレーションにおいて、『静止湖』および『乾き湖』の両定常状態を正確に保持するバランスの取れた有限体積スキームの開発。
• 特に数値オシレートが深さを負に駆り立てる危険性がある場合に、湿り/乾きセルにおける水深の正性を保証すること。
• 乾き領域および衝撃波が乾き領域に衝突する状況下でも、2次精度と頑健性を維持すること。
• 時間ステップを短縮せずに、湿り/乾き前線の安定的かつ正確なシミュレーションを可能にすること。
• 連続的で区分線形の床地形を扱えるように中央上流スキームを拡張し、バランスの取れ方および正性を保つこと。

提案手法

• 湿り/乾きセルにおける流れ変数のための新規再構成技術を提案し、『静止湖』および『乾き湖』の両定常状態を正確に保持するように設計されている。
• 数値的安定性および精度を向上させるために、床地形の連続的で区分線形離散化を用いる。
• 局所的排水時間がグローバル時間ステップより短い場合に、流出フラックスを制限することで水深の正性を強制する。
• Kurganov & Petrova (2007) の2次半離散中央上流有限体積法の枠組み内に本スキームを構築する。
• 再構成により、平衡状態（湿り/乾きの混合状態を含む）において、数値フラックスとソース項が正確にバランスをとる。
• 湿り/乾き前線における特異性を回避するため、Manning摩擦項に対して部分的陰解法を適用し、スキームの不安定化を防ぐ。

実験結果

リサーチクエスチョン

• RQ1 連続的で区分線形の床地形再構成を用いた、バランスの取れた正性を保証するスキームを、浅水域方程式に対して構築可能か？
• RQ2 湿り/乾き前線をどのように再構成すれば、『静止湖』および『乾き湖』の両定常状態を同時に保持できるか？
• RQ3 時間ステップ制限を課さずに、水深の正性を保証することは可能か？
• RQ4 衝撃波が乾き領域に衝突する場合や、強い摩擦項が存在する場合、本スキームはどのように性能を発揮するか？
• RQ5 提案された再構成は、非平衡状態の非定常流れにおいても2次精度およびバランスの取れ方を維持するか？

主な発見

• 特別に設計された再構成のおかげで、混合湿り/乾き領域においても『静止湖』および『乾き湖』の両定常状態を正確に保持する。
• 本スキームは空間的に2次精度を維持し、すべてのテスト例で実験的収束次数を達成する。
• 局所的排水時間におけるフラックス制限により、湿り/乾き前線全域で水深の正性が保証され、時間ステップ短縮の必要がない。
• Manning摩擦を伴う実験的ダムブレイク（三角形の盛り上がり上を通過）のシミュレーションにおいて、7か所の測定ゲージ点での水深が高精度で測定値と一致した。
• 乾き領域に衝撃波が衝突する場合でも本スキームは頑健であり、Manning摩擦項の特異性に対しても効果的に対処できる。
• 事例研究により、従来版のスキームが平衡状態に収束しないことが示され、新規再構成の必要性が裏付けられた。