[論文レビュー] A Second-Order Crank-Nicolson Method for Subdiffusion

Bangti Jin, Buyang Li+1 more|arXiv (Cornell University)|Jul 23, 2016

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

ひとこと要約

本稿では、$\alpha \in (0,1)$ の順序をもつキャプートー微分を含むサブディュージョン方程式に対して、2次精度のクランク・ニコルソン型時間離散化を提案する。後退オイラー畳み込みクワドラチャと $\theta = \alpha/2$ の $\theta$-法を組み合わせ、初期補正を適用することにより、空間的にガラーキン有限要素法を用いる場合、滑らかでないデータに対しても時間方向で2次収束を達成する。

ABSTRACT

In this work, we analyze a Crank-Nicolson type time stepping scheme for the subdiffusion equation, which involves a Caputo fractional derivative of order $\alpha\in (0,1)$ in time. It combines the backward Euler convolution quadrature with a $ heta$-type method, with the parameter $ heta$ dependent on the fractional order $\alpha$ by $ heta=\alpha/2$, and it naturally generalizes the classical Crank-Nicolson method. We develop essential initial corrections for the fractional Crank-Nicolson scheme, and together with the Galerkin finite element method in space, obtain a fully discrete scheme. A complete error analysis of the fully discrete scheme is provided, and a second-order accuracy in time is established for both smooth and nonsmooth problem data. Extensive numerical experiments are provided to illustrate its accuracy, efficiency and robustness, and further, the comparative study indicates that the scheme is competitive with existing schemes.

研究の動機と目的

• サブディュージョン方程式にキャプートー分数マージン（順序 $\alpha \in (0,1)$）を含む場合に、時間方向で2次精度を達成する時間離散化スキームの開発。
• 古典的クランク・ニコルソン法を分数マージンサブディュージョンに一般化するため、$\theta = \alpha/2$ に依存するパラメータ $\theta$ を導入する。
• 滑らかでない初期データに対して収束次数が低下する問題を解消するため、重要な初期補正を導入する。
• 提案された時間スキームを空間方向のガラーキン有限要素法と組み合わせ、完全離散スキームを構築する。
• 広範な数値実験を通じて、スキームの頑健性と効率性を示す包括的な誤差解析を提供する。

提案手法

• 時間離散化は、$\theta = \alpha/2$ の $\theta$-法に基づくバックワード・オイラー畳み込みクワドラチャを用いたクランク・ニコルソン型スキームであり、古典的手法を分数マージンに一般化する。
• 初期補正は導出され、滑らかでない問題データに対して2次収束を回復させるために適用され、そうでない場合に収束次数が低下するのを防ぐ。
• 空間離散化にはガラーキン有限要素法が用いられ、最適な安定性と収束性を有する完全離散スキームが得られる。
• 誤差解析は $L^2$ 範囲およびエネルギー範囲の両方で実施され、滑らかでない初期データに対しても時間方向で2次収束が確立される。
• 広範な数値実験によりスキームが検証され、既存手法との比較を通じて、その正確性、効率性、頑健性が示される。

実験結果

リサーチクエスチョン

• RQ1 キャプートー分数マージンを含むサブディュージョンに対して、時間方向で2次精度を達成するクランク・ニコルソン型スキームを構築可能か？
• RQ2 滑らかでない初期データに対して、分数マージンサブディュージョンスキームの2次収束を回復させるために、初期補正をどのように設計できるか？
• RQ3 提案されたスキームは、空間方向にガラーキン有限要素法を組み合わせた場合でも、2次収束を維持するか？
• RQ4 既存のサブディュージョン数値スキームと比較して、提案スキームの正確性と効率性はどの程度か？
• RQ5 分数マージンの順序 $\alpha$ に応じて $\theta$ パラメータを調整することで、古典的クランク・ニコルソン法を分数マージンサブディュージョン設定に一般化可能か？

主な発見

• 提案スキームは、滑らかでない初期データに対しても時間方向で2次収束を達成し、分数マージンサブディュージョンで一般的に見られる次数低下を克服する。
• 初期補正は、特に滑らかでないデータに対して2次精度を維持するために不可欠であり、解析で厳密に導出されている。
• 時間スキームとガラーキン有限要素法を組み合わせた完全離散スキームは、無条件安定かつ最適次数の収束性を有する。
• 数値実験により理論的誤差推定が確認され、さまざまなテストケースにおいてスキームの頑健性が示された。
• 比較的検証により、本スキームは正確性と計算効率の面で既存手法と同等か、それ以上の性能を示した。
• $\theta = \alpha/2$ のパラメータ選択は、古典的クランク・ニコルソン法を分数マージンサブディュージョン設定に効果的に一般化する。