[論文レビュー] Compact ADI method for solving two-dimensional Riesz space fractional diffusion equation

Sohrab Valizadeh, Alaeddin Malek+1 more|arXiv (Cornell University)|Feb 6, 2018

Fractional Differential Equations Solutions参考文献 33被引用数 3

ひとこと要約

本稿は、2次元リーズ空間分数拡散方程式を解くためのコンactな交替方向隠し（ADI）法を提案する。空間方向で2次精度を達成し、無条件安定性を保証する。コンパクト有限差分スキームとクランク・ニコルソン法、および作用素分割技術を組み合わせることで、高い精度と安定性を実現し、理論的分析と数値実験により検証されている。

ABSTRACT

In this paper, a compact alternating direction implicit (ADI) method has been developed for solving two-dimensional Riesz space fractional diffusion equation. The precision of the discretization method used in spatial directions is twice the order of the corresponding fractional derivatives. It is proved that the proposed method is unconditionally stable via the matrix analysis method and the maximum error in achieving convergence is discussed. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed technique.

研究の動機と目的

• 2次元リーズ空間分数拡散方程式を解くための高次精度かつ無条件安定な数値解法を開発すること。
• リーズ微分を伴う分数拡散問題において、空間方向で2次精度を達成する課題に対処すること。
• コンパクト差分スキームと交替方向隠し（ADI）法を組み合わせることで、計算効率と安定性を向上させること。
• 行列解析を用いて、提案スキームの無条件安定性と収束次数を理論的に分析すること。
• 定量的誤差解析を伴う数値例を通じて、手法の妥当性と適用可能性を示すこと。

提案手法

• 空間方向のリーズ空間分数微分を、コンパクト有限差分スキームで離散化し、空間方向で2次精度を達成する。
• 時間方向にクランク・ニコルソンスキームを適用することで、2次精度の時間積分と無条件安定性を確保する。
• AD I法による作用素分割技術を用いて、2次元問題を逐次1次元問題に分解し、計算複雑度を低減する。
• 三重対角行列 $ A_x, B_x, A_y, B_y $ を用いて、空間離散化作用素を行列形式で表現し、対角化による効率的解法を可能にする。
• 空間作用素の固有値分解を用い、分数微分作用素行列を $ S_x = \frac{k_t}{2}\mathcal{C}_x P \mathbf{D}_x^{\alpha/2} P^{-1} $ および $ T_y = \frac{k_t}{2}\mathcal{C}_y Q \mathbf{D}_y^{\beta/2} Q^{-1} $ と表現する。
• 行列 $ S_x $ と $ T_y $ が対称的かつ可換であることを保証することで、各時間ステップにおける線形方程式系の安定的かつ効率的な解法が可能になる。

実験結果

リサーチクエスチョン

• RQ1 コンパクトAD I法は、2次元リーズ空間分数拡散方程式に対して、空間方向で2次精度を達成できるか？
• RQ2 提示された定式化のもとで、提案されたコンパクトAD Iスキームは無条件安定か？
• RQ3 数値スキームの空間および時間方向の収束次数は何か？
• RQ4 既存手法と比較して、最大誤差と計算効率の観点から、本手法の性能はいかがなものか？
• RQ5 非局所的性質を持つリーズ分数微分を効果的に扱いながら、高精度と安定性を維持できるか？

主な発見

• 提案されたコンパクトAD I法は、コンパクト差分スキームにより切り捨て誤差を $ \mathcal{O}(h^4) $ に低減することで、空間方向で2次精度を達成する。
• 行列解析を用いた理論的分析により、空間作用素の正定値性と可換性に基づき、無条件安定性が証明された。
• 理論的解析により、時間および空間方向で2次収束が確認され、最大誤差が $ \mathcal{O}(k_t^2 + h_x^4 + h_y^4) $ で有界であることが示された。
• 数値実験により、手法の妥当性と適用可能性が確認され、一貫した収束率と低誤差レベルを示した。
• 空間作用素 $ -A_x^{-1}B_x $ および $ -A_y^{-1}B_y $ の固有値は正であり、かつ相異なるため、安定性と対角化可能性が保証された。
• 行列 $ S_x $ と $ T_y $ は対称的かつ可換であるため、各時間ステップにおける線形方程式系が対角化により安定的かつ効率的に解ける。