A fully implicit finite element method for bidomain models of cardiac electrophysiology

This work introduces a novel, unconditionally stable and fully coupled finite element method for the bidomain system of equations of cardiac electrophysiology. The transmembrane potential and the extracellular potential are treated as independent variables. To this end, the respective reaction-diffusion equations are recast into weak forms via a conventional isoparametric Galerkin approach. The resultant nonlinear set of residual equations is consistently linearised. The method results in a symmetric set of equations, which reduces the computational time significantly compared to the conventional solution algorithms. The proposed method is inherently modular and can be combined with phenomenological or ionic models across the cell membrane. The efficiency of the method and the comparison of its computational cost with respect to the simplified monodomain models are demonstrated through representative numerical examples.

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

[2]  K Skouibine,et al.  A numerically efficient model for simulation of defibrillation in an active bidomain sheet of myocardium. , 2000, Mathematical biosciences.

[3]  Lawrence W. Solomon,et al.  Yale University School of Medicine Heart Book , 1992, The Yale Journal of Biology and Medicine.

[4]  B. Roth,et al.  Approximate analytical solutions of the Bidomain equations for electrical stimulation of cardiac tissue with curving fibers. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  Barry L. Zaret,et al.  Yale University School of Medicine Heart Book , 1992 .

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

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

[8]  S. Göktepe,et al.  Atrial and ventricular fibrillation: computational simulation of spiral waves in cardiac tissue , 2010 .

[9]  D. Geselowitz,et al.  Simulation Studies of the Electrocardiogram: I. The Normal Heart , 1978, Circulation research.

[10]  S. Göktepe,et al.  Computational modeling of electrocardiograms: A finite element approach toward cardiac excitation , 2010 .

[11]  Alan Garny,et al.  A numerical guide to the solution of the bi-domain equations of cardiac electrophysiology. , 2010, Progress in biophysics and molecular biology.

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

[13]  C. Henriquez,et al.  Cardiac propagation simulation. , 1992, Critical reviews in biomedical engineering.

[14]  R. Aliev,et al.  A simple two-variable model of cardiac excitation , 1996 .

[15]  S. Göktepe,et al.  Computational modeling of cardiac electrophysiology: A novel finite element approach , 2009 .

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

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

[18]  Natalia A. Trayanova,et al.  Concepts of ventricular defibrillation , 2001, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

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

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

[21]  A. Tveito,et al.  An operator splitting method for solving the bidomain equations coupled to a volume conductor model for the torso. , 2005, Mathematical biosciences.

[22]  G Plank,et al.  Computational tools for modeling electrical activity in cardiac tissue. , 2003, Journal of electrocardiology.

[23]  J. Keener,et al.  A numerical method for the solution of the bidomain equations in cardiac tissue. , 1998, Chaos.

[24]  Peter R Johnston,et al.  A finite volume method solution for the bidomain equations and their application to modelling cardiac ischaemia , 2010, Computer methods in biomechanics and biomedical engineering.