KAUST DepartmentExtreme Computing Research Center
Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
Online Publication Date2019-01-05
Print Publication Date2018
Permanent link to this recordhttp://hdl.handle.net/10754/631257
MetadataShow full item record
AbstractWe 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.
CitationFranzone PC, Pavarino LF, Scacchi S, Zampini S (2018) Scalable Cardiac Electro-Mechanical Solvers and Reentry Dynamics. Domain Decomposition Methods in Science and Engineering XXIV: 31–43. Available: http://dx.doi.org/10.1007/978-3-319-93873-8_3.