An a posteriori error estimator for model adaptivity in electrocardiology

We introduce an a posteriori modeling error estimator for the effective computation of electric potential propagation in the heart. Starting from the Bidomain problem and an extended formulation of the simplified Monodomain system, we build a hybrid model, called Hybridomain, which is dynamically adapted to be either Bi- or Monodomain ones in different regions of the computational domain according to the error estimator. We show that accurate results can be obtained with the adaptive Hybridomain model with a reduced computational cost compared to the full Bidomain model. We discuss the effectivity of the estimator and the reliability of the results on simulations performed on real human left ventricle geometries retrieved from healthy subjects.

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

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

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

[4]  J. Oden,et al.  Estimation of local modeling error and goal-oriented adaptive modeling of heterogeneous materials: Part II: a computational environment for adaptive modeling of heterogeneous elastic solids , 2001 .

[5]  J. Tinsley Oden,et al.  Estimation of modeling error in computational mechanics , 2002 .

[6]  S. Perotto.  Adaptive modeling for free-surface flows , 2006 .

[7]  J. Tinsley Oden,et al.  Estimation of local modeling error and goal-oriented adaptive modeling of heterogeneous materials: I. Error estimates and adaptive algorithms , 2000 .

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

[9]  D. Farina,et al.  Model-based approach to the localization of infarction , 2007, 2007 Computers in Cardiology.

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

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

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

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

[14]  Raducanu Razvan,et al.  MATHEMATICAL MODELS and METHODS in APPLIED SCIENCES , 2012 .

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

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

[17]  Alexandre Ern,et al.  A Posteriori Control of Modeling Errors and Discretization Errors , 2003, Multiscale Model. Simul..

[18]  Paul T. Bauman,et al.  Computational analysis of modeling error for the coupling of particle and continuum models by the Arlequin method , 2008 .

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

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

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

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

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

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

[25]  Luca F. Pavarino,et al.  Multilevel Schwarz and Multigrid Preconditioners for the Bidomain System , 2008 .

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

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

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

[29]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.

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

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

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

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

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

[35]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .