[論文レビュー] A fourth-order compact solver for fractional-in-time fourth-order diffusion equations

Jialing Zhong, Hong-lin Liao+2 more|arXiv (Cornell University)|Jul 3, 2019

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

ひとこと要約

本稿では、1次Dirichlet境界条件を満たす時間分数階4次サブディフュージョン方程式に対して、4次精度のコンパクト有限差分スキームを提案する。問題を2階系に再定式化し、空間近似に平均化されたコンパクト作用素を用いることで、4次空間精度を達成するとともに、Caputo微分の初期特異性に対処するため、階層的非一様L1公式を用いる。離散分数階Grönwall不等式と誤差畳み込み構造を用いた厳密な安定性および収束解析により、最適な時間収束率 $ O(\tau^{\min\{\gamma\sigma,2-\alpha\}}) $ が確認された。

ABSTRACT

A fourth-order compact scheme is proposed for a fourth-order subdiffusion equation with the first Dirichlet boundary conditions. The fourth-order problem is firstly reduced into a couple of spatially second-order system and we use an averaged operator to construct a fourth-order spatial approximation. This averaged operator is compact since it involves only two grid points for the derivative boundary conditions. The L1 formula on irregular mesh is considered for the Caputo fractional derivative, so we can resolve the initial singularity of solution by putting more grid points near the initial time. The stability and convergence are established by using three theoretical tools: a complementary discrete convolution kernel, a discrete fractional Gronwall inequality and an error convolution structure. Some numerical experiments are reported to demonstrate the accuracy and efficiency of our method.

研究の動機と目的

• 時間分数階4次サブディフュージョン方程式と1次Dirichlet境界条件を対象とする高次精度数値解法の開発。
• Caputo分数階微分に起因する解の初期特異性に対処するため、階層的非一様時間メッシュを用いる。
• 平均化作用素に基づくコンパクト有限差分スキームを用いて、4次空間精度を達成する。
• 離散分数階Grönwall不等式や誤差畳み込み構造といった理論的道具を用いて、スキームの厳密な安定性および収束性を確立する。
• さまざまな分数階数およびメッシュパラメータを用いた広範な数値実験を通じて理論的予測の妥当性を検証する。

提案手法

• 補助変数 $ v = \partial_x^2 u $ を用いて、4次サブディフュージョン方程式を2階系の連立系に再定式化する。
• 境界条件に2グリッドポイントのみを含むコンパクトスタencilを維持しつつ、4次空間精度を達成するための平均化コンパクト作用素を構築する。
• 初期特異性に対処するため、階層的非一様メッシュ $ t_k = (k/N)^\gamma T $（$ \gamma \geq 1 $）上でのCaputo微分に非一様L1公式を適用する。
• 補完的離散畳み込みカーネル、離散分数階Grönwall不等式、誤差畳み込み構造の3つの主要な道具を用いて、安定性および収束性を証明する。
• 空間では一様メッシュ、時間では階層的メッシュを用いた実装を行い、収束率は $ \text{Order}(M) \approx \log_2(e(M,N)/e(2M,N)) $ および $ \text{Order}(N) \approx \log_2(e(M,N)/e(M,2N)) $ で計算する。
• さまざまな状況下で $ \alpha $, $ \sigma $, $ \gamma $, メッシュサイズを変化させた数値実験を実施し、精度および効率性を検証する。

実験結果

リサーチクエスチョン

• RQ1 時間分数階4次サブディフュージョン方程式と1次Dirichlet境界条件に対して、4次精度コンパクトスキームを構築可能か？
• RQ2 非一様時間メッシュを用いて、解の初期特異性を効果的に解消できるか？
• RQ3 スキームの最適時間収束率は何か？ また、階層パラメータ $ \gamma $ および解の正則性パラメータ $ \sigma $ にどのように依存するか？
• RQ4 離散分数階Grönwall不等式および誤差畳み込み構造を用いて、スキームの安定性および収束性を厳密に証明できるか？
• RQ5 さまざまな分数階数およびメッシュパラメータ下でも、4次空間精度および最適時間収束性を達成できるか？

主な発見

• 数値実験により、すべてのテストケースで空間収束率が4.00に近く、4次空間精度が確認された。
• 時間収束率は理論的予測 $ O(\tau^{\min\{\gamma\sigma,2-\alpha\}}) $ と一致し、$ \gamma > (2-\alpha)/\sigma $ の場合に最適な $ O(\tau^{2-\alpha}) $ が観測された。
• 一様メッシュ（$ \gamma = 1 $）の場合、時間収束率は $ O(\tau^{\sigma}) $ となり、理論的解析および表3の観測結果と整合的であった。
• 高正則性解（$ \sigma = 1.9 $）に対しても、スキームは安定かつ収束性を維持し、$ \alpha = 0.9 $ の場合に時間収束率 $ \approx 1.10 $ を達成した。
• 表4〜6の数値結果により、$ L^\infty $ 誤差推定 $ \|u - U^n\|_{\infty} = O(\tau^{\min\{\gamma\sigma,2-\alpha\}}) $ が鋭く、観測された収束率が理論的上限と一致した。
• 階層的メッシュによる処理により初期層の挙動が効果的に解消され、強い初期特異性下でも一貫した収束率が得られた。