×

心脏电生理学的梯度流和变分原理:实现心脏电活动的高效和稳健的数值模拟。 (英语) Zbl 1296.76186号

摘要:心脏电活动的计算机模拟在过去十年中取得了巨大进展。然而,计算方法在医学界的接受程度将在很大程度上取决于其可靠性、效率和鲁棒性。在这项工作中,我们提出了心脏电生理方程的梯度流格式,并针对时间离散的电生理问题提出了一个极大极小变分原理。基于变分分析的结果,我们导出了保证鞍点存在唯一性的时间步长界,进而导出了电生理增量问题的弱解。我们还展示了极大极小问题等价于有效极小原理的条件,该原理适用于Rayleigh-Ritz有限元分析。导出的时间步长边界保证了空间离散化目标函数的严格凸性,从而保证了梯度衰减方法的收敛性。将该理论应用于广泛应用的FitzHugh-Nagumo模型,结果表明该模型符合本文提出的变分框架。该方法的适用性及其对时间积分鲁棒性的影响通过对单细胞、三维楔形和双心室几何结构中电行为的数值模拟得到了证明。我们预计,拟议的框架将为开发健壮高效的电生理模型和模拟打开大门。

MSC公司:

76Z05个 生理流
92C30型 生理学(一般)
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
74升15 生物力学固体力学
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 阿利耶夫,R。;潘菲洛夫,A.V.,心脏兴奋的简单双变量模型,混沌孤子分形,7,3,293-301(1996)
[2] Ambrosio,L。;Gigli,N。;Savaré,G.,度量空间和概率测度空间中的梯度流(2008),Birkhäuser,Verlag AG:Birkháuser·Zbl 1145.35001号
[3] 比勒,B.Y.G.W。;Reuter,H.,心室肌纤维动作电位的重建,《生理学杂志》。,268, 177-210 (1977)
[4] 博伊德,S。;Vandenberghe,L.,《凸优化》,第25卷(2004),剑桥大学出版社·兹比尔1058.90049
[5] Chen,R.,紧黎曼流形上的Neumann特征值估计,Proc。美国数学。《社会学杂志》,108,4961-970(1990)·Zbl 0688.53017号
[7] Courant,R.,解决平衡和振动问题的变分方法,Bull。美国数学。Soc.,49,1,1-23(1943年)·Zbl 0810.65100号
[8] Courtemanche,M。;斯卡格斯,W。;Winfree,A.,《FitzHugh-Nagumo模型中的稳定三维动作电位循环》,Physica D,41,173-182(1990)·Zbl 0711.35068号
[9] Dacorogna,B.,《变分法导论》(2009),帝国理工学院出版社:帝国理工大学出版社伦敦
[10] De Giorgi,E。;Franzoni,T.,《阿提·阿卡德汇聚之路》(Sur un tipo di convergenza variazionale)。纳粹。临北人。Cl.科学。财政部。材料自然。,8 (1975) ·Zbl 0339.49005号
[11] 埃克兰,I。;Témam,R.,凸分析和变分问题(1999),SIAM:费城SIAM·Zbl 0939.49002号
[12] El Sayed,T。;Mota,A。;F.兄弟会。;Ortiz,M.,《软生物组织的变分本构模型》,J.Biomech。,41, 7, 1458-1466 (2008)
[13] Evans,L.C.,偏微分方程(1998),美国数学学会·Zbl 0902.35002号
[14] 芬顿,F.H。;Cherry,E.M.,《心脏细胞模型》,Scholarpedia出版社,第3期,第8期,1868页(2008年)
[15] 芬顿,F.H。;Karma,A.,纤维旋转三维连续心肌的涡旋动力学:纤维不稳定性和纤颤,Chaos,8,1,20-47(1998)·Zbl 1069.92503号
[16] FitzHugh,R.,神经膜理论模型中的冲动和生理状态,生物物理学。J.,1,6,445-466(1961)
[17] Friesecke,G。;Dolzmann,G.,具有非凸能量的拟线性发展方程的隐式时间离散和全局存在性,SIAM J.Math。分析。,28, 2, 363 (1997) ·兹比尔0872.35026
[18] 加维尼,V。;克纳普,J。;巴塔查亚,K。;Ortiz,M.,无轨道密度泛函理论的非周期有限元公式,J.Mech。物理学。固体,55,4669-696(2007)·兹比尔1162.74461
[19] 哥克特佩,S。;Kuhl,E.,《心脏电生理的计算建模:一种新的有限元方法》,国际期刊《数值》。方法工程,79,2,156-178(2009)·Zbl 1171.92310号
[20] (Guccione,J.;Kassab,G.;Ratcliffe,M.,《计算心血管力学:心力衰竭中的建模和应用》(2010),Springer:Springer New York,Dordrecht,Heidelberg,London)
[21] Harrild,D.M。;Henriquez,C.S.,心脏传播的有限体积模型,Ann.Biomed。工程师,25,2,315-334(1997)
[22] 霍奇金。;Huxley,A.,《膜电流的定量描述及其在神经传导和兴奋中的应用》,《生理学杂志》。,117, 500-544 (1952)
[24] Hurtado,D.E。;Ortiz,M.,《表面效应和镍微柱的尺寸依赖性硬化和强化》,J.Mech。物理学。固体,60,8,1432-1446(2012)
[25] Hurtado,D.E。;Ortiz,M.,晶体塑性中几何必要位错的有限元分析,国际J。数值。方法工程,93,1,66-79(2013)·Zbl 1352.74079号
[26] Johnston,P.R.,bidomain方程的有限体积法解及其在心脏缺血建模中的应用,计算。方法生物技术。生物识别。工程,13,2,157-170(2010)
[27] Keener,J。;Sneyd,J.,《数学生理学:细胞生理学》,跨学科应用数学,第1卷(1998年),Springer Verlag:Springer Verlag纽约,NY·Zbl 0913.92009号
[28] Langwallner,B。;奥特纳,C。;Süli,E.,电子密度泛函Galerkin近似的存在性和收敛性结果,数学。模型方法应用。科学。,20, 12, 2237-2265 (2010) ·Zbl 1208.82063号
[29] 线路,G.T。;布伊斯特,M.L。;Grottum,P。;普兰,A.J。;Sundnes,J。;Tveito,A.,心脏电生理正向问题的数学模型和数值方法,计算机。目视检查。科学。,5, 215-239 (2003) ·Zbl 1050.92017年
[30] 罗,C.H。;Rudy,Y.,《心室动作电位的动态模型》。I.离子电流和浓度变化的模拟,Circ。第74、6、1071-1096号决议(1994年)
[31] 莫斯勒,J。;Ortiz,M.,关于变分任意拉格朗日-欧拉(VALE)公式的数值实现,Int.J.Numer。方法工程,67,9,1272-1289(2006)·Zbl 1113.74073号
[32] Nagumo,J。;Arimoto,S。;Yoshizawa,S.,模拟神经轴突的主动脉冲传输线,Proc。IRE,50,10,2061-2070(1962)
[33] 尼德勒,S.A。;Kerfoot,E。;本森,美联社。;伯纳乌,M.O。;伯纳斯,O。;布拉德利,C。;Cherry,E.M。;克莱顿,R。;芬顿,F.H。;Garny,A。;Heidenreich,E。;Land,S。;Maleckar,M。;Pathmanathan,P。;普朗克,G。;罗德里格斯,J.F。;罗伊,我。;Sachse,F.B。;Seemann,G。;O.斯卡夫豪格。;Smith,N.P.,使用N版基准验证心脏组织电生理模拟器,Philos。事务处理。序列号。数学。物理学。工程科学。,369, 4331-4351 (2011)
[34] Noble,D.,适用于Purkinje纤维动作和起搏器电位的Hodgkin-Huxley方程的修正,《生理学杂志》。,160, 317-352 (1962)
[35] 奥尔蒂斯,M。;Stainier,L.,粘塑性本构更新的变分公式,计算。方法应用。机械。工程,171,419-444(1999)·Zbl 0938.74016号
[36] 潘菲洛夫,A.V。;Keener,J.,《心脏解剖模型中的再实验》,混沌孤子分形,5,3-4,681-689(1995)·Zbl 0925.92046号
[37] Pennacchio,M。;萨瓦雷,G。;Colli Franzone,P.,心脏生物电活动的多尺度建模,SIAM J.Math。分析。,37, 4, 1333 (2005) ·Zbl 1113.35019号
[38] 普兰,A.J。;布伊斯特,M.L。;Cheng,L.K.,《心脏电活动的数学建模:从细胞到体表再回来》(2005),世界科学出版社:新加坡世界科学出版社·Zbl 1120.92015年
[39] 拉多维茨基,R。;Ortiz,M.,强非线性动力学问题中的误差估计和自适应网格划分,计算机。方法应用。机械。工程,172,1-4,203-240(1999)·Zbl 0957.74058号
[40] 罗杰斯,J.M。;McCulloch,A.D.,心脏动作电位传播的配置-伽辽金有限元模型,IEEE Trans。生物识别。工程,41,8,743-757(1994)
[41] Sanfelici,S.,心电图退化进化问题的Galerkin近似的收敛性,数值。方法部分差异。Equ.、。,18, 2, 218-240 (2002) ·Zbl 1002.65100号
[42] 夏普,G.H。;Joyner,R.W.,《心脏动作电位的模拟传播》,生物物理学。J.,31,3,403-423(1980)
[43] Suryanarayana,P。;加维尼,V。;布莱森,T。;巴塔查亚,K。;Ortiz,M.,Kohn-Sham密度泛函理论的非周期有限元公式,J.Mech。物理学。固体,58,2,256-280(2010)·Zbl 1193.81006号
[44] 10 Tusscher,K.H.W.J。;Noble,D。;Noble,P.J。;潘菲洛夫,A.V.,《人类心室组织模型》,美国物理杂志。心脏循环。生理学。,286,4,H1573-H1589(2004年)
[45] Trayanova,N.,《整体心脏建模:心脏电生理和电机学的应用》,Circ。研究,108,1,113-128(2011)
[46] Veeser,A。;Verfurth,R.,有限元恒星的庞加莱常数,IMA J.Numer。分析。,32, 1, 30-47 (2011) ·Zbl 1235.65127号
[47] Veneroni,M.,心脏电场宏观双畴模型的反应扩散系统,非线性分析。真实世界应用。,10, 2, 849-868 (2009) ·Zbl 1167.35403号
[48] 维格蒙德,E.J。;Aguel,F。;Trayanova,N.,求解三维双域方程的计算技术,IEEE Trans。生物识别。工程师,49,11,1260-1269(2002)
[49] 温伯格,K。;Ortiz,M.,冲击波对肾脏组织的损伤,计算机。马特。科学。,32, 3-4, 588-593 (2005)
[50] Wong,J。;哥克特佩,S。;Kuhl,E.,《电化学耦合的计算建模:心脏电生理离子模型的新型有限元方法》,计算。方法应用。机械。工程,200,45-46,3139-3158(2011)·Zbl 1230.92013年
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。