[论文解读] The Riemann problem for the shallow water equations with discontinuous topography

Philippe G. LeFloch, Mai Duc Thanh|ArXiv.org|Dec 21, 2007

Navier-Stokes equation solutions参考文献 12被引用 6

一句话总结

本论文为具有不连续地形的浅水方程构造了黎曼问题的显式解，通过引入三种类型的可接受波——膨胀波、激波和静止波，解决了非严格双曲性与非守恒形式的挑战。关键贡献是基于波相互作用分析与（MC）熵准则，提出了一套完全显式且可数值实现的解框架，能够准确逼近变地形条件下的广义初值问题。

ABSTRACT

We construct the solution of the Riemann problem for the shallow water equations with discontinuous topography. The system under consideration is non-strictly hyperbolic and does not admit a fully conservative form, and we establish the existence of two-parameter wave sets, rather than wave curves. The selection of admissible waves is particularly challenging. Our construction is fully explicit, and leads to formulas that can be implemented numerically for the approximation of the general initial-value problem.

研究动机与目标

• 求解非严格双曲性且无完全守恒形式的具有不连续地形的浅水方程黎曼问题。
• 在广义Rankine-Hugoniot条件框架下，定义并表征可接受波——膨胀波、激波与静止波。
• 为所有可能的初始状态建立完整且显式的解结构，包括涉及同一族多波的复杂波序列。
• 基于波相互作用与（MC）熵准则，提供一种完全可数值实现的解框架，适用于有限体积格式。

提出的方法

• 通过分析三种特征族（1-波：膨胀/激波，2-波：膨胀/激波，3-波：静止）构造非严格双曲系统的波集（非波曲线）。
• 应用（MC）熵准则选择可接受的波序列，确保波序正确且物理一致。
• 利用时间无关的光滑解推导静止波关系，满足 [hu] = 0 与 [u²/2 + g(h+a)] = 0，且激波速度为零。
• 通过自相似解与守恒律推导波速与状态的显式公式，实现数值可实施性。
• 根据波行为将初始状态划分为区域 A₁、A₂、A₃，并应用几何与分析条件（如 h ≤ h̄₁）判断波的存在性。
• 通过波相互作用分析验证解，包括1-波或2-波重复出现的情况，以及三股零速波连续出现的情形。

实验结果

研究问题

• RQ1 当系统为非严格双曲性且非守恒形式时，如何求解具有不连续地形的浅水方程黎曼问题？
• RQ2 在广义Rankine-Hugoniot条件下，该系统的可接受波类型（膨胀波、激波、静止波）是什么？
• RQ3 如何构建并验证涉及同一族多波（如两个1-波或两个2-波）的波序列？
• RQ4 在多个波具有零速时，如何确保波序正确且波序列可接受？
• RQ5 如何使解完全显式化并实现于一般初始条件的数值计算中？

主要发现

• 通过自相似性与守恒律确定波速，显式构造了三种波型的黎曼解：1-波（膨胀或激波）、2-波（膨胀或激波）与3-波（静止波）。
• 由于非严格双曲性，存在四波解（如 R₁ ⊕ W₃ ⊕ W₁ ⊕ W₂），展示了复杂的波相互作用。
• 条件 h₂ ≤ h̄₁ 确保序列中2-波速度非正，从而实现解中波序的正确排列。
• 当静止波从 a_L 跳跃至中间水平 a_m，随后为零速 k-激波（k=1 或 3），再跃迁至 a_R 时，存在三股零速波的解。
• 集合 SW(U_L, [a_L, a_R]) 及其变换形式 ŨSW(U_L, [a_L, a_R]) 可用于构造含两股静止波与一股零速1-激波的解。
• 当 ŨSW(U_L, [a_L, a_R]) ∩ W₂^B(U_R) ≠ ∅ 时，解结构 W₃(U_L, U₁) ⊕ S₁(U₁, Ũ₁) ⊕ W₃(Ũ₁, U₂) ⊕ W₂(U₂, U_R) 有效，为 U_L ∈ A₁ 提供完整的波序列。