[论文解读] Error estimates of finite difference schemes for the Korteweg-de Vries equation

Clémentine Courtès, Fré́dé́ric Lagoutière+1 more|arXiv (Cornell University)|Dec 6, 2017

Advanced Mathematical Physics Problems参考文献 44被引用 24

一句话总结

该论文为求解Korteweg-de Vries（KdV）方程的有限差分格式建立了收敛速率，采用显式Rusanov格式处理非线性通量项，采用四点θ-格式处理色散项。当 $ \theta \geq \frac{1}{2} $ 时，在双曲型CFL条件下，对 $ s \geq 6 $ 的初值数据，证明了在 $ H^s(\mathbb{R}) $ 范数下的一阶收敛率；当 $ \theta < \frac{1}{2} $ 时，采用“Airy”型CFL条件，对 $ s \geq \frac{3}{4} $ 的粗糙初值数据，收敛阶数降低。数值结果表明，当 $ s \geq 3 $ 时收敛阶数达到最优。

ABSTRACT

This article deals with the numerical analysis of the Cauchy problem for the Korteweg-de Vries equation with a finite difference scheme. We consider the Rusanov scheme for the hyperbolic flux term and a 4-points $\ heta$-scheme for the dispersive term. We prove the convergence under a hyperbolic Courant-Friedrichs-Lewy condition when $\ heta\\geq \\frac{1}{2}$ and under an "Airy" Courant-Friedrichs-Lewy condition when $\ heta&lt;\\frac{1}{2}$. More precisely, we get the first order convergence rate for strong solutions in the Sobolev space $H^s(\\mathbb{R})$, $s \\geq 6$ and extend this result to the non-smooth case for initial data in $H^s(\\mathbb{R})$, with $s\\geq \\frac{3}{4}$ , to the price of a loss in the convergence order. Numerical simulations indicate that the orders of convergence may be optimal when $s\\geq3$.

研究动机与目标

• 填补有限差分格式应用于KdV方程时缺乏对低正则性初值的严格收敛速率分析的空白。
• 建立依赖于初值正则性的收敛速率，拓展至经典光滑解之外的情形。
• 分析 $ \theta $-格式中时间步长参数 $ \theta $ 对稳定性与收敛性的影响，区分 $ \theta \geq \frac{1}{2} $ 与 $ \theta < \frac{1}{2} $ 的情况。
• 通过在误差估计中利用色散PDE理论中的色散增强平滑效应，弥合数值分析与色散PDE理论之间的鸿沟。
• 为将该方法推广至三阶色散系统（如 $ abcd $-系统）奠定基础。

提出的方法

• 采用显式Rusanov格式对双曲通量项 $ \partial_x(u^2/2) $ 进行离散化，确保在双曲型CFL条件下保持稳定性。
• 采用四点 $ \theta $-格式近似三阶色散项 $ \partial_x^3 u $，其稳定性与精度取决于 $ \theta $ 的取值。
• 通过求和-by-parts技巧与离散分部积分法推导离散能量估计，以控制非线性和色散项。
• 利用离散的Gagliardo-Nirenberg型不等式，将高阶离散导数控制在低阶导数范围内。
• 结合稳定性估计与连续问题的良好适定性结果，以及KdV方程的色散增强平滑效应，控制全局误差。
• 利用离散的 $ L^2 $-守恒律与离散 $ \ell^2 $-范数估计，控制数值解的演化过程。

实验结果

研究问题

• RQ1 当初值数据正则性有限，具体为 $ H^s(\mathbb{R}) $ 且 $ s \geq \frac{3}{4} $ 时，有限差分格式求解KdV方程的收敛速率是多少？
• RQ2 在色散项的 $ \theta $-格式中，$ \theta $ 的选择如何影响格式的稳定性与收敛性？
• RQ3 当 $ \theta < \frac{1}{2} $ 时，收敛所需的CFL条件为何种类型？它与标准双曲型CFL条件有何不同？
• RQ4 能否通过在数值分析中利用色散增强平滑效应，改善或保持粗糙初值数据下的收敛速率？
• RQ5 理论收敛速率是否在数值模拟中得到验证，特别是在 $ s \geq 3 $ 的情况下？

主要发现

• 当 $ \theta \geq \frac{1}{2} $ 且在双曲型CFL条件下，对 $ s \geq 6 $ 的强解，该格式在 $ H^s(\mathbb{R}) $ 范数下实现一阶收敛。
• 对于 $ H^s(\mathbb{R}) $ 中初值数据满足 $ s \geq \frac{3}{4} $ 的情形，由于正则性损失，收敛阶数降低，但收敛性仍可证明。
• 当 $ \theta < \frac{1}{2} $ 时，收敛性在“Airy”型CFL条件下得到证明，该条件比双曲型CFL条件更宽松，但仍能保证稳定性。
• 数值模拟表明，当 $ s \geq 3 $ 时，收敛阶数可能达到最优，表明该理论界在该参数区间内是紧致的。
• 色散增强平滑效应在控制低正则性初值数据的误差中起关键作用，使得尽管初值光滑性有限，仍能实现收敛。
• 该分析具有足够的鲁棒性，可推广至其他三阶色散系统（如 $ abcd $-系统），但此类推广需要更强的连续问题良好适定性结果。