[论文解读] Hamiltonian interpolation of splitting approximations for nonlinear PDEs

Erwan Faou, Benoît Grébert|arXiv (Cornell University)|Dec 15, 2009

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

一句话总结

该论文为具有多项式非线性的非线性哈密顿PDE的分裂格式建立了哈密顿插值，表明经过修正的哈密顿量的数值流在指数时间尺度内与精确流高度吻合。通过利用反向误差分析构造一个考虑谱滤波的修正哈密顿量，该方法在非共振和可积性条件下确保了能量和正则性的长期保持。

ABSTRACT

We consider a wide class of semi linear Hamiltonian partial differential equa- tions and their approximation by time splitting methods. We assume that the nonlinearity is polynomial, and that the numerical tra jectory remains at least uni- formly integrable with respect to an eigenbasis of the linear operator (typically the Fourier basis). We show the existence of a modified interpolated Hamiltonian equation whose exact solution coincides with the discrete flow at each time step over a long time depending on a non resonance condition satisfied by the stepsize. We introduce a class of modified splitting schemes fulfilling this condition at a high order and prove for them that the numerical flow and the continuous flow remain close over exponentially long time with respect to the step size. For stan- dard splitting or implicit-explicit scheme, such a backward error analysis result holds true on a time depending on a cut-off condition in the high frequencies (CFL condition). This analysis is valid in the case where the linear operator has a discrete (bounded domain) or continuous (the whole space) spectrum.

研究动机与目标

• 将反向误差分析扩展到具有无界线性算子的半线性哈密顿PDE，克服有限维估计中常数发散的问题。
• 证明存在一个修正哈密顿量，其精确流在指数时间尺度内与每个时间步的离散分裂格式完全一致。
• 分析在谱滤波下对非线性薛定谔方程和波动方程的分裂格式（如Lie-Trotter格式），以确保长期稳定性。
• 证明在傅里叶空间中的非共振和一致可积性条件下，数值解能保持与连续流的接近性。
• 推导出能量保持和正则性在依赖于时间步长的指数时间尺度内的定量界。

提出的方法

• 形式反向误差分析将修正哈密顿量 $ H_h $ 构造为 $ h $ 的幂级数，利用哈密顿量的线性部分与非线性部分之间的迭代泊松括号。
• 引入谱滤波函数 $ \alpha_h(x) $ 以正则化线性算子，确保修正哈密顿量的有界性并避免出现无界特征值。
• 该方法适用于离散（环面）和连续（全空间）谱，采用傅里叶表示和伪谱方法。
• 关键估计依赖于零动量条件以及解在傅里叶空间中的 $ L^1 $-有界性，以控制高频模态的增长。
• 通过频率配置上的多指标积分，利用李-泊松括号形式递归求解修正哈密顿量方程。
• 施加非共振条件 $ |\Lambda(\boldsymbol{\tau})| \geq 2\pi - \delta $ 以确保级数收敛及误差项的指数小性。

实验结果

研究问题

• RQ1 能否构造一个修正哈密顿量，使得 $ H_h $ 的精确流在每个时间步与离散分裂映射完全一致？
• RQ2 在何种条件下，分裂格式的数值解能在指数时间尺度内保持与连续流的接近性？
• RQ3 谱滤波如何用于稳定具有无界线性算子的PDE的分裂方法？
• RQ4 傅里叶空间中的一致可积性在确保长期精度中起什么作用？
• RQ5 在该修正哈密顿量框架下，能量和正则性能在多大程度上在指数时间尺度内保持？

主要发现

• 分裂格式的数值流与修正哈密顿量 $ H_h $ 的精确流之间的误差在 $ nh \leq \exp(c h^{-1/2}) $ 时达到 $ \mathcal{O}(\exp(-c h^{-1/2})) $ 的指数小量。
• 对于立方非线性薛定谔方程，当 $ \alpha_h(x) = \sqrt{h} \arctan(x/\sqrt{h}) $ 时，修正哈密顿量 $ H_h^{(1)} $ 在指数时间尺度内将初始能量保持在 $ \mathcal{O}(h) $ 以内。
• 若初始数据属于 $ H^1(\mathbb{R}^3) $，解在 $ H^1 $-范数下在指数时间尺度内保持有界，且通过滤波器控制高频分量。
• 若解保持在 $ B^s_M $ 内且 $ s \geq 3 $，则原始哈密顿量 $ H $ 在指数时间尺度内保持在 $ \mathcal{O}(h) $ 以内。
• 非共振条件 $ |\Lambda(\boldsymbol{\tau})| \leq 2\pi - \delta $ 确保了形式级数的收敛性以及修正系统的指数稳定性。
• 该方法在离散和连续谱（包括全空间 $ \mathbb{R}^d $）下均适用，前提是采用适当的谱滤波。