[论文解读] Solutions of certain fractional kinetic equations and a fractional diffusion equation

R. K. Saxena, A. M. Mathai+1 more|arXiv (Cornell University)|Apr 15, 2007

Fractional Differential Equations Solutions被引用 6

一句话总结

本文利用拉普拉斯变换和广义Mittag-Leffler函数，推导出统一分数阶动力学方程和分数阶扩散方程的显式解。主要贡献在于以H-函数和广义Mittag-Leffler函数表示的紧凑且可解析处理的解，扩展了先前结果，并为具有长记忆效应的异常扩散与动力学过程提供了统一框架。

ABSTRACT

In view of the usefulness and importance of the kinetic equation in certain physical problems, the authors derive the explicit solution of a fractional kinetic equation of general character, that unifies and extends earlier results. Further, an alternative shorter method based on a result developed by the authors is given to derive the solution of a fractional diffusion equation.

研究动机与目标

• 将现有的分数阶动力学方程解统一并扩展为一个单一的广义框架。
• 为Kochubei（1990）先前推导出的分数阶扩散方程提供一种替代且更简短的证明。
• 推导解的级数表示和渐近展开，以支持数值计算与解析分析。
• 作为解过程的副产品，建立单边Lévy稳定密度的H-函数表示。
• 使系统在长程相互作用和Lévy噪声下的应用成为可能，例如在异常扩散和复杂随机过程中。

提出的方法

• 使用拉普拉斯变换方法求解涉及多个Riemann-Liouville分数阶积分的统一分数阶动力学方程。
• 应用卷积定理和广义Mittag-Leffler函数的级数展开，推导出以 $ E^{(l+1)}_{\nu_1, \nu} $ 为核的积分形式的解。
• 利用反拉普拉斯变换公式 $ L^{-1}\{s^{-\gamma}(1 - a s^{-\beta})^{-\delta}; t\} = t^{\gamma-1} E^{\delta}_{\beta,\gamma}(a t^{\beta}) $ 恢复时域解。
• 利用特殊函数的H-函数表示，以紧凑形式表达解，适用于渐近分析与级数分析。
• 应用已知的H-函数恒等式，对解进行简化与重表达，包括幂律标度和标度不变性。
• 利用已知的积分表示和伽马函数恒等式，推导解的级数展开与渐近行为。

实验结果

研究问题

• RQ1 如何利用积分变换以闭式表达包含多个分数阶导数的一般分数阶动力学方程的解？
• RQ2 分数阶扩散方程的解在特殊函数中的解析结构是什么？
• RQ3 当空间或时间变量较大或较小时，解的渐近行为如何？
• RQ4 能否以与解框架一致的方式，通过H-函数表示单边Lévy稳定密度？
• RQ5 广义Mittag-Leffler函数在统一不同分数阶动力学与扩散模型的解中起什么作用？

主要发现

• 统一分数阶动力学方程的解以涉及广义Mittag-Leffler函数的双重级数与积分形式表达，且给出了收敛条件。
• 对于分数阶扩散方程，解以两种等价形式推导得出：一种涉及 $ H^{2,0}_{1,2} $ H-函数，另一种为带缩放参数的 $ H^{2,0}_{1,2} $ 形式。
• 当 $ n=1 $ 时，解由包含 $ \Gamma(1 - \alpha(l+1)/2) $ 的级数给出，表现出时间上的幂律衰减。
• 当 $ n=2 $ 时，解表现出对数发散 $ \sim \ln(t^{\alpha/2}/x) $，与二维情况下的已知结果一致。
• 当 $ n=3 $ 时，解在大 $ x $ 处衰减为 $ 1/x $，表明存在长程空间依赖性。
• 作为解过程的副产品，推导出单边Lévy稳定密度的H-函数表示：$ \Phi_\rho(t) = \frac{1}{\rho} H^{1,0}_{1,1}\left[ \frac{1}{t} \middle| {}^{(1,1)}_{(1/\rho,1/\rho)} \right] $，与先前工作一致。