[论文解读] A numerical scheme for space-time fractional advection-dispersion equation

Shahnam Javadi, M. Jani+1 more|arXiv (Cornell University)|Dec 21, 2015

Fractional Differential Equations Solutions参考文献 18被引用 6

一句话总结

本文提出了一种新颖的谱伽辽金方法，结合乘积积分与伯恩斯坦多项式，以求解时空分数阶对流-弥散方程。通过首先对时间进行离散化，并对空间分数阶导数应用乘积积分，该方法将问题转化为线性代数系统，实现了二阶收敛，并在数值实验中表现出高精度与高效率。

ABSTRACT

In this paper, we develop a numerical resolution of the space-time fractional advection-dispersion equation. After time discretization, we utilize collocation technique and implement a product integration method in order to simplify the evaluation of the terms involving spatial fractional order derivatives. Then utilizing Bernstein polynomials as basis, the problem is transformed into a linear system of algebraic equations. Error analysis and order of convergence for the proposed method are also discussed. Some numerical experiments are presented to demonstrate the effectiveness of the proposed method and to confirm the analytic results.

研究动机与目标

• 开发一种高效数值方法，用于求解具有Caputo导数的时空分数阶对流-弥散方程。
• 解决时间与空间分数阶PDEs耦合求解的挑战，尽管其具有重要的物理意义，但目前仍研究不足。
• 将乘积积分与伯恩斯坦配置法结合，以简化空间分数阶导数的计算。
• 与传统的有限差分或有限元方法相比，降低计算成本与矩阵规模。
• 通过数值实验提供严格的误差分析，并验证收敛速率。

提出的方法

• 首先应用时间离散化，将初值问题转化为适合配置方法求解的边值问题序列。
• 采用乘积积分规则处理空间分数阶导数，实现对各项的显式且稳定的计算。
• 使用伯恩斯坦多项式作为基函数来近似解，从而得到一组线性代数方程组。
• 该方法利用等距配置点的谱配置法，确保高阶精度。
• 通过逼近理论推导出代数系统的误差界，并进行数值求解。
• 该方法避免求解复杂的非线性系统，从而提升了计算效率。

实验结果

研究问题

• RQ1 能否有效结合乘积积分技术与伯恩斯坦多项式，以求解时空分数阶对流-弥散方程？
• RQ2 所提出的方法在空间方向是否如理论分析所预测的那样实现二阶收敛？
• RQ3 与现有的分数阶PDE有限差分或有限元格式相比，该方法在精度与计算成本方面表现如何？
• RQ4 该方法在保持稳定性和收敛性的同时，能否处理可变的分数阶阶数 α, β, 和 γ？
• RQ5 对于不同的分数阶阶数与网格细化，实验收敛速率是多少？

主要发现

• 空间离散化的实验收敛率为约 2.0，与理论预测一致。
• 当 α = 0.2，β ∈ (1,2)，γ ∈ (0,1) 时，L2 误差范数从 h = 1/4 时的 3.315×10⁻³ 降低至 h = 1/16 时的 2.199×10⁻⁴，收敛率保持一致。
• 当 α 变化而 β = 1.5、γ = 0.5 时，L∞ 误差范数从 4.476×10⁻³ 降低至 4.43×10⁻⁴，收敛率接近 1.9。
• 当 γ 变化而 α = 0.5、β = 1.5 时，L2 误差范数随 γ 增加从 3.927×10⁻³ 降至 2.899×10⁻³，收敛率约为 1.9。
• 当时间步长 τ 从 1/10 减小至 1/40 时，CPU 时间从 0.16 秒增加至 6.266 秒，表明计算成本可控。
• 即使在复杂分数阶阶数下，该方法仍保持高精度，所有测试配置下的 L2 误差均低于 5×10⁻³。