[论文解读] Numerical solutions to boundary value problem for anomalous diffusion equation with Riesz-Feller fractional operator

Mariusz Ciesielski, Jacek Leszczyński|arXiv (Cornell University)|Jul 5, 2006

Fractional Differential Equations Solutions参考文献 8被引用 28

一句话总结

本文提出了一种分数阶有限差分法（FFDM），用于求解一维空间中涉及Riesz-Feller分数阶导数的两点边值问题，实现了异常扩散的数值模拟。该方法推广了经典有限差分格式，关键结果表明当 α < 2 时呈现非线性温度分布，且通过 α = 0.35 和 θ = -0.055 可准确拟合实验测得的纳米管热学数据。

ABSTRACT

In this paper, we present a numerical solution to an ordinary differential equation of a fractional order in one-dimensional space. The solution to this equation can describe a steady state of the process of anomalous diffusion. The process arises from interactions within complex and non-homogeneous background. We present a numerical method which is based on the finite differences method. We consider a boundary value problem (Dirichlet conditions) for an equation with the Riesz-Feller fractional derivative. In the final part of this paper, same simulation results are shown. We present an example of non-linear temperature profiles in nanotubes which can be approximated by a solution to the fractional differential equation.

研究动机与目标

• 开发一种用于求解涉及Riesz-Feller分数阶导数的边值问题的数值方法，以描述异常扩散现象。
• 在标准扩散模型失效的复杂非均匀介质（如纳米管）中，建立稳态温度分布的模型。
• 将经典有限差分方法扩展至分数阶导数，以捕捉长程空间相关性及重尾粒子跳跃行为。
• 通过将数值解与纳米管实验数据对比，验证该方法的有效性。
• 研究偏度参数 θ 和阶数 α 对解对称性与形状的影响。

提出的方法

• 使用有限差分法（FDM）对空间中的Riesz-Feller分数阶导数算子进行离散化。
• 通过域内函数值的加权和推导Riesz-Feller导数的数值逼近，其中权重系数依赖于 α 和 θ。
• 构建线性系统 A·T = B，其中矩阵 A 编码分数阶导数的模板，向量 B 包含狄利克雷边界条件。
• 当 α = 2 且 θ = 0 时，该格式退化为二阶导数的经典中心差分法。
• 该方法考虑了非局部行为：每个点的导数依赖于整个域内的值，而不仅邻近点。
• 通过求解线性系统获得解，由于分数阶导数的非局部性，边界条件会影响所有内部节点。

实验结果

研究问题

• RQ1 如何在边值问题中利用有限差分格式准确逼近Riesz-Feller分数阶导数？
• RQ2 分数阶阶数 α 和偏度参数 θ 的变化对稳态温度分布形状有何影响？
• RQ3 所提出的分数阶有限差分法能否准确再现标准扩散模型失效时纳米管中的实验温度分布？
• RQ4 当 α 趋近于 1 且 θ 趋近于 ±1 时，解的行为如何变化？其物理意义是什么？
• RQ5 分数阶导数参数与异常扩散中重尾概率分布出现之间的关系是什么？

主要发现

• 当 α = 2 且 θ = 0 时，解为线性，与经典热传导方程一致。
• 当 α < 2 时，解呈现非线性分布，表明存在异常扩散行为。
• 当 α → 1+ 且 θ → ±1+ 时，解趋近于一阶波动方程的稳态解。
• 偏度参数 θ 引入了解的非对称性，当 θ ∈ (0,1) 时产生非对称分布。
• 采用 α = 0.35 和 θ = -0.055 的模型与 Zhang 和 Li（2005）的实验纳米管温度数据拟合最佳。
• 分数阶导数能够捕捉重尾粒子跳跃行为，从而实现对高斯统计无法描述的罕见但极端事件的建模。