[论文解读] Uniformly accurate exponential-type integrators for klein-gordon equations with asymptotic convergence to classical splitting schemes in the nonlinear schroedinger limit

Simon Baumstark, Erwan Faou+1 more|arXiv (Cornell University)|Jun 15, 2016

Numerical methods for differential equations参考文献 12被引用 4

一句话总结

本论文提出了一类用于克莱因-戈登方程的统一精确指数型积分格式，可在所有参数范围内保持一阶与二阶收敛性——包括相对论性（$c=1$）和高度振荡的非相对论性（$c\gg1$）情形——且无需渐近展开或依赖于$c$的步长限制。当$c\to\infty$时，这些格式在非线性薛定谔方程极限下收敛于经典的李与施朗格分裂格式，且在最小正则性假设下实现统一精度。

ABSTRACT

We introduce efficient and robust exponential-type integrators for Klein-Gordon equations which resolve the solution in the relativistic regime as well as in the highly-oscillatory non-relativistic regime without any step-size restriction, and under the same regularity assumptions on the initial data required for the integration of the corresponding limit system. In contrast to previous works we do not employ any asymptotic/multiscale expansion of the solution. This allows us derive uniform convergent schemes under far weaker regularity assumptions on the exact solution. In particular, the newly derived exponential-type integrators of first-, respectively, second-order converge in the non-relativistic limit to the classical Lie, respectively, Strang splitting in the nonlinear Schr{ö}dinger limit.

研究动机与目标

• 解决在高度振荡的非相对论性参数范围（$c\gg1$）内求解克莱因-戈登方程的数值挑战，且无需步长限制。
• 克服现有方法依赖渐近展开或多尺度展开并受$c$相关稳定性约束的局限性。
• 开发在时间上统一收敛的积分格式，并在非线性薛定谔方程极限下收敛于经典分裂格式。
• 在弱于以往统一精度格式所要求的正则性假设下，实现最优收敛速率。
• 确保在$c\to\infty$极限下，一阶与二阶格式均与经典分裂格式（李与施朗格）保持一致。

提出的方法

• 提出克莱因-戈登方程在二阶形式下的指数型积分格式，避免对解进行渐近展开。
• 基于方程结构提出一种新颖的分裂策略，通过指数积分格式直接处理高度振荡项。
• 推导出一阶与二阶格式（公式28与81），在$H^7$与$H^6$正则性假设下，对所有$c\geq1$实现时间上的统一收敛。
• 证明当$c\to\infty$时，这些格式收敛于非线性薛定谔极限系统（2）的经典李与施朗格分裂格式。
• 采用傅里叶伪谱方法进行空间半离散化，使用$K=2^{10}$个模式以保证空间谱精度。
• 采用指数积分器的改进形式，通过$c$相关的色散关系$\langle\nabla\rangle_c = (I - c^{-2}\Delta)^{1/2}$显式引入$c$依赖性。

实验结果

研究问题

• RQ1 能否为克莱因-戈登方程构造出在所有$c$值（包括非相对论极限）下均保持统一收敛性的指数型积分格式？
• RQ2 这些格式在$c\to\infty$时是否收敛于非线性薛定谔方程极限下的经典李与施朗格分裂格式？
• RQ3 能否在弱于多尺度或渐近展开方法所需正则性假设下实现统一收敛？
• RQ4 在$c\to\infty$参数范围内，指数积分格式与经典分裂格式之间存在何种关系？
• RQ5 在不同参数范围内，这些格式在数值上表现出怎样的收敛阶次与$c$无关性？

主要发现

• 一阶指数积分格式（28）在$H^1$范数下对所有$c$实现$\mathcal{O}(\tau)$的统一收敛性，适用于$H^7$正则解，且满足$\sup_t \|z(t)\|_{H^7} + c^{-2}\|\partial_t z(t)\|_{H^6} \leq 1$。
• 二阶指数积分格式（81）在足够光滑的$H^7$解下，对所有$c$实现$\mathcal{O}(\tau^2)$的统一收敛性。
• 二阶格式（81）在$c\to\infty$时收敛于非线性薛定谔极限系统（2）的经典施朗格分裂格式，误差满足$\|u_*^n - u_{*,\infty}^n\|_r \leq c^{-1}k_r(M_3,T)$，条件为$\tau c \geq 1$。
• 数值实验验证了在$c = 1, 5, 10, \dots, 10000$下$H^1$范数中的一阶与二阶收敛性，斜率与理论阶数一致。
• 即使在$c=10000$时，格式仍保持稳定与精确，表明在高度振荡参数范围内具有鲁棒性，且无需步长限制。
• 该方法避免了高斯奇型格式的$c^4\tau^2$全局误差，后者要求$c^2\tau < 1$，从而克服了先前方法的关键局限。