[論文レビュー] 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

ひとこと要約

本稿は、リーマン＝フェラー空間分数級微分とキャプート時間分数級微分を有する空間時間分数級拡散方程式の基本解（グリーン関数）を導出する。フーリエ＝ラプラス変換およびメリン＝バーナーズ積分表現を用いて、スケーリング性質を確立し、還元されたグリーン関数に対する収束級数および漸近展開を明示的に提供することで、非対称的かつ非マルコフ的状況を含む多様なパrameter範囲における確率密度解釈を可能にする。

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.

研究の動機と目的

• 一般パrameter α ∈ (0,2], β ∈ (0,2], および歪度 θ を有する空間時間分数級拡散方程式の基本解を導出すること。
• 標準的でない状況（例：空間または時間分数級拡散）を超えて、より広範なパrameter範囲におけるグリーン関数の確率密度解釈を拡張すること。
• 還元グリーン関数 Kα,βθ(x) のための、一般化された計算フレームワークをメリン＝バーナーズ積分を用いて構築すること。
• 物理的に関連するすべてのパrameter範囲にわたるグリーン関数の正確な数値計算を可能にする収束級数および漸近展開を提供すること。
• 代表的なパrameter値に対するグリーン関数の可視化を提示し、α, β, および θ の変化に伴う挙動を示すこと。

提案手法

• 空間時間分数級拡散方程式のフーリエ＝ラプラス変換を用いてグリーン関数を導出する。
• 相似変数 x/tβ/α を通じたスケーリング不変性を確立し、解を Gα,βθ(x,t) = t−β/α Kα,βθ(x/tβ/α) の形で表現する。
• グリーン関数をパrameter範囲にわたり解析接続するためにメリン＝バーナーズ積分表現を用いる。
• メリン＝バーナーズ表現を用いて、すべての α ∈ (0,2], β ∈ (0,2], |θ| ≤ min{α,2−α} に対して Kα,βθ(x) の収束級数展開を導出する。
• x → ∞ における Kα,βθ(x) の漸近展開を導出し、べき則および伸びた指数的減衰を含む。α > 1 および β > 1 の場合に明示的な係数を提供する。
• 収束級数と漸近展開の間の一致戦略を適用し、x 全体にわたる数値的正確性を保証する。

実験結果

リサーチクエスチョン

• RQ1 空間時間分数級拡散方程式の基本解は、特殊関数および積分変換の観点からどのように表現できるか？
• RQ2 グリーン関数が確率密度関数として解釈可能なパrameter範囲はどのようになるか？
• RQ3 一般の α, β, θ に対して還元グリーン関数 Kα,βθ(x) の解析的構造はどのようなものか？
• RQ4 グリーン関数の漸近的挙動は、特に重尾および伸びた指数的領域において、パrameter α, β, θ にどのように依存するか？
• RQ5 物理的に関連するすべてのパrameter値にわたるグリーン関数を正確に計算するための統一的計算フレームワークを構築できるか？

主な発見

• グリーン関数は、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 であるため強い局在化を示す。
• プロットは、グリーン関数が対称的レヴィ安定分布形（α < 2, β = 1）からガウス型（α = 2, β → 2）に移行し、θ ≠ 0 の場合に歪度を示すことを確認している。