[论文解读] Resonances in long time integration of semi linear Hamiltonian PDEs

Erwan Faou, Benoît Grébert|ArXiv.org|Apr 9, 2009

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

一句话总结

本文表明，对于半线性哈密顿PDE的标准辛积分器，若不截断线性算子的高频分量，可能因数值共振而失效。通过引入频率截断并采用分裂格式，作者在严格满足CFL型条件的前提下，证明了作用量与Sobolev范数的长时间保持性，从而在所有时间步长下避免了虚假共振。

ABSTRACT

We consider a class of Hamiltonian PDEs that can be split into a linear unbounded operator and a regular non linear part, and we analyze their numerical discretizations by symplectic methods when the initial value is small in Sobolev norms. The goal of this work is twofold: First we show how standard approximation methods cannot in general avoid resonances issues, and we give numerical examples of pathological behavior for the midpoint rule and implicit-explicit integrators. Such phenomena can be avoided by suitable truncations of the linear unbounded operator combined with classical splitting methods. We then give a sharp bound for the cut-off depending on the time step. Using a new normal form result, we show the long time preservation of the actions for such schemes for all values of the time step, provided the initial continuous system does not exhibit resonant frequencies.

研究动机与目标

• 识别为何标准辛积分器（如中点法则和IMEX格式）在半线性哈密顿PDE的长时间模拟中可能因数值共振而失效。
• 表明数值共振源于时间步长与本征频率之间的非线性关系，特别是当 $ h ilde{ heta} o 2 ilde{ heta} ilde{ heta} ilde{ heta} $ 且 $ ilde{ heta} eq 0 $ 时，导致守恒性质的丧失。
• 建立在应用分裂方法前对线性算子的高频模态进行截断，可消除这些虚假共振。
• 推导出截断频率的精确边界，其依赖于时间步长 $ h $，以确保长时间稳定性与作用量保持。
• 证明在原始系统满足非共振条件时，修正后的格式可保持物理守恒律，且不引入人工共振。

提出的方法

• 在Sobolev空间框架下，将哈密顿PDE形式化为 $ H = H_0 + P $，其中 $ H_0 $ 为无界线性算子，$ P $ 为光滑非线性项。
• 在线性算子 $ H_0 $ 的谱上引入频率截断 $ K $，令 $ ilde{ heta}_j = 0 $ 对于 $ |j| > Kh^{-1} $，以消除高频贡献。
• 应用形式为 $ ilde{ heta}_K(hH_0) imes ext{flow of } P $ 的分裂格式，其中 $ ilde{ heta}_K $ 为线性流的截断版本。
• 利用正规形变换控制形如 $ |e^{ih ilde{ heta}^h(oldsymbol{j})} - 1| $ 的小分母，通过非共振条件确保有界性。
• 推导同态方程 $ (e^{ih ilde{ heta}^h(oldsymbol{j})} - 1) ilde{ heta}_{oldsymbol{j}} + hZ_{oldsymbol{j}} = hG_{oldsymbol{j}} $ 的估计，其界依赖于 $ ilde{ heta}^h(oldsymbol{j}) $ 及 $ oldsymbol{j} $ 的指标结构。
• 证明正规形中的余项为 $ ilde{ heta}^{r+1} $ 阶，从而确保初始数据大小为 $ ilde{ heta} $ 时，系统在 $ ilde{ heta}^{-r} $ 时间尺度内保持长期稳定。

实验结果

研究问题

• RQ1 标准辛积分器（如中点法则或IMEX格式）能否在半线性哈密顿PDE的长时间积分中避免数值共振？
• RQ2 当时间步长 $ h $ 未受CFL条件约束时，为何辛格式中会出现数值共振？
• RQ3 对线性算子 $ H_0 $ 的高频模态进行截断如何影响分裂格式的长时间行为？
• RQ4 为确保作用量与Sobolev范数的长时间保持，截断频率 $ K $ 的精确边界是什么（以时间步长 $ h $ 表示）？
• RQ5 能否使用正规形论证来控制辛积分器中由非线性相位关系引发的小分母？

主要发现

• 当 $ h ilde{ heta}^h(oldsymbol{j}) o 2 ilde{ heta} ilde{ heta} ilde{ heta} $ 且 $ ilde{ heta} eq 0 $ 时，即使原始系统为非共振，标准辛积分器（如中点法则和IMEX格式）仍可能发生数值共振。
• 对完整系统应用中点法则可能破坏物理守恒律，或因人工共振而意外保持作用量，具体取决于 $ h $。
• 在频率 $ |j| > Kh^{-1} $ 处对线性算子 $ H_0 $ 进行截断，可消除虚假共振并恢复长时间稳定性。
• 推导出精确边界 $ K hicksim h^{-1} $，使得当 $ K hicksim h^{-1} $ 时，该格式可在初始数据大小为 $ ilde{ heta} $ 的 $ ilde{ heta}^{-r} $ 时间尺度内保持作用量与Sobolev范数。
• 当 $ |e^{ih ilde{ heta}^h(oldsymbol{j})} - 1| hicksim h ilde{ heta}^h(oldsymbol{j}) $ 时，正规形的同态方可求解，该条件在截断条件与非共振假设下成立。
• 正规形中的余项为 $ ilde{ heta}^{r+1} $ 阶，且变换保持在 $ ilde{ heta}_{oldsymbol{j}} otin ilde{ heta}_{oldsymbol{j}} $ 类别中，从而确保长时间控制。