Alya Red CCM: HPC-Based Cardiac Computational Modelling

This paper describes Alya Red CCM, a cardiac computational modelling tool for supercomputers. It is based on Alya, a parallel simulation code for multiphysics and multiscale problems, which can deal with all the complexity of biological systems simulations. The final goal is to simulate the pumping action of the heart: the model includes the electrical propagation, the mechanical contraction and relaxation and the blood flow in the heart cavities and main vessels. All sub-problems are treated as fully transient and solved in a staggered fashion. Electrophysiology and mechanical deformation are solved on the same mesh, with no interpolation. Fluid flow is simulated on a moving mesh using an Arbitrary Lagrangian-Eulerian (ALE) strategy, being the mesh deformation computed through an anisotropic Laplacian equation. The parallel strategy is based on an automatic mesh partition using Metis and MPI tasks. When required and in order to take profit of multicore clusters, an additional OpenMP parallelization layer is added. The paper describes the tool and its different parts. Alya’s flexibility allows to easily program a large variety of physiological models for each of the sub-problems, including the mutual coupling. This flexibility, added to the parallel efficiency to solve multiphysics problems in complex geometries render Alya Red CCM a well suited tool for cardiac biomedical research at either industrial or academic environments.

[1]  Arif Masud,et al.  A multiscale stabilized ALE formulation for incompressible flows with moving boundaries , 2010 .

[2]  Guillaume Houzeaux,et al.  A massively parallel fractional step solver for incompressible flows , 2009, J. Comput. Phys..

[3]  Guillaume Houzeaux,et al.  A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling , 2013 .

[4]  Alejandro F. Frangi,et al.  Characterization and Modeling of the Peripheral Cardiac Conduction System , 2013, IEEE Transactions on Medical Imaging.

[5]  Rainald Löhner,et al.  A linelet preconditioner for incompressible flow solvers , 2003 .

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

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

[8]  Gerald D Buckberg,et al.  'The electrical spiral of the heart': its role in the helical continuum. The hypothesis of the anisotropic conducting matrix. , 2006, European journal of cardio-thoracic surgery : official journal of the European Association for Cardio-thoracic Surgery.

[9]  Takumi Washio,et al.  A Multi-Scale Heart Simulation on Massively Parallel Computers , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.

[10]  Guillaume Houzeaux,et al.  Hybrid MPI-OpenMP performance in massively parallel computational fluid dynamics , 2010 .

[11]  Guillaume Houzeaux,et al.  Alya: Towards Exascale for Engineering Simulation Codes , 2014, 1404.4881.

[12]  D. Noble,et al.  A model for human ventricular tissue. , 2004, American journal of physiology. Heart and circulatory physiology.

[13]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[14]  Guillaume Houzeaux,et al.  Deflated preconditioned conjugate gradient solvers for the pressure‐Poisson equation: Extensions and improvements , 2011 .

[15]  Guillaume Houzeaux,et al.  Extension of fractional step techniques for incompressible flows: The preconditioned Orthomin(1) for the pressure Schur complement , 2011 .

[16]  Ruth Arís Sánchez Electromechanical large scale computational models of the ventricular myocardium , 2014 .

[17]  Guillaume Houzeaux,et al.  Coupled electromechanical model of the heart: Parallel finite element formulation , 2012, International journal for numerical methods in biomedical engineering.

[18]  P. Hunter,et al.  A quantitative analysis of cardiac myocyte relaxation: a simulation study. , 2006, Biophysical journal.

[19]  Guillaume Houzeaux,et al.  Parallel uniform mesh multiplication applied to a Navier–Stokes solver , 2013 .

[20]  Gerhard A Holzapfel,et al.  Constitutive modelling of passive myocardium: a structurally based framework for material characterization , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

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

[22]  A. Huerta,et al.  Arbitrary Lagrangian–Eulerian Methods , 2004 .

[23]  Ekkehard Ramm,et al.  A strong coupling partitioned approach for fluid–structure interaction with free surfaces , 2007 .

[24]  G. Houzeaux,et al.  A variational subgrid scale model for transient incompressible flows , 2008 .

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