Comparing the bidomain and monodomain models in electro-cardiology through convergence analysis

The monodomain and bidomain models are widely used in electro cardiology to simulate spreading of excitation potential waves in the myocardium. The bidomain model is quite popular for its physiological foundation and relevance whereas the monodomain model simply is a heuristic approximation of the previous one, lacking this physiological foundation but providing computational facilities. The purpose of the present article is to numerically compare these two models using a method of (numerical) convergence analysis. This method enables to reach two different objectives of first importance in biomedical engineering. Firstly it provides the discrepancy between the models at the continuous level (and not between the discretised equations only) by getting rid of the discretisation errors. Secondly, it allows to estimate the discretisation error so providing necessary grid resolution in order to run accurate enough simulation. The comparison is held in terms of activation times, a quantity of major physiological importance. Two test cases are considered, both including enhanced cell membrane kinetic description, tissue anisotropy and realistic macroscopic tissue parameters. The first test case is based on an academic (a unit square) geometry. The second one involves a 2d cut of a segmented human heart. It has been built from 3D medical data of the two ventricles and incorporates anisotropy occurring from muscular fibre rotation around the ventricles. Two conclusions are drawn from this study. In terms of activation time relative error, the discrepancy between the two models is quite small: of order 1\% or even below. Moreover this error is smaller than the discretisation error resulting from commonly used mesh size in biomedical engineering.

[1]  Zhiqiang Cai,et al.  The finite volume element method for diffusion equations on general triangulations , 1991 .

[2]  L. Clerc.  Directional differences of impulse spread in trabecular muscle from mammalian heart. , 1976, The Journal of physiology.

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

[4]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[5]  Yves Coudière,et al.  2D/3D Discrete Duality Finite Volume (DDFV) scheme for anisotropic- heterogeneous elliptic equations, application to the electrocardiogram simulation. , 2008 .

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

[7]  R C Barr,et al.  Electrophysiological effects of remodeling cardiac gap junctions and cell size: experimental and model studies of normal cardiac growth. , 2000, Circulation research.

[8]  Luca Gerardo-Giorda,et al.  OPTIMIZED SCHWARZ COUPLING OF BIDOMAIN AND MONODOMAIN MODELS IN ELECTROCARDIOLOGY , 2011 .

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

[10]  Kenneth H. Karlsen,et al.  Convergence of discrete duality finite volume schemes for the cardiac bidomain model , 2010, Networks Heterog. Media.

[11]  D. Durrer,et al.  Total Excitation of the Isolated Human Heart , 1970, Circulation.

[12]  Marc Ethier,et al.  Semi-Implicit Time-Discretization Schemes for the Bidomain Model , 2008, SIAM J. Numer. Anal..

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

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

[15]  L. Ambrosio,et al.  On the asymptotic behaviour of anisotropic energies arising in the cardiac bidomain model , 2000 .

[16]  Olivier Rousseau,et al.  Geometrical modeling of the heart , 2010 .

[17]  R C Barr,et al.  Extracellular discontinuities in cardiac muscle: evidence for capillary effects on the action potential foot. , 1998, Circulation research.

[18]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

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

[20]  Andrew J. Pullan,et al.  The effect of torso impedance on epicardial and body surface potentials: a modeling study , 2003, IEEE Transactions on Biomedical Engineering.

[21]  P. Savard,et al.  Extracellular Measurement of Anisotropic Bidomain Myocardial Conductivities. I. Theoretical Analysis , 2001, Annals of Biomedical Engineering.

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

[23]  Nicolas P Smith,et al.  Altered T Wave Dynamics in a Contracting Cardiac Model , 2003, Journal of cardiovascular electrophysiology.

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

[25]  Marie-Claude Trudel,et al.  Simulation of QRST integral maps with a membrane-based computer heart model employing parallel processing , 2004, IEEE Transactions on Biomedical Engineering.

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

[27]  W. Krassowska,et al.  Homogenization of syncytial tissues. , 1993, Critical reviews in biomedical engineering.

[28]  Xing Cai,et al.  On the Computational Complexity of the Bidomain and the Monodomain Models of Electrophysiology , 2006, Annals of Biomedical Engineering.