[論文レビュー] A Finite Difference Scheme based on Cubic Trigonometric B-splines for Time Fractional Diffusion-wave Equation

Muhammad Abbas|arXiv (Cornell University)|May 23, 2017

Fractional Differential Equations Solutions被引用数 3

ひとこと要約

本稿では、反応項を含む時間分数拡散波方程式を解くために、立方体三角的Bスプラインを用いた新規な有限差分スキームを提案する。コッホの時間微分（1 < γ ≤ 2）に対して標準的な2次精度の有限差分近似を用い、空間微分は立方体三角的Bスプラインを用いて離散化することで、高精度（最大10⁻¹¹）かつ無条件安定性を達成し、既存のヘルミート法および径路基底関数法よりも精度と計算効率の両面で優れている。

ABSTRACT

In this paper, we propose an efficient numerical scheme for the approximate solution of the time fractional diffusion-wave equation with reaction term based on cubic trigonometric basis functions. The time fractional derivative is approximated by the usual finite difference formulation and the derivative in space is discretized using cubic trigonometric B-spline functions. A stability analysis of the scheme is conducted to confirm that the scheme does not amplify errors. Computational experiments are also performed to further establish the accuracy and validity of the proposed scheme. The results obtained are compared with a finite difference schemes based on the Hermite formula and radial basis functions. It is found that our numerical approach performs superior to the existing methods due to its simple implementation, straight forward interpolation and very less computational cost.

研究の動機と目的

• 時間分数拡散波方程式に反応項を含む高精度な数値スキームの開発を目的とする。
• 解析的解が存在しない分数階微分方程式を解く課題に取り組む。
• 精度を向上させるとともに計算コストを低減することで、既存の有限差分法およびメッシュフリー法を改善すること。
• 特に長期シミュレーションにおいても数値的安定性を保証するため、厳密な解析を実施すること。
• ヘルミート式および径路基底関数に基づく既存スキームと比較して、優れた性能を示すことを示すこと。

提案手法

• コッホ時間分数微分（1 < γ ≤ 2）の導関数は、標準的な2次精度有限差分公式を用いて近似する。
• 空間微分は、C²連続性と最小サポートを有する立方体三角的Bスプライン基底関数を用いて離散化する。
• 得られる方程式系は、各時間ステップで反復的に解き、Bスプラインの局所近似特性を活用する。
• 安定性は、ヴァン・ニューマン型の手法を用いて解析され、スキームの無条件安定性が確認された。
• 解は立方体三角的Bスプライン関数の線形結合として表現されるコロケーションフレームワークで実装される。
• 精度と収束性を評価するために、既知の解析解を有する3つのベンチマーク例を用いた数値的検証が実施された。

実験結果

リサーチクエスチョン

• RQ1 立方体三角的Bスプラインに基づく有限差分スキームは、時間分数拡散波方程式に対して、既存手法よりも高い精度を達成できるか？
• RQ2 提案スキームは無条件安定であるか？ また、ヘルミート法のような条件付き安定スキームと比較して、安定性はどのように異なるか？
• RQ3 分数階数γの異なる値において、計算コストと収束速度の観点から、本手法はどのように性能を発揮するか？
• RQ4 標準的な有限差分法と比較して、三角関数的Bスプラインを用いることで、補間精度と解の滑らかさはどの程度向上するか？
• RQ5 本スキームは、γや境界条件が異なる複数のテストケースにおいても、高精度（例：10⁻¹¹）を維持できるか？

主な発見

• 例1では、提案スキームが最大絶対誤差10⁻¹¹を達成し、ヘリミート法（精度10⁻⁵）を著しく上回った。
• ヴァン・ニューマン安定性解析により、本スキームは無条件安定であることが確認された。これに対して[14]に示されたスキームは条件付き安定であった。
• 例2では、時間tと空間xに対してu(x,t) = t²x(1−x)の解析解と数値解が、γの値に関係なく非常に近い一致を示した。N=60、Δt=0.001の条件下で、絶対誤差は10⁻¹⁰未満に保たれた。
• 例3では、正確な解u(x,t) = t²sinh(x)を正確に捉えており、誤差プロファイルは時間ステップにわたり一貫した収束を示した。
• 例2および例3における数値解と解析解の3次元プロットは、ほぼ完全に一致しており、空間的・時間的精度の高さが確認された。
• 表7～表9における径路基底関数法との比較により、本スキームはより低い絶対誤差と優れた収束特性を示した。