[论文解读] Spectrally accurate fully discrete schemes for some nonlocal and nonlinear integrable PDEs via explicit formulas

Yvonne Alama Bronsard, Xi Chen+1 more|arXiv (Cornell University)|Jan 1, 2024

Differential Equations and Numerical Methods被引用 1

一句话总结

本文提出了一类谱精度的全离散数值格式，用于环面上的 Benjamin–Ono 方程、Calogero–Sutherland DNLS 方程和立方 Szegő 方程，利用其 Lax 对结构导出的显式解公式。该格式在空间上实现谱收敛，对于 $H^s(\mathbb{T})$ 初始数据（$s > 1$），$L^2$ 范数下的误差以 $K^{-s+1}$ 速率衰减，误差常数随最终时间线性增长，从而实现计算成本恒定的长时间模拟。

ABSTRACT

We construct fully-discrete schemes for the Benjamin-Ono, Calogero-Sutherland DNLS, and cubic Szegő equations on the torus, which are $ extit{exact in time}$ with $ extit{spectral accuracy}$ in space. We prove spectral convergence for the first two equations, of order $K^{-s+1}$ in $L^2$ norm for initial data in $H^s(\mathbb T)$, $s>1$, with an error constant depending $ extit{linearly}$ on the final time instead of exponentially. These schemes are based on $ extit{explicit formulas}$, which have recently emerged in the theory of nonlinear integrable equations. Numerical simulations show the strength of the newly designed methods both at short and long time scales, thanks to the remarkable fact that the computational cost of the method is independent of the final time. These schemes open doors for the understanding of the long-time dynamics of integrable equations.

研究动机与目标

• 为三个关键的非局部、非线性可积 PDE 建立全离散数值格式：Benjamin–Ono 方程、Calogero–Sutherland DNLS 方程和立方 Szegő 方程。
• 利用源自 Lax 对结构的显式解公式，在空间上实现谱精度，在时间上实现精确性。
• 证明严格的收敛性界，误差随最终时间线性增长而非指数增长，从而支持长时间模拟。
• 弥合可积系统解析解公式的最新进展与长期动力学计算方法之间的鸿沟。

提出的方法

• 利用基于环面上 Benjamin–Ono、CS-DNLS 和立方 Szegő 方程 Lax 对的显式解公式。
• 通过将谱空间离散化与基于显式公式的精确时间积分相结合，构造全离散格式。
• 使用截断的 Lax 算子 $L_{u_0,K}$ 定义离散演化算子 $e^{itA_K}$，确保谱精度。
• 基于全算子与截断算子之间的差异进行误差分析，将误差以 $K^{-s}$ 的形式有界。
• 利用酉演化和投影算子控制范数，并在 $H^s$ 空间中推导稳定性界。
• 通过归纳法和涉及 $\|u_0\|_{H^s}$ 及时间相关常数的范数估计建立收敛性。

实验结果

研究问题

• RQ1 能否基于显式解公式，为非局部和非线性可积 PDE 构建谱精度的全离散格式？
• RQ2 对于 $H^s(\mathbb{T})$ 初始数据（$s > 1$），此类格式在 $L^2$ 范数下的收敛速率如何？误差如何依赖于最终时间？
• RQ3 误差常数能否随时间线性增长而非指数增长，从而实现计算成本恒定的长时间模拟？
• RQ4 在时间独立成本下，这些格式在短时间和长时间尺度上的数值表现如何？

主要发现

• 所提格式在空间上实现谱收敛，对于 $H^s(\mathbb{T})$ 初始数据（$s > 1$），$L^2$ 范数下的误差以 $K^{-s+1}$ 速率衰减，其中 $K$ 为谱截断参数。
• 收敛界中的误差常数随最终时间 $t$ 线性增长，而非指数增长，相较于标准方法有显著改进。
• 该格式的计算成本与最终时间无关，从而可高效模拟可积 PDE 的长时间动力学。
• 数值模拟验证了格式在短时间和长时间尺度下的鲁棒性与高精度。
• 对于小初始数据，误差界 $\|u - u_K\|_{H^r} \leq C_5(1 + tK)^{2K^{2r-2s}}$ 中的常数 $C_5$ 随 $\|u_0\|_{H^s}$ 线性趋于零。
• 在短时间 $t = O(K^{-1})$ 下，该方法实现了最优衰减速率 $K^{-s+r}$，表明在早期动力学中具有高精度。