CardioModel - New Software for Cardiac Electrophysiology Simulation

The rise of supercomputing technologies during the last decade has enabled significant progress towards the invention of a personal biologically relevant computer model of a human heart. In this paper we present a new code for numerical simulation of cardiac electrophysiology on supercomputers. Having constructed a personal segmented tetrahedral grid of the human heart based on a tomogram, we solve the bidomain equations of cardiac electrophysiology using the finite element method thus achieving the goal of modeling of the action potential propagation in heart. Flexible object-oriented architecture of the software allows us to expand its capabilities by using relevant cell models, preconditioners and numerical methods for solving SLAEs. The results of numerical modeling of heart under normal conditions as well as a number of simulated pathologies are in a good agreement with theoretical expectations. The software achieves at least 75% scaling efficiency on the 120 ranks on the Lobachevsky supercomputer.

[1]  Victor P. Gergel,et al.  High Performance Computing in Biomedical Applications , 2013, ICCS.

[2]  A. Tveito,et al.  Numerical solution of the bidomain equations , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

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

[5]  Alexander G. Fletcher,et al.  Chaste: An Open Source C++ Library for Computational Physiology and Biology , 2013, PLoS Comput. Biol..

[7]  Yuri V. Vassilevski,et al.  Personalized Anatomical Meshing of the Human Body with Applications , 2015 .

[8]  Peter J. Hunter,et al.  The CellML Model Repository , 2008, Bioinform..

[9]  N. Trayanova.  Whole-heart modeling: applications to cardiac electrophysiology and electromechanics. , 2011, Circulation research.

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

[11]  Yuri V. Vassilevski,et al.  Parallel software platform INMOST: a framework for numerical modeling , 2015, Supercomput. Front. Innov..

[12]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

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

[14]  Otto H. Schmitt,et al.  Biological Information Processing Using the Concept of Interpenetrating Domains , 1969 .

[15]  Viatcheslav Gurev,et al.  Towards real-time simulation of cardiac electrophysiology in a human heart at high resolution , 2013, Computer methods in biomechanics and biomedical engineering.

[16]  J. C. Cajas,et al.  Alya Red CCM: HPC-Based Cardiac Computational Modelling , 2015 .

[17]  Gernot Plank,et al.  Development of an anatomically detailed MRI-derived rabbit ventricular model and assessment of its impact on simulations of electrophysiological function , 2009, American journal of physiology. Heart and circulatory physiology.

[18]  Vipin Kumar,et al.  Parallel Multilevel series k-Way Partitioning Scheme for Irregular Graphs , 1999, SIAM Rev..

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

[20]  Nicholas Wilson,et al.  Chaste: A case study of parallelisation of an open source finite-element solver with applications to computational cardiac electrophysiology simulation , 2014, Int. J. High Perform. Comput. Appl..

[21]  Mateo Valero,et al.  Alya: Multiphysics engineering simulation toward exascale , 2016, J. Comput. Sci..

[22]  Konstantin Ushenin,et al.  Human Heart Simulation Software for Parallel Computing Systems , 2015 .

[23]  Iosif Meyerov,et al.  Dynamic Parallelization Strategies for Multifrontal Sparse Cholesky Factorization , 2015, PaCT.

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

[25]  P. Hunter,et al.  Computational physiology and the physiome project , 2004, Experimental physiology.

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

[27]  Vipin Kumar,et al.  Parallel Multilevel k-way Partitioning Scheme for Irregular Graphs , 1996, Proceedings of the 1996 ACM/IEEE Conference on Supercomputing.

[28]  Katherine C. Wu,et al.  Arrhythmia risk stratification of patients after myocardial infarction using personalized heart models , 2016, Nature Communications.