Boundary-Consistent L1–Legendre Spectral Approximation of Nonlinear Caputo Reaction–Diffusion Problems with Integral Constraints

技术摘要：具有非局部积分约束的非线性 Caputo 反映-扩散问题之边界一致性 L1–Legendre 谱近似

问题陈述
本文研究了涉及 Caputo 导数（$0 < \alpha < 1$）且带有非局部积分端点约束的非线性时间分数阶反应-扩散方程的数值解。控制方程为：
$${}^C_0D^\alpha_t u(x, t) = \kappa u_{xx}(x, t) + f(x, t, u(x, t))$$
定义在 $(a, b) \times (0, T]$ 上，并耦合了将端点值与全局空间积分联系起来的边界条件：
$$u(a, t) + \int_a^b \gamma_1(x)u(x, t) \, dx = \phi_1(t)$$
$$u(b, t) + \int_a^b \gamma_2(x)u(x, t) \, dx = \phi_2(t)$$
这些条件与标准的 Dirichlet 或 Robin 条件有本质区别，因为它们并不强制解在边界处为零。作者指出，虽然经典的 L1 时间离散化和谱空间方法对于标准边界条件已非常成熟，但将其应用于非局部积分约束时，需要对空间算子进行仔细处理，以确保与原始强形式微分方程的一致性。

方法论
所提方法结合了边界自适应 Legendre 谱空间与 L1 时间离散化，重点在于构建一个一致的离散算子。

1. 边界提升与空间构建：
   基于前人的工作（Alhejaili, Alsa'di, 和 Salas [11]），作者采用仿射提升函数 $R(x,t)$ 将问题转化为具有齐次积分约束的问题。解分解为 $u = R + y$，其中 $y$ 满足齐次约束。空间逼近空间 $V_N^\gamma$ 使用经过修正的移位 Legendre 多项式构建，以精确满足齐次积分约束。至关重要的是，该空间中的基函数并不一定在端点（$x=a, b$）处为零。

2. 强形式离散化：
   核心的方法论贡献是组装强形式下的空间算子矩阵。作者定义了刚度矩阵 $K$，其分量为 $K_{ij} = (-\psi_j'', \psi_i)_{L^2}$。
   他们强调了针对此类约束的标准谱方法中存在的一个关键不一致性：通常的导数-刚度形式 $S_{ij} = (\psi_j', \psi_i')_{L^2}$ 在由于分部积分产生的非零边界项情况下，与强形式是不等价的：
   $$(-\psi_j'', \psi_i) = (\psi_j', \psi_i') - [\psi_j' \psi_i]_a^b$$
   由于基函数在边界处不为零，项 $[\psi_j' \psi_i]_a^b$ 不为零。论文认为，使用 $S$ 代替 $K$ 实际上是在求解一个不同的离散边界值问题。

3. 时间离散化与非线性求解器：
   Caputo 导数在均匀网格上使用标准的 L1 方案进行近似。生成的全离散系统是非线性的（由于 $f$ 的存在），并通过残差阻尼不动点迭代进行求解。作者在涉及非线性 Lipschitz 常数和逆离散算子范数的显式收缩条件下，对该迭代的收敛性进行了分析。

4. 稳定性与收敛性：
   在**离散增量性假设（discrete accretivity hypothesis）**下推导了稳定性与误差估计。该假设假定强形式矩阵 $K$ 的对称部分相对于质量矩阵是一致正定的。由于对于积分约束，Poincaré 不等式并不像在 Dirichlet 问题中那样能自动保证强制性，因此这一性质在研究中通过数值验证。

主要贡献
本文明确将自身的贡献界定为对现有提升技术的验证与改进，而非基函数的全新推导。具体贡献包括：

• (C1) 推导并使用了边界一致的强形式矩阵 $K$，明确识别了标准导数-刚度近似中所丢失的端点项。
• (C2) 构建了包含非齐次提升项的完全离散 L1 方案，并对残差阻尼非线性求解器进行了分析。
• (C3) 推导了非对称强形式离散的条件稳定性与收敛界限。
• (C4) 在多种分数阶阶数下验证了时间收敛率（$2-\alpha$），避免了单一数值测试。
• (C5) 引入了一个平滑的非多项式基准函数，实现了空间误差与时间误差的分离，证明了导数-刚度矩阵 ($S$) 会导致误差平台（error plateau），而强矩阵 ($K$) 则能收敛至机器精度。
• (C6) 提供可复现的 Python 代码以重新生成所有数值结果。

结果
数值研究通过四个独立的检查验证了该方法：

1. 时间收敛性： 对于具有多项式空间形状的线性基准测试，该方法在不同 $\alpha$ 值（0.25, 0.5, 0.75）下均恢复了预期的 L1 收敛阶 $2-\alpha$。
2. 空间收敛性与算子一致性： 使用平滑非多项式解时，强形式方法表现出谱收敛（误差随多项式阶数 $N$ 指数级下降）。相比之下，若将 $K$ 替换为标准刚度矩阵 $S$，则无论如何增加 $N$，误差都会在 $2.8 \times 10^{-2}$ 附近出现平台期。这证实了 $S$ 求解的是错误的问题。
3. 离散增量性： 通过对 $K$ 的对称部分相对于质量矩阵的广义特征值进行扫描，确认了在测试的核函数范围内，该性质在广泛的阶数范围（$N=4$ 到 $30$）内是正定的。
4. 非线性性能： 非线性测试（使用二次反应项）保持了预测的时间阶，并将积分约束强制执行至机器精度（舍入水平）。

意义与主张
本文谦逊地声称解决了将谱方法应用于非局部积分约束时的一个特定不一致问题。文章强调，虽然边界自适应基函数可以强制执行约束，但它们并不能自动确保空间算子与强微分形式保持一致。其主要意义在于证明了：默默地用导数-刚度算子替换强算子会导致一个虽然收敛但错误的解（即误差平台）。

作者指出，其分析依赖于离散增量性假设，该假设已针对所测试的核函数进行了数值验证，但在一般核函数下的解析表征仍是一个开放性问题。该工作并不声称解决了一般核函数的解析特征问题，而是为满足该条件的此类问题提供了一个稳健且经过验证的数值框架。作者建议的未来扩展方向包括处理初始奇异性的分级网格、快速历史压缩以及多维约束。