A numerical approach for solving the time-Fractional Mobile-Immobile Transport Equation

技术摘要：求解时间分数阶移动-不动传输方程的一种数值方法

问题陈述
本文研究了一维时间分数阶移动-不动（Mobile-Immobile, MIM）传输方程的数值解，该模型用于描述多孔介质中的异常扩散现象。与经典的 Fick 扩散模型不同，MIM 模型通过阶数为 $\alpha \in (0, 1)$ 的 Caputo 分数阶导数来描述由移动区域与不动区域之间的质量传递导致的溶质传输延迟。所考虑的具体初值-边界值问题（IBVP）为：
$$\lambda_1 \frac{\partial u}{\partial t} + \lambda_2 {}^C D^\alpha_{0,t} u + Lu = f(x,t),$$
其中 $L$ 是空间椭圆算子，系统满足零 Dirichlet 边界条件。主要的挑战在于分数阶导数的非局部性质，以及需要空间和时间上的高阶精度，以高效捕捉记忆效应和异常传输行为。

方法论
作者开发了两种将空间离散化与时间离散化结合在统一框架内的全离散数值格式：

1. 空间离散化： 采用了非对称内罚间断有限元法（Non-Symmetric Interior Penalty Discontinuous Galerkin, NIPG）。该方法利用破碎 Sobolev 空间，并允许单元边界处存在不连续性。空间算子通过包含惩罚项的双线性形式 $B_h(\cdot, \cdot)$ 进行近似，以确保稳定性与强制性。
2. 时间离散化： 提出了两种不同的策略来近似时间导数（包括一阶导数和 Caputo 分数阶导数）：
   • 方案 1 (CN-L1-NIPG)： 将 NIPG 方法与应用于标准 L1 公式（针对 Caputo 导数）的 Crank-Nicolson (CN) 框架相结合。该方案的目标是实现 $(2-\alpha)$ 阶的时间精度。
   • 方案 2 (L2-1$\sigma$-NIPG)： 将 NIPG 方法与 L2-1$\sigma$ 公式（由 Alikhanov 引入，用于 Caputo 导数）相结合，同样是在 CN 框架内。该方案旨在实现二阶时间精度。

分析过程依赖于使用能量范数的严格稳定性证明，以及通过椭圆投影和离散 Grönwall 型不等式推导出的误差估计。

主要贡献

• 新颖的格式： 本文引入了两种针对时间分数阶 MIM 模型设计的不同全离散格式（CN-L1-NIPG 和 L2-1$\sigma$-NIPG），与仅关注简单分数阶扩散方程的现有文献相比，这种组合受到的关注较少。
• 理论分析：
  • 稳定性： 为两种格式建立了无条件稳定性。分析表明，数值解受初始数据和源项限制，且与网格尺寸无关。
  • 收敛性： 推导出了最优误差估计。CN-L1-NIPG 方案实现了 $O(\tau^{2-\alpha} + h^{k+1})$ 的误差界，而 L2-1$\sigma$-NIPG 方案实现了 $O(\tau^2 + h^{k+1})$ 的误差界，其中 $\tau$ 是时间步长，$h$ 是空间网格尺寸，$k$ 是多项式阶数。这些收敛阶数取决于精确解的正则性。
• 综合误差界： 误差分析明确考虑了解的正则性，为观察到的收敛速率提供了理论基础。

结果
通过一个具有已知精确解的测试问题进行了数值实验，以验证理论预测。

• CN-L1-NIPG： 数值结果证实了在不同 $\alpha \in (0, 1)$ 取值下，其在时间方向上的理论收敛阶为 $(2-\alpha)$。
• L2-1$\sigma$-NIPG： 实验表明，该方案实现了二阶时间精度，在时间导数的收敛速率方面优于基于 L1 的方案。
• 鲁棒性： 两种方案在不同分数阶和网格加密下均表现出高效性和鲁棒性，其中 L2-1$\sigma$ 方法展现出更优越的时间精度。

意义与主张
本文声称为模拟多孔介质中复杂的分数阶传输过程提供了一个可靠且精确的计算框架。通过将 NIPG 方法的灵活性与高阶时间离散化相结合，所提出的方案为模拟非 Fick 传输现象提供了一种统一的方法。作者断言，其工作填补了关于 MIM 模型高阶数值方法的文献空白，为水文学、污染物传输和沉积物弥散等应用领域提供了一个既有理论依据又具实际效率的工具。研究结论指出，这些方法为未来向高维及变阶分数阶模型的扩展铺平了道路。