[论文解读] The fundamental solution of the space-time fractional diffusion equation

Francesco Mainardi, Yuri Luchko+1 more|ArXiv.org|Feb 18, 2007

Fractional Differential Equations Solutions参考文献 43被引用 728

一句话总结

本文推导了具有Riesz-Feller空间分数阶导数和Caputo时间分数阶导数的时空分数阶扩散方程的基本解（格林函数）。通过傅里叶-拉普拉斯变换和梅林-巴恩斯积分表示，建立了其标度性质，并为约化格林函数提供了显式收敛级数和渐近展开式，从而在包括非对称和非马尔可夫情形在内的广泛参数范围内实现了概率密度解释。

ABSTRACT

We deal with the Cauchy problem for the space-time fractional diffusion-wave equation, which is obtained from the standard diffusion equation by replacing the second-order space derivative with a Riesz-Feller derivative of order alpha in (0,2] and skewness theta, and the first-order time derivative with a Caputo derivative of order beta in (0,2]. The fundamental solution is investigated with respect to its scaling and similarity properties, starting from its Fourier-Laplace representation. By using the Mellin transform, we provide a general representation of the solution in terms of Mellin-Barnes integrals in the complex plane, which allows us to extend the probability interpretation known for the standard diffusion equation to suitable ranges of the relevant parameters alpha and beta. We derive explicit formulae (convergent series and asymptotic expansions), which enable us to plot the corresponding spatial probability densities.

研究动机与目标

• 推导具有通用参数 α ∈ (0,2], β ∈ (0,2], 和偏度 θ 的时空分数阶扩散方程的基本解。
• 将格林函数的概率密度解释从标准情形（如空间或时间分数阶扩散）扩展至更广泛的参数范围。
• 利用梅林-巴恩斯积分开发一个针对约化格林函数 Kα,βθ(x) 的通用计算框架。
• 为在所有物理相关参数范围内准确进行数值计算，提供收敛级数和渐近展开式。
• 展示代表性参数值下格林函数的可视化图像，以说明其在不同 α, β, 和 θ 下的行为。

提出的方法

• 通过时空分数阶扩散方程的傅里叶-拉普拉斯变换推导格林函数。
• 利用相似变量 x/tβ/α 建立标度不变性，将解表示为 Gα,βθ(x,t) = t−β/α Kα,βθ(x/tβ/α)。
• 利用梅林-巴恩斯积分表示，将格林函数解析延拓至其仍可解释为概率密度的参数范围。
• 利用梅林-巴恩斯表示推导 Kα,βθ(x) 的收敛级数展开式，适用于所有 α ∈ (0,2], β ∈ (0,2], |θ| ≤ min{α,2−α}。
• 推导 Kα,βθ(x) 在 x → ∞ 时的渐近展开式，包括幂律衰减和拉伸指数衰减，对 α > 1 且 β > 1 的情形给出显式系数。
• 采用收敛展开式与渐近展开式之间的匹配策略，确保在整个 x 定义域内保持数值精度。

实验结果

研究问题

• RQ1 时空分数阶扩散方程的基本解如何用特殊函数和积分变换表示？
• RQ2 格林函数在哪些参数范围内可被解释为概率密度函数？
• RQ3 对于一般 α, β, 和 θ，约化格林函数 Kα,βθ(x) 的解析结构是什么？
• RQ4 格林函数的渐近行为如何依赖于参数 α, β, 和 θ，特别是在重尾和拉伸指数衰减区域？
• RQ5 能否开发一个统一的计算框架，以在所有物理相关参数值下准确计算格林函数？

主要发现

• 格林函数表现出通用标度形式 Gα,βθ(x,t) = t−β/α Kα,βθ(x/tβ/α)，其中约化格林函数 Kα,βθ(x) 仅依赖于相似变量。
• 通过梅林-巴恩斯积分表示，严格将概率密度解释扩展至参数范围 {0 < α ≤ 2} ∩ {0 < β ≤ 1} 和 {1 < β ≤ α ≤ 2}。
• 当 α = 0.5, β = 0.5, θ = 0 时，约化格林函数在 |x| 较大时按幂律衰减，指数为 −1.5，与稳定分布行为一致。
• 当 1 < α < 2 且 β = 1 时，渐近形式 Kα,βθ(x) ∼ A x^a e^{−b x^c} 成立，其中显式系数 A, a, b, c 依赖于 α 和 β。
• 当 α = 1.5, β = 1.25, θ = −0.50 时，渐近衰减形式为 x^{−0.5} e^{−1.25 x^{3}}，其中 c = 3，表明强局域化。
• 图像结果表明，格林函数从对称的Lévy稳定形式（α < 2, β = 1）过渡到高斯型（α = 2, β → 2），且当 θ ≠ 0 时表现出偏度。