[論文レビュー] A symmetric fractional-order reduction method for direct nonuniform approximations of semilinear diffusion-wave equations

Pin Lyu, Seakweng Vong|arXiv (Cornell University)|Jan 24, 2021

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

ひとこと要約

本稿では、一般非一様時刻メッシュ上での半線形時間分数拡散波方程式に対する高次精度数値スキームを構築するため、対称的分数階還元（SFOR）法を提案する。元の問題を分数階 $α/2$ の導関数を有する連立系に還元することで、著者らは線形化L1スキームおよびアリカノフスキームを開発し、弱いメッシュ条件および正則性仮定のもとで最適な $ϴ(\tau^2)$ 時間収束を達成した。この結果は、適応的時間ステッピングおよび数値実験によって検証された。

ABSTRACT

We introduce a symmetric fractional-order reduction (SFOR) method to construct numerical algorithms on general nonuniform temporal meshes for semilinear fractional diffusion-wave equations. By using the novel order reduction method, the governing problem is transformed to an equivalent coupled system, where the explicit orders of time-fractional derivatives involved are all $α/2$ $(1

研究の動機と目的

• 非一様時刻メッシュ上での初期特異性が弱い半線形拡散波方程式に対して最適収束を達成する課題に対処すること。
• 初期特異性による解の低正則性に対しても高次精度を保つ強固な数値スキームを開発すること。
• L1およびアリカノフスキームに基づく適応的時間ステッピング戦略を設計し、計算効率と精度を向上させること。
• 最小限の正則性およびメッシュ制約のもとで、$H^2$エネルギー法を用いて厳密な収束解析を確立すること。
• 支配方程式をバランスの取れた分数階導関数を有する連立系に簡略化する対称的分数階還元技術を提供すること。

提案手法

• 元の $\alpha$ 階時間分数拡散波方程式を $\alpha/2$ 階導関数を有する等価な連立系に変換するため、対称的分数階還元（SFOR）法を提案する。
• 一般非一様時刻メッシュ上に線形化L1スキームおよびアリカノフスキームを適用して、還元系を離散化する。
• 弱いメッシュグレーディングおよび正則性仮定のもとで、最適な $\mathcal{O}(\tau^2)$ 収束率を導出するために $H^2$ エネルギー法を用いる。
• 動的時間ステップを調整可能な、（高速化された）線形化L1およびアリカノフスキームに基づく適応的時間ステッピング戦略を導入する。
• 両スキームのためのグローバル一貫性誤差推定を用いて、時間離散化誤差を定義および解析する。
• 離散畳み込みカーネル $P_{n-j}^{(n)}$ および $A_{j-k}^{(j)}$ を含む和不等式を用いて誤差バインドを確立する。

実験結果

リサーチクエスチョン

• RQ1 非一様メッシュ上での半線形拡散波方程式に対して、精度および安定性を向上させる対称的分数階還元法を開発することは可能か？
• RQ2 弱い正則性およびメッシュ制約のもとで、線形化L1およびアリカノフスキームを適用した場合に達成可能な最適収束率は何か？
• RQ3 計算コストを低減しつつ高精度を維持するため、適応的時間ステッピング戦略をどのように設計できるか？
• RQ4 時間誤差推定値が、メッシュグレーディングパラメータ $\gamma$ および解の正則性パラメータ $\sigma_1, \sigma_2$ にどの程度依存するか？
• RQ5 非一様で非一様グレーディングなメッシュに対し、$H^2$ エネルギー法を効果的に適用して、分数自己拡散問題における最適収束結果を導出できるか？

主な発見

• 線形化L1スキームは、条件 $\gamma\sigma_1 \geq 2$ のもとで $\mathcal{O}(\tau^2)$ 収束率を達成し、誤差バインド $\sum_{j=1}^{n}P_{n-j}^{(n)}\|\Delta_h({\cal T}_f)^{n-\theta}\| \leq C\tau^{\min\{2\gamma\sigma_1,2\}}$ を満たす。
• アリカノフスキームは、$\gamma\sigma_1 \geq 3 - \beta$ のもとで $\mathcal{O}(\tau^2)$ 収束率を達成し、誤差バインド $\sum_{j=1}^{n}P_{n-j}^{(n)}\|\Delta_h({\cal T}_u)^{j-\theta}\| \leq C\tau^{\min\{3-\beta,\gamma\sigma_1\}}$ を満たす。
• 空間誤差は $\|{\cal S}^n\| \leq C_u h^2$ で有界であり、2次空間収束を確認した。
• L1スキームの場合、正則性仮定 $|g''(t)| \leq C_g(1 + t^{\sigma-2})$ のもとで、時間誤差は $C\tau^{\min\{2-\beta,\gamma\sigma_2\}}$ で有界である。
• アリカノフスキームの場合、仮定 $|g^{\prime\prime\prime}(t)| \leq C_g(1 + t^{\sigma-3})$ のもとで、誤差は $C\tau^{\min\{3-\beta,\gamma\sigma_1\}}$ で有界であり、追加的に $\tau_1^{\sigma}$ および $\tau_k^3$ を含む項が存在する。
• L1およびアリカノフスキームに基づく適応的時間ステッピング戦略は、数値例によって確認されたように、精度と効率を両立させた。