Error analysis of the Strang splitting for the 3D semilinear wave equation with finite-energy data

技术摘要：有限能量数据下三维半线性波动方程的 Strang 分裂误差分析

问题陈述
本文研究了有限能量条件下三维环面 $\mathbb{T}^3$ 上半线性波动方程的时间积分。该方程由下式给出：
$$\partial_{tt} u - \Delta u + \mu u^\alpha = 0, \quad (t, x) \in [0, T] \times \mathbb{T}^3,$$
初始数据为 $(u_0, v_0) \in H^1(\mathbb{T}^3) \times L^2(\mathbb{T}^3)$。研究聚焦于幂次非线性项，其中 $\alpha \in \{3, 4\}$（三次和四次情形）以及符号 $\mu \in \{-1, 1\}$。所解决的主要挑战是：当解仅具有有限能量（正则性为 $H^1 \times L^2$）时，推导数值格式的最优收敛速率。在此情形下，由于 $H^1$ 不构成代数，标准的 Sobolev 嵌入工具在维度 $d \ge 2$ 时不足以使用。

方法论
作者分析了一种结合特定频率滤波器的 Strang 分裂格式（一种二阶三角积分器）的变体。半离散格式通过以下方式计算近似解 $U_n \approx (u(n\tau), \partial_t u(n\tau))$：
$$U_{n+1/2} = e^{\tau A}[U_n + \frac{\tau}{2}G(\Pi_{\tau^{-1}}U_n)], \quad U_{n+1} = U_{n+1/2} + \frac{\tau}{2}G(\Pi_{\tau^{-1}}U_{n+1/2}),$$
其中 $A$ 是线性波动算子，$G$ 代表非线性项，$\Pi_{\tau^{-1}}$ 是一个傅里叶乘子滤波器，将频率限制在 $[-\tau^{-1}, \tau^{-1}]^3$ 范围内。

关键的方法论组成部分包括：

1. 离散时间 Strichartz 估计：本文将离散时间 Strichartz 估计从全空间 $\mathbb{R}^3$ 推广到环面 $\mathbb{T}^3$。这依赖于波动方程的有限传播速度，使得环面上的局部时间估计能够与 $\mathbb{R}^3$ 上的估计相匹配，这与 Schrödinger 方程形成鲜明对比，后者的环面估计受到限制。
2. 分部求和：为了处理由局部截断误差引起的误差项，作者结合滤波器算子使用了分部求和（离散分部积分）。该技术利用了误差项中的相消效应，这一策略受非线性 Schrödinger 方程先前工作的启发，但在此处针对波动方程的特定结构进行了调整。
3. 连续与离散估计：概念上的新颖之处在于，在三次情形（$\alpha=3$）中，使用连续时间 Strichartz 估计来界定局部误差项，从而避免了通常将离散 Strichartz 估计应用于局部误差时导致的阶数损失。对于四次情形（$\alpha=4$），离散 Strichartz 估计对于稳定性是必要的，但仅应用于滤波后的近似解，而非原始数值解。
4. 全离散化：分析扩展到了全离散设置，在空间上采用傅里叶伪谱方法，利用基于 $L^q$ 的插值误差估计来处理负阶 Sobolev 空间。

主要贡献与结果
本文建立了有限能量假设 $(u, \partial_t u) \in H^1 \times L^2$ 下 Strang 分裂的误差界。结果分别针对半离散（仅时间）和全离散（时间和空间）设置进行了展示。

• 三次非线性（$\alpha = 3$）：

  • $L^2 \times H^{-1}$ 范数：几乎二阶收敛，具体为 $O(\tau^2 |\log \tau|)$。
  • $H^1 \times L^2$（能量）范数：几乎一阶收敛，具体为 $O(\tau |\log \tau|)$。
  • 对数因子的出现源于使用了端点 Strichartz 估计（$L^2_T L^\infty$），这需要一种对数修正。

• 四次非线性（$\alpha = 4$）：

  • $L^2 \times H^{-1}$ 范数：收敛阶为 $O(\tau^{3/2})$。
  • $H^1 \times L^2$（能量）范数：收敛阶为 $O(\tau^{1/2})$。
  • 本文指出，与 $\alpha=3$ 相比，$\alpha=4$ 存在“阶数降低”现象。这归因于一阶形式中的非线性项形式上获得了一个导数，而利用 Strichartz 估计仅能对较小的 $\alpha$ 最优地利用这一特性。

• 全离散格式：

  • 对于 $\alpha=3$，误差界为 $C(\tau^2 |\log \tau| + K^{-1})$。
  • 对于 $\alpha=4$，误差界为 $C(\tau^{3/2} + K^{-1})$。
  • 此处 $K$ 是空间离散化参数。鉴于解的 $H^1$ 正则性，$O(K^{-1})$ 的空间收敛速率是最优的。

意义与主张
作者声称，这些是有限能量条件下三维三次和四次波动方程目前可获得的最佳收敛结果。具体而言：

• 结果恢复了三次情形下的最优时间收敛速率（至多相差对数因子），与已知的 1D 结果一致，但将其推广到了 $H^1$ 不构成代数的 3D 情形。
• 分析表明，$\alpha=4$ 时的阶数降低是有限能量下三维波动方程的固有特征，与非线性 Schrödinger 方程的行为不同。
• 本文提供了一个严格的框架，用于在 3D 中使用带滤波器的 Strang 分裂，而无需更高正则性（例如 $H^2$）假设，而此前的高维分析需要这些假设。
• 数值实验证实了理论收敛速率及结果的尖锐性，包括 $\alpha=4$ 时的阶数降低以及临界幂次 $\alpha=5$ 的行为。

该工作与“低正则性积分器”（如修正的 Lie 分裂）的区别在于，它表明经典的二阶格式若结合适当的滤波和 Strichartz 估计，在有限能量下也能实现高收敛速率，尽管本文指出，对于 3D 情形，目前尚无法证明低正则性积分器在此特定设置下优于 Strang 分裂。