A model-based block-triangular preconditioner for the Bidomain system in electrocardiology

We introduce a preconditioner for the solution of the Bidomain system governing the propagation of action potentials in the myocardial tissue. The Bidomain model is a degenerate parabolic set of nonlinear reaction-diffusion equations. The nonlinear term describes the ion flux at the cellular level. The degenerate nature of the problem results in a severe ill conditioning of its discretization. Our preconditioning strategy is based on a suitable adaptation of the Monodomain model, a simplified version of the Bidomain one, which is by far simpler to solve, nevertheless is unable to capture significant features of the action potential propagation. The Monodomain preconditioner application to a non-symmetric formulation of the Bidomain system results at the algebraic level in a lower block-triangular preconditioner. We prove optimality of the preconditioner with respect to the mesh size, and corroborate our theoretical results with 3D numerical simulations both on idealized and real ventricle geometries.

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

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

[3]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[4]  P. Hunter,et al.  Laminar structure of the heart: a mathematical model. , 1997, The American journal of physiology.

[5]  G. Huiskamp,et al.  Simulation of depolarization in a membrane-equations-based model of the anisotropic ventricle , 1998, IEEE Transactions on Biomedical Engineering.

[6]  A. McCulloch,et al.  A collocation-Galerkin finite element model of cardiac action potential propagation , 1994, IEEE Transactions on Biomedical Engineering.

[7]  C. Henriquez.  Simulating the electrical behavior of cardiac tissue using the bidomain model. , 1993, Critical reviews in biomedical engineering.

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

[9]  S. Scacchi.  Mathematics and Statistics Multilevel Schwarz preconditioners for the Bidomain system and applications to electrocardiology , 2007 .

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

[11]  Andrew J. Pullan,et al.  Mathematically Modelling the Electrical Activity of the Heart: From Cell to Body Surface and Back Again , 2005 .

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

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

[14]  Pavel B. Bochev,et al.  On the Finite Element Solution of the Pure Neumann Problem , 2005, SIAM Rev..

[15]  Marco Veneroni,et al.  Reaction–diffusion systems for the microscopic cellular model of the cardiac electric field , 2006 .

[16]  Isaac Fried.  Bounds on the extremal eigenvalues of the finite element stiffness and mass matrices and their spectral condition number , 1972 .

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

[18]  Russell D Folks,et al.  Onset of left ventricular mechanical contraction as determined by phase analysis of ECG-gated myocardial perfusion SPECT imaging: Development of a diagnostic tool for assessment of cardiac mechanical dyssynchrony , 2005, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

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

[20]  B. Taccardi,et al.  Spread of excitation in 3-D models of the anisotropic cardiac tissue. III. Effects of ventricular geometry and fiber structure on the potential distribution. , 1997, Mathematical biosciences.

[21]  Frank B. Sachse,et al.  Computational Cardiology , 2004, Lecture Notes in Computer Science.

[22]  Luca Gerardo-Giorda,et al.  A Robin?Robin preconditioner for advection?diffusion equations with discontinuous coefficients , 2004 .

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

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

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

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

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

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

[29]  R. Winslow,et al.  Mechanisms of altered excitation-contraction coupling in canine tachycardia-induced heart failure, II: model studies. , 1999, Circulation research.

[30]  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.

[31]  R. FitzHugh.  Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

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

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

[34]  B. Roth,et al.  Action potential propagation in a thick strand of cardiac muscle. , 1991, Circulation research.

[35]  Jeroen J. Bax,et al.  Left ventricular dyssynchrony assessed by two three-dimensional imaging modalities: phase analysis of gated myocardial perfusion SPECT and tri-plane tissue Doppler imaging , 2007, European Journal of Nuclear Medicine and Molecular Imaging.

[36]  Andrew D. McCulloch,et al.  Computational Methods for Cardiac Electrophysiology , 2004 .

[37]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[38]  Giuseppe Savaré,et al.  Computational electrocardiology: mathematical and numerical modeling , 2006 .

[39]  Yoram Rudy,et al.  Computational biology in the study of cardiac ion channels and cell electrophysiology , 2006, Quarterly Reviews of Biophysics.