[論文レビュー] Numerical Solutions of a Boundary Value Problem for the Anomalous Diffusion Equation with the Riesz Fractional Derivative

Mariusz Ciesielski, Jacek Leszczyński|ArXiv.org|Jun 27, 2005

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

ひとこと要約

本稿では、Riesz-Feller分数階微分を伴う1次元境界値問題の数値的有限差分法を提示する。この手法は、分数階微分の非局所的性質に対処するために陰解法と陽解法を併用し、シミュレーションによりα < 2のとき、初期段階では拡散が遅くなるが、長尾分布を示し、Lévy飛行の挙動を捉える。また、歪度パラメータθを変化させることで非対称な確率密度関数が得られる。

ABSTRACT

In this paper we present in one-dimensional space a numerical solution of a partial differential equation of fractional order. This equation describes a process of anomalous diffusion. The process arises from the interactions within the complex and non-homogeneous background. We presented a numerical method which bases on the finite differences method. We considered pure initial and boundary-initial value problems for the equation with the Riesz-Feller fractional derivative. In the final part of this paper sample results of simulation were shown.

研究の動機と目的

• 有界領域におけるRiesz-Feller分数階微分を伴う拡散方程式を解くための数値的手法を開発すること。
• 標準的なブラウン運動から逸脱する異常拡散過程、特に長尾を持つLévy飛行をモデル化すること。
• 分数階微分の非局所的かつ歴史依存性を扱えるように、古典的な有限差分法を拡張すること。
• 歪度パラメータθおよび次数αが解の挙動に与える影響を分析すること、特に波動的ダイナミクスへの遷移を含む。

提案手法

• 分数階拡散方程式に含まれる時間および空間微分を、Riesz-Feller作用素を用いて有限差分法で離散化する。
• Riesz-Feller微分は、αおよびθに依存する係数を有する左および右のRiemann-Liouville分数階微分の重み付き和として近似される。
• 陰解法（σ = 0.5）と陽解法（σ = 1）が定式化され、陰解法は安定性を保証する。
• 境界条件は、gLおよびgRを含む源項を介して組み込まれ、分数階微分の非局所性のため、各時刻における解は、すべての過去の空間点に依存する。
• α = 2およびθ = 0のとき、このスキームは古典的なFTCS法に還元され、標準的拡散との一貫性が検証される。
• 歪度パラメータθの調整により非対称な確率密度関数を扱えるようにし、指向的輸送のモデル化が可能になる。

実験結果

リサーチクエスチョン

• RQ1 有界領域におけるDirichlet境界条件を満たすRiesz-Feller分数階微分を伴う方程式に対して、安定的かつ高精度な有限差分スキームを構築する方法は何か？
• RQ2 標準的拡散（α = 2）と比較して、分数階次数α ∈ (1,2] が拡散ダイナミクスに与える影響は何か？
• RQ3 歪度パラメータθが解の形状および非対称性に与える影響は何か、特にα → 1+およびθ → ±1+の極限においては？
• RQ4 提案されたスキームは、異常拡散に特徴的なLévy飛行の長尾挙動を正確に捉えることができるか？
• RQ5 分数階拡散方程式が、双曲型（波動的）な挙動に遷移する条件は何か？

主な発見

• α < 2のとき、解は標準的ブラウン運動（α = 2）と比較して初期段階で拡散が遅く、これは準拡散的挙動を示す。
• 長時間における確率密度関数は、Lévy飛行に特徴的な長尾を示し、極端な粒子のジャンプが希に発生する。
• α → 1+およびθ → ±1+のとき、解は1階微分の波動方程式に収束し、拡散的輸送から球形輸送への遷移を示す。
• 歪度パラメータθにより非対称な確率密度関数が生成され、θ > 0では右方向への伝播が促進され、θ < 0では左方向への伝播が促進される。
• 陰解法は安定性を保証するが、α → 1+では陽解法は無条件に不安定であることが確認され、理論的限界が裏付けられる。
• α = 2およびθ = 0のとき、数値スキームは古典的なFTCS法に正確に還元され、標準的拡散との一貫性が検証される。