×

A monolithic algorithm for the simulation of cardiac electromechanics in the human left ventricle. (English) Zbl 1442.92006

Summary: In this paper, we propose a monolithic algorithm for the numerical solution of the electromechanics model of the left ventricle in the human heart. Our coupled model integrates the monodomain equation with the Bueno-Orovio minimal model for electrophysiology and the Holzapfel-Ogden constitutive law for the passive mechanics of the myocardium; a distinguishing feature of our electromechanics model is the use of the active strain formulation for muscle contraction, which we exploit – for the first time in this context – by means of a transmurally variable active strain formulation. We use the finite element method for space discretization and backward differentiation formulas for time discretization, which we consider for both implicit and semi-implicit schemes. We compare and discuss the two schemes in terms of computational efficiency as the semi-implicit scheme poses significant restrictions on the time step size due to stability considerations, while the implicit scheme yields instead a nonlinear problem, which we solve by means of the Newton method. Emphasis is laid on preconditioning strategy of the linear solver, which we perform by factorizing a block Gauss-Seidel preconditioner in combination with combination with parallel preconditioners for each of the single core models composing the integrated electromechanics model. We carry out several numerical simulations in the high performance computing framework for both idealized and patient-specific left ventricle geometries and meshes, which we obtain by segmenting medical MRI images. We produce personalized pressure-volume loops by means of the computational procedure, which we use to synthetically interpret and analyze the outputs of the electromechanics model.

MSC:

92C10 Biomechanics
92C30 Physiology (general)
35Q92 PDEs in connection with biology, chemistry and other natural sciences
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
PDFBibTeX XMLCite
Full Text: DOI

References:

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.