[论文解读] Unconditional Stability Of A Two-Step Fourth-Order Modified Explicit Euler/Crank-Nicolson Approach For Solving Time-Variable Fractional Mobile-Immobile Advection-Dispersion Equation

Eric Ngondiep|arXiv (Cornell University)|May 6, 2022

Fractional Differential Equations Solutions被引用 5

一句话总结

本文提出了一种用于求解具有变阶Caputo导数的时间可变分数阶移动-非移动对流-弥散方程的两步四阶改进显式欧拉/克兰克-尼科尔森格式。该方法实现了无条件稳定性和 $O(k + h^4)$ 阶的收敛性，相较于文献中现有的单步格式表现出更优的效率和精度。

ABSTRACT

This paper considers a two-step fourth-order modified explicit Euler/Crank-Nicolson numerical method for solving the time-variable fractional mobile-immobile advection-dispersion model subjects to suitable initial and boundary conditions. Both stability and error estimates of the new approach are deeply analyzed in the $L^{\infty}(0,T;L^{2})$-norm. The theoretical studies show that the proposed technique is unconditionally stable with convergence of order $O(k+h^{4})$, where $h$ and $k$ are space step and time step, respectively. This result indicate that the two-step fourth-order formulation is more efficient than a broad range of numerical schemes widely studied in the literature for the considered problem. Numerical experiments are performed to verify the unconditional stability and convergence rate of the developed algorithm.

研究动机与目标

• 开发一种针对具有变阶Caputo导数的时间可变分数阶移动-非移动对流-弥散方程的高阶、无条件稳定数值方法。
• 克服现有单步格式常存在的时间精度较低或条件稳定性等局限性。
• 在保持计算效率的同时，实现四阶空间精度和一阶时间收敛性。
• 针对所提出的格式，在 $L^{ty}(0,T;L^{2})$-范数下提供严格的稳定性与误差分析。
• 通过具有精确解的数值实验和收敛率分析，验证理论结果。

提出的方法

• 采用两步四阶格式，结合改进的显式欧拉与克兰克-尼科尔森格式，以平衡精度与稳定性。
• 时间导数采用二阶精度的两步格式近似，空间导数则使用四阶紧凑差分进行离散化。
• 基于包含奇异核的定义，利用 quadrature 规则对变阶Caputo分数阶导数进行近似。
• 通过能量方法在 $L^{ty}(0,T;L^{2})$-范数下进行稳定性分析，证明该方法的无条件稳定性与时间步长和空间步长无关。
• 该格式被表述为线性方程组，可通过高效的迭代求解器求解，未来工作可考虑预条件技术。
• 对分数阶导数项采用半隐式处理，以确保稳定性并允许采用更大的时间步长。

实验结果

研究问题

• RQ1 两步四阶数值格式能否对时间可变分数阶移动-非移动对流-弥散方程实现无条件稳定性？
• RQ2 所提方法在时间和空间上的收敛阶是多少？是否优于现有的单步格式？
• RQ3 当分数阶在时间和空间上变化时，该方法在稳定性和精度方面表现如何？
• RQ4 在 $L^{ty}(0,T;L^{2})$-范数下的理论误差界是否可通过具有精确解的数值实验得到验证？
• RQ5 与该类问题中广泛使用的文献格式相比，所提格式在计算效率上是否更具优势？

主要发现

• 所提格式在 $L^{ty}(0,T;L^{2})$-范数下具有无条件稳定性，意味着其稳定性对所有时间步长和空间步长均成立。
• 该方法在时间方向实现一阶收敛（$O(k)$），在空间方向实现四阶收敛（$O(h^4)$），整体收敛阶为 $O(k + h^4)$。
• 数值实验验证了理论收敛速率，多个测试案例中观测到的空间收敛阶约为 4.0，时间收敛阶为 3.8–4.1。
• 在示例1中，当 $h = 1/16$ 时，最大误差为 $1.1488 \times 10^{-6}$，收敛率为 4.1029，证实了四阶空间精度。
• 在示例2中，当 $\alpha = 0.25$ 时，随着 $h$ 减小，误差从 $1.0918 \times 10^{-2}$ 降低至 $2.5995 \times 10^{-6}$，收敛率为 4.1328，进一步证实了高阶精度。
• 该方法在效率和速度方面优于广泛存在的现有格式，理论分析与数值结果均表明了这一点。