[論文レビュー] Numerical solutions to boundary value problem for anomalous diffusion equation with Riesz-Feller fractional operator

Mariusz Ciesielski, Jacek Leszczyński|arXiv (Cornell University)|Jul 5, 2006

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

ひとこと要約

本稿では、1次元空間におけるリーマン＝フィラー分数階微分を含む境界値問題を解くための分数階有限差分法（FFDM）を提示する。この手法により、異常拡散の数値シミュレーションが可能となる。本手法は古典的有限差分スキームの一般化であり、主な結果として、α < 2 の場合に非線形温度プロファイルが得られ、α = 0.35 および θ = -0.055 を用いて実験的ナノチューブ熱伝導データに良好にフィットすることが示された。

ABSTRACT

In this paper, we present a numerical solution to an ordinary differential equation of a fractional order in one-dimensional space. The solution to this equation can describe a steady state of the process of anomalous diffusion. The process arises from interactions within complex and non-homogeneous background. We present a numerical method which is based on the finite differences method. We consider a boundary value problem (Dirichlet conditions) for an equation with the Riesz-Feller fractional derivative. In the final part of this paper, same simulation results are shown. We present an example of non-linear temperature profiles in nanotubes which can be approximated by a solution to the fractional differential equation.

研究の動機と目的

• 異常拡散を含むリーマン＝フィラー分数階微分を含む境界値問題を解くための数値的手法を開発すること。
• 標準的な拡散法では失敗する、ナノチューブのような複雑で非均質な媒体における定常状態温度プロファイルをモデル化すること。
• 長距離空間的相関や重い尾を持つ粒子のジャンプを捉えるために、長距離依存性を持つ分数階微分に一般化された古典的有限差分法を拡張すること。
• ナノチューブからの実験データと数値解を比較することで、手法の妥当性を検証すること。
• 歪度パラメータ θ および次数 α が解の対称性と形状に与える影響を調査すること。

提案手法

• 空間的リーマン＝フィラー分数階微分作用素を離散化するために、有限差分法（FDM）を用いる。
• 関数値の重み付き和を用いてリーマン＝フィラー微分の数値近似を導出する。重み係数は α および θ に依存する。
• 行列 A が分数階微分ステンシルを記述する線形方程式系 A·T = B を構築する。ベクトル B はディリクレ境界条件を組み込む。
• α = 2 および θ = 0 の場合、本スキームは2階微分の古典的中央差分法に帰着する。
• 非局所的挙動を考慮：各点の微分値は、局所的な近傍値だけでなく、領域全体の値に依存する。
• 線形方程式系を解くことで解が得られ、分数階微分の非局所的性質のため、境界条件がすべての内部ノードに影響を与える。

実験結果

リサーチクエスチョン

• RQ1 境界値問題において、リーマン＝フィラー分数階微分を有限差分スキームで正確に近似する方法は何か？
• RQ2 分数階の次数 α および歪度パラメータ θ を変化させた場合、定常状態温度プロファイルの形状にどのような影響があるか？
• RQ3 提案された分数階有限差分法は、標準的拡散法では失敗するナノチューブにおける実験的温度プロファイルを正確に再現できるか？
• RQ4 α → 1+ および θ → ±1+ に近づけると解の挙動はどのように変化し、これはどのような物理的解釈に対応するか？
• RQ5 分数階微分パラメータと異常拡散における長尾確率分布の出現との間にはどのような関係があるか？

主な発見

• α = 2 および θ = 0 の場合、解は線形であり、古典的熱方程式との一貫性が確認された。
• α < 2 の場合、解は非線形プロファイルを示し、異常拡散の挙動を示している。
• α → 1+ および θ → ±1+ に近づくと、解は1階波動方程式の定常状態に近づく。
• 歪度パラメータ θ は解の非対称性を導入し、θ ∈ (0,1) の場合に非対称なプロファイルが得られる。
• α = 0.35 および θ = -0.055 のモデルが、Zhang と Li (2005) の実験的ナノチューブ温度データに最も良好にフィットした。
• 分数階微分は長尾粒子ジャンプを捉え、ガウス統計では表現できないが、まれに発生する極端な出来事のモデル化を可能にする。