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

ひとこと要約

本稿では、ラプラス変換および一般化されたミタグ・レフラー関数を用いて、統一的分数階動的方程式および分数階拡散方程式の明示的解を導出する。主な貢献は、H関数および一般化されたミタグ・レフラー関数を用いた、コンactで解析的に取り扱いやすい解の形式であり、既存の結果を拡張し、長記憶効果を有する異常拡散および動的過程を統合的に扱うフレームワークを提供する。

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) が以前に導出した分数階拡散方程式の解に対して、代替的で短い証明を提供すること。
• 数値計算および解析的分析のための、解の級数表現および漸近展開を導出すること。
• 解の導出過程の副産物として、片側レビ安定分布のH関数表現を確立すること。
• 長距離相互作用およびレビノイズを有する系、例えば異常拡散および複雑な確率過程への応用を可能にすること。

提案手法

• 複数のリーマン＝リウヴィル分数階積分を含む統一的分数階動的方程式を解くために、ラプラス変換法を用いる。
• 畳み込み定理および一般化されたミタグ・レフラー関数の級数展開を適用し、$ 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関数表現を活用し、漸近的および級数的解析に適したコンactな形で解を表現する。
• 既知のH関数の恒等式を適用して、解を代替形式に簡略化・再表現し、べき則スケーリングおよびスケーリング不変性を含む。
• 既知の積分表現およびガンマ関数の恒等式を用いて、解の級数展開および漸近的挙動を導出する。

実験結果

リサーチクエスチョン

• RQ1 複数の分数階微分を含む一般化された分数階動的方程式は、積分変換を用いて閉形式でどのように解けるか？
• RQ2 分数階拡散方程式の解の解析的構造は、特殊関数の観点からどのように記述できるか？
• RQ3 空間的または時間的変数が非常に大きいか非常に小さい場合、解はどのように漸近的に振る舞うか？
• RQ4 片側レビ安定分布は、解の枠組みと整合する形でH関数を用いて表現可能か？
• RQ5 一般化されたミタグ・レフラー関数は、異なる分数階動的および拡散モデル間の解を統合する役割を果たすか？

主な発見

• 統一的分数階動的方程式の解は、一般化されたミタグ・レフラー関数を含む二重級数および積分形式で表され、収束条件が明示されている。
• 分数階拡散方程式に関しては、二通りの同等な形で解が得られた：一つは $ 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) $ のように対数発散を示し、2次元における既知の結果と整合的である。
• n=3 の場合、解は大きな $ x $ において $ 1/x $ のように減衰し、長距離空間的依存性を示している。
• 片側レビ安定分布のH関数表現として、$ \Phi_\rho(t) = \frac{1}{\rho} H^{1,0}_{1,1}\left[ \frac{1}{t} \middle| {}^{(1,1)}_{(1/\rho,1/\rho)} \right] $ が導出され、先行研究と整合的であることが確認された。