高精度多不变量保持的隐式-显式Runge-Kutta方法：基于拉格朗日乘子框架在守恒非线性演化方程中的应用

《Communications in Nonlinear Science and Numerical Simulation》：High-order structure-preserving implicit-explicit Runge–Kutta methods with Lagrange multiplier approach for conservative nonlinear evolution equations

本研究综述了一类新型结构保持数值方法，其结合了隐式-显式龙格-库塔（IERK）方法与拉格朗日乘子（LM）框架，旨在为非线性演化方程构建可保持多个高阶不变量（如能量、质量）的高阶、线性隐式数值格式。该方法克服了经典投影法无法保持线性不变量等局限，在实现任意高阶时间精度的同时，保证了计算效率，在非线性光学、量子力学及水波理论等领域具有重要应用潜力。文中通过理论分析与大量数值算例验证了方法的有效性和长期模拟的稳定性。

　　

Highlight（研究亮点）

• •
  一旦拉格朗日乘子被确定，得到的格式保持线性隐式特性。
• •
  所提出的框架能够保持多个原始的高阶离散非线性不变量。
• •
  所提出的框架可以实现任意高阶的时间精度。

高精度多不变量保持格式

在本节中，我们将基于隐式-显式龙格-库塔（IERK）方法，为等价系统（2.11）开发一类能够保持多个不变量的高阶隐式-显式格式。

我们首先回顾s级q阶IERK（IERKq）方法。令 {a_ij, b_i, a?_ij, b?_i, c_i} (i, j ∈ Z₁^s) 为实数，并定义 A = (a_ij)_s×s, b = (b₁, …, b_s)^T, A? = (a?_ij)_s×s, b? = (b?₁, …, b?_s)^T, c = (c₁, …, c_s)^T。则IERK方法可以使用以下Butcher表来表示：

c	A
bT

=

c₁| a₁₁

c₂| a₂₁a₂₂

? | ? ? ?

c_s| a_s1a_s2? a_ss

| b₁b₂? b_s

,

c	A?
b?T

=

c₁| 0

c₂| a?₂₁0

? | ? ? ?

c_s| a?_s1a?_s2? a?_ss

| b?₁b?₂? b?_s

主要结果

LM-IERK格式（3.2）的适定性取决于乘子λ^k的存在唯一性，这源于IERK方法的线性隐式结构。由于当 ?^k+1= u^k时 λ^k= 0，后续分析仅关注非平凡情况 ?^k+1≠ u^k。

定理 2 (λ^k+1的局部存在唯一性) 对于LM-IERK格式（3.2），假设IERK方法 {A, b, A?, b?} 是q阶的（q ≥ 1）且 I_i[u^k] = I_i[φ], i ∈ Z₁^p, k ∈ Z₁ⁿ。如果 (1/τ)B(0,0) 是非奇异的，则存在一个常数 τ₀> 0，使得对于 0 < τ < τ₀，代数系统... [定理内容因文档截断而未完整呈现]。

数值结果

本节中，我们将通过几个数值算例来展示所提出方法的有效性。

结论与展望

本文为保守非线性演化方程的数值求解开发了一类新颖的IERK方法。所提出的框架融合了LM方法，以确保多个离散不变量的保持。详细的适定性分析确立了相关拉格朗日乘子的存在性和唯一性，并证实了所得格式能达到高阶时间精度。大量的数值实验进一步证明了该方法在... [文档截断，后续内容未完整提供]。