A fully discrete GL-ADI scheme for 2D time-fractional reaction-subdiffusion equation

IF 3.5 2区数学 Q1 MATHEMATICS, APPLIED

Applied Mathematics and Computation Pub Date : 2024-10-24 DOI:10.1016/j.amc.2024.129147

Yubing Jiang , Hu Chen , Chaobao Huang , Jian Wang

{"title":"A fully discrete GL-ADI scheme for 2D time-fractional reaction-subdiffusion equation","authors":"Yubing Jiang , Hu Chen , Chaobao Huang , Jian Wang","doi":"10.1016/j.amc.2024.129147","DOIUrl":null,"url":null,"abstract":"<div><div>Alternating direction implicit (ADI) difference method for solving a 2D reaction-subdiffusion equation whose solution behaves a weak singularity at <span><math><mi>t</mi><mo>=</mo><mn>0</mn></math></span> is studied in this paper. A Grünwald-Letnikov (GL) approximation is used for the discretization of Caputo fractional derivative (of order <em>α</em>, with <span><math><mn>0</mn><mo><</mo><mi>α</mi><mo><</mo><mn>1</mn></math></span>) on a uniform mesh. Stability and convergence of the fully discrete ADI scheme are rigorously established. With the help of a discrete fractional Gronwall inequality, we get the sharp error estimate. The stability in <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> norm and the convergence of the GL-ADI scheme are strictly proved, where the convergent order is <span><math><mi>O</mi><mo>(</mo><mi>τ</mi><msubsup><mrow><mi>t</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>α</mi><mo>−</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msup><mrow><mi>τ</mi></mrow><mrow><mn>2</mn><mi>α</mi></mrow></msup><mo>+</mo><msubsup><mrow><mi>h</mi></mrow><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>h</mi></mrow><mrow><mn>2</mn></mrow><mrow><mn>2</mn></mrow></msubsup><mo>)</mo></math></span>. Numerical experiments are given to verify the theoretical analysis.</div></div>","PeriodicalId":55496,"journal":{"name":"Applied Mathematics and Computation","volume":"488 ","pages":"Article 129147"},"PeriodicalIF":3.5000,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Applied Mathematics and Computation","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0096300324006088","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}

引用次数: 0

Abstract

Alternating direction implicit (ADI) difference method for solving a 2D reaction-subdiffusion equation whose solution behaves a weak singularity at

t = 0

is studied in this paper. A Grünwald-Letnikov (GL) approximation is used for the discretization of Caputo fractional derivative (of order α, with

0 < α < 1

) on a uniform mesh. Stability and convergence of the fully discrete ADI scheme are rigorously established. With the help of a discrete fractional Gronwall inequality, we get the sharp error estimate. The stability in

L^{2}

norm and the convergence of the GL-ADI scheme are strictly proved, where the convergent order is

O (τ t_{s}^{α - 1} + τ^{2 α} + h_{1}^{2} + h_{2}^{2})

. Numerical experiments are given to verify the theoretical analysis.

查看原文本刊更多论文

二维时间分数反应-次扩散方程的完全离散 GL-ADI 方案

本文研究了交替方向隐式（ADI）差分法求解二维反应-次扩散方程，该方程的解在 t=0 处表现为弱奇点。在均匀网格上对 Caputo 分数导数（阶数 α，0<α<1）的离散化采用了 Grünwald-Letnikov (GL) 近似方法。严格确定了完全离散 ADI 方案的稳定性和收敛性。在离散分式 Gronwall 不等式的帮助下，我们得到了尖锐的误差估计。严格证明了 GL-ADI 方案在 L2 规范下的稳定性和收敛性，其中收敛阶数为 O(τtsα-1+τ2α+h12+h22)。数值实验验证了理论分析。

本文章由计算机程序翻译，如有差异，请以英文原文为准。

求助全文

约1分钟内获得全文求助全文

来源期刊

Applied Mathematics and Computation 数学-应用数学

CiteScore

7.90

自引率

10.00%

发文量

755

审稿时长

36 days

期刊介绍： Applied Mathematics and Computation addresses work at the interface between applied mathematics, numerical computation, and applications of systems – oriented ideas to the physical, biological, social, and behavioral sciences, and emphasizes papers of a computational nature focusing on new algorithms, their analysis and numerical results. In addition to presenting research papers, Applied Mathematics and Computation publishes review articles and single–topics issues.