[論文レビュー] Solutions of fractional reaction-diffusion equations in terms of the H-function

H. J. Haubold, A. M. Mathai+1 more|ArXiv.org|Apr 3, 2007

Fractional Differential Equations Solutions参考文献 14被引用数 11

ひとこと要約

本稿では、ラプラス変換およびフーリエ変換を用いて、Caputo時間分数導関数とRiesz-Feller空間分数導関数を含む統一された分数マクロ反応拡散方程式の解析的解を導出する。解はH関数を用いた閉形式で表され、空間時間、時間分数、空間分数拡散の既知の結果を一般化し、Lévy安定密度やガウス拡散を極限形として含む。

ABSTRACT

This paper deals with the investigation of the solution of an unified fractional reaction-diffusion equation associated with the Caputo derivative as the time-derivative and Riesz-Feller fractional derivative as the space-derivative. The solution is derived by the application of the Laplace and Fourier transforms in closed form in terms of the H-function. The results derived are of general nature and include the results investigated earlier by many authors, notably by Mainardi et al. (2001, 2005) for the fundamental solution of the space-time fractional diffusion equation, and Saxena et al. (2006a, b) for fractional reaction- diffusion equations. The advantage of using Riesz-Feller derivative lies in the fact that the solution of the fractional reaction-diffusion equation containing this derivative includes the fundamental solution for space-time fractional diffusion, which itself is a generalization of neutral fractional diffusion, space-fractional diffusion, and time-fractional diffusion. These specialized types of diffusion can be interpreted as spatial probability density functions evolving in time and are expressible in terms of the H-functions in compact form.

研究の動機と目的

• Caputo時間導関数とRiesz-Feller空間導関数を有する分数反応拡散方程式の統一的解析的解を導出すること。
• Riesz-Feller導関数を組み込むことにより、空間時間、時間分数、空間分数拡散に関する既存の結果を一般化すること。
• 解をH関数を用いて表現することで、複雑な拡散過程を簡潔に記述すること。
• 解がガウス拡散やLévy安定密度といった既知の特殊ケースを極限形として回復できること。

提案手法

• 初期条件および分数階動的特性を扱うために、Caputo時間分数導関数にラプラス変換を適用する。
• Riesz-Feller空間分数導関数を、その特性関数 $ \Psi^\theta_\alpha(k) = |k|^\alpha \exp[i(\text{sign}\,k)\frac{\theta\pi}{2}] $ を用いて表現する。
• 変換領域における解の導出の後、H関数の既知の積分公式を用いて逆変換を行う。
• H関数のコサイン変換公式を用いて、フーリエ逆変換をH関数の形で評価する。
• H関数のスケーリング性質 $ H^{m,n}_{p,q}[x^\delta] = \frac{1}{\delta}H^{m,n}_{p,q}[x] $ （$ \delta > 0 $）を用いて式を簡略化する。
• ガウス、Lévy、中立的分数拡散密度といった既知の特殊ケースを回復することで、解の正当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 CaputoおよびRiesz-Feller導関数を有する統一的分数反応拡散方程式を、閉形式で解析的に解く方法は何か？
• RQ2 H関数は、異なる分数階数においてこのような方程式の基本解を表現する上で果たす役割は何か？
• RQ3 Riesz-Feller導関数の組み込みが、古典的および分数拡散モデルをどのように一般化するか？
• RQ4 解が時間分数、空間分数、標準拡散を極限形として回復できるか？
• RQ5 反応項が非線形または時間依存である場合、解の解析的構造はどのようになるか？

主な発見

• 統一的分数反応拡散方程式の解は、$ N(x,t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} f^*(k) E_{1/2,1}(-\eta t^\beta \Psi^\theta_\alpha(k)) \exp(-ikx) dk + \cdots $ の形で閉形式で表され、H関数とMittag-Leffler関数を含む。
• 空間時間分数拡散の基本解は、$ \alpha = \beta $、$ \theta \leq \min(\alpha, 2 - \alpha) $、$ \beta = 1 $ の特別な場合に回復され、H関数表現が得られる。
• $ \beta = 1 $、$ 0 < \alpha < 2 $、$ \theta = 0 $ の場合、解は $ \alpha $-安定Lévy密度 $ L^\theta_\alpha(x) = \frac{1}{\alpha (\eta t)^{1/\alpha}} H^{1,1}_{2,2}\left[ \frac{(|x|)}{(\eta t)^{1/\alpha}} \middle| \cdots \right] $ に簡略化される。
• $ \alpha = 2 $、$ \beta = 1 $、$ \theta \to 0 $ のとき、解は古典的ガウス密度 $ N(x,t) = (4\pi\eta t)^{-1/2} \exp\left(-\frac{|x|^2}{4\eta t}\right) $ に還元され、標準拡散と整合性が確認される。
• $ \beta = 1/2 $ の場合、Mittag-Leffler関数 $ E_{1/2,1/2} $ を含む新しい形の解が得られ、半階時間分数反応拡散モデルに適用可能である。
• H関数表現により、中立的分数拡散、空間分数拡散、時間分数拡散といったすべての特殊ケースが、統一的な数学的枠組みで簡潔に表現可能である。