[論文レビュー] Homotopy perturbation method for fractional-order Burgers-Poisson equation

Caibin Zeng, Qigui Yang+1 more|arXiv (Cornell University)|Mar 9, 2010

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

ひとこと要約

本稿では、時間微分を Caputo の分数階微分（$\alpha \in (0,1]$）に置き換えることで、分数階 Burgers-Poisson 方程式を導入する。ホモトピー摂動法（HPM）を用いて、正確かつ近似された明示的解が導出され、この手法は非線形分数階微分方程式を解く際に極めて効果的かつ正確であることが示され、特に $\alpha = 1$ の場合に既知の正確解が回復されることを確認する。

ABSTRACT

In this paper, the fractional-order Burgers-Poisson equation is introduced by replacing the first-order time derivative by fractional derivative of order $α$. Both exact and approximate explicit solutions are obtained by employing homotopy perturbation method. The numerical results reveal that the proposed method is very effective and simple for handling fractional-order differential equations.

研究の動機と目的

• 古典的 Burgers-Poisson 方程式を Caputo 時間分数階微分を用いて分数階形式に拡張すること。
• 得られた分数階系を解くためにホモトピー摂動法（HPM）を構築し、適用すること。
• 分数階 Burgers-Poisson 方程式に対して正確解および近似解を明示的に導出すること。
• 正確解との数値的比較を通じて、HPM の精度と効率を検証すること。

提案手法

• 時間微分を $\alpha \in (0,1]$ の Caputo 分数階微分に置き換えることで、分数階 Burgers-Poisson 方程式を定式化する。
• 埋め込みパラメータ $p \in [0,1]$ を用いてホモトピーを構築し、元の式をより単純な部分問題の系列へ変形する。
• 解 $u(x,t)$ を $p$ のべき級数として展開する：$u = u_0 + p u_1 + p^2 u_2 + \cdots$。
• 各摂動級数の項に分数階積分作用素 $J^\alpha$ を適用し、補題1を用いて $u_j(x,t)$ を再帰的に解く。
• 高次項 $u_j(x,t)$（$j \geq 4$）を求めるために記号計算ソフトウェア（Mathematica）を用いる。
• 級数の最初の4項を組み合わせ、$u(x,t)$ に対する4項近似明示解を構成する。

実験結果

リサーチクエスチョン

• RQ1 ホモトピー摂動法は、非線形分散項を含む分数階 Burgers-Poisson 方程式を効果的に解くことができるか？
• RQ2 \alpha = 1 の場合に、HPM に基づく近似解は正確解とどの程度一致するか？
• RQ3 \alpha を変化させた際の解の挙動はどのように変化するか？また、\alpha \to 1 の極限で、古典的解が保存されるか？
• RQ4 分数階 Burgers-Poisson 方程式に対して、HPM の級数解の収束性と安定性はいかがなっているか？

主な発見

• HPM により、分数階 Burgers-Poisson 方程式に対する4項近似明示解が得られた：$u(x,t) = x - \frac{1+x}{\Gamma(1+\alpha)}t^{\alpha} + \frac{2(1+x)}{\Gamma(2+\alpha)}t^{2\alpha} - \frac{(1+x)(4\Gamma(1+\alpha)^2 + \Gamma(2+\alpha))}{\Gamma(1+\alpha)^2\Gamma(3+\alpha)}t^{3\alpha}$。
• $\alpha = 1$ の場合に、HPM は正確解 $u(x,t) = \frac{1+x}{1+t} - 1$ を回復し、古典的結果と整合性があることが確認された。
• 表1における数値的比較では、近似解が正確解と良好に一致しており、テストされた $x$ および $t$ 値に対して絶対誤差は 0.002462 から 0.083333 の範囲に収まっている。
• 誤差は $t$ が小さくなるほど減少しており、時間間隔が小さいほど精度が向上することを示している。
• 非線形項 $3u_x u_{xx} + u u_{xxx}$ を含む複雑な非線形項を有する微分方程式に対しても、本手法は極めて効果的かつ簡便であることが示された。
• 結果から、HPM は特に完全非線形分散項を有する分数階非線形偏微分方程式を解く強力で信頼性の高い手法であると示唆される。