[论文解读] On the Polysymplectic Integrator for the Short Pulse Equation

Monika E. Pietrzyk, Igor V. Kanatchikov|arXiv (Cornell University)|Dec 30, 2015

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

一句话总结

本论文提出了一种基于 De Donder-Weyl 哈密顿形式的几何多辛积分方法，用于求解短脉冲方程（SPE），在空间和时间方向上均采用中点格式进行离散化。该方法在长传播时间下，精度比伪谱方法高出一个数量级，计算速度提升2.5倍，同时保持能量守恒并确保数值鲁棒性。

ABSTRACT

The polysymplectic analysis of the Short Pulse Equation known in nonlinear optics is used in order to construct a geometric polysymplectic integrator for it. The proposed scheme turns out to be much more effective than other standard integration schemes for nonlinear PDEs, such as the pseudo-spectral integrator. In our numerical experiments the polysymplectic integrator appears to be an order of magnitude more precise and approximately 2.5 times faster at long propagation times than the pseudo-spectral method.

研究动机与目标

• 开发一种用于短脉冲方程（SPE）的几何数值积分方法，以保持其内在的多辛结构。
• 在非线性薛定谔方程因慢包络近似失效而失效的超短光脉冲模拟中，提升数值稳定性和精度。
• 在精度、速度和能量守恒方面，将所提出的多辛积分方法与标准方法（如伪谱积分）进行比较。
• 建立一种鲁棒且高效的长期模拟SPE及其在非线性光学中推广形式的方案。

提出的方法

• 使用一阶拉格朗日函数表述SPE，并通过多动量 $ p^t $ 和 $ p^x $ 推导出 De Donder-Weyl（DW）哈密顿形式。
• 利用 Duffin-Kemmer-Petiau 矩阵 $ \beta^x $ 和 $ \beta^t $ 将DW哈密顿方程表示为矩阵形式，从而实现多辛结构。
• 在空间 $ x $ 和时间 $ t $ 上均采用二阶中点格式对DW方程进行离散化，使用交错网格点表示变量，如 $ \phi_{i,j+1/2} $ 和 $ p^t_{i+1/2,j} $。
• 在每个网格点上，通过求解关于 $ p^t_{i,j+1} $ 的三次方程，随后显式更新 $ \phi_{i,j+1} $ 和 $ p^x_{i,j+1} + p^x_{i,j} $，求解得到的非线性代数方程组。
• 通过 $ u_{i,j+1} = 2p^t_{i,j+1} $ 重构电场振幅，确保与原始SPE中 $ u_{xt} = u + \frac{1}{6}(u^3)_{xx} $ 的一致性。
• 应用周期性或吸收边界条件，初始条件基于精确的Sakovich解获得，并强制能量泛函守恒。

实验结果

研究问题

• RQ1 能否基于 De Donder-Weyl 形式，为短脉冲方程构造一种几何多辛积分方法？
• RQ2 在长时间模拟中，所提出的多辛积分方法在精度和计算效率方面与伪谱方法相比表现如何？
• RQ3 多辛格式是否能保持能量泛函并维持SPE数值解的长期稳定性？
• RQ4 该离散化方案在多大程度上保持了其内在的多辛结构及其相关的守恒律？

主要发现

• 该多辛积分方法的精度约为伪谱方法的10倍，在 $ N = 2^{17} $ 和 $ \Delta t = 0.0001 $ 条件下，$ t = 100 $ 时的误差为 $ \sigma \approx 6.5 \times 10^{-6} $。
• 该方法计算速度约为2.5倍快，相同条件下仅需40,000秒CPU时间，而伪谱方法需100,000秒。
• 该格式在长时间传播过程中保持了能量泛函守恒，相较于标准积分格式展现出更优的鲁棒性与稳定性。
• 关于 $ p^t_{i,j+1} $ 的三次方程（18a）确保了解的连续性，且根的选择过程维持了数值一致性。
• 数值结果证实，多辛积分方法在长时间模拟中仍保持精度与稳定性，经与精确Sakovich解对比验证。
• 该方法在多种离散化步长下均表现有效，随着 $ \Delta x $ 和 $ \Delta t $ 的减小，误差持续减小。