[论文解读] Hodograph Method and Numerical Integration of Two Hyperbolic Quasilinear Equations. Part I. The Shallow Water Equations

E. V. Shiryaeva, M. Yu. Zhukov|arXiv (Cornell University)|Oct 10, 2014

Nonlinear Waves and Solitons参考文献 6被引用 3

一句话总结

本文提出了一种基于速度势图的数值方法，用于求解两类双曲型拟线性偏微分方程，特别是浅水方程，通过将原始的柯西问题转化为常微分方程（ODE）系统。该方法利用守恒律和显式的黎曼-格林函数，实现精确的非近似数值积分，从而能够构建多值解，并以高精度捕捉波破碎现象。

ABSTRACT

In paper [S.I. Senashov, A. Yakhno. 2012. SIGMA. Vol.8. 071] the variant of the hodograph method based on the conservation laws for two hyperbolic quasilinear equations of the first order is described. Using these results we propose a method which allows to reduce the Cauchy problem for the two quasilinear PDE's to the Cauchy problem for ODE's. The proposed method is actually some similar method of characteristics for a system of two hyperbolic quasilinear equations. The method can be used effectively in all cases, when the linear hyperbolic equation in partial derivatives of the second order with variable coefficients, resulting from the application of the hodograph method, has an explicit expression for the Riemann-Green function. One of the method's features is the possibility to construct a multi-valued solutions. In this paper we present examples of method application for solving the classical shallow water equations.

研究动机与目标

• 开发一种无需近似的数值方法，用于求解双曲型拟线性PDE，利用守恒律。
• 实现多值解的构建，尤其适用于涉及波破碎的问题。
• 通过将PDE系统简化为ODE的柯西问题，提供传统有限差分或有限体积方法的替代方案。
• 在具有周期性初始条件的典型浅水方程上验证该方法。
• 展示该方法在捕捉复杂波动力学（如波破碎和液滴形成）方面的有效性。

提出的方法

• 将速度势图方法应用于两个在黎曼不变量下的双曲型拟线性一阶PDE，利用守恒律推导出一个二阶线性双曲型PDE。
• 该方法需要所推导的二阶PDE的黎曼-格林函数具有显式解析表达式，以实现隐式解的表示。
• 将原始PDE的柯西问题重新表述为速度势图变量中的ODE系统，随后进行数值求解。
• 该变换依赖于由守恒律导出的特征速度和通量密度，边界条件在特征线上指定。
• 使用标准方法（如龙格-库塔法）对ODE进行数值求解，避免了直接求解非线性超越方程。
• 采用浅水方程的周期性初始数据对方法进行验证，结果在多个时间步长下可视化。

实验结果

研究问题

• RQ1 该速度势图方法能否有效应用于将双曲型拟线性PDE系统转化为可解的ODE系统以实现数值积分？
• RQ2 该方法在浅水流动中捕捉波破碎和多值解的精度如何？
• RQ3 与从隐式解中导出的非线性超越方程相比，求解ODE在计算上具有何种优势？
• RQ4 该方法能否应用于其他具有已知黎曼-格林函数的双曲系统，如气体动力学或孤子方程？
• RQ5 该方法在模拟浅水系统中周期性初始条件的长时间演化时表现如何？

主要发现

• 该方法成功地将浅水方程转化为ODE的柯西问题，实现了无近似、高精度的数值积分。
• 在 t = 5.855 时，数值结果中清晰观察到波破碎，表明该方法具备捕捉间断性的能力。
• 在 t = 7.633 和 t = 8.008 时，水表面形成气泡或液滴的视觉迹象明显，表明该方法能捕捉复杂的自由表面动力学。
• 通过使用超几何函数和椭圆函数，即使在负自变量情况下也能通过复数运算或实变量恒等式实现黎曼-格林函数的精确计算。
• 与直接求解超越方程相比，该方法更具优势，因为ODE求解器在鲁棒性和效率上优于牛顿型方法。
• 该方法可推广至其他具有显式黎曼-格林函数的系统，包括气体动力学、孤子方程和色谱模型。