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

IF 4.3 3区材料科学 Q1 ENGINEERING, ELECTRICAL & ELECTRONIC

ACS Applied Electronic Materials 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":3,"journal":{"name":"ACS Applied Electronic Materials","volume":null,"pages":null},"PeriodicalIF":4.3000,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"ACS Applied Electronic Materials","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0096300324006088","RegionNum":3,"RegionCategory":"材料科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"ENGINEERING, ELECTRICAL & ELECTRONIC","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分钟内获得全文求助全文

来源期刊

ACS Applied Electronic Materials Multiple-

CiteScore

7.20

自引率

4.30%

发文量

567