[论文解读] A weighted finite difference method for the fractional diffusion equation based on the Riemann-Liouville derivative

Ercı́lia Sousa, Can Li|arXiv (Cornell University)|Sep 11, 2011

Fractional Differential Equations Solutions参考文献 17被引用 13

一句话总结

本文提出了一种针对一维分数阶扩散方程（具有阶数 $1 < \alpha < 2$ 的Riemann-Liouville空间导数）的二阶精度加权有限差分方法。该方法提出了一种新颖的Riemann-Liouville导数的二阶离散化方案，并构建了一个无条件稳定的 $\tau$-格式（包含Crank-Nicolson和后向Euler变体），通过von Neumann分析证明了其稳定性，并通过数值实验验证了在 $\tau = 1/2$ 时空间和时间方向均呈现二阶收敛。该方法能够以高精度有效捕捉异常超扩散的动力学行为。

ABSTRACT

A one dimensional fractional diffusion model with the Riemann-Liouville fractional derivative is studied. First, a second order discretization for this derivative is presented and then an unconditionally stable weighted average finite difference method is derived. The stability of this scheme is established by von Neumann analysis. Some numerical results are shown, which demonstrate the efficiency and convergence of the method. Additionally, some physical properties of this fractional diffusion system are simulated, which further confirm the effectiveness of our method.

研究动机与目标

• 为具有无界域中Riemann-Liouville空间导数的分数阶扩散方程开发一种二阶精度的数值方法。
• 解决非局部Riemann-Liouville导数在实现二阶空间精度方面所面临的挑战，该导数通常以一阶Grünwald-Letnikov格式近似。
• 构建一个加权平均 $\tau$-格式（$\tau \in [1/2,1]$），确保对所有 $\tau \in [1/2,1]$ 实现无条件稳定，并在 $\tau = 1/2$ 时实现二阶精度。
• 通过数值实验验证该方法，包括收敛性研究和异常扩散的物理模拟。

提出的方法

• 通过基于级数的离散化方法（而非有限和）推导出阶数 $\alpha \in (1,2)$ 的Riemann-Liouville分数阶导数的二阶有限差分近似。
• 该方法采用加权平均时间离散化格式，其中 $\tau \in [1/2,1]$，涵盖Crank-Nicolson（$\tau = 1/2$）和后向Euler（$\tau = 1$）方法。
• 通过von Neumann分析严格证明了该格式的稳定性，建立了对所有 $\tau \in [1/2,1]$ 的无条件稳定性。
• 通过仔细逼近Riemann-Liouville导数的非局部积分形式，构建空间离散化以保持二阶精度。
• 该格式被应用于具有可变扩散系数 $d(x)$ 和源项 $p(x,t)$ 的一维分数阶扩散方程，满足初始条件和无穷远处衰减的边界条件。
• 数值实现采用空间和时间上的均匀网格，通过最终时刻的 $l_\infty$ 范数计算收敛率和误差。

实验结果

研究问题

• RQ1 能否在无界域中为阶数 $\alpha \in (1,2)$ 的Riemann-Liouville分数阶导数构造出二阶空间离散化？
• RQ2 当与该二阶空间离散化结合时，加权平均 $\tau$-格式在时间积分中是否能产生无条件稳定且二阶精度的方法？
• RQ3 与基于移位Grünwald-Letnikov公式的现有方法相比，该方法在精度和收敛性方面表现如何？
• RQ4 该数值格式能够多准确地捕捉异常超扩散的动力学行为，特别是当 $\alpha \to 2$ 时？
• RQ5 当 $\tau \neq 1/2$ 时，该格式的收敛率是多少？其是否与理论预期的 $O(\Delta t + \Delta x^2)$ 一致？

主要发现

• 所提出的Riemann-Liouville导数的二阶空间离散化在理论上被证明具有二阶一致性。
• 加权 $\tau$-格式对所有 $\tau \in [1/2,1]$ 均为无条件稳定，该结论经von Neumann稳定性分析确认。
• 当 $\tau = 1/2$ 时，该方法在时间和空间方向均达到二阶精度，数值实验中观测到的收敛率约为2.0。
• 该方法在性能上优于基于移位Grünwald-Letnikov公式的Crank-Nicolson格式，在相同网格分辨率下绝对误差降低了数量级或更多。
• 对具有高斯脉冲初始值问题的数值模拟表明，解表现出超扩散所特有的长尾行为，且随着 $\alpha \to 2$，非Fickian特性逐渐减弱，与物理预期一致。
• 当 $\tau \neq 1/2$ 时，收敛率符合理论预测 $O(\Delta t + \Delta x^2)$，当 $\Delta t = \Delta x^2$ 时，观测到的收敛率趋近于2.0。