Benjamin–Bona–Mahony–Burgers 方程的四阶守恒紧致有限差分方法分析

《Mathematics》：Analysis of a Fourth-Order Conservative Compact Finite Difference Method for Benjamin–Bona–Mahony–Burgers Equation

　　

本研究针对Benjamin–Bona–Mahony–Burgers (BBMB)方程（一种描述多种波动现象动力学的非线性长波方程）提出了一种四阶隐式紧致有限差分格式。通过采用辅助变量的降阶框架，研究人员构建了一种紧致差分格式，该格式产生一个具有窄带状结构的非线性代数系统。由于连续BBMB模型受内在守恒律支配，所提出的数值方法旨在离散意义上保持这一结构性质。该格式的离散守恒性、有界性和唯一可解性得到了严格确立，并推导出了最优的离散最大模误差估计。最后，开展了全面的数值实验，以验证所提出格式的保守性质和理论精度阶。

**论文解读：Benjamin–Bona–Mahony–Burgers 方程的四阶守恒紧致有限差分方法**

**研究背景与问题**

Benjamin–Bona–Mahony–Burgers (BBMB) 方程是描述非线性色散介质中小振幅长波单向传播的基本模型，广泛应用于浅水波动力学、分层流体中的内波传播、等离子体物理以及地球物理流体力学等领域。该方程结合了Benjamin–Bona–Mahony (BBM) 方程与Burgers型耗散项，在能量耗散显著的场景中提供更全面的物理框架。然而，现有数值方法在同时实现高精度、结构守恒和计算效率方面存在不足。尽管已有多种有限差分、有限元及样条方法被提出，但大多仅达到二阶空间精度。目前唯一达到四阶空间精度的紧致格式是Zhang和Liu提出的三层线性化格式，但该格式存在三个关键缺陷：第一，它需要同时存储两个完整时间层（uⁿ和uⁿ⁺¹）的数组，在细网格上显著增加内存占用；第二，它需要专用的第一步求解器生成u¹，增加了实现复杂度且可能导致起始近似阶数不匹配；第三，由于线性化，它不能直接继承连续问题的守恒律，仅保持一个两层平均不变量，而非严格的单层不变量。因此，开发一种同时满足四阶空间精度、严格单层守恒、仅需一个历史时间层的两层紧致格式成为亟待解决的问题。

**研究内容与结论**

研究人员针对带周期边界条件的BBMB方程，通过引入辅助变量v = u_x的降阶框架，构造了一个两层非线性隐式紧致有限差分格式。该格式在时间上达到二阶精度，在空间上达到四阶精度，且严格遵守离散守恒律。研究人员运用Browder不动点定理证明了每个时间层离散解的存在性，并在步长条件τ < 2/α下证明了唯一性。通过离散内积和Gronwall不等式，推导出最优的离散最大模误差估计O(τ² + h⁴)。数值实验（包括高斯型、双曲正割型、严格周期正弦型初始条件，以及非齐次源项和广义非线性势扩展）全面验证了理论精度和守恒性质。该工作发表在《Mathematics》。

**主要关键技术方法**

1. **降阶框架**：引入辅助变量v = u_x，将原三阶BBMB方程（2）转化为一阶耦合系统（10）–（11），从而允许使用紧致算子进行空间离散。
2. **紧致差分格式**：在时空网格点上应用泰勒展开，对线性项采用四阶Padé型紧致算子，对非线性项利用引理4中的线性化技巧，得到局部截断误差为O(τ² + h⁴)的离散格式（30）–（31）。
3. **非线性求解器**：每个时间步求解一个由块三对角雅可比矩阵构成的非线性系统，采用牛顿迭代法，初始猜测为上一时间层解，迭代终止条件为最大修正小于10^-12。
4. **误差分析**：利用离散Sobolev不等式、Young不等式、引理6和Gronwall不等式，在假设精确解足够光滑（u ∈ C^4,3）且步长满足τ < 1/(2C₁)的条件下，推导出最大模误差估计。

**研究结果**

**3.1 离散守恒性质**
通过取离散内积并应用引理3和引理5，研究人员证明了定理1：离散解满足守恒量Iⁿ = ||uⁿ||² + hΣ_i=1^M (v_i-1/2ⁿ)²在全时间层上严格守恒，即Iⁿ = I⁰。推论1进一步给出解的有界性：||uⁿ|| ≤ C，||u_xⁿ|| ≤ C，以及||uⁿ||_∞ ≤ C。

**3.2 可解性**
研究人员利用Browder不动点定理（定理2）证明了定理3：对于每个时间层，存在解uⁿ⁺¹和v^n+1/2。定理4在条件τ < 2/α下证明了唯一性：通过假设两组解，推导误差方程并取内积，得到误差能量不等式，进而得出误差为零。

**4 误差分析**
定理5建立了误差估计：设eⁿ = Uⁿ - uⁿ，f^n+1/2 = V^n+1/2 - v^n+1/2，在τ < 1/(2C₁)条件下，存在常数C使得||eⁿ|| ≤ C(τ² + h⁴)，||f^n+1/2|| ≤ C(τ² + h⁴)。证明过程通过误差方程（68）–（69），取内积，利用引理6、Young不等式和Gronwall不等式逐步推导。

**5.2 数值实验**
研究人员设计了5个数值例子，采用离散L²范数和L^∞范数计算误差，通过网格细化法获得实验收敛阶。对于齐次方程（例1–3），分别采用高斯、双曲正割和正弦初始条件，在固定时间步长下改变空间网格数M以测试空间精度，在固定空间网格下改变时间步数N以测试时间精度。结果明确显示：空间收敛阶接近4，时间收敛阶接近2，且守恒量Iⁿ的绝对偏差保持在机器精度水平。对于非齐次方程（例4），添加已知解析解的强制项，直接计算全局误差，验证了四阶空间和二阶时间精度。对于广义势项扩展（例5），修改牛顿迭代的右端项和雅可比矩阵，同样获得一致的精度和守恒性。

**结论**

本研究设计、分析并验证了一种针对Benjamin–Bona–Mahony–Burgers (BBMB)方程的新型两层四阶结构保持紧致有限差分格式。通过策略性地采用辅助变量v = u_x的降阶框架，研究人员成功将高阶控制方程转化为耦合低阶系统，该格式保持紧凑窄带宽度，显著减少了边界闭合复杂性并降低了整体计算模板尺寸。该格式在时间上严格保持离散不变量Iⁿ，确保数值解忠实模仿连续模型的几何和物理性质。利用Browder不动点定理严格证明了每个时间层离散解的存在性，并正式证明了其条件唯一性。研究人员还推导了最优阶误差估计，表明该格式在离散L²网格范数下达到O(τ² + h⁴)收敛。广泛的数值实验（涵盖高斯、双曲正割和严格周期正弦初始条件，以及非齐次源框架和广义非线性势扩展）完全证实了理论四阶空间精度和二阶时间精度。关键的是，实证结果确认绝对不变量偏差稳定保持在机器精度附近，凸显了所提出框架模拟复杂非线性波现象的鲁棒可靠性、几何保真度和长期稳定性。