×

双域和单域模型的高阶算子分裂。 (英语) Zbl 1385.92005年

摘要:双畴和单畴模型是描述心肌组织电活动的多尺度反应扩散模型。由于生成有意义数据所需的系统的大小和特定结构,通常通过应用运算符拆分(OS)方法来找到数值解。一阶和二阶OS方法已成功用于求解双域和单域模型的数值解。众所周知,阶数大于2的操作系统方法需要在每个步骤中进行向后时间集成。因此,可以得出结论,阶数高于2的分裂方法不适用于包含确定性抛物方程的模型,因为必要的向后时间积分会导致不稳定性。在本文中,我们证明了应用于双畴和单畴模型的三阶OS方法确实可以获得稳定的结果。此外,我们还证明,与低阶方法相比,高阶方法在一些示例性问题上的效率提高显著提高了三阶方法的性能。

MSC公司:

92C05型 生物物理学
92C30型 生理学(一般)
65平方米 涉及偏微分方程的初值和初边值问题的线方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] U.M.Ascher、S.J.Ruuth和R.J.Spiteri,{含时偏微分方程的隐式显式Runge-Kutta方法},应用。数字。数学。,25(1997),第151-167页·Zbl 0896.65061号
[2] A.A.Cuellar、C.M.Lloyd、P.F.Nielsen、D.P.Bullivant、D.P.Nickerson和P.J.Honter,《CellML概述》(1.1),生物模型描述语言,《模拟》,79(2003),第740-747页。
[3] W.Auzinger和W.Herfort,{指数分裂格式的李元素的局部误差结构和阶条件},Opuscula Math。,34 (2014). ·兹比尔1333.65055
[4] W.Auzinger、W.Herfort、H.Hofsta¨tter和O.Koch,{分裂方法的序条件设置},《科学计算中的计算机代数学报:第18届国际研讨会》,纽约斯普林格,2016年,第30-42页·Zbl 1453.65115号
[5] W.Auzinger,H.Hofstaötter,D.Ketcheson,and O.Koch,{非线性发展方程自适应积分的实用分裂方法。第一部分:优化方案和方案对的构建},BIT,57(2017),第55-74页·Zbl 1359.65080号
[6] W.Bao,S.Jin,和P.A.Markowich,{关于半经典Schro¨dinger方程的分时谱逼近},J.Compute。物理。,175(2002),第487-524页·Zbl 1006.65112号
[7] S.Blanes、F.Casas、A.Farres、J.Laskar、J.Makazaga和A.Murua,{动力天文学数值积分辛分裂方法的新家族},应用。数字。数学。,68(2013),第58-72页·Zbl 1263.85007号
[8] C.Cantwell等人,{\ it Nektar++:开源光谱/hp元素框架},计算。物理学。社区。,192(2015),第205-219页·兹比尔1380.65465
[9] E.M.Cherry、H.S.Greenside和C.S.Henriquez,《模拟复杂心脏动力学的时空自适应方法》,Phys。修订稿。,84(2000),第1343-1346页。
[10] S.A.Chin,{量子统计计算和辛校正算法},Phys。E版,69(2004),046118。
[11] M.Ethier和Y.Bourgault,{双域模型的半隐式时间离散化方案},SIAM J.Numer。分析。,46(2008),第2443-2468页·Zbl 1182.92009年
[12] I.Faragoí,P.G.Thomsen,and Z.Zlatev,{关于加法分裂过程及其计算机实现},Appl。数学。型号。,32(2008),第1552-1569页·Zbl 1176.65065号
[13] R.FitzHugh,《神经兴奋和传播的数学模型》,摘自《生物工程》,麦格劳-希尔出版社,纽约,1966年,第1-85页。
[14] J.Geiser,《微分方程的分解方法:理论与应用》,Chapman&Hall/CRC Numer。分析。科学。计算。,CRC出版社,佛罗里达州博卡拉顿,2009年·Zbl 1188.65132号
[15] L.Geraldo-Giorda,{生物系统中的非线性动力学},SEMA SIMAI Springer 7,Springer,Cham,2016年。
[16] A.Gerisch和J.G.Verwer,《出租车扩散反应模型的算子分裂和近似因子分解》,应用。数字。数学。,42(2002),第159-176页·Zbl 0998.65102号
[18] R.Glowinski,T.W.Pan,T.I.Hesla,D.D.Joseph,and J.Peíriaux,{\it《运动刚体不可压缩粘性流直接数值模拟的虚拟域方法:颗粒流应用》,J.Compute。物理。,169(2001),第363-426页·Zbl 1047.76097号
[19] S.K.Godunov,{流体动力学方程间断解数值计算的差分方法},Mat.Sb.,89(1959),第271-306页·Zbl 0171.46204号
[20] G.Goldman和T.J.Kaper,{N阶算子分裂方案和不可逆系统},SIAM J.Numer。分析。,33(1996年),第349-367页·兹比尔0849.65070
[21] W.Groébner,{\it Die Lie-Reihen und Ihre Anwendungen},数学。单声道。,德国维森沙芬出版社,柏林,1960年·Zbl 0141.08502号
[22] J.L.Guermond、P.Minev和J.Shen,《不可压缩流投影方法概述》,《计算》。方法应用。机械。工程,195(2006),第6011-6045页·Zbl 1122.76072号
[23] E.Hairer、C.Lubich和G.Wanner,《几何-数值积分:常微分方程的结构保持算法》,Springer Ser。计算。数学。31,Springer,纽约,2006年·Zbl 1094.65125号
[24] L.Hatton和A.Roberts,科学软件的准确性如何。柔和。《工程》,20(1994),第785-797页。
[25] W.Hundsdorfer和J.Verwer,{依赖时间的对流扩散反应方程的数值解},Springer Ser。计算。数学。33,Springer,纽约,2003年·兹比尔1030.65100
[26] S.Jin,P.A.Markowich,C.Zheng,{广义Zakharov系统的数值模拟},J.Compute。物理。,201(2004),第376-395页·Zbl 1079.65101号
[27] J.P.Keener和K.Bogar,《心脏组织中bidomain方程解的数值方法》,Chaos,8(1998),第234-241页·Zbl 1069.92507号
[28] O.Koch,C.Neuhauser和M.Thalhammer,《非线性发展方程时间积分的嵌入指数算子分裂方法》,应用。数字。数学。,63(2013),第14-24页·Zbl 1255.65102号
[29] A.Kriston、G.Inzelt、I.Faragaí和T.Szabo³,{通过实时应用的操作员分裂技术模拟燃料电池的瞬态行为},计算。《化学工程》,34(2010),第339-348页。
[30] R.J.LeVeque,{守恒定律的数值方法},Birkhauser,巴塞尔,1992年·Zbl 0847.65053号
[31] S.Lie和F.Engel,《转型理论》,第一卷,美国社会,普罗维登斯,1970年。
[32] C.H.Luo和Y.Rudy,{心室心脏动作电位的模型。去极化、复极及其相互作用},《循环研究》,68(1991),第1501-1526页。
[33] G.I.Marchuk,《分裂方法》,瑙卡,莫斯科,1988年。
[34] G.I.Marchuk,{分裂和交替方向方法},《数值分析手册》,第一卷,北荷兰,阿姆斯特丹,1990年,第197-462页·Zbl 0875.65049号
[35] M.E.Marsh、S.T.Ziaratgahi和R.J.Spiteri,《Rush-Larsen方法及其推广的成功秘诀》,IEEE Trans。生物识别。工程,59(2012),第2506-2515页。
[36] G.R.Mirams等人,《{it Chaste:计算生理学和生物学的开源(C++)库》,《公共科学图书馆·计算》。《生物学》,9(2013)。
[37] S.A.Niederer等人,{使用n版本基准验证心脏组织电生理模拟器},Philos。事务处理。A、 369(2011),第4331-4351页。
[38] E.Oran和J.Boris,《反应流数值模拟》,剑桥大学出版社,英国剑桥,2005年·Zbl 1092.76002号
[39] P.Pathmanathan、G.R.Mirams、J.Southern和J.P.Whiteley,《选择离子电流积分方法在心脏电生理模拟中的显著影响》,Int.J.Numer。方法生物识别。《工程》,27(2011),第1751-1770页·Zbl 1242.92017年
[40] Z.Qu和A.Garfinkel,{一种求解心脏传导偏微分方程的高级算法},IEEE Trans。生物识别。工程,46(1999),第1166-1168页。
[41] L.F.Richardson,{微分方程物理问题的有限差分近似算术解,及其在砌石坝应力中的应用},Philos。事务处理。A、 210(1911),第307-357页·JFM 42.0873.02号
[42] T.Roy、Y.Bourgault和C.Pierre,《单域模型时间步长方法分析》,2017年·Zbl 1463.65243号
[43] R.D.Ruth,《规范积分技术》,IEEE Trans。核科学。,30(1983年),第2669-2671页。
[44] Q.Sheng,{用指数分裂法求解线性偏微分方程},IMA J.Numer。分析。,9(1989),第199-212页·Zbl 0676.65116号
[45] A.T.Sornborger,{确定性抛物方程的高阶算子分裂方法},国际计算杂志。数学。,84(2007),第887-893页·Zbl 1169.65332号
[46] A.T.Sornborger和E.D.Stewart,《量子计算机模拟的高阶方法》,Phys。A版,60(1999),第765-789页。
[47] J.A.Southern、G.Plank、E.J.Vigmond和J.P.Whiteley,《求解耦合系统可提高双域方程的计算效率》,IEEE Trans。生物识别。工程,56(2009),第2404-2412页。
[48] R.J.Spiteri和R.C.Dean,心脏电生理模型的刚度分析,Ann.Biomed。工程,38(2010),第3592-3604页。
[49] R.J.Spiteri和S.Torabi Ziaratgahi,《重访双域模型的算子分裂》,J.Compute。申请。数学。,296(2016),第550-563页·Zbl 1341.78025号
[50] G.Strang,{关于差分格式的构造和比较},SIAM J.Numer。分析。,5(1968年),第506-517页·Zbl 0184.38503号
[51] J.Sundnes、G.Lines和X.Cai,《计算心脏电活动》,纽约斯普林格,2006年·Zbl 1182.92020年
[52] J.Sundnes,G.T.Lines和A.Tveito,{一种求解耦合到躯干体积导体模型的双畴方程的算子分裂方法},数学。生物科学。,194(2005),第233-248页·Zbl 1063.92018年
[53] 铃木,《分形路径积分的一般理论及其在多体理论和统计物理中的应用》,J.Math。物理。,32(1991),第400-407页·Zbl 0735.47009号
[54] K.H.W.J.Ten Tuscher和A.V.Panfilov,《人类心室组织模型中的Alternans和螺旋破裂》,Amer。《生理学-心脏循环生理学》,291(2006),第H1088-H1100页。
[55] M.Thalhammer,{含时薛定谔方程的高阶指数算子分裂方法},SIAM J.Numer。分析。,46(2008),第2022-2038页·Zbl 1170.65061号
[56] M.Thalhammer、M.Caliari和C.Neuhauser,{高阶时间分裂Hermite和Fourier谱方法},J.Compute。物理。,228(2009),第822-832页·Zbl 1158.65340号
[57] H.Trotter,{\it关于算子半群的乘积},Proc。阿默尔。数学。《社会学杂志》,第10卷(1959年),第545-551页·Zbl 0099.10401号
[58] L.Tung,{描述缺血心肌D-C电位的双域模型},麻省理工学院博士论文,1978年。
[59] R.Tyson、L.G.Stern和R.J.LeVeque,《应用于趋化模型的分步方法》,J.Math。生物学,41(2000),第455-475页·Zbl 1002.92003年
[60] N.N.Yanenko,{分数步法。多变量数学物理问题的求解},Springer,纽约,1971·Zbl 0209.47103号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。