[论文解读] On the Rate of Error Growth in Time for Numerical Solutions of Nonlinear Dispersive Wave Equations

Hendrik Ranocha, Manuel Quezada de Luna+1 more|arXiv (Cornell University)|Feb 15, 2021

Advanced Mathematical Physics Problems参考文献 83被引用 13

一句话总结

本文研究了非线性色散波方程数值解中的时变误差增长行为，通过大量数值实验表明，保守数值格式——尤其是保持非线性不变量的格式——表现出线性误差增长，而非保守格式则表现出二次误差增长。该研究将此行为扩展至一大类方程，包括非色散双曲系统和二维浅水方程，其中保守方法在精度和长期稳定性方面始终优于非保守方法。

ABSTRACT

We study the numerical error in solitary wave solutions of nonlinear dispersive wave equations. A number of existing results for discretizations of solitary wave solutions of particular equations indicate that the error grows quadratically in time for numerical methods that do not conserve energy, but grows only linearly for conservative methods. We provide numerical experiments suggesting that this result extends to a very broad class of equations and numerical methods.

研究动机与目标

• 研究在孤立波解中观察到的保守数值方法的有利线性误差增长是否可推广至更广泛的非线性色散波方程类。
• 检验当守恒不变量为非二次形式时，保守格式是否仍保持线性误差增长。
• 探究该行为是否存在于不具有经典孤立波解与平移对称性的非色散一阶双曲系统中。
• 评估保守方法在二维系统（如具有变底地形的浅水方程）中的性能。
• 证明保守方法不仅在渐近意义下，甚至在短时间尺度内，也能产生显著更小的全局误差，即使理论框架中未严格强制能量或质量守恒。

提出的方法

• 采用求和按部分（SBP）算子结合同时近似项（SATs）构造可证明稳定且保守的空间离散化。
• 使用松弛Runge-Kutta时间积分方法，实现显式或低成本隐式时间推进，同时保持离散不变量。
• 基于傅里叶谱近似应用谱配点方法以实现高阶空间精度。
• 为各类方程（包括Fornberg-Whitham、Camassa-Holm、Degasperis-Procesi、Holm-Hone、BBM-BBM和浅水方程）实现保守与非保守数值格式的变体。
• 进行长时间模拟，比较保守与非保守格式之间的时间误差增长速率。
• 通过数值实验评估误差增长，即随时间计算数值解与参考解之间L2范数的差异。

实验结果

研究问题

• RQ1 已知色散方程（如KdV、NLS）中保守格式所表现出的随时间线性误差增长，是否可推广至更广泛的非线性色散波方程类？
• RQ2 线性误差增长是否依赖于守恒不变量的具体形式，还是即使在非二次不变量下也成立？
• RQ3 在不具有经典孤立波解的非色散一阶双曲系统中，保守数值方法是否仍能实现更优的误差行为？
• RQ4 在二维系统（如具有变底地形的浅水方程）中，保守方法的优势是否依然存在？
• RQ5 空间分辨率和数值误差对复杂现实系统（如浅水方程）中观测到的误差增长速率有何影响？

主要发现

• Fornberg-Whitham、Camassa-Holm、Degasperis-Procesi、Holm-Hone和BBM-BBM方程的保守数值格式表现出随时间的线性误差增长，而非保守格式则表现出二次误差增长。
• 误差增长改善不仅限于二次不变量；本研究证实，即使守恒量为非线性形式，线性误差增长依然成立。
• 在一维变系数p-系统中，该系统缺乏标准孤立波对称性，保守格式仍显著降低误差增长，表明该现象可超越平移不变解的范围。
• 在具有变底地形的二维浅水模拟中，保守方法在计算成本相近时，精度比非保守方法高出一个数量级，尽管不存在精确孤立波解。
• 即使在短时间尺度内，保守格式产生的全局误差也显著小于非保守格式，表明其在全时间尺度上均具有一致优势。
• 本研究提供了计算证据，表明质量守恒是观察到线性误差增长的必要条件，如一例违反质量守恒的情况导致误差增长转为二次。