A computational model applied to myocardial perfusion in the human heart: From large coronaries to microvasculature

Published on Jan 1, 2021in Journal of Computational Physics2.985
· DOI :10.1016/J.JCP.2020.109836
Simone Di Gregorio3
Estimated H-index: 3
(Polytechnic University of Milan),
Marco Fedele5
Estimated H-index: 5
(Polytechnic University of Milan)
+ 4 AuthorsAlfio Quarteroni83
Estimated H-index: 83
(EPFL: École Polytechnique Fédérale de Lausanne)
Abstract 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.
📖 Papers frequently viewed together
47 Citations
9 Authors (John Lee, ..., Nicolas P. Smith)
11 Citations
#1Ettore VidottoH-Index: 4
#2Timo KochH-Index: 6
Last. Barbara WohlmuthH-Index: 51
view all 5 authors...
In this paper, we are concerned with the simulation of blood flow in microvascular networks and the surrounding tissue. To reduce the computational complexity of this issue, the network structures ...
14 CitationsSource
#1Clément Daviller (UCBL: Claude Bernard University Lyon 1)H-Index: 1
#2Thomas Grenier (UCBL: Claude Bernard University Lyon 1)H-Index: 9
Last. Magalie Viallon (French Institute of Health and Medical Research)H-Index: 33
view all 6 authors...
Abstract Even if cardiovascular magnetic resonance (CMR) perfusion imaging has proven its relevance for visual detection of ischemia, myocardial blood flow (MBF) quantification at the voxel observation scale remains challenging. Integration of an automated segmentation step, prior to perfusion index estimation, might be a significant reconstruction component that could allow sustainable assumptions and constraint enlargement prior to advanced modeling. Current clustering techniques, such as bull...
1 CitationsSource
#1Alfio QuarteroniH-Index: 83
#2Luca DedèH-Index: 21
Last. Christian VergaraH-Index: 24
view all 4 authors...
Mathematical and numerical modelling of the human cardiovascular system has attracted remarkable research interest due to its intrinsic mathematical difficulty and the increasing impact of cardiovascular diseases worldwide. This book addresses the two principal components of the cardiovascular system: arterial circulation and heart function. It systematically describes all aspects of the problem, stating the basic physical principles, analysing the associated mathematical models that comprise PD...
13 Citations
#1Liza ThomasH-Index: 30
#2Thomas H. Marwick (Baker IDI Heart and Diabetes Institute)H-Index: 137
Last. Luigi P. Badano (University of Milano-Bicocca)H-Index: 69
view all 5 authors...
Highlights •LA function, measured as LA reservoir strain, is an important metric in diagnosing LVDD, grading its severity, and monitoring the effects of treatment. •LA volume has been used as a surrogate of the chronicity and severity of LVDD. However, volume is an insensitive biomarker of the early phases of LVDD. •LA phasic function can be assessed both by volumetric analysis, using 3-dimensional echocardiography, and by strain/strain-rate analysis, using speckle-tracking echocardiography. The...
91 CitationsSource
#1Clara Jaquet (University of Paris)H-Index: 3
#2Laurent Najman (Université Paris-Saclay)H-Index: 35
Last. Charles TaylorH-Index: 162
view all 9 authors...
Objective : In this paper, we propose an algorithm for the generation of a patient-specific cardiac vascular network starting from segmented epicardial vessels down to the arterioles. Method : We extend a tree generation method based on satisfaction of functional principles, named constrained constructive optimization, to account for multiple, competing vascular trees. The algorithm simulates angiogenesis under vascular volume minimization with flow-related and geometrical constraints adapting t...
12 CitationsSource
#1Ruben Doste (UPF: Pompeu Fabra University)H-Index: 2
#2David Soto-Iglesias (UPF: Pompeu Fabra University)H-Index: 12
Last. Oscar Camara (UPF: Pompeu Fabra University)H-Index: 25
view all 10 authors...
Rule-based methods are often used for assigning fiber orientation to cardiac anatomical models. However, existing methods have been developed using data mostly from the left ventricle. As a consequence, fiber information obtained from rule-based methods often does not match histological data in other areas of the heart such as the right ventricle, having a negative impact in cardiac simulations beyond the left ventricle. In this work, we present a rule-based method where fiber orientation is sep...
27 CitationsSource
#1Ravi NamaniH-Index: 7
#2Ghassan S. KassabH-Index: 34
Last. Yoram Lanir (Technion – Israel Institute of Technology)H-Index: 29
view all 3 authors...
Experimental limitations in measurements of coronary flow in the beating heart have led to the development of in silico models of reconstructed coronary trees. Previous coronary reconstructions relied primarily on anatomical data, including statistical morphometry (e.g., diameters, length, connectivity, longitudinal position). Such reconstructions are non-unique, however, often leading to unrealistic predicted flow features. Thus, it is necessary to impose physiological flow constraints to ensur...
6 CitationsSource
#1Yonggyun Bae (KU: Korea University)H-Index: 5
#2Seung Tae Hwang (KU: Korea University)H-Index: 2
Last. Jin-Ho Choi (SMC: Samsung Medical Center)H-Index: 43
view all 9 authors...
Invasive procedure is a prerequisite for studying coronary physiology. We established the measurement of non-invasive physiological parameters including coronary blood flow (CBF), flow velocity, and microvascular resistance using coronary computed tomography angiography (CCTA). Vessel-specific CBF was derived from transluminal attenuation flow encoding (TAFE) and then tested using three separate datasets consisted of computational simulation, human perfusion CT, and human CCTA. TAFE-derived CBF ...
5 CitationsSource
#1Gianluca PontoneH-Index: 54
#2Daniele Andreini (University of Milan)H-Index: 46
Last. Mauro PepiH-Index: 45
view all 25 authors...
Abstract Objectives The goal of this study was to evaluate the diagnostic accuracy of stress computed tomography myocardial perfusion (CTP) for the detection of functionally significant coronary artery disease (CAD) by using invasive coronary angiography (ICA) plus invasive fractional flow reserve (FFR) as the reference standard in consecutive intermediate- to high-risk symptomatic patients. Background Stress CTP recently emerged as a potential strategy to combine the anatomic and functional eva...
45 CitationsSource
#1Eduard Rohan (University of West Bohemia)H-Index: 15
#2Vladimír Lukeš (University of West Bohemia)H-Index: 7
Last. Alena Jonášová (University of West Bohemia)H-Index: 6
view all 3 authors...
The paper deals with modeling the liver perfusion intended to improve quantitative analysis of the tissue scans provided by the contrast-enhanced computed tomography (CT). For this purpose, we developed a model of dynamic transport of the contrast fluid through the hierarchies of the perfusion trees. Conceptually, computed time-space distributions of the so-called tissue density can be compared with the measured data obtained from CT; such a modeling feedback can be used for model parameter iden...
8 CitationsSource
Cited By5
#1Hana Kolesová (First Faculty of Medicine, Charles University in Prague)H-Index: 4
#1Hana Kolesová (First Faculty of Medicine, Charles University in Prague)
Last. David Sedmera (First Faculty of Medicine, Charles University in Prague)H-Index: 38
view all 5 authors...
Summary Tissue imaging in 3D using visible light is limited and various clearing techniques were developed to increase imaging depth, but none provides universal solution for all tissues at all developmental stages. In this review, we focus on different tissue clearing methods for 3D imaging of heart and vasculature, based on chemical composition (solvent-based, simple immersion, hyperhydration, and hydrogel embedding techniques). We discuss in detail compatibility of various tissue clearing tec...
#1Gonzalo D. Maso Talou (University of Auckland)H-Index: 4
#2Soroush Safaei (University of Auckland)H-Index: 7
Last. Pablo J. Blanco (NIST: National Institute of Standards and Technology)H-Index: 24
view all 4 authors...
Mimicking angiogenetic processes in vascular territories acquires importance in the analysis of the multi-scale circulatory cascade and the coupling between blood flow and cell function. The present work extends, in several aspects, the Constrained Constructive Optimisation (CCO) algorithm to tackle complex automatic vascularisation tasks. The main extensions are based on the integration of adaptive optimisation criteria and multi-staged space-filling strategies which enhance the modelling capab...
#1Marco Fedele (Polytechnic University of Milan)H-Index: 5
#2Alfio Quarteroni (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 83
In order to simulate the cardiac function for a patient-specific geometry, the generation of the computational mesh is crucially important. In practice, the input is typically a set of unprocessed polygonal surfaces coming either from a template geometry or from medical images. These surfaces need ad-hoc processing to be suitable for a volumetric mesh generation. In this work we propose a set of new algorithms and tools aiming to facilitate the mesh generation process. In particular, we focus on...
#1Luca DedèH-Index: 21
#2Francesco Regazzoni (Polytechnic University of Milan)H-Index: 28
Last. Alfio QuarteroniH-Index: 83
view all 10 authors...
Emerging studies address how COVID-19 infection can impact the human cardiovascular system This relates particularly to the development of myocardial injury, acute coronary syndrome, myocarditis, arrhythmia, and heart failure Prospective treatment approach is advised for these patients To study the interplay between local changes (reduced contractility), global variables (peripheral resistances, heart rate) and the cardiac function, we considered a lumped parameters computational model of the ca...
#1Jakub Wiktor BothH-Index: 5
#2Nicolás BarnafiH-Index: 1
Last. Alfio QuarteroniH-Index: 83
view all 5 authors...
We address numerical solvers for a poromechanics model particularly adapted for soft materials, as it generally respects thermodynamics principles and energy balance. Considering the multi-physics nature of the problem, which involves solid and fluid species, interacting on the basis of mass balance and momentum conservation, we decide to adopt a solution strategy of the discrete problem based on iterative splitting schemes. As the model is similar (but not equivalent to) the Biot poromechanics ...