Abstract
We present a scalable solver for the three-dimensional cardiac electro-mechanical coupling (EMC) model, which represents, currently, the most complete mathematical description of the interplay between the electrical and mechanical phenomena occurring during a heartbeat. The most computational demanding parts of the EMC model are: the electrical current flow model of the cardiac tissue, called Bidomain model, consisting of two non-linear partial differential equations of reaction-diffusion type; the quasi-static finite elasticity model for the deformation of the cardiac tissue. Our finite element parallel solver is based on: Block Jacobi and Multilevel Additive Schwarz preconditioners for the solution of the linear systems deriving from the discretization of the Bidomain equations; Newton-Krylov-Algebraic-Multigrid or Newton-Krylov-BDDC algorithms for the solution of the non-linear algebraic system deriving from the discretization of the finite elasticity equations. Three-dimensional numerical test on two linux clusters show the effectiveness and scalability of the EMC solver in simulating both physiological and pathological cardiac dynamics.
This is a preview of subscription content, log in via an institution to check access.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
S. Balay et al., PETSc users manual. Tech. Rep. ANL-95/11 - Revision 3.3, Argonne National Laboratory, 2012.
P. Colli Franzone, L.F. Pavarino, S. Scacchi, Mathematical Cardiac Electrophysiology, MSA, vol. 13 (Springer, New York, 2014)
P. Colli Franzone, L.F. Pavarino, S. Scacchi, Bioelectrical effects of mechanical feedbacks in a strongly coupled cardiac electro-mechanical model. Math. Models Methods Appl. Sci. 26, 27–57 (2016)
P. Colli Franzone, L.F. Pavarino, S. Scacchi, Effects of mechanical feedback on the stability of cardiac scroll waves: a bidomain electro-mechanical simulation study. Chaos 27, 093905 (2017)
C.R. Dohrmann, A preconditioner for substructuring based on constrained energy minimization. SIAM J. Sci. Comput. 25, 246–258 (2003)
T.S.E. Eriksson et al., Influence of myocardial fiber/sheet orientations on left ventricular mechanical contraction. Math. Mech. Solids 18, 592–606 (2013)
V. Gurev et al., Models of cardiac electromechanics based on individual hearts imaging data: image-based electromechanical models of the heart. Biomech. Model Mechanobiol. 10, 295–306 (2011)
A. Klawonn, O. Rheinbach, Highly scalable parallel domain decomposition methods with an application to biomechanics. ZAMM-Z. Angew. Math. Mech. 90, 5–32 (2010)
A. Klawonn, O.B. Widlund, Dual-primal FETI methods for linear elasticity. Commun. Pure Appl. Math., 59, 1523–1572 (2006)
S. Land et al., An analysis of deformation-dependent electromechanical coupling in the mouse heart. J. Physiol. 590, 4553–4569 (2012)
J. Mandel, C.R. Dohrmann, Convergence of a balancing domain decomposition by constraints and energy minimization. Linear Algebra Appl. 10, 639–659 (2003)
S.A. Niederer, N.P. Smith, A mathematical model of the slow force response to stretch in rat ventricular myocites. Biophys. J. 92, 4030–4044 (2007)
L.F. Pavarino, S. Scacchi, Multilevel additive Schwarz preconditioners for the Bidomain reaction-diffusion system. SIAM J. Sci. Comput. 31, 420–443 (2008)
L.F. Pavarino, S. Zampini, O.B. Widlund, BDDC preconditioners for spectral element discretizations of almost incompressible elasticity in three dimensions. SIAM J. Sci. Comput. 32(6), 3604–3626 (2010)
L.F. Pavarino, S. Scacchi, S. Zampini, Newton-krylov-BDDC solvers for non-linear cardiac mechanics. Comput. Methods Appl. Mech. Eng. 295, 562–580 (2015)
S. Rossi et al., Orthotropic active strain models for the numerical simulation of cardiac biomechanics. Int. J. Numer. Methods Biomed. Eng. 28, 761–788 (2012)
J. Sundnes et al., Improved discretisation and linearisation of active tension in strongly coupled cardiac electro-mechanics simulations. Comput. Methods Biomech. Biomed. Eng. 17, 604–615 (2014)
K.H.W.J. ten Tusscher et al., A model for human ventricular tissue. Am. J. Phys. Heart. Circ. Physiol. 286, H1573–H1589 (2004)
K.H.W.J. ten Tusscher, A.V. Panfilov, Alternans and spiral breakup in a human ventricular tissue model. Am. J. Physiol. Heart Circ. Physiol. 291, H1088–H1100 (2006)
A. Toselli, O.B. Widlund, Domain Decomposition Methods: Algorithms and Theory (Springer, Berlin, 2004)
F.J. Vetter, A.D. McCulloch, Three-dimensional stress and strain in passive rabbit left ventricle: a model study. Ann. Biomed. Eng. 28, 781–792 (2000)
S. Zampini, Dual-primal methods for the cardiac bidomain model. Math. Models Methods Appl. Sci. 24, 667–696 (2014)
S. Zampini, PCBDDC: a class of robust dual-primal preconditioners in PETSc. SIAM J. Sci. Comput. 38(5), S282–S306 (2016)
Author information
Authors and Affiliations
Corresponding authors
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2018 Springer International Publishing AG, part of Springer Nature
About this paper
Cite this paper
Franzone, P.C., Pavarino, L.F., Scacchi, S., Zampini, S. (2018). Scalable Cardiac Electro-Mechanical Solvers and Reentry Dynamics. In: Bjørstad, P., et al. Domain Decomposition Methods in Science and Engineering XXIV . DD 2017. Lecture Notes in Computational Science and Engineering, vol 125. Springer, Cham. https://doi.org/10.1007/978-3-319-93873-8_3
Download citation
DOI: https://doi.org/10.1007/978-3-319-93873-8_3
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-93872-1
Online ISBN: 978-3-319-93873-8
eBook Packages: Mathematics and StatisticsMathematics and Statistics (R0)