Semi-implicit Non-conforming Finite-Element Schemes for Cardiac Electrophysiology: A Framework for Mesh-Coarsening Heart Simulations

The field of computational cardiology has steadily progressed toward reliable and accurate simulations of the heart, showing great potential in clinical applications such as the optimization of cardiac interventions and the study of pro-arrhythmic effects of drugs in humans, among others. However, the computational effort demanded by in-silico studies of the heart remains challenging, highlighting the need of novel numerical methods that can improve the efficiency of simulations while targeting an acceptable accuracy. In this work, we propose a semi-implicit non-conforming finite-element scheme (SINCFES) suitable for cardiac electrophysiology simulations. The accuracy and efficiency of the proposed scheme are assessed by means of numerical simulations of the electrical excitation and propagation in regular and biventricular geometries. We show that the SINCFES allows for coarse-mesh simulations that reduce the computation time when compared to fine-mesh models while delivering wavefront shapes and conduction velocities that are more accurate than those predicted by traditional finite-element formulations based on the same coarse mesh, thus improving the accuracy-efficiency trade-off of cardiac simulations.

[1]  Gernot Plank,et al.  Image-based models of cardiac structure with applications in arrhythmia and defibrillation studies. , 2009, Journal of electrocardiology.

[2]  S. Göktepe,et al.  Electromechanics of the heart: a unified approach to the strongly coupled excitation–contraction problem , 2010 .

[3]  Lawrence Mitchell,et al.  Simulating Human Cardiac Electrophysiology on Clinical Time-Scales , 2011, Front. Physio..

[4]  D. Hurtado,et al.  Gradient flows and variational principles for cardiac electrophysiology: Toward efficient and robust numerical simulations of the electrical activity of the heart , 2014 .

[5]  Daniel E Hurtado,et al.  Uncertainty quantification of 2 models of cardiac electromechanics , 2017, International journal for numerical methods in biomedical engineering.

[6]  S. Pezzuto,et al.  Space‐discretization error analysis and stabilization schemes for conduction velocity in cardiac electrophysiology , 2016, International journal for numerical methods in biomedical engineering.

[7]  F. Fenton,et al.  Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation. , 1998, Chaos.

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

[9]  D. D. Streeter,et al.  Engineering Mechanics for Successive States in Canine Left Ventricular Myocardium: II. Fiber Angle and Sarcomere Length , 1973, Circulation research.

[10]  I. Tezaur,et al.  Uncertainty Quantification , 2011, Encyclopedia of Parallel Computing.

[11]  Ellen Kuhl,et al.  Computational modelling of electrocardiograms: repolarisation and T-wave polarity in the human heart , 2014, Computer methods in biomechanics and biomedical engineering.

[12]  Daniel B. Ennis,et al.  Electrophysiology of Heart Failure Using a Rabbit Model: From the Failing Myocyte to Ventricular Fibrillation , 2016, PLoS Comput. Biol..

[13]  Guillaume Houzeaux,et al.  A massively parallel computational electrophysiology model of the heart , 2011 .

[14]  Jiang Yao,et al.  Predicting drug‐induced arrhythmias by multiscale modeling , 2018, International journal for numerical methods in biomedical engineering.

[15]  Natalia A Trayanova,et al.  Personalized Imaging and Modeling Strategies for Arrhythmia Prevention and Therapy. , 2018, Current opinion in biomedical engineering.

[16]  D. Hurtado,et al.  Computational modeling of non-linear diffusion in cardiac electrophysiology: A novel porous-medium approach , 2016 .

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

[18]  E. Wilson,et al.  A non-conforming element for stress analysis , 1976 .

[19]  Robert Michael Kirby,et al.  High-order spectral/hp element discretisation for reaction–diffusion problems on surfaces: Application to cardiac electrophysiology , 2014, J. Comput. Phys..

[20]  Ellen Kuhl,et al.  The importance of mechano-electrical feedback and inertia in cardiac electromechanics. , 2017, Computer methods in applied mechanics and engineering.

[21]  A. Quarteroni,et al.  Finite element and finite volume‐element simulation of pseudo‐ECGs and cardiac alternans , 2015 .

[22]  Alessio Gizzi,et al.  Effects of Pacing Site and Stimulation History on Alternans Dynamics and the Development of Complex Spatiotemporal Patterns in Cardiac Tissue , 2013, Front. Physiol..

[23]  Eric Kerfoot,et al.  Verification of cardiac tissue electrophysiology simulators using an N-version benchmark , 2011, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[24]  Ellen Kuhl,et al.  Generating Purkinje networks in the human heart. , 2016, Journal of biomechanics.

[25]  Daniel E. Hurtado,et al.  Non-conforming finite-element formulation for cardiac electrophysiology: an effective approach to reduce the computation time of heart simulations without compromising accuracy , 2018 .

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

[27]  Yves Coudière,et al.  Very high order finite volume methods for cardiac electrophysiology , 2017, Comput. Math. Appl..

[28]  David Kay,et al.  Efficient simulation of cardiac electrical propagation using high-order finite elements II: Adaptive p-version , 2013, J. Comput. Phys..

[29]  Robert L. Taylor,et al.  FEAP - - A Finite Element Analysis Program , 2011 .

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

[31]  Shankarjee Krishnamoorthi,et al.  Numerical quadrature and operator splitting in finite element methods for cardiac electrophysiology , 2013, International journal for numerical methods in biomedical engineering.

[32]  Jonathan P. Whiteley,et al.  An Efficient Numerical Technique for the Solution of the Monodomain and Bidomain Equations , 2006, IEEE Transactions on Biomedical Engineering.

[33]  Wenjun Ying,et al.  Efficient Fully Implicit Time Integration Methods for Modeling Cardiac Dynamics , 2008, IEEE Transactions on Biomedical Engineering.

[34]  Pablo Lamata,et al.  An accurate, fast and robust method to generate patient-specific cubic Hermite meshes , 2011, Medical Image Anal..

[35]  D. D. Streeter.  Engineering Mechanics for Successive States in Canine Left Ventricular Myocardium: I. CAVITY AND WALL GEOMETRY , 1973, Circulation research.

[36]  Joakim Sundnes,et al.  Computing the electrical activity in the heart , 2006 .

[37]  D. Nordsletten,et al.  A displacement-based finite element formulation for incompressible and nearly-incompressible cardiac mechanics , 2014, Computer methods in applied mechanics and engineering.

[38]  Edward L. Wilson,et al.  Incompatible Displacement Models , 1973 .

[39]  M. Nash,et al.  Electromechanical model of excitable tissue to study reentrant cardiac arrhythmias. , 2004, Progress in biophysics and molecular biology.