[論文レビュー] Fractional Calculus in Wave Propagation Problems

Francesco Mainardi|arXiv (Cornell University)|Feb 1, 2012

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

ひとこと要約

本稿では、線形粘弾性媒体における拡散と古典的波動運動の中間的波動伝播現象をモデル化する枠組みとして、分数階微積分を提案する。波動方程式における整数階時間微分を 0 < β ≤ 2 の分数量微分に置き換えることで、べき則的クリープ挙動を捉え、M-Wright関数に従う解が得られ、拡散的および振動的ダイナミクスの間を滑らかに補間する。

ABSTRACT

Fractional calculus, in allowing integrals and derivatives of any positive order (the term "fractional" kept only for historical reasons), can be considered a branch of mathematical physics which mainly deals with integro-differential equations, where integrals are of convolution form with weakly singular kernels of power law type. In recent decades fractional calculus has won more and more interest in applications in several fields of applied sciences. In this lecture we devote our attention to wave propagation problems in linear viscoelastic media. Our purpose is to outline the role of fractional calculus in providing simplest evolution processes which are intermediate between diffusion and wave propagation. The present treatment mainly reflects the research activity and style of the author in the related scientific areas during the last decades.

研究の動機と目的

• 拡散と波動伝播の中間的過程をモデル化する分数量微積分の役割を調査すること。
• べき則的クリープを示す線形粘弾性材料における波動伝播を分析すること。これは、粘性流体と弾性固体の中間の挙動を示す。
• 時間領域における分数量微分の次数 β（0 < β ≤ 2）が、標準の拡散方程式および波動方程式を一般化する積分微分方程式を導く仕組みを示すこと。
• 初期値問題における物理的応用において、標準的な初期条件と整合性を持つことから、Caputo分数量微分が最も適したツールであることを確立すること。
• 分数量波動方程式の基本解としてのM-Wright関数の解析的および漸近的性質を調査すること。

提案手法

• 標準の波動方程式または拡散方程式における1階または2階時間微分を、Caputo分数量微分（0 < β ≤ 2）に置き換えることで、一般化された発展方程式を構築する。
• 整数階微分を含む物理的初期条件と整合性を持つように、分数量微分のCaputo定義を用いる。
• 分数量積分を適用することで、分数量波動方程式の等価な積分微分形式を導出し、初期条件を明示的に組み込む。
• 分数量微分方程式にラプラス変換を適用し、周波数領域における解表現を導出する。
• M-Wright関数 Mν(x) を用いて基本解を特徴づけ、解析的漸近展開および数値プロットを用いてその振る舞いを分析する。
• 線形および対数スケールを用いて、ν 値（0 から 1 の間）における M-Wright 関数の挙動を比較し、ν → 1⁻ がディラックのデルタ関数に対応する極限を検討する。

実験結果

リサーチクエスチョン

• RQ1 分数量微積分を用いて、拡散と古典的波動運動の中間的波動伝播をどのようにモデル化できるか？
• RQ2 波動伝播問題において、1階または2階の時間微分を 0 < β ≤ 2 の分数量微分に置き換える物理的意味は何か？
• RQ3 粘弾性波動モデルの初期値問題において、Riemann-Liouville微分よりもなぜ Caputo微分が好まれるのか？
• RQ4 M-Wright関数 Mν(x) は、分数量波動方程式の基本解をどのように記述するか？また、β の変化に伴いその形状はどのように変化するか？
• RQ5 特に ν → 1⁻ のとき、原点付近および無限遠における M-Wright 関数の漸近的および数値的挙動は何か？

主な発見

• Caputo微分を用いた分数量波動方程式（β（0 < β ≤ 2））の解は、β < 1 のとき拡散的挙動、β > 1 のとき振動的挙動の間を滑らかに補間する。
• β = 1 のとき、方程式は標準の拡散方程式に簡略化される。β = 2 のとき、古典的波動方程式に還元される。
• 分数量波動方程式の基本解は、ν = β/2 − 1/2 として定義される M-Wright 関数 Mν(x) で表され、ν = 0 のとき指数的減衰、ν = 1/2 のときガウス型、ν = 1 のときディラックのデルタ関数へと変化する。
• β → 2⁻（ν → 1⁻）に伴い、M-Wright 関数は x ≈ ±1 付近に鋭いピークを形成し、有限の前線速度を伴う波動的伝播の兆候が現れる。
• ν が 1 に近い（例：ν = 0.99 または ν = 0.999）場合の Mν(x) の数値プロットは、x ≈ 1 の周辺に高振幅で狭いピークを示し、波動的挙動への移行を確認する。
• Caputo微分のラプラス変換は、初期条件の物理的形を保ち続けるため、測定可能な初期状態を持つ現実の粘弾性系のモデル化に最適である。