Numerical simulation of electrocardiograms for full cardiac cycles in healthy and pathological conditions

This work is dedicated to the simulation of full cycles of the electrical activity of the heart and the corresponding body surface potential. The model is based on a realistic torso and heart anatomy, including ventricles and atria. One of the specificities of our approach is to model the atria as a surface, which is the kind of data typically provided by medical imaging for thin volumes. The bidomain equations are considered in their usual formulation in the ventricles, and in a surface formulation on the atria. Two ionic models are used: the Courtemanche-Ramirez-Nattel model on the atria and the 'minimal model for human ventricular action potentials' by Bueno-Orovio, Cherry, and Fenton in the ventricles. The heart is weakly coupled to the torso by a Robin boundary condition based on a resistor-capacitor transmission condition. Various electrocardiograms (ECGs) are simulated in healthy and pathological conditions (left and right bundle branch blocks, Bachmann's bundle block, and Wolff-Parkinson-White syndrome). To assess the numerical ECGs, we use several qualitative and quantitative criteria found in the medical literature. Our simulator can also be used to generate the signals measured by a vest of electrodes. This capability is illustrated at the end of the article. Copyright © 2015 John Wiley & Sons, Ltd.

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

[2]  Casey J Kohen,et al.  7. Electrocardiogram Interpretation , 2014 .

[3]  Daniel B. Ennis,et al.  Simulation Methods and Validation Criteria for Modeling Cardiac Ventricular Electrophysiology , 2014, PloS one.

[4]  Dominique Chapelle,et al.  A SURFACE-BASED ELECTROPHYSIOLOGY MODEL RELYING ON ASYMPTOTIC ANALYSIS AND MOTIVATED BY CARDIAC ATRIA MODELING , 2013 .

[5]  Bin He,et al.  Noninvasive reconstruction of three-dimensional ventricular activation sequence from the inverse solution of distributed equivalent current density , 2006, IEEE Transactions on Medical Imaging.

[6]  M. Courtemanche,et al.  Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model. , 1998, The American journal of physiology.

[7]  A. Bayés‐Genís,et al.  Third-degree inter-atrial block and supraventricular tachyarrhythmias. , 1999, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[8]  Kanu Chatterjee,et al.  Systolic and diastolic heart failure: differences and similarities. , 2007, Journal of cardiac failure.

[9]  Pascal Frey,et al.  YAMS A fully Automatic Adaptive Isotropic Surface Remeshing Procedure , 2001 .

[10]  Kawal S. Rhode,et al.  Personalization of Atrial Anatomy and Electrophysiology as a Basis for Clinical Modeling of Radio-Frequency Ablation of Atrial Fibrillation , 2013, IEEE Transactions on Medical Imaging.

[11]  M Boulakia,et al.  Reduced-order modeling for cardiac electrophysiology. Application to parameter identification. , 2012, International journal for numerical methods in biomedical engineering.

[12]  Robert Plonsey,et al.  Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields , 1995 .

[13]  Nigel H. Lovell,et al.  A Simplified 3D Model of Whole Heart Electrical Activity and 12-Lead ECG Generation , 2013, Comput. Math. Methods Medicine.

[14]  A. Moorman,et al.  Atrioventricular junctional tissue. Discrepancy between histological and electrophysiological characteristics. , 1996, Circulation.

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

[16]  M. Mansour,et al.  Left Atrial Wall Thickness Variability Measured by CT Scans 
in Patients Undergoing Pulmonary Vein Isolation , 2011, Journal of cardiovascular electrophysiology.

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

[18]  Daming Wei,et al.  Comparative simulation of excitation and body surface electrocardiogram with isotropic and anisotropic computer heart models , 1995, IEEE Transactions on Biomedical Engineering.

[19]  Cristian Lorenz,et al.  Modeling Atrial Fiber Orientation in Patient-Specific Geometries: A Semi-automatic Rule-Based Approach , 2011, FIMH.

[20]  G. Huiskamp,et al.  A new method for myocardial activation imaging , 1997, IEEE Transactions on Biomedical Engineering.

[21]  Dd. Streeter,et al.  Gross morphology and fiber geometry of the heart , 1979 .

[22]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[23]  Mark Potse,et al.  Cardiac anisotropy in boundary-element models for the electrocardiogram , 2009, Medical & Biological Engineering & Computing.

[24]  F. Fenton,et al.  Minimal model for human ventricular action potentials in tissue. , 2008, Journal of theoretical biology.

[25]  M. Lysaker,et al.  On the possibility for computing the transmembrane potential in the heart with a one shot method: an inverse problem. , 2007, Mathematical biosciences.

[26]  F. Liu,et al.  Two hybrid regularization frameworks for solving the electrocardiography inverse problem. , 2008, Physics in medicine and biology.

[27]  Leslie Tung,et al.  A bi-domain model for describing ischemic myocardial d-c potentials , 1978 .

[28]  Olaf Dössel,et al.  Influence of ${I_{Ks}}$ Heterogeneities on the Genesis of the T-wave: A Computational Evaluation , 2012, IEEE Transactions on Biomedical Engineering.

[29]  D. Schaeffer,et al.  A two-current model for the dynamics of cardiac membrane , 2003, Bulletin of mathematical biology.

[30]  Stefan Neubauer,et al.  Normal human left and right ventricular and left atrial dimensions using steady state free precession magnetic resonance imaging. , 2005, Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance.

[31]  F. Prinzen,et al.  Patient-specific modelling of cardiac electrophysiology in heart-failure patients , 2014, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[32]  Robert Michael Kirby,et al.  Inverse electrocardiographic source localization of ischemia: An optimization framework and finite element solution , 2013, J. Comput. Phys..

[33]  Miguel A. Fernández,et al.  A Coupled System of PDEs and ODEs Arising in Electrocardiograms Modeling , 2010 .

[34]  Jean-Frédéric Gerbeau,et al.  MagnetoHemoDynamics in the aorta and electrocardiograms , 2012, Physics in medicine and biology.

[35]  Miguel A. Fernández,et al.  Mathematical Modeling of Electrocardiograms: A Numerical Study , 2010, Annals of Biomedical Engineering.

[36]  Daniel B Ennis,et al.  Changes in Mitral Annular Geometry and Dynamics With &bgr;-Blockade in Patients With Degenerative Mitral Valve Disease , 2010, Circulation. Cardiovascular imaging.

[37]  O. Dossel,et al.  Computer based modeling of the congenital long-QT 2 syndrome in the Visible Man torso: From genes to ECG , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

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

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

[40]  A. J. Pullan,et al.  Mathematical models and numerical methods for the forward problem in cardiac electrophysiology , 2002 .

[41]  P. Hunter,et al.  Computational Mechanics of the Heart , 2000 .

[42]  M. Potse,et al.  ECG simulations with realistic human membrane, heart, and torso models , 2003, Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE Cat. No.03CH37439).

[43]  C. Tei,et al.  Size and Motion of the Mitral Valve Annulus in Man: II. Abnormalities in Mitral Valve Prolapse , 1982, Circulation.

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

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

[46]  W. Brady,et al.  Electrocardiography in the patient with the Wolff-Parkinson-White syndrome: diagnostic and initial therapeutic issues. , 1999, The American journal of emergency medicine.

[47]  Dominique Chapelle,et al.  Surface-Based Electrophysiology Modeling and Assessment of Physiological Simulations in Atria , 2013, FIMH.

[48]  Ling Xia,et al.  Simulation of biatrial conduction via different pathways during sinus rhythm with a detailed human atrial model , 2012, Journal of Zhejiang University SCIENCE B.

[49]  C. Henriquez,et al.  A computer model of normal conduction in the human atria. , 2000, Circulation research.

[50]  D. Sánchez-Quintana,et al.  The importance of atrial structure and fibers , 2009, Clinical anatomy.

[51]  P. Hunter,et al.  Computational mechanics of the heart : From tissue structure to ventricular function , 2000 .

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

[53]  Miguel A. Fernández,et al.  Numerical Simulation of the Electromechanical Activity of the Heart , 2009, FIMH.

[54]  B. Khandheria,et al.  Three methods for evaluation of left atrial volume. , 2007, European journal of echocardiography : the journal of the Working Group on Echocardiography of the European Society of Cardiology.

[55]  Iwona Cygankiewicz,et al.  Interatrial blocks. A separate entity from left atrial enlargement: a consensus report. , 2012, Journal of electrocardiology.

[56]  Robert H. Anderson,et al.  Atrial structure and fibres: morphologic bases of atrial conduction. , 2002, Cardiovascular research.

[57]  W. Stewart,et al.  Reference values for normal adult transesophageal echocardiographic measurements. , 1995, Journal of the American Society of Echocardiography : official publication of the American Society of Echocardiography.

[58]  R Hoekema,et al.  On selecting a body surface mapping procedure. , 1999, Journal of Electrocardiology.

[59]  Pyotr G. Platonov,et al.  Clinical predictors of termination and clinical outcome of catheter ablation for persistent atrial fibrillation. , 2009, Journal of the American College of Cardiology.

[60]  Gan-XinYan,et al.  Cellular Basis for the Normal T Wave and the Electrocardiographic Manifestations of the Long-QT Syndrome , 1998 .

[61]  C. Otto,et al.  Changes in Aortic Annulus Diameter During the Cardiac Cycle and its Effect on Predicting Aortic Valve Prosthesis Size , 1994 .

[62]  Gernot Plank,et al.  Bidomain ECG Simulations Using an Augmented Monodomain Model for the Cardiac Source , 2011, IEEE Transactions on Biomedical Engineering.

[63]  Cesare Corrado,et al.  Identification of weakly coupled multiphysics problems. Application to the inverse problem of electrocardiography , 2015, J. Comput. Phys..

[64]  Richard A. Johnson,et al.  Applied Multivariate Statistical Analysis , 1983 .