ABSTRACT
Accurate modeling of blood dynamics in the coronary microcirculation is a crucial step toward the clinical application of in silico methods for the diagnosis of coronary artery disease. In this work, we present a new mathematical model of microcirculatory hemodynamics accounting for microvasculature compliance and cardiac contraction; we also present its application to a full simulation of hyperemic coronary blood flow and 3D myocardial perfusion in real clinical cases. Microvasculature hemodynamics is modeled with a compliant multi-compartment Darcy formulation, with the new compliance terms depending on the local intramyocardial pressure generated by cardiac contraction. Nonlinear analytical relationships for vessels distensibility are included based on experimental data, and all the parameters of the model are reformulated based on histologically relevant quantities, allowing a deeper model personalization. Phasic flow patterns of high arterial inflow in diastole and venous outflow in systole are obtained, with flow waveforms morphology and pressure distribution along the microcirculation reproduced in accordance with experimental and in vivo measures. Phasic diameter change for arterioles and capillaries is also obtained with relevant differences depending on the depth location. Coronary blood dynamics exhibits a disturbed flow at the systolic onset, while the obtained 3D perfusion maps reproduce the systolic impediment effect and show relevant regional and transmural heterogeneities in myocardial blood flow (MBF). The proposed model successfully reproduces microvasculature hemodynamics over the whole heartbeat and along the entire intramural vessels. Quantification of phasic flow patterns, diameter changes, regional and transmural heterogeneities in MBF represent key steps ahead in the direction of the predictive simulation of cardiac perfusion.
ABSTRACT
Predictive modeling of hyperemic coronary and myocardial blood flow (MBF) greatly supports diagnosis and prognostic stratification of patients suffering from coronary artery disease (CAD). In this work, we propose a novel strategy, using only readily available clinical data, to build personalized inlet conditions for coronary and MBF models and to achieve an effective calibration for their predictive application to real clinical cases. Experimental data are used to build personalized pressure waveforms at the aortic root, representative of the hyperemic state and adapted to surrogate the systolic contraction, to be used in computational fluid-dynamics analyses. Model calibration to simulate hyperemic flow is performed in a "blinded" way, not requiring any additional exam. Coronary and myocardial flow simulations are performed in eight patients with different clinical conditions to predict FFR and MBF. Realistic pressure waveforms are recovered for all the patients. Consistent pressure distribution, blood velocities in the large arteries, and distribution of MBF in the healthy myocardium are obtained. FFR results show great accuracy with a per-vessel sensitivity and specificity of 100% according to clinical threshold values. Mean MBF shows good agreement with values from stress-CTP, with lower values in patients with diagnosed perfusion defects. The proposed methodology allows us to quantitatively predict FFR and MBF, by the exclusive use of standard measures easily obtainable in a clinical context. This represents a fundamental step to avoid catheter-based exams and stress tests in CAD diagnosis.
