[论文解读] Semi-explicit solutions to the water-wave dispersion relation and their role in the nonlinear Hamiltonian coupled-mode theory

Theodosios K. Papathanasiou, Christos Papoutsellis+1 more|arXiv (Cornell University)|Feb 22, 2018

Coastal and Marine Dynamics参考文献 85被引用 6

一句话总结

本文提出了一套高度精确的半显式解析公式，用于计算水波色散关系的根，特别是倏逝模态的根，采用混合Picard-Householder迭代算法。这些公式可在三轮以内实现机器精度的根求解，显著提升了哈密顿耦合模理论（HCMT）在任意水深条件下对完全非线性水波模拟的效率与鲁棒性。

ABSTRACT

The Hamiltonian coupled-mode theory (HCMT), recently derived by Athanassoulis and Papoutsellis [1], provides an efficient new approach for solving fully nonlinear water-wave problems over arbitrary bathymetry. In HCMT, heavy use is made of the roots of a local, water-wave dispersion relation with varying parameter, which have to be calculated at every horizontal position and every time instant. Thus, fast and accurate calculation of these roots, valid for all possible values of the varying parameter, are of fundamental importance. In this paper, new, semi-explicit and highly accurate root-finding formulae are derived, especially for the roots corresponding to evanescent modes. The derivation is based on the successive application of a Picard-type iteration and the Householder's root finding method. Explicit approximate formulae of very good accuracy are obtained, which are adequate to support HCMT for many types of applications. In most demanding cases, e.g. very steep, deep-water waves, machine-accurate determination of the required roots is achieved by no more than three iterations, using the explicit forms as initial values. Exploiting this root-finding procedure in the HCMT, results in an efficient, numerical solver able to treat fully nonlinear water waves over arbitrary bathymetry. Applications to demanding nonlinear problems demonstrate the efficiency and the robustness of the present approach.

研究动机与目标

• 解决哈密顿耦合模理论（HCMT）中因需在每个空间和时间点求解水波色散关系而产生的计算瓶颈。
• 克服标准根求解方法（如牛顿-拉夫森法）在深水条件下初始猜测不佳时的不稳定与低效问题。
• 开发一种鲁棒的高阶迭代方法，确保在所有水深条件下仅用极少迭代次数即可实现机器精度的高精度求解。
• 通过提供准确的闭式初始近似值来实现对复杂水深地形下完全非线性水波的稳定、长时间模拟。
• 通过提供可重用、可靠的根求解公式，促进HCMT及类似多模态方法的广泛应用。

提出的方法

• 采用两步法推导局部水波色散关系根的半显式近似：首先使用Picard型迭代生成良好的初始猜测；其次应用Householder的高阶根求解方法。
• 将Householder迭代应用于色散关系方程，该方程在波数 $k_n$ 上为非线性，且频率参数 $\mu(x,t)$ 随空间和时间变化。
• 构建 $k_n(\mathbf{x},t)$ 初始近似的显式闭式表达式，适用于所有水深条件，包括倏逝模态。
• 证明复合算法（Picard后接Householder）保持了Householder的收敛阶，但误差常数显著降低。
• 将所得根求解过程集成到HCMT框架中，其中耦合模系统中的系数 $A$、$B$ 和 $C$ 使用这些精确的 $k_n$ 值计算。
• 在HCMT对完全非线性波问题的模拟中验证该方法，包括孤立波演化、波的散射以及在水深突变区域的波破碎现象。

实验结果

研究问题

• RQ1 如何在所有水深条件下（包括倏逝模态）高效且精确地计算水波色散关系的根？
• RQ2 何种迭代策略可确保在牛顿-拉夫森法在深水中失效时仍能实现收敛与高精度？
• RQ3 结合低阶与高阶迭代方法的混合方法是否相比标准方法具有更优的收敛特性？
• RQ4 半显式初始近似在多大程度上减少了达到机器精度根求解所需的迭代次数？
• RQ5 所提出的根求解方法在多大程度上提升了HCMT在长时间完全非线性水波模拟中的稳定性和效率？

主要发现

• 所提出的二阶与三阶半显式方法在波数计算中可实现机器精度，且最多仅需三轮迭代。
• 尽管三阶方法具有更高的理论收敛阶，但二阶方法在计算效率上表现更优。
• 该方法在所有水深条件下均具有鲁棒性：在牛顿-拉夫森法因初始猜测不当而发散的深水区域，本方法仍能稳定收敛。
• 对 $k_n$ 的闭式初始近似足够精确，可作为迭代求解器的可靠初始点，从而实现稳定长期模拟。
• 采用新根求解流程的HCMT模拟成功再现了复杂非线性现象，包括在水深突变区域的孤立波破碎以及在正弦波纹地形上的布雷达散射。
• 结果表明，该方法可实现高效、降维且高度精确的完全非线性水波模拟，适用于任意水深地形，其精度与直接数值方法相当，但计算成本显著降低。