[论文解读] From Power Laws to Fractional Diffusion: the Direct Way

Rudolf Gorenflo, E. A. Abdel-Rehim|ArXiv.org|Dec 30, 2007

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

一句话总结

本文從具有幂律尾部等待时间与跳跃的连续时间随机游走（CTRWs）出發，直接推導出时空分数阶扩散方程。透過應用良好尺度的極限與分佈函數的Tauber型引理，證明了停留時間密度弱收斂至一個具有時間導數階β ∈ (0,1) 與空間導數階α ∈ (0,2) 的分數階擴散方程的解，該方程由Mittag-Leffler函數與分數階微分所控制。

ABSTRACT

Starting from the model of continuous time random walk, we focus our interest on random walks in which the probability distributions of the waiting times and jumps have fat tails characterized by power laws with exponent between 0 and 1 for the waiting times, between 0 and 2 for the jumps. By stating the relevant lemmata (of Tauber type) for the distribution functions we need not distinguish between continuous and discrete space and time. We will see that, by a well-scaled passage to the diffusion limit, generalized diffusion processes, fractional in time as well as in space, are obtained. The corresponding equation of evolution is a linear partial pseudo-differential equation with fractional derivatives in time and in space, the orders being equal to the above exponents. Such processes are well approximated and visualized by simulation via various types of random walks. For their explicit solutions there are available integral representations that allow to investigate their detailed structure.

研究动机与目标

• 建立具有重尾等待時間與跳躍的連續時間隨機游走（CTRWs）與分數階擴散過程之間的直接聯繫。
• 證明此類CTRWs的良好尺度極限會產生一個具有時間與空間分數階導數的時空分數階擴散方程的解。
• 利用分佈函數與拉普拉斯與傅里葉變換的漸近分析，統一處理連續與離散的時空過程。
• 為分數階擴散的出現提供嚴謹的理論基礎，無需依賴次屬化或中間尺度步驟。

提出的方法

• 以CTRWs的Montroll-Weiss方程為起點，透過停留時間密度的拉普拉斯與傅里葉變換推導得出。
• 將Tauber型引理應用於等待時間與跳躍的累積分佈函數，以統一處理連續與離散分佈。
• 在空間步長h與等待時間τ之間施加一個尺度關係h ∝ τ^β/α，以確保獲得良好尺度的擴散極限。
• 分析等待時間與跳躍分佈在幂律尾部（指數β ∈ (0,1) 與α ∈ (0,2)）下的拉普拉斯與傅里葉變換的漸近行為。
• 推導出極限演化方程為一個線性偏偽微分方程，其時間導數階為β，空間導數為α階的Riesz分數階微分。
• 利用包含Mittag-Leffler函數的積分表達式來表徵解，並驗證其收斂至分數階擴散過程。

实验结果

研究问题

• RQ1 具有幂律尾部等待時間與跳躍的連續時間隨機游走的良態極限是否能產生分數階擴散過程？
• RQ2 空間與時間增量之間的何種尺度關係可確保收斂至時空分數階擴散方程？
• RQ3 如何利用分佈函數與拉普拉斯與傅里葉變換的漸近分析，嚴謹地建立收斂性？
• RQ4 Mittag-Leffler函數在表徵分數階擴散方程極限解的過程中發揮何種作用？
• RQ5 為何標準的CTRW框架在具有幂律等待時間時，除非經過適當尺度化，否則無法產生分數階擴散？

主要发现

• 經尺度化的CTRW之停留時間密度弱收斂至具有時間導數階β ∈ (0,1) 與空間導數階α ∈ (0,2) 的時空分數階擴散方程的解。
• 極限方程由Caputo型時間分數階導數與Riesz空間分數階導數給出：$ {}_0^C D_t^eta u(x,t) = R^eta u(x,t) $，初始條件為$ u(x,0) = ho(x) $。
• 解透過Mittag-Leffler函數明確表達，此函數自然來自具有幂律尾部的等待時間分佈的拉普拉斯變換。
• 收斂性透過分佈函數的Tauber引理建立，使連續與離散分佈得以在統一框架下處理。
• 該方法避免了次屬化的需要，直接分析CTRW轉移機率的尺度極限。
• 結果確認分數階擴散的出現並非來自任意的幂律，而是僅在空間與時間增量之間存在特定良好尺度關係時才會發生。