用于非线性波动方程的能量守恒型逐次割线差分方法

《Applied Numerical Mathematics》：Energy-conserving successive secant-difference method for nonlinear wave equations

　　

研究人员提出了一种针对非线性波动方程的高阶时间离散方法，该方法能够保持能量守恒。该方法建立在割线差分方法（secant-difference method）基础之上，而该方法本身具有内禀的能量守恒性质。为在保持能量守恒所需结构性质的同时实现更高时间精度，研究人员提出了一种基于逐次策略的新型高阶多阶段方法。尽管逐次割线差分方法（successive secant-difference, SSD）偏离了标准 Runge–Kutta（RK）框架，研究人员仍利用经典 Taylor 展开技术推导了实现高阶精度所需的阶条件。研究进一步给出了直到六阶的完整阶条件，并通过考察其线性相关性确定了最小独立约束组。该推导不仅提供了细致的理论认识，也引入了区别于传统组合方法（composition methods）的新视角。研究人员通过数值算例验证了所提方法的时间精度与无条件能量稳定性，这些算例包括 sine-Gordon 方程和 Boussinesq 方程。研究还给出了与四阶平均向量场方法（averaged vector field methods, AVF）的数值比较。

该文发表于《Applied Numerical Mathematics》，围绕非线性波动方程的保结构时间离散问题展开，核心目标是在高阶时间精度与严格能量守恒之间建立兼容的数值框架。研究背景在于，双曲型偏微分方程（hyperbolic partial differential equations, PDEs）广泛用于描述声学、弹性力学、医学成像、流体力学、量子力学与电磁学中的波传播过程。无论在线性还是非线性情形，能量守恒都是模型物理可信性与长时间数值稳定性的关键要求。现有几何数值积分方法（geometric numerical integrators）虽然在保持辛结构、能量或离散守恒律方面已取得大量进展，但对于非线性波动方程，同时实现高阶精度、较少阶段数以及无条件能量守恒仍然具有明显挑战。特别是，传统组合方法通常依赖对称二阶格式的递归组合，不仅阶段数增长迅速，而且天然局限于偶数阶方法；而若采用适用于线性可加结构的逐次多阶段方法，则在半线性波动方程中又不能保证能量守恒。因此，有必要发展一种既继承割线差分格式能量守恒特征、又能系统提升时间阶数的新方法。

研究人员考察的方程可统一写为^?2u/_?t2=?grad_VE(u)=?G_VδE(u)/δu，其中 E(u) 为定义于特定内积空间 V 上的能量泛函（energy functional），G_V 为与 V 相关的对称正定线性算子。文中给出了 L² 与 H^?1 两类内积空间下的典型情形，并采用周期边界条件或零 Neumann 边界条件。相应 Hamilton 泛函（Hamiltonian functional）写为 H(u)=1/2‖?u/?t‖_V²+E(u)，连续解满足 dH/dt=0。研究以此为基础，选取 sine-Gordon 方程与 Boussinesq 方程作为代表性非线性波动模型，说明所构造能量守恒方法的普适意义与应用价值。

在方法上，研究人员首先将原二阶波动方程通过引入辅助变量 v=u_t 改写为一阶规范系统，从而得到关于位移变量 u 与速度变量 v 的 Hamilton 形式。随后，论文从二阶割线差分方法出发构造逐次割线差分方法（SSD）。该方法的关键思想并不是将非线性项作可加分裂并嵌入经典 Runge–Kutta 方法，而是沿各个中间阶段依次应用割线差分结构，使每一步更新都继承能量平衡恒等式。由于割线差分形式与中点展开之间存在天然联系，研究人员得以围绕阶段中点建立逐次 Taylor 展开，从而推导高阶 SSD 格式的代数阶条件。论文特别指出，虽然 SSD 不属于标准 RK 方法，但其高阶分析仍可借助经典 Taylor 定理完成，并可系统整理出直至六阶的完整条件组，再通过线性相关性分析提炼最小独立约束。这一点构成了本文相较于传统组合方法的重要理论推进。

方法概括而言，作者主要使用了四类关键技术：其一，将非线性波动方程转化为一阶 Hamilton 系统，并在 L² 或 H^?1 内积框架下表述能量守恒律；其二，以割线差分（secant-difference）构造无条件能量守恒的基础离散格式；其三，采用逐次多阶段策略（successive multistage strategy）与阶段中点 Taylor 展开推导高阶 SSD 方法及其阶条件；其四，通过数值实验在 sine-Gordon 方程与 Boussinesq 方程上验证收敛阶、能量稳定性，并与四阶 AVF 方法进行比较。

在“Successive secant-difference methods”部分，研究人员首先定义辅助变量 v=u_t，将原方程改写为 ?u/?t=v，?v/?t=?G_VδE(u)/δu，对应 Hamilton 泛函为 H(u,v)=1/2‖v‖_V²+E(u)。这一改写使问题适于构造成单步多阶段时间离散格式。该部分的核心结论是：基于割线差分构造的 SSD 方法天然保持离散能量守恒，且这一性质并不依赖于时间步长限制，因此具有无条件能量稳定性。文中强调，正是割线差分结构本身简化了能量守恒证明，并使其区别于可加型 RK 离散思路。

在“Order Conditions for the SSD method up to fourth order”部分，研究人员围绕四阶以内高阶格式建立理论基础。通过定义阶段变量 w_i=(u_i,v_i) 及其中点变量 w￣_i=1/2(w_i+w_i?1)，SSD 方法被重写为围绕中点的统一形式，再据此开展逐次 Taylor 展开。该部分得到的主要结论是：尽管 SSD 格式不属于标准 RK 体系，但仍可明确写出保证高阶一致性的代数阶条件；特别是，基于中点表达和阶段递推，可系统获得四阶 SSD 方法的构造条件。这为后续更高阶推广提供了严格基础。

在“Numerical results up to fourth-order accuracy”部分，研究人员对半线性波动方程 u_tt+G_V(Lu+Φ′(u))=0 进行数值检验，其对应 Hamilton 泛函为 H(u,v)=1/2‖v‖_V²+∫_Ω(1/2uLu+Φ(u))dx。文中选取两个代表性算例：一是非均匀介质中的 sine-Gordon 方程，其中 G_V=I，Lu=??·(M(x)?u)，Φ(u)=1?cosu；二是 Boussinesq 方程，其中 G_V=?Δ，Lu=?Δu，Φ(u)=1/2u²?u³。该部分通过计算表明，所提 SSD 方法在这些模型上达到了理论预期的时间收敛阶，同时在计算过程中稳定保持离散能量不变，验证了方法的高精度与保结构特性。

在“Higher-order SSD methods and numerical examples”部分，研究人员进一步研究六阶以内 SSD 高阶方法的构造与性能。该部分的主要工作包括：给出直到六阶精度为止的完整阶条件集合，分析这些条件之间的线性相关关系，并据此识别每一阶所需满足的最小独立约束数。由此得到的结论是，SSD 框架不仅可以构造偶数阶方法，也能够构造奇数阶方法，这一点突破了经典对称组合方法通常仅产生偶数阶格式的限制。随后，论文通过 sine-Gordon 方程上的数值实验展示了相应高阶格式的实际精度，说明理论阶条件与数值表现一致。

在“Comparison with higher-order averaged vector field methods”部分，研究人员将 SSD 方法与高阶平均向量场方法（AVF）进行关系辨析与数值比较。文中指出，单阶段二阶 SSD 方法可以解释为离散梯度方法（discrete gradient, DG）框架中的一种形式，因此更高阶 SSD 格式与基于组合的 DG 方法存在自然联系。另一方面，高阶 AVF 方法可通过不同途径扩展到更高精度。该部分的结论在于：SSD 与 AVF 在低阶保能量思想上存在联系，但在高阶构造机理、代数条件组织方式以及阶段效率方面具有明显差异。论文尤其说明，本文提出的五阶与六阶系数表虽然借鉴了辛 RK 方法中的某些系数结构，但 SSD 方法本身具有内禀能量守恒性质，而传统辛 RK 方法并不自动保证这一点。因此，相较于经典组合型保结构方法，SSD 有望以更少阶段数实现高阶能量守恒积分。

在讨论层面，论文突出强调了两个方面。其一，非线性项的存在使割线差分方法无法像线性可加问题那样化为标准可加 RK 形式，因此必须建立不同于经典 RK 分析的高阶理论；其二，通过围绕中点变量组织 Taylor 展开并梳理线性相关约束，SSD 方法的阶条件呈现出清晰的代数结构，这为系统设计高阶保能量格式提供了新的分析视角。论文还表明，SSD 方法与经典组合方法并非简单替代关系，而是在结构保持、可达阶数及阶段效率上展现出互补优势。

研究结论部分可译为：研究人员提出了一种用于非线性波动方程的高阶时间离散方法，并且该方法能够保持能量守恒。研究以能量守恒的割线差分方法为起点，发展出逐次割线差分方法（SSD），以获得更高的时间精度。由于方程的非线性，割线差分方法不能简化为可加型格式，而逐次策略也不会产生 Runge–Kutta 方法。然而，通过在各阶段中点应用 Taylor 定理，研究人员仍建立了实现高阶精度所需的阶条件体系。总体而言，该研究为非线性波动方程的高阶、保能量时间积分提供了系统的理论构造与数值验证，对长时间稳定模拟与保结构科学计算具有重要意义。