[論文レビュー] A MacCormack Method for Complete Shallow Water Equations with Source Terms

Eric Ngondiep, Alqahtani T. Rubayyi+1 more|arXiv (Cornell University)|Mar 20, 2019

Computational Fluid Dynamics and Aerodynamics参考文献 25被引用数 5

ひとこと要約

本稿では、源項を含む1次元完全な浅い水方程式に対して、2次精度のMacCormack差分スキームを提示する。分数ステップ法を用いて剛性のある摩擦勾配を扱い、振動を抑えるために上流化を適用する。この手法は、適切な時間刻みとメッシュサイズの制約下で2次収束を達成し、解析的および数値的比較によって検証された安定性と精度を示す。

ABSTRACT

In the last decades, more or less complex physically-based hydrological models, have been developed to solve the shallow water equations or their approximations using various numerical methods. The MacCormack method was developed for simulating overland flow with spatially variable infiltration and microtopography using the hydrodynamic flow equations. The basic MacCormack scheme is enhanced when it uses the method of fractional steps to treat the friction slope or a stiff source term and to upwind the convection term in order to control the numerical oscillations and stability. In this paper we describe, the MacCormack scheme for 1D complete shallow water equations with source terms, analyze the stability condition of the method and we provide the convergence rate of the algorithm. This work improves some well known results deeply studied in the literature which concern the Saint-Venant problem and it represents an extension of the time dependent shallow water equations without source terms. The numerical evidences consider the rate of convergence of the method and compares the numerical solution respect to the analytical one.

研究の動機と目的

• 源項（摩擦および底面勾配効果を含む）を含む1次元完全な浅い水方程式に対する堅牢な数値スキームの開発。
• 基本的なMacCormackスキームに分数ステップ法と上流化を組み合わせ、安定性を向上させ、数値的振動を低減する。
• フーリエ解析を用いて修正されたスキームの安定性条件を分析し、数値的に収束率を検証する。
• 時間依存降雨超過を伴う河川流の解析解に対して、手法の妥当性を検証する。
• 微小な水深および急勾配が存在する状況下でも、スキームがバランスの取れた状態を維持し、正値性を保つことを保証する。

提案手法

• スキームは、時間方向に2次テイラー級数、空間方向に4次精度を用いた2段階の予測・修正アプローチを採用する。
• 分数ステップ法を用いて源項（特に摩擦勾配）を分離し、剛性を扱う。
• 対流項に上流化を適用して数値的振動を抑制し、安定性を向上させる。
• 系を保存則形式に再定式化する：∂W/∂t + ∂F/∂x = S、ここで W = (A, Q)、F = (Q, gA²/(2T) + Q²/A)、S = (r, gA(S₀ - S_f))。
• 安定性はフーリエ解析により分析され、Δt と Δx に対するCFLに類似した制限が得られる。
• 収束性は、メッシュを細分化した際の断面積 A および流量 Q の L² 誤差を用いて数値的に評価する。

実験結果

リサーチクエスチョン

• RQ1 修正されたMacCormackスキームは、源項を含む1次元浅い水方程式に対して2次収束を維持するか？
• RQ2 分数ステップ法と上流化の導入が、MacCormackスキームの安定性および精度に与える影響は何か？
• RQ3 メッシュサイズおよび時間刻みが、数値解の安定性および収束性に与える影響は何か？
• RQ4 微小な水深領域において、スキームが正値性を保ち、不自然な振動を回避できるか？
• RQ5 時間依存降雨超過ケースにおける数値解は、解析解とどの程度一致するか？

主な発見

• MacCormackスキームは、表1に示すL²誤差ノルムによって確認されたように、断面積 A および流量 Q の両方で2次収束を達成する。
• A の収束率は約2.0であり、連続するメッシュの細分化において、2.0000、2.0645、1.9787 の値を示す。
• テストケースにおいて、流量 Q の数値解は解析解と完全に一致し、L²誤差はゼロである。
• 安定性は条件的である：Δt と Δx が導出された安定性制限式 (101) を満たす場合にのみ、スキームは安定する。
• 数値解は時間経過とともに増大せず、解析解に収束するため、適切なパrameter選択下で堅牢であることが示された。
• 上流化と分数ステップ法を適用することで、微小な水深領域においても、数値的振動が効果的に制御され、正値性が維持される。