OPTIMIZED SCHWARZ COUPLING OF BIDOMAIN AND MONODOMAIN MODELS IN ELECTROCARDIOLOGY

The Bidomain model is nowadays one of the most accurate mathematical descriptions of the action potential propagation in the heart. However, its numerical approximation is in general fairly expensive as a consequence of the mathematical features of this system. For this reason, a simplification of this model, called Monodomain problem is quite often adopted in order to reduce computational costs. Reliability of this model is however questionable, in particular in the presence of applied currents and in the regions where the upstroke or the late recovery of the action potential is occurring. In this paper we investigate a domain decomposition approach for this problem, where the entire computational domain is suitably split and the two models are solved in different subdomains. Since the mathematical features of the two problems are rather different, the heterogeneous coupling is non trivial. Here we investigate appropriate interface matching conditions for the coupling on non overlapping domains. Moreover, we pursue an Optimized Schwarz approach for the numerical solution of the heterogeneous problem. Convergence of the iterative method is analyzed by means of a Fourier analysis. We investigate the parameters to be selected in the matching radiation-type conditions to accelerate the convergence. Numerical results both in two and three dimensions illustrate the effectiveness of the coupling strategy.

[1]  Zi-Cai Li,et al.  Schwarz Alternating Method , 1998 .

[2]  Alessandro Veneziani,et al.  An a posteriori error estimator for model adaptivity in electrocardiology , 2011 .

[4]  Alfio Quarteroni,et al.  Domain Decomposition Methods for Partial Differential Equations , 1999 .

[5]  J. Nenonen,et al.  Activation Dynamics in Anisotropic Cardiac Tissue via Decoupling , 2004, Annals of Biomedical Engineering.

[6]  Kenneth H. Karlsen,et al.  Analysis of a class of degenerate reaction-diffusion systems and the bidomain model of cardiac tissue , 2006, Networks Heterog. Media.

[7]  J. Clark,et al.  Mathematical model of an adult human atrial cell: the role of K+ currents in repolarization. , 1998, Circulation research.

[8]  Natalia Trayanova,et al.  Defibrillation of the heart: insights into mechanisms from modelling studies , 2006, Experimental physiology.

[9]  Giuseppe Savaré,et al.  Degenerate Evolution Systems Modeling the Cardiac Electric Field at Micro- and Macroscopic Level , 2002 .

[10]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[11]  P. Hunter,et al.  Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. , 1995, The American journal of physiology.

[12]  B. Roth,et al.  Electrical stimulation of cardiac tissue by a bipolar electrode in a conductive bath , 1998, IEEE Transactions on Biomedical Engineering.

[13]  Dd. Streeter,et al.  Gross morphology and fiber geometry of the heart , 1979 .

[14]  Marco Veneroni,et al.  Reaction–diffusion systems for the macroscopic bidomain model of the cardiac electric field , 2009 .

[15]  Natalia A. Trayanova,et al.  Computational techniques for solving the bidomain equations in three dimensions , 2002, IEEE Transactions on Biomedical Engineering.

[16]  Simone Scacchi,et al.  A hybrid multilevel Schwarz method for the bidomain model , 2008 .

[17]  B. Roth.  A comparison of two boundary conditions used with the bidomain model of cardiac tissue , 2006, Annals of Biomedical Engineering.

[18]  P. C. Franzone,et al.  A PARALLEL SOLVER FOR REACTION-DIFFUSION SYSTEMS IN COMPUTATIONAL ELECTROCARDIOLOGY , 2004 .

[19]  Martin J. Gander,et al.  Optimized Schwarz Methods for Maxwell's Equations , 2006, SIAM J. Sci. Comput..

[20]  V. Simoncini,et al.  Efficient algebraic solution of reaction-diffusion systems for the cardiac excitation process , 2002 .

[21]  Martin J. Gander,et al.  Optimized Schwarz Methods , 2006, SIAM J. Numer. Anal..

[22]  Rodrigo Weber dos Santos,et al.  Parallel multigrid preconditioner for the cardiac bidomain model , 2004, IEEE Transactions on Biomedical Engineering.

[23]  H Zhang,et al.  Models of cardiac tissue electrophysiology: progress, challenges and open questions. , 2011, Progress in biophysics and molecular biology.

[24]  G Plank,et al.  Solvers for the cardiac bidomain equations. , 2008, Progress in biophysics and molecular biology.

[25]  Alfio Quarteroni,et al.  Complex Systems in Biomedicine , 2006 .

[26]  Mark Potse,et al.  A Comparison of Monodomain and Bidomain Reaction-Diffusion Models for Action Potential Propagation in the Human Heart , 2006, IEEE Transactions on Biomedical Engineering.

[27]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[28]  Martin J. Gander,et al.  Optimized Schwarz Methods without Overlap for the Helmholtz Equation , 2002, SIAM J. Sci. Comput..

[29]  Luca F. Pavarino,et al.  Multilevel Additive Schwarz Preconditioners for the Bidomain Reaction-Diffusion System , 2008, SIAM J. Sci. Comput..

[30]  Aslak Tveito,et al.  Optimal monodomain approximations of the bidomain equations , 2007, Appl. Math. Comput..

[31]  Luca Gerardo-Giorda,et al.  A model-based block-triangular preconditioner for the Bidomain system in electrocardiology , 2009, J. Comput. Phys..

[32]  J. Keener,et al.  Direct activation and defibrillation of cardiac tissue. , 1996, Journal of theoretical biology.

[33]  Editors , 1986, Brain Research Bulletin.

[34]  Luca Gerardo-Giorda,et al.  New Nonoverlapping Domain Decomposition Methods for the Harmonic Maxwell System , 2006, SIAM J. Sci. Comput..

[35]  Jens Lang,et al.  Konrad-zuse-zentrum F ¨ Ur Informationstechnik Berlin Adaptivity in Space and Time for Reaction-diffusion Systems in Electrocardiology Adaptivity in Space and Time for Reaction-diffusion Systems in Electrocardiology , 2022 .

[36]  A V Panfilov,et al.  A guide to modelling cardiac electrical activity in anatomically detailed ventricles. , 2008, Progress in biophysics and molecular biology.

[37]  C. Henriquez,et al.  Anisotropy, Fiber Curvature, and Bath Loading Effects on Activation in Thin and Thick Cardiac Tissue Preparations: , 1996, Journal of cardiovascular electrophysiology.

[38]  Robert F Gilmour,et al.  Ionic mechanism of electrical alternans. , 2002, American journal of physiology. Heart and circulatory physiology.

[39]  Simona Sanfelici.  Convergence of the Galerkin approximation of a degenerate evolution problem in electrocardiology , 2002 .

[40]  Y. Bourgault,et al.  Existence and uniqueness of the solution for the bidomain model used in cardiac electrophysiology , 2009 .

[41]  James P. Keener,et al.  Mathematical physiology , 1998 .

[42]  S. Nepomnyaschikh,et al.  Domain Decomposition Methods , 2007 .

[43]  C. Luo,et al.  A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction. , 1991, Circulation research.

[44]  B. Taccardi,et al.  Simulating patterns of excitation, repolarization and action potential duration with cardiac Bidomain and Monodomain models. , 2005, Mathematical biosciences.

[45]  A. Veneziani,et al.  AN EFFICIENT GENERALIZATION OF THE RUSH-LARSEN METHOD FOR SOLVING ELECTRO-PHYSIOLOGY MEMBRANE EQUATIONS ∗ , 2009 .