[论文解读] Symmetric high order Gautschi-type exponential wave integrators pseudospectral method for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime

Yan Wang, Xiaofei Zhao|arXiv (Cornell University)|Nov 4, 2016

Numerical methods for differential equations参考文献 21被引用 10

一句话总结

该论文提出了一种针对非相对论极限 regime 下非线性 Klein-Gordon 方程的对称、高阶 Gautschi 类指数波积分（EWI）傅里叶伪谱方法。该方法实现了任意偶数阶的时间精度和谱空间精度，采用最优网格策略 τ = O(ε²) 和 h = O(1)，在相同 ε 可扩展性条件下，其精度和能量守恒性能优于现有方法如 GIFP 和 RK4FP。

ABSTRACT

A group of high order Gautschi-type exponential wave integrators (EWIs) Fourier pseudospectral method are proposed and analyzed for solving the nonlinear Klein-Gordon equation (KGE) in the nonrelativistic limit regime, where a parameter $0

研究动机与目标

• 解决非线性 Klein-Gordon 方程在非相对论极限 regime 下的数值挑战，其中解表现出波长为 O(ε²) 的高频时间振荡。
• 克服现有指数波积分器的局限性，后者最多只能实现二阶时间精度。
• 开发一种完全显式、时间对称、高阶的数值格式，保持最优的 ε 可扩展性和谱精度。
• 建立严格的误差估计，证明在最优网格策略 τ = O(ε²)，h = O(1) 下方法的收敛性。
• 在精度和长时间能量守恒方面，相较于经典方法如 GIFP 和 RK4FP 展现出优越性能。

提出的方法

• 为指数波积分器中的时间积分提出一族高阶 Gautschi 类求积规则，实现时间上任意偶数阶精度。
• 将这些求积规则与傅里叶伪谱空间离散化相结合，实现空间上的谱精度。
• 确保格式的时间对称性和显式性，以保持 PDE 的结构特性。
• 推导并分析严格的误差界，确认在 0 < ε ≪ 1 条件下，最优网格策略 τ = O(ε²) 和 h = O(1) 的有效性。
• 构建方法时保证在解具有光滑性假设下仍保持高阶精度。
• 利用多尺度分析和渐近分析，证明设计的合理性与在非相对论极限下的收敛性。

实验结果

研究问题

• RQ1 能否为非相对论极限 regime 下的非线性 Klein-Gordon 方程构造出具有任意偶数阶时间精度的高阶指数波积分器？
• RQ2 所提方法是否在 τ = O(ε²) 和 h = O(1) 条件下保持最优的 ε 可扩展性，与时间上的物理波长相匹配？
• RQ3 在相同网格策略下，所提方法与现有经典方法如 GIFP 和 RK4FP 相比，在精度和能量守恒方面表现如何？
• RQ4 在实现高阶精度的同时，能否保持格式的时间对称性和显式结构？
• RQ5 该方法在高度振荡 regime 下的收敛速率和长期能量行为如何？

主要发现

• 4阶-GIFP 和 6阶-GIFP 格式分别实现了 4 阶和 6 阶的时间收敛率，且具有谱空间精度，验证了理论误差估计。
• 该方法在非相对论极限 regime 下保持了最优网格策略 τ = O(ε²) 和 h = O(1)，这是经典方法中最佳的。
• 数值实验表明，在相同 ε 可扩展性条件下，所提方法的计算误差显著小于 GIFP 和 RK4FP。
• 4阶-GIFP 和 6阶-GIFP 的能量误差在整个时间范围内保持有界且极小，而 RK4FP 的能量误差随时间持续增长。
• 该方法表现出优异的长时间能量守恒特性，能量波动保持在机器精度水平。
• 在相同网格条件下，所提方法在精度和稳定性方面均优于经典 Runge-Kutta 方法和现有 EWI 方法。