[論文レビュー] Efficient multistep methods for tempered fractional calculus: Algorithms and Simulations

Ling Guo, Fanhai Zeng+3 more|arXiv (Cornell University)|Dec 3, 2018

Fractional Differential Equations Solutions被引用数 14

ひとこと要約

本稿では、生成関数から導かれる畳み込み・クワドラチャ・ウェイトを用いて、トゥーメンテッド分数階微積分問題を解くための、高速でメモリ効率の良い2段階法—Fast Method I および Fast Method II—を提案する。Fast Method II は、均一な実数直線トラペゾイド近似に基づくもので、計算コスト $O(Q)$、メモリコスト $O(Qn_T)$ であり、分数階反応拡散モデル（ブリュッセルァー系を含む）の長時間シミュレーションにより検証された。Fast Method I よりも優れた精度、効率、コーディングの簡潔さを達成している。

ABSTRACT

In this work, we extend the fractional linear multistep methods in [C. Lubich, SIAM J. Math. Anal., 17 (1986), pp.704--719] to the tempered fractional integral and derivative operators in the sense that the tempered fractional derivative operator is interpreted in terms of the Hadamard finite-part integral. We develop two fast methods, Fast Method I and Fast Method II, with linear complexity to calculate the discrete convolution for the approximation of the (tempered) fractional operator. Fast Method I is based on a local approximation for the contour integral that represents the convolution weight. Fast Method II is based on a globally uniform approximation of the trapezoidal rule for the integral on the real line. Both methods are efficient, but numerical experimentation reveals that Fast Method II outperforms Fast Method I in terms of accuracy, efficiency, and coding simplicity. The memory requirement and computational cost of Fast Method II are $O(Q)$ and $O(Qn_T)$, respectively, where $n_T$ is the number of the final time steps and $Q$ is the number of quadrature points used in the trapezoidal rule. The effectiveness of the fast methods is verified through a series of numerical examples for long-time integration, including a numerical study of a fractional reaction-diffusion model.

研究の動機と目的

• 異常拡散過程に現れるトゥーメンテッド分数階微分方程式を解くための、効率的かつメモリ効率の良い数値的手法の開発。
• ハダマール有限部積分の解釈を用いて、ルビッチの分数階線形多段法をトゥーメンテッド積分および微分作用素へと拡張すること。
• 離散畳み込みの計算コストを $O(n_T^2)$ から $O(Qn_T)$ に、メモリを $O(n_T)$ から $O(Q)$ に削減する高速アルゴリズムの設計。
• 2つの新規高速法—Fast Method I（積分路に基づく）および Fast Method II（実数直線トラペゾイドに基づく）—の、精度、効率、実装の複雑さという観点での性能比較。
• 分数階反応拡散モデルの長時間シミュレーションを通じた手法の検証、特に分数階ブリュッセルァー系および活性化子・抑制子系を含む。

提案手法

• トゥーメンテッド分数階微分作用素をハダマール有限部積分として解釈することで、ルビッチの分数階線形多段法（FLMMs）を直接的にトゥーメンテッド作用素へと拡張可能となる。
• 畳み込み重み $\omega^{(\alpha,\sigma)}_n$ は生成関数から導かれ、積分として表現される：Fast Method I は複素平面における積分路を用い、Fast Method II は実数直線における積分を用いる。
• Fast Method I は、各 $n$ に対して局所的積分路四角法（例：タロプ法や双曲線法）を用いて積分路積分を近似する。複素数演算を要し、不安定な係数を有する常微分方程式を解く必要がある。
• Fast Method II は、実数直線上で一様なトラペゾイド則近似を採用し、安定な再帰関係を可能とし、実数演算のみを用いる。
• 四角則による指数関数の和近似を活用することで、時間ステップ数に対して線形の複雑度を達成し、$Q$ 個の四角則点が精度とコストを決定する。
• 数値的シミュレーションでは、さまざまな初期条件を用いた分数階ブリュッセルァー系およびギエラー＝マインハルトモデルを用い、長時間安定性と精度を検証した。

実験結果

リサーチクエスチョン

• RQ1 ハダマール有限部積分の解釈を用いて、ルビッチの分数階線形多段法をトゥーメンテッド分数階微積分に効果的に拡張できるか？
• RQ2 トゥーメンテッド分数階作用素のための離散畳み込みを、線形複雑度と低メモリで効率的に計算する方法は何か？
• RQ3 積分路に基づく（Fast Method I）と実数直線トラペゾイドに基づく（Fast Method II）アプローチの間で、精度、安定性、実装の観点から相違点は何か？
• RQ4 Fast Method II は、計算コスト、丸め誤差、コーディングの簡潔さという観点で、Fast Method I よりも優れた性能を発揮できるか？
• RQ5 提案された高速法は、ブリュッセルァー系や活性化子・抑制子系を含む分数階反応拡散モデルの長時間シミュレーションにおいて、どのように性能を発揮するか？

主な発見

• Fast Method II は $O(Q)$ のメモリと $O(Qn_T)$ の計算コストを達成し、直接計算の $O(n_T^2)$ のコストを顕著に削減している。
• 分数階ブリュッセルァー系モデルにおける高速解と直接解との最大差は $10^{-9}$ 未満であり、256個の四角則点を用いることで高い精度を達成していることが示された。
• 長時間シミュレーションにおいて、高速法は約 4000 秒（約1時間）で実行されたが、直接法は約 87,000 秒（約1日2時間）を要した。これにより、20倍の高速化が確認された。
• Fast Method II は、安定な実数演算と一様近似のおかげで、精度、効率、実装の簡潔さにおいて Fast Method I を上回っている。
• Fast Method II では、常微分方程式の解法において係数の正の実部が生じる不安定性を避けるために、安定な再帰関係が用いられている。これは Fast Method I に影響を与える可能性がある。
• 数値的結果から、Fast Method II が分数階偏微分方程式の長時間積分において、強固で正確であることが確認された。特に、分数階ブリュッセルァー系およびギエラー＝マインハルトモデルにおいて有効である。