TY - JOUR
T1 - A computational model applied to myocardial perfusion in the human heart
T2 - From large coronaries to microvasculature
AU - Di Gregorio, Simone
AU - Fedele, Marco
AU - Pontone, Gianluca
AU - Corno, Antonio F.
AU - Zunino, Paolo
AU - Vergara, Christian
AU - Quarteroni, Alfio
N1 - Funding Information:
This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 740132 , iHEART - An Integrated Heart Model for the simulation of the cardiac function, P.I. Prof. A. Quarteroni).
Publisher Copyright:
© 2020 Elsevier Inc.
Copyright:
Copyright 2020 Elsevier B.V., All rights reserved.
PY - 2021/1/1
Y1 - 2021/1/1
N2 - In this paper we present a mathematical and numerical model for human cardiac perfusion which accounts for the different length scales of the vessels in the coronary tree. Epicardial vessels are represented with fully three-dimensional (3D) fluid-dynamics, whereas intramural vessels are modeled as a multi-compartment porous medium. The coupling of these models takes place through interface conditions based on the continuity of mass and momentum. Instead, is neglected in this first preliminary model the myocardium deformation. To estimate the physical parameters of the multi-compartment model, a virtual intramural vascular network is generated using a novel algorithm which works in non-convex domains. Modeling epicardial vessels with a 3D model and intramural ones with a porous medium approach makes it possible to apply the proposed strategy to patient-specific heart geometries reconstructed from clinical imaging data. We also address the derivation of numerical solvers for the coupled problem. In particular, we propose a splitting algorithm for the monolithic problem, with the corresponding convergence analysis performed in a simplified linearized case, and a suitable preconditioner for the multi-compartment porous sub-model. Finally, we test the computational framework in a realistic human heart, obtaining results that fall in the physiological range for both pressures and local myocardial flows.
AB - In this paper we present a mathematical and numerical model for human cardiac perfusion which accounts for the different length scales of the vessels in the coronary tree. Epicardial vessels are represented with fully three-dimensional (3D) fluid-dynamics, whereas intramural vessels are modeled as a multi-compartment porous medium. The coupling of these models takes place through interface conditions based on the continuity of mass and momentum. Instead, is neglected in this first preliminary model the myocardium deformation. To estimate the physical parameters of the multi-compartment model, a virtual intramural vascular network is generated using a novel algorithm which works in non-convex domains. Modeling epicardial vessels with a 3D model and intramural ones with a porous medium approach makes it possible to apply the proposed strategy to patient-specific heart geometries reconstructed from clinical imaging data. We also address the derivation of numerical solvers for the coupled problem. In particular, we propose a splitting algorithm for the monolithic problem, with the corresponding convergence analysis performed in a simplified linearized case, and a suitable preconditioner for the multi-compartment porous sub-model. Finally, we test the computational framework in a realistic human heart, obtaining results that fall in the physiological range for both pressures and local myocardial flows.
KW - Cardiac perfusion
KW - Finite elements
KW - Intramural vessel network
KW - Iterative numerical scheme
KW - Multi-compartment Darcy model
KW - Perfusion regions
UR - http://www.scopus.com/inward/record.url?scp=85091902432&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85091902432&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2020.109836
DO - 10.1016/j.jcp.2020.109836
M3 - Article
AN - SCOPUS:85091902432
VL - 424
JO - Journal of Computational Physics
JF - Journal of Computational Physics
SN - 0021-9991
M1 - 109836
ER -