一种用于非线性波动方程的高阶多辛有限元方法

《Physica C: Superconductivity and its Applications》：A high-order multi-symplectic finite element method for nonlinear wave equations

　　彭城程 摘要：为了实现非线性波动方程长期可靠模拟，保持几何结构与守恒行为至关重要。本文针对一类一维多辛非线性波动方程，结合基于有限元外微分法的相容空间离散化、对非线性势项的几何一致处理方式以及四阶高斯-勒让德辛时间积分器，提出了一种高阶结构保持方法。研究表明，该完全离散方案满足离散多辛守恒律，能够实现离散能量与动量的有限时间近似守恒，并在适当光滑性假设下具有O(hk+1+Δt4)的误差界。此外，在固定网格上还可得到有界时间区间内的非线性稳定性估计。数值实验验证了理论结果：构造解测试表明收敛阶符合预期，长时Klein-Gordon方程模拟展现了良好的能量行为，且数值诊断结果与混沌动力学特征一致；另外还通过二维正弦-戈登呼吸子问题作为高维场景下的对比示例，说明在所测试参数下，该方法相比二阶参考方案具有更小的能量误差和更低的数值辐射，但每步计算成本更高。这些结果表明，基于有限元外微分法的高阶结构保持框架，是那些对几何兼容性与守恒行为有较高要求的非线性波动方程模拟的可行方法。

引言：非线性波动方程是一类用于描述从流体力学、等离子体物理到非线性光学、凝聚态物理以及基本粒子理论等广泛物理现象的核心数学模型[1][2]。诸如Klein-Gordon方程和正弦-戈登方程这样的典型模型，体现了波色散与非线性作用之间的本质关联，往往会导致孤子、呼吸子等强健的、空间局域化的相干结构的产生与相互作用[2][3]。准确预测这类系统的长期演化过程，尤其是在涉及复杂相互作用、模式形成或向混沌状态转变的情况下，是一项极具挑战性的任务，这需要高度可靠的数值模拟工具[4]。传统的数值方法，包括标准有限差分或有限元离散化结合通用时间步进方案，往往难以在长时间模拟过程中保持精度[5][6]。这些方法可能会引入人为的数值耗散或色散效应，不仅可能违反能量、动量、质量等基本物理守恒定律，还可能破坏控制方程中固有的几何结构。对于哈密顿系统而言，这种结构通常为辛结构或多辛结构[7]。结构退化带来的后果可能十分严重，会导致模拟结果出现定性错误，如解的人为阻尼、波速不准确、虚假的数值不稳定性，以及无法捕捉相干结构的长期稳定性或混沌现象的起始等[7][8]。

为了解决这些问题，几何数值积分领域逐渐兴起，其核心目标是开发专门用于保持连续动力系统几何特性的数值算法[7][8][9]。对于哈密顿常微分方程，辛积分器因其能够保持相空间体积并展现出优异的长期能量行为（通常表现为围绕真实值的有限波动，而非长期漂移）而备受推崇[7][10]。变分积分器则源自哈密顿原理的离散化形式，为构建既能保持系统对称性对应的动量映射、又能实现结构保持的辛方法提供了系统化的途径[11][12]。在处理哈密顿偏微分方程时，则出现了多辛性概念，它体现为涉及时空结构双方的守恒律[13][14]。多辛积分器的目标便是保持这种局部时空守恒律的离散形式[15]，人们认为这有助于提升全局守恒性质及整体模拟精度[14][16]。目前已有多种多辛方案被提出，包括类似普雷斯曼盒法的有限差分方法及其能量特性研究[14][17]，以及多辛龙格-库塔方法[18][19][20]。这类方案的长期行为一直是研究热点，尤其是针对非线性Klein-Gordon方程等特定非线性波动方程[21]。虽然早期研究主要集中于有限差分方案，但在广义多辛积分器方面也取得了显著进展：胡、邓等人[22]建立了适用于哈密顿非线性波动偏微分方程的广义多辛积分器框架，洪和李则[23]将多辛龙格-库塔方法扩展到了非线性狄拉克方程领域。此外，孙[24]对分块龙格-库塔方法中的二次不变量与多辛性进行了严格分析，为这类方案奠定了坚实的理论基础。

就本研究所使用的特定方程而言，王等人做出了重要贡献：他们为正弦-戈登方程开发了高阶辛格式[25]，并对Klein-Gordon-薛定谔方程的傅里叶伪谱守恒格式进行了分析[26]，从而证明了结构保持方法在长期模拟中的优势。至于在我们的工作中处于核心地位的有限元离散化方法，甄等人[27]研究了基于简单有限元方法的辛格式与多辛格式。最近，塞莱多尼和杰克曼[28]则对多辛偏微分方程的有限元离散化形式的离散守恒律进行了详细分析。这些研究成果，再加上金等人[29]在结构保持型神经网络方面的最新进展，都凸显出几何积分策略的持续重要性与生命力。然而，要构建能够在时空两方面都实现高精度、可灵活嵌入有限元法等现有框架中、且能有效处理潜在非线性的多辛积分器，仍然是一项挑战。许多现有的多辛方法都以有限差分为基础，且精度通常较低。此外，标准的有限元投影操作可能会破坏几何结构。有限元外微分法[30]则为在尊重偏微分方程的几何与拓扑结构的同时对其进行离散化提供了一种系统化方法，因此很适合用于几何数值积分中的空间离散化操作[9]。要确保非线性项，尤其是势能项的离散化形式与守恒量（如动量）以及底层几何结构相兼容，需要仔细设计，这往往涉及到保持对称性、采用特定的投影方法[7][11]，或运用离散梯度等概念来确保能量守恒或其他结构特性[31]。此外，像高斯-勒让德格式这类高阶隐式方法的实现也会带来一系列挑战，尤其是在处理受约束或结构复杂的系统时[32]。为了利用其在高阶场景下的结构保持特性，人们还研究了分裂法等先进的 time integration策略[33]。

虽然有限元外微分法[34]最初是为了解决混合变分问题中的稳定性问题而开发的——例如通过确保inf-sup条件来避免锁模现象或虚假模态的出现[30]——但它对同调代数的依赖为几何数值积分带来了独特优势。有限元外微分法的核心创新在于构建了构成连续de Rham复形微分子复形的离散空间。通过为这些空间配备与外导数可交换的有界上链投影（即可交换图性质），有限元外微分法能够确保微分算子的代数性质，尤其是d2=0和斯托克斯定理，在离散层面也严格成立。在多辛积分的背景下，这种结构保持并非仅仅是稳定性要求，更是几何上的必要条件。它能够保证离散辛形式保持封闭性（dhωh=0），并且确保从变分原理导出的离散守恒律是离散化过程本身固有的，而非渐进近似结果。

本文介绍并分析了一类具有多辛哈密顿表述的一维非线性波动方程的高阶结构保持数值格式。我们的方法将以下几部分有机结合：首先，利用基于有限元外微分法原理的空间离散化方法，保持与空间导数相关的几何结构；其次，以符合系统哈密顿本质及相应守恒律的方式处理非线性势项；最后，采用四阶辛高斯-勒让德龙格-库塔方法进行时间积分，从而确保时间维度上的高精度与结构保持[7]。本文的主要贡献如下：首先，我们为这类多辛非线性波动方程构建了一个统一的高阶完全离散框架，该框架结合了基于有限元外微分法原理的相容空间离散化、对非线性项的几何一致处理方式，以及四阶辛高斯-勒让德时间积分器。其次，我们在当前研究框架下对该方案进行了理论分析，具体包括推导出离散多辛守恒律，建立离散能量与动量的有限时间近似守恒结果，得到固定网格上有限时间区间内的非线性稳定性估计，并在适当光滑性假设下证明了O(hk+1+Δt4)的收敛速率。第三，我们通过三项数值测试来展示该方法的表现：构造解实验验证了预期的收敛阶；Klein-Gordon方程的长时模拟既保持了良好的能量行为，其数值诊断结果也与混沌动力学特征相符；此外，还通过二维正弦-戈登呼吸子问题的对比实验作为高维场景下的示例，说明在所测试参数下，该方法相比二阶参考方案具有更小的能量误差和更低的数值辐射，但每步的计算成本更高。

本文结构如下：第2节详细介绍了模型问题、其多辛哈密顿结构以及相关的守恒定律；第3节阐述了数值方法，包括空间离散化（有限元外微分法与非线性项处理）、时间积分（高斯-勒让德方法）以及实现细节；第4节给出了支撑该方法各项性能（守恒性、稳定性、收敛性）的理论分析；第5节报告了数值实验的结果，用以验证该方案的性能；最后，第6节总结了研究结论，并提出了未来研究的方向。

模型问题与哈密顿结构：本节介绍了模型问题，即一类非线性波动方程，并阐明了其背后的哈密顿结构。理解这一结构对于开发能够保持系统定性特征（如守恒定律）的数值方法至关重要。

数值方法：本节详细描述了用于离散化非线性波动方程（2.1式）同时保持其固有几何结构的数值方法。该方法结合了基于有限元外微分法的结构保持型空间离散化、动量映射处理方式以及辛时间积分方案。尽管目前的空间离散化在理念上与离散外微分法存在几何关联，但实际上所采用的结构是基于……的。

离散多辛守恒律：定理1 离散多辛守恒律：该完全离散方案在每个时空单元[tn,tn+1]×[xm,xm+1]上满足以下离散局部守恒律：
1Δt(ωhn+1,m?ωhn,m)+1Δx(κˉhn,m+1?κˉhn,m)=0，
其中ωhn,m表示在空间元素上积分后的离散辛形式（时间部分），κˉhn,m则表示时间平均后的空间通量形式。具体表达式为：
ωhn,m=∫xmxm+1(δuhn∧δvhn)dx，
κˉhn,j=Δt∑i=1sbi(δUi(xj)∧δWi(xj))，
其中s=2为积分阶数，bi为高斯-勒让德权重，(Ui, Wi)为……

数值实验：本节通过三项独立的数值实验来验证所提出的结构保持算法：首先使用构造解方法验证理论上的收敛速率，其次研究Klein-Gordon方程中的混沌动力学现象，最后将该方法与其他用于模拟正弦-戈登方程中二维呼吸子动力学的结构保持方法进行对比研究。

结论：在本研究中，我们为具有多辛哈密顿结构的非线性波动方程开发并分析了一种高阶结构保持方法。该方法结合了基于有限元外微分法原理的相容空间离散化、对非线性势项的几何一致处理方式，以及四阶辛高斯-勒让德时间积分器。从理论层面来看，我们在当前分析范围内确立了该完全离散方案的若干性质。在……

利益冲突声明：作者声明自己不存在任何可能影响本文研究结果的已知财务利益或个人关系。