×

多维分数阶Allen-Cahn方程的分维指数时间差分格式。 (英文) Zbl 1466.65099号

本工作的目的是为求解多维空间分数阶Allen-Cahn(FAC)方程提供有效的数值方法。由于空间离散矩阵的维数分裂,提出了一种新的ETD(指数时间差分)型多维FAC方程积分器。结果表明,所提出的分维ETD-Runge-Kutta(DSETDRK2)格式保留了最大值原理,具有二阶精度。此外,给出了最大范数下的误差估计。该方案导致线性系统和仅与具有Toeplitz结构的1D离散化矩阵相关的矩阵指数。在该方案的实现中,使用了Toeplitz逆和Toeplitz-矩阵指数的快速算法。文中给出了二维和三维数值算例,以说明该方法的有效性。

MSC公司:

65平方米 含偏微分方程初值和初边值问题离散方程的数值解
65号06 含偏微分方程边值问题的有限差分方法
65层10 线性系统的迭代数值方法
2015财年65 矩阵特征值和特征向量的数值计算
65升05 常微分方程初值问题的数值方法
65层60 矩阵指数和相似矩阵函数的数值计算
65岁15岁 涉及PDE的初值和初边值问题的误差界
15个B05 Toeplitz、Cauchy和相关矩阵
35克53 KdV方程(Korteweg-de-Vries方程)
35兰特 分数阶偏微分方程
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Allen,SM公司;Cahn,JW,反相边界运动微观理论及其在反相畴粗化中的应用,金属学报。,27, 1085-1095 (1979) ·doi:10.1016/0001-6160(79)90196-2
[2] Al-Mohy,AH;Higham,NJ,计算矩阵指数的作用,并应用于指数积分器,SIAM J.Sci。计算。,33, 488-511 (2011) ·12346.5028兹比尔 ·doi:10.1137/100788860
[3] 布埃诺·奥罗维亚,A。;凯·D·。;Burrage,K.,分数维空间反应扩散方程的傅里叶谱方法,BIT,54,937-954(2014)·Zbl 1306.65265号 ·doi:10.1007/s10543-014-0484-2
[4] Burrage,K。;黑尔,N。;Kay,D.,分数空间反应扩散方程的高效隐式有限元格式,SIAM J.Sci。计算。,34,A2145-A2172(2012)·Zbl 1253.65146号 ·doi:10.1137/110847007
[5] Chan,R。;Jin,X.,《迭代Toeplitz解算器简介》(2007),费城:SIAM,费城·Zbl 1146.65028号 ·数字对象标识代码:10.1137/1.9780898718850
[6] Chan,R。;Ng,M.,Toeplitz系统的共轭梯度法,SIAM Rev.,38,427-482(1992)·Zbl 0863.65013号 ·doi:10.1137/S0036144594276474
[7] 考克斯,SM;Matthews,PC,刚性系统的指数时间差分,J.Compute。物理。,176, 430-455 (2002) ·Zbl 1005.65069号 ·doi:10.1006/jcph.2002.6995
[8] 杜琪。;朱伟,指数时间差分格式的稳定性分析和应用,J.Compute。数学。,22, 200-209 (2004) ·Zbl 1052.65081号
[9] 杜琪。;朱伟,指数时间差分格式及其轮廓积分修正的分析与应用,BIT,45,307-328(2005)·Zbl 1080.65074号 ·doi:10.1007/s10543-005-7141-8
[10] 杜琪。;Yang,J.,非局部Allen-Cahn方程的渐近兼容傅里叶谱近似,SIAM J.Numer。分析。,54, 1899-1919 (2016) ·Zbl 1342.65198号 ·doi:10.1137/15M1039857
[11] 杜琪。;Ju,L。;李,X。;Qiao,Z.,非局部Allen-Cahn方程的保最大原理指数时间差分格式,SIAM J.Numer。分析。,57, 875-898 (2019) ·Zbl 1419.65018号 ·doi:10.1137/18M118236X
[12] Du,Q.,Ju,L.,Li,X.,Qiao,Z.:一类半线性抛物方程和指数时间差分格式的最大界原理。SIAM版本接受,2020年
[13] X·冯。;Prohl,A.,Allen-Cahn方程的数值分析和平均曲率流的近似,Numer。数学。,94, 33-65 (2003) ·兹比尔1029.65093 ·doi:10.1007/s00211-002-0413-1
[14] X·冯。;宋,H。;Tang,T。;Yang,J.,Allen-Cahn方程隐式-显式方法的非线性稳定性,逆Probl。成像,7679-695(2013)·Zbl 1273.65111号 ·doi:10.3934/ip.2013.7.679
[15] X·冯。;Tang,T。;Yang,J.,相场模型的稳定曲柄-尼科尔森/阿达姆-巴霍斯方案,东亚应用杂志。数学。,3, 59-80 (2013) ·Zbl 1302.65224号 ·doi:10.4208/eajam.20113.220213a
[16] Gaudreault,S。;雨水,G。;Tokman,M.,KIOPS:指数积分器的快速自适应Krylov子空间解算器,J.Compute。物理。,372, 236-255 (2018) ·Zbl 1418.65074号 ·doi:10.1016/j.jcp.2018.06.026
[17] 戈伯格,I。;Olshevsky,V.,矩阵的循环、位移和分解,积分。等于。操作。理论,15730-743(1992)·Zbl 0772.15010号 ·doi:10.1007/BF01200697
[18] 他,D。;潘,K。;Hu,H.,多维分数阶Allen-Cahn方程的空间四阶保最大值原理算子分裂格式,应用。数字。数学。,151, 44-63 (2020) ·Zbl 1434.65117号 ·doi:10.1016/j.apnum.2019.12.018
[19] 新泽西州海姆,《矩阵函数:理论与计算》(2008),费城:SIAM,费城·Zbl 1167.15001号 ·数字对象标识代码:10.1137/1.9780898717778
[20] 新泽西州海姆;Al-Mohy,AH,计算矩阵函数,Acta Numer。,19, 159-208 (2010) ·Zbl 1242.65090号 ·doi:10.1017/S0962492910000036
[21] Hochbruck,M。;Lubich,C.,关于矩阵指数算子的Krylov子空间近似,SIAM J.Numer。分析。,34, 1911-1925 (1997) ·Zbl 0888.65032号 ·doi:10.137/S0036142995280572
[22] Hochbruck,M。;Ostermann,A.,指数积分器,《数值学报》。,1920年至286年(2010年)·Zbl 1242.65109号 ·doi:10.1017/S0962492910000048
[23] 霍恩,R。;Johnson,C.,《矩阵分析专题》(1991),剑桥:剑桥大学出版社,剑桥·Zbl 0729.15001号 ·doi:10.1017/CBO9780511840371
[24] Hou,T。;Tang,T。;Yang,J.,分数维空间Allen-Cahn方程的全离散Crank-Nicolson格式的数值分析,J.Sci。计算。,72, 1214-1231 (2017) ·Zbl 1379.65063号 ·doi:10.1007/s10915-017-0396-9
[25] Jin,X.,块Toeplitz迭代求解器的开发和应用(2002),Dordrecht:Kluwer,Dordracht
[26] 阿拉斯加州卡萨姆;Trefethen,LN,刚性PDE的四阶时间步进,SIAM J.Sci。计算。,26, 1214-1233 (2005) ·Zbl 1077.65105号 ·doi:10.1137/S1064827502410633
[27] Lee,S。;刘,X。;Sun,H-W,带跳跃期权定价的快速指数时间积分方案,Numer。线性代数应用。,19, 87-101 (2012) ·兹比尔1274.91481 ·doi:10.1002/nla.749
[28] Lee,S。;Pang,H.等人。;Sun,H-W,Toeplitz矩阵指数的Shift-in-vert-Arnoldi近似,SIAM J.Sci。计算。,32, 774-792 (2010) ·Zbl 1210.65079号 ·doi:10.1137/090758064
[29] Liao,H-L;Tang,T。;Zhou,T.,关于Allen-Cahn方程的能量稳定、最大值原理保持的二阶变步长BDF格式,SIAM J.Numer。分析。,58, 2294-2314 (2020) ·兹比尔1447.65083 ·doi:10.1137/19M1289157
[30] Minchev,B.V.,Wright,W.M.:一阶半线性问题指数积分器综述。预印本,NTNU特隆赫姆,特隆赫姆挪威,2005
[31] 尼森,J。;Wright,WM,算法919:用于评估指数积分器中出现的函数的Krylov子空间算法,ACM Trans。数学。软质。,38, 22 (2012) ·Zbl 1365.65185号 ·doi:10.1145/2168773.2168781
[32] Pang,H.等人。;Sun,H-W,对称半正定Toeplitz矩阵指数的Shift-in-vert-Lanczos方法,Numer。线性代数应用。,18, 603-614 (2011) ·Zbl 1249.65104号 ·doi:10.1002/nla.747
[33] Pang,H.等人。;Sun,H-W,随机波动率跳扩散模型中定价期权的快速指数时间积分,东亚应用杂志。数学。,4, 53-68 (2014) ·Zbl 1292.91186号 ·doi:10.4208/eajam.280313.061013a
[34] Podlubny,I.,分数微分方程(1999),纽约:学术出版社,纽约·Zbl 0924.34008号
[35] Schmelzer,T。;Trefethen,LN,通过Carathéodory-Fejér近似和轮廓积分评估指数积分器的矩阵函数,Electron。变速器。数字。分析。,29, 1-18 (2007) ·Zbl 1186.65092号
[36] 沈,J。;Tang,T。;Yang,J.,关于广义Allen-Cahn方程的最大原理保持格式,Commun。数学。科学。,14, 1517-1534 (2016) ·Zbl 1361.65059号 ·doi:10.4310/CMS.2016.v14.n6.a3
[37] 沈杰。;Yang,X.,Allen-Cahn和Cahn-Hilliard方程的数值近似,离散。Contin公司。动态。系统。,28, 1669-1691 (2010) ·兹比尔1201.65184 ·doi:10.3934/cds.2010.28.1669
[38] Tang,T。;Yang,J.,Allen-Cahn方程的隐式显式格式保留了最大值原理,J.Compute。数学。,34, 471-481 (2016) ·Zbl 1374.65154号
[39] 田伟。;周,H。;邓,W.,解空间分数阶扩散方程的一类二阶差分逼近,数学。计算。,841703-1727(2015年)·Zbl 1318.65058号 ·doi:10.1090/S0025-5718-2015-02917-2
[40] 杨,X.,Allen-Cahn方程稳定半隐式方法的误差分析,离散Contin。动态。系统。序列号。B、 11057-1070(2009)·Zbl 1201.65170号
[41] 张杰。;Du,Q.,尖锐界面极限下Allen-Cahn方程离散近似的数值研究,SIAM J.Sci。计算。,31, 3042-3063 (2009) ·Zbl 1198.82045号 ·doi:10.1137/080738398
[42] 张,L。;Sun,H。;Pang,H.,分数阶扩散方程的指数求积快速数值解,J.Compute。物理。,299, 130-143 (2015) ·Zbl 1352.65304号 ·doi:10.1016/j.jcp.2015.07.001
[43] 张,L。;Sun,H.,用指数Runge-Kutta方法求解多维Riesz分数阶非线性反应扩散方程,J.Appl。数学。计算。,62, 449-472 (2020) ·Zbl 1475.65054号 ·doi:10.1007/s12190-019-01291-w
[44] 张,L。;张,Q。;Sun,H.,二维非线性分数阶复Ginzburg-Landau方程的指数Runge-Kutta方法,J.Sci。计算。,83, 59 (2020) ·Zbl 1442.65340号 ·doi:10.1007/s10915-020-01240-x
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。