[论文解读] Fractional Diffusion Processes: Probability Distributions and Continuous Time Random Walk

Rudolf Gorenflo, Francesco Mainardi|ArXiv.org|Sep 25, 2007

Fractional Differential Equations Solutions参考文献 47被引用 7

一句话总结

本文提出了一种基于Riesz-Feller和Caputo导数的时空分数阶扩散方程，用于建模异常扩散，表明其基本解对应于自相似的概率密度函数。其主要贡献在于证明了该方程作为连续时间随机行走（CTRW）主方程在等待时间与跳跃宽度被压缩时的标度极限而出现，从而将分数阶动力学与具有重尾分布的随机过程联系起来。

ABSTRACT

A physical-mathematical approach to anomalous diffusion may be based on fractional diffusion equations and related random walk models. The fundamental solutions of these equations can be interpreted as probability densities evolving in time of peculiar self-similar stochastic processes: an integral representation of these solutions is here presented. A more general approach to anomalous diffusion is known to be provided by the master equation for a continuous time random walk (CTRW). We show how this equation reduces to our fractional diffusion equation by a properly scaled passage to the limit of compressed waiting times and jump widths. Finally, we describe a method of simulation and display (via graphics) results of a few numerical case studies.

研究动机与目标

• 建立一种使用空间和时间分数阶导数的广义扩散模型，以描述异常输运。
• 将时空分数阶扩散方程的基本解解释为随时间演化的概率密度函数。
• 通过标度极限，证明分数阶扩散方程与连续时间随机行走（CTRW）模型之间的联系。
• 提供一种用于模拟分数阶扩散过程的方法，并通过数值案例研究进行验证。
• 通过引入幂律尾部的等待时间与跳跃分布，将异常扩散的理解扩展至布朗运动之外。

提出的方法

• 利用阶数β ∈ (0,1]的Caputo时间导数和阶数α ∈ (0,2]、偏度为θ的Riesz-Feller空间导数，构建时空分数阶扩散方程。
• 通过Mellin-Barnes积分表示法推导基本解（格林函数），其描述了分数阶扩散过程中随时间演化的概率密度。
• 在等待时间与跳跃宽度被压缩的标度极限下，表明CTRW主方程可简化为分数阶扩散方程，其关系由λτ^β = μh^α决定。
• 使用Mittag-Leffler函数，以及一种更简单的可逆替代形式Ψ*(t) = 1/(1 + Γ(1−β)t^β)，来建模具有幂律衰减的等待时间分布。
• 采用完全单调性与拉普拉斯变换技术，确保等待时间与跳跃宽度分布的物理一致性。
• 基于CTRW结构开发模拟框架，对选定参数情形的数值结果以图形方式展示。

实验结果

研究问题

• RQ1 如何利用空间和时间的分数阶导数来建模超越标准布朗运动的异常扩散？
• RQ2 时空分数阶扩散方程的基本解具有何种概率解释？
• RQ3 在适当的标度下，连续时间随机行走（CTRW）模型如何收敛至分数阶扩散方程？
• RQ4 导致分数阶扩散的等待时间与跳跃宽度分布的渐近特性是什么？
• RQ5 能否开发一种数值上可行且物理解释一致的分数阶扩散过程模拟方法？

主要发现

• 时空分数阶扩散方程的基本解是一种自相似的概率密度函数，其推广了布朗运动的高斯分布。
• 当β = 1时，等待时间分布呈指数衰减（平均等待时间有限），恢复标准扩散；当0 < β < 1时，其衰减形式为t^{-(β+1)}，表明平均等待时间无限，具有重尾行为。
• Mittag-Leffler函数Ψ(t) = E_β(−t^β)与替代形式Ψ*(t) = 1/(1 + Γ(1−β)t^β)在t较大时均表现出相同的渐近幂律衰减t^{-(β+1)}，从而确保长尾等待时间。
• 时间（τ）与空间（h）压缩因子之间的标度关系λτ^β = μh^α确保了CTRW向分数阶扩散方程的收敛。
• 数值模拟证实了异常扩散行为的出现，β=1时表现为布朗运动特征，而β<1时则表现出超扩散或亚扩散动力学。
• 基于CTRW的模拟方法能够高效研究分数阶扩散过程，尤其适用于长时间或自相似状态。