MATHICSE Technical Report : Integrated Heart – Coupling multiscale and multiphysics models for the simulation of the cardiac function

Mathematical modelling of the human heart and its function can expand our understanding of various cardiac diseases, which remain the most common cause of death in the developed world. Like other physiological systems, the heart can be understood as a complex multiscale system involving interacting phenomena at the molecular, cellular, tissue, and organ levels. This article addresses the numerical modelling of many aspects of heart function, including the interaction of the cardiac electrophysiology system with contractile muscle tissue, the sub-cellular activation-contraction mechanisms, as well as the haemodynamics inside the heart chambers. Resolution of each of these sub-systems requires separate mathematical analysis and specially developed numerical algorithms, which we review in detail. By using specific sub-systems as examples, we also look at systemic stability, and explain for example how physiological concepts such as microscopic force generation in cardiac muscle cells, translate to coupled systems of differential equations, and how their stability properties influence the choice of numerical coupling algorithms. Several numerical examples illustrate three fundamental challenges of developing multiphysics and multiscale numerical models for simulating heart function, namely: (i) the correct upscaling from single-cell models to the entire cardiac muscle, (ii) the proper coupling of electrophysiology and tissue mechanics to simulate electromechanical feedback, and (iii) the stable simulation of ventricular haemodynamics during rapid valve opening and closure.

Integrated Heart-Coupling multiscale and multiphysics models for the simulation of the cardiac function

Toni Mikael Lassila, Alfio Quarteroni, Simone Rossi, Ricardo Esteban Ruiz Baier

Mathematical modeling of the human heart and its function can expand our understanding of various cardiac diseases, which remain the most common cause of death in the developed world. Like other physiological systems, the heart can be understood as a complex multiscale system involving interacting phenomena at the molecular, cellular, tissue, and organ levels. This article addresses the numerical modeling of many aspects of heart function, including the interaction of the cardiac electrophysiology system with contractile muscle tissue, the sub-cellular activation contraction mechanisms, as well as the hemodynamics inside the heart chambers. Resolution of each of these sub-systems requires separate mathematical analysis and specially developed numerical algorithms, which we review in detail. By using specific sub-systems as examples, we also look at systemic stability, and explain for example how physiological concepts such as microscopic force generation in cardiac muscle cells, translate to coupled systems of differential equations, and how their stability properties influence the choice of numerical coupling algorithms. Several numerical examples illustrate three fundamental challenges of developing multiphysics and multiscale numerical models for simulating heart function, namely: (i) the correct upscaling from single-cell models to the entire cardiac muscle, (ii) the proper coupling of electrophysiology and tissue mechanics to simulate electromechanical feedback, and (iii) the stable simulation of ventricular hemodynamics during rapid valve opening and closure. (C) 2016 Elsevier B.V. All rights reserved.

Elsevier Science Sa2017

Phenomenological and physiological approaches in the study of cardiac alternans: Theoretical background, mathematical modeling and numerical simulations

Marie Dupraz

This project presents the theoretical background in electrophysiology that is used as a basis in the development of numerical methods for the simulation of heart’s electrical activity. We discuss the mathematical models currently used in cardiac electrophysiology research to simulate the propagation of an action potential wave, and then focus on the description of the ionic currents. This study allows us to anticipate the numerical model sensitivity with regard to the discretization, with the aim of recovering the physiological behaviors of the tissue (i.e. conduction velocity, electrical restitution, etc). Thereafter, we introduce the pseudo-ECG signal reconstruction method from the action potential map. This leads us to examine the approximation of such a signal using the finite element basis coming from the discretization of the monodomain equations. We discuss then a particular cardiac rhythm disorder: action potential alternans. In a second step, we perform an extended simulation study in order to validate the previous numerical methods. We carefully check the convergence of the conduction velocity with respect to the mesh size, and we bring out the advantage of phenomenological approach in terms of electrical restitutions and computational cost. We reproduce typical rhythm disorders and show that the pseudo-ECG signal is able to detect them. Finally, we discuss the influence of anisotropy in cardiac wave propagation and alternans development, characterizing the effect of fibers direction and pacing site in the resulting spatio-temporal distribution.

2013

Parallel Algorithms for the Solution of Large-Scale Fluid-Structure Interaction Problems in Hemodynamics

Davide Forti

This thesis addresses the development and implementation of efficient and parallel algorithms for the numerical simulation of Fluid-Structure Interaction (FSI) problems in hemodynamics. Indeed, hemodynamic conditions in large arteries are significantly affected by the interaction of the pulsatile blood flow with the arterial wall. The simulation of fluid-structure interaction problems requires the approximation of a coupled system of Partial Differential Equations (PDEs) and the set up of efficient numerical solution strategies. Blood is modeled as an incompressible Newtonian fluid whose dynamics is governed by the Navier-Stokes equations. Different constituive models are used to describe the mechanical response of the arterial wall; specifically, we rely on hyperelastic isotropic and anistotropic material laws. The finite element method is used for the space discretization of both the fluid and structure problems. In particular, for the Navier-Stokes equations we consider a semi-discrete formulation based on the Variational Multiscale (VMS) method. Among a wide range of possible solution strategies for the FSI problem, here we focus on strongly coupled monolithic approaches wherein the nonlinearities are treated in a fully implicit mode. To cope with the high computational complexity of the three dimensional FSI problem, a parallel solution framework is often mandatory. To this end, we develop a new block parallel preconditioner for the coupled linearized FSI system obtained after space and time discretization. The proposed preconditioner, named FaCSI, exploits the factorized form of the FSI Jacobian matrix, the use of static condensation to formally eliminate the interface degrees of freedom of the fluid equations, and the use of a SIMPLE preconditioner for unsteady Navier-Stokes equations. In FSI problems, the different resolution requirements in the fluid and structure physical domains, as well as the presence of complex interface geometries make the use of matching fluid and structure meshes problematic. In such situations, it is much simpler to deal with discretizations that are nonconforming at the interface, provided however that the matching conditions at the interface are properly fulfilled. In this thesis we develop a novel interpolation-based method, named INTERNODES, for numerically solving partial differential equations by Galerkin methods on computational domains that are split into two (or several) subdomains featuring nonconforming interfaces. By this we mean that either a priori independent grids and/or local polynomial degrees are used to discretize each subdomain. INTERNODES can be regarded as an alternative to the mortar element method: it combines the accuracy of the latter with the easiness of implementation in a numerical code. The aforementioned techniques have been applied for the numerical simulation of large-scale fluid-structure interaction problems in the context of biomechanics. The parallel algorithms developed showed scalability up to thousands of cores utilized on high performance computing machines.

EPFL2016

Parallel Algorithms for the Solution of Large-Scale Fluid-Structure Interaction Problems in Hemodynamics

Davide Forti

EPFL2016