Functional and anatomical measures for outflow boundary conditions in atherosclerotic coronary bifurcations

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary anatomies combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. We propose an optimization-based method to personalize multiscale coronary flow simulations by assimilating clinical CT myocardial perfusion imaging and cardiac function measurements to yield patient-specific flow distributions and model parameters. Using this proof-of-concept study on a cohort of six patients, we reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based purely on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary models combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. Personalized flow distributions and model parameters are informed by clinical CT myocardial perfusion imaging and cardiac function using surrogate-based optimization. We reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

An in silico trials platform for the evaluation of stent design effect in post-implantation outcomes

Bioresorbable Vascular Scaffolds (BVS), developed to allow drug deliver and mechanical support, followed by complete resorption, have revolutionized atherosclerosis treatment. InSilc is a Cloud platform for in silico clinical trials (ISCT) used in the design, development and evaluation pipeline of stents. The platform integrates beyond the state-of-the-art multi-disciplinary and multiscale models, which predict the scaffold's performance in the short/acute and medium/long term. In this study, a use case scenario of two Bioabsorbable Vascular Stents (BVSs) implanted in the same arterial anatomy is presented, allowing the whole InSilc in silico pipeline to be applied and predict how the different aspects of this intervention affect the success of stenting process.

Reduced order methods for parametric optimal flow control in coronary bypass grafts, toward patient-specific data assimilation

Coronary artery bypass grafts (CABG) surgery is an invasive procedure performed to circumvent partial or complete blood flow blockage in coronary artery disease. In this work, we apply a numerical optimal flow control model to patient-specific geometries of CABG, reconstructed from clinical images of real-life surgical cases, in parameterized settings. The aim of these applications is to match known physiological data with numerical hemodynamics corresponding to different scenarios, arisen by tuning some parameters. Such applications are an initial step toward matching patient-specific physiological data in patient-specific vascular geometries as best as possible. Two critical challenges that reportedly arise in such problems are: (a) lack of robust quantification of meaningful boundary conditions required to match known data as best as possible and (b) high computational cost. In this work, we utilize unknown control variables in the optimal flow control problems to take care of the first challenge. Moreover, to address the second challenge, we propose a time-efficient and reliable computational environment for such parameterized problems by projecting them onto a low-dimensional solution manifold through proper orthogonal decomposition-Galerkin.

Coronary CT angiography-based estimation of myocardial perfusion territories for coronary artery FFR and wall shear stress simulation

This study aims to apply a CCTA-derived territory-based patient-specific estimation of boundary conditions for coronary artery fractional flow reserve (FFR) and wall shear stress (WSS) simulation. The non-invasive simulation can help diagnose the significance of coronary stenosis and the likelihood of myocardial ischemia. FFR is often regarded as the gold standard to evaluate the functional significance of stenosis in coronary arteries. In another aspect, proximal wall shear stress (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{{WSS}_{prox}}$$\end{document}WSSprox) can also be an indicator of plaque vulnerability. During the simulation process, the mass flow rate of the blood in coronary arteries is one of the most important boundary conditions. This study utilized the myocardium territory to estimate and allocate the mass flow rate. 20 patients are included in this study. From the knowledge of anatomical information of coronary arteries and the myocardium, the territory-based FFR and the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{{WSS}_{prox}}$$\end{document}WSSprox can both be derived from fluid dynamics simulations. Applying the threshold of distinguishing between significant and non-significant stenosis, the territory-based method can reach the accuracy, sensitivity, and specificity of 0.88, 0.90, and 0.80, respectively. For significantly stenotic cases (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{FFR}_{m}$$\end{document}FFRm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le$$\end{document}≤ 0.80), the vessels usually have higher wall shear stress in the proximal region of the lesion.

Comparison of LBM and FVM in the estimation of LAD stenosis

The finite volume method (FVM)-based computational fluid dynamics (CFD) technology has been applied in the non-invasive diagnosis of coronary artery stenosis. Nonetheless, FVM is a time-consuming process. In addition to FVM, the lattice Boltzmann method (LBM) is used in fluid flow simulation. Unlike FVM solving the Navier-Stokes equations, LBM directly solves the simplified Boltzmann equation, thus saving computational time. In this study, 12 patients with left anterior descending (LAD) stenosis, diagnosed by CTA, are analysed using FVM and LBM. The velocities, pressures, and wall shear stress (WSS) predicted using FVM and LBM for each patient is compared. In particular, the ratio of the average and maximum speed at the stenotic part characterising the degree of stenosis is compared. Finally, the golden standard of LAD stenosis, invasive fractional flow reserve (FFR), is applied to justify the simulation results. Our results show that LBM and FVM are consistent in blood flow simulation. In the region with a high degree of stenosis, the local flow patterns in those two solvers are slightly different, resulting in minor differences in local WSS estimation and blood speed ratio estimation. Notably, these differences do not result in an inconsistent estimation. Comparison with invasive FFR shows that, in most cases, the non-invasive diagnosis is consistent with FFR measurements. However, in some cases, the non-invasive diagnosis either underestimates or overestimates the degree of stenosis. This deviation is caused by the difference between physiological and simulation conditions that remains the biggest challenge faced by all CFD-based non-invasive diagnostic methods.

The improvement of the shear stress and oscillatory shear index of coronary arteries during Enhanced External Counterpulsation in patients with coronary heart disease

Background

Enhanced External Counterpulsation (EECP) can chronically relieve ischemic chest pain and improve the prognosis of coronary heart disease (CHD). Despite its role in mitigating heart complications, EECP and the mechanisms behind its therapeutic nature, such as its effects on blood flow hemodynamics, are still not fully understood. This study aims to elucidate the effect of EECP on significant hemodynamic parameters in the coronary arterial tree.

Methods

A finite volume method was used in conjunction with the inlet pressure wave (surrogated by the measured aortic pressure) before and during EECP and outlet flow resistance, assuming the blood as newtonian fluid. The time-average wall shear stress (TAWSS) and oscillatory shear index (OSI) were determined from the flow field.

Results

Regardless of the degree of vascular stenosis, hemodynamic conditions and flow patterns could be improved during EECP. In comparison with the original tree, the tree with a severe stenosis (75% area stenosis) demonstrated more significant improvement in hemodynamic conditions and flow patterns during EECP, with surface area ratio of TAWSS risk area (SAR-TAWSS) reduced from 12.3% to 6.7% (vs. SAR-TAWSS reduced from 7.2% to 5.6% in the original tree) and surface area ratio of OSI risk area (SAR-OSI) reduced from 6.8% to 2.5% (vs. SAR-OSI of both 0% before and during EECP in the original tree because of mild stenosis). Moreover, it was also shown that small ratio of diastolic pressure (D) and systolic pressure (S) (D/S) could only improve the hemodynamic condition mildly. The SAR-TAWSS reduction ratio significantly increased as D/S became larger.

Conclusions

A key finding of the study was that the improvement of hemodynamic conditions along the LMCA trees during EECP became more significant with the increase of D/S and the severity degree of stenoses at the bifurcation site. These findings have important implications on EECP as adjuvant therapy before or after percutaneous coronary intervention (PCI) in patients with diffuse atherosclerosis.

Assessment of boundary conditions for CFD simulation in human carotid artery

Computational fluid dynamics (CFD) is an increasingly used method for investigation of hemodynamic parameters and their alterations under pathological conditions, which are important indicators for diagnosis of cardiovascular disease. In hemodynamic simulation models, the employment of appropriate boundary conditions (BCs) determines the computational accuracy of the CFD simulation in comparison with pressure and velocity measurements. In this study, we have first assessed the influence of inlet boundary conditions on hemodynamic CFD simulations. We selected two typical patients suspected of carotid artery disease, with mild stenosis and severe stenosis. Both patients underwent digital subtraction angiography (DSA), magnetic resonance angiography, and the invasive pressure guide wire measured pressure profile. We have performed computational experiments to (1) study the hemodynamic simulation outcomes of distributions of wall shear stress, pressure, pressure gradient and (2) determine the differences in hemodynamic performances caused by inlet BCs derived from DSA and Womersley analytical solution. Our study has found that the difference is related to the severity of the stenosis; the greater the stenosis, the more the difference ensues. Further, in our study, the two typical subjects with invasively measured pressure profile and thirty subjects with ultrasound Doppler velocimeter (UDV) measurement served as the criteria to evaluate the hemodynamic outcomes of wall shear stress, pressure, pressure gradient and velocity due to different outlet BCs based on the Windkessel model, structured-tree model, and fully developed flow model. According to the pressure profiles, the fully developed model appeared to have more fluctuations compared with the other two models. The Windkessel model had more singularities before convergence. The three outlet BCs models also showed good correlation with the UDV measurement, while the Windkessel model appeared to be slightly better ([Formula: see text]). The structured-tree model was seen to have the best performance in terms of available computational cost and accuracy. The results of our numerical simulation and the good correlation with the computed pressure and velocity with their measurements have highlighted the effectiveness of CFD simulation in patient-specific human carotid artery with suspected stenosis.

Computational Fluid Dynamics as an Engineering Tool for the Reconstruction of Hemodynamics after Carotid Artery Stenosis Operation: A Case Study

Background and objectives: Brain ischemic stroke is caused by impaired or absolutely blocked blood flow into the brain regions. Despite the large number of possible origins, there is no general strategy for preventive treatment. In this paper, we aimed to predict the hemodynamics in a patient who experienced a critical stenosis operation in the carotid artery. This is a unique study where we used medical data together with the computational fluid (CFD) technique not to plan the surgery, but to predict its outcome. Materials and Methods: AngioCT data and blood perfusion of brain tissue (CT-perfusion) together with CFD technique were applied for stroke formation reconstruction in different clinical conditions. With the use of self-made semiautomatic algorithm for image processing and 3DDoctror software, 3D-vascular geometries before and after surgical intervention were reconstructed. As the paper is focused on the analysis of stroke appearance, apparent stroke was simulated as higher and lower pressure values in the cranial part due to different outcomes of the surgical intervention. This allowed to investigate the influence of spatial configuration and pressure values on blood perfusion in the analyzed circulatory system. Results: Application of CFD simulations for blood flow reconstruction for clinical conditions in the circulatory system accomplished on average 98.5% and 98.7% accuracy for CFD results compared to US-Doppler before and after surgical intervention, respectively. Meanwhile, CFD results compared to CT-perfusion indicated an average 89.7% and 92.8% accuracy before and after surgical intervention, respectively. Thus, the CFD is a reliable approach for predicting the patient hemodynamics, as it was confirmed by postoperative data. Conclusions: Our study indicated that the application of CFD simulations for blood flow reconstruction for clinical conditions in circulatory system reached 98% and 90% accuracy for US-Doppler and CT-perfusion, respectively. Therefore, the proposed method might be used as a tool for reconstruction of specific patients' hemodynamics after operation of critical stenosis in the carotid artery. However, further studies are necessary to confirm its usefulness in clinical practice.

Patient-Specific Modeling of Stented Coronary Arteries Reconstructed from Optical Coherence Tomography: Towards a Widespread Clinical Use of Fluid Dynamics Analyses

The recent widespread application of optical coherence tomography (OCT) in interventional cardiology has improved patient-specific modeling of stented coronary arteries for the investigation of local hemodynamics. In this review, the workflow for the creation of fluid dynamics models of stented coronary arteries from OCT images is presented. The algorithms for lumen contours and stent strut detection from OCT as well as the reconstruction methods of stented geometries are discussed. Furthermore, the state of the art of studies that investigate the hemodynamics of OCT-based stented coronary artery geometries is reported. Although those studies analyzed few patient-specific cases, the application of the current reconstruction methods of stented geometries to large populations is possible. However, the improvement of these methods and the reduction of the time needed for the entire modeling process are crucial for a widespread clinical use of the OCT-based models and future in silico clinical trials.

Fusion of fibrous cap thickness and wall shear stress to assess plaque vulnerability in coronary arteries: a pilot study

Purpose

Identification of rupture-prone plaques in coronary arteries is a major clinical challenge. Fibrous cap thickness and wall shear stress are two relevant image-based risk factors, but these two parameters are generally computed and analyzed separately. Accordingly, combining these two parameters can potentially improve the identification of at-risk regions. Therefore, the purpose of this study is to investigate the feasibility of the fusion of wall shear stress and fibrous cap thickness of coronary arteries in patient data.

Methods

Fourteen patients were included in this pilot study. Imaging of the coronary arteries was performed with optical coherence tomography and with angiography. Fibrous cap thickness was automatically quantified from optical coherence tomography pullbacks using a contour segmentation approach based on fast marching. Wall shear stress was computed by applying computational fluid dynamics on the 3D volume reconstructed from two angiograms. The two parameters then were co-registered using anatomical landmarks such as side branches.

Results

The two image modalities were successfully co-registered, with a mean (±SD) error corresponding to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8.6\,\pm \,6.7\,\%$$\end{document}8.6±6.7% of the length of the analyzed region. For all the analyzed participants, the average thinnest portion of each fibrous cap was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$129\,\pm \,69\,\upmu \text {m}$$\end{document}129±69μm, and the average WSS value at the location of the fibrous cap was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.46\,\pm \,1.16\,\text {Pa}$$\end{document}1.46±1.16Pa. A unique index was finally generated for each patient via the fusion of fibrous cap thickness and wall shear stress measurements, to translate all the measured parameters into a single risk map.

Conclusion

The introduced risk map integrates two complementary parameters and has potential to provide valuable information about plaque vulnerability.

Influence of the Accuracy of Angiography-Based Reconstructions on Velocity and Wall Shear Stress Computations in Coronary Bifurcations: A Phantom Study

INTRODUCTION

Wall shear stress (WSS) plays a key role in the onset and progression of atherosclerosis in human coronary arteries. Especially sites with low and oscillating WSS near bifurcations have a higher propensity to develop atherosclerosis. WSS computations in coronary bifurcations can be performed in angiography-based 3D reconstructions. It is essential to evaluate how reconstruction errors influence WSS computations in mildly-diseased coronary bifurcations. In mildly-diseased lesions WSS could potentially provide more insight in plaque progression.

MATERIALS METHODS

Four Plexiglas phantom models of coronary bifurcations were imaged with bi-plane angiography. The lumens were segmented by two clinically experienced readers. Based on the segmentations 3D models were generated. This resulted in three models per phantom: one gold-standard from the phantom model itself, and one from each reader. Steady-state and transient simulations were performed with computational fluid dynamics to compute the WSS. A similarity index and a noninferiority test were used to compare the WSS in the phantoms and their reconstructions. The margin for this test was based on the resolution constraints of angiography.

RESULTS

The reconstruction errors were similar to previously reported data; in seven out of eight reconstructions less than 0.10 mm. WSS in the regions proximal and far distal of the stenosis showed a good agreement. However, the low WSS areas directly distal of the stenosis showed some disagreement between the phantoms and the readers. This was due to small deviations in the reconstruction of the stenosis that caused differences in the resulting jet, and consequently the size and location of the low WSS area.

DISCUSSION

This study showed that WSS can accurately be computed within angiography-based 3D reconstructions of coronary arteries with early stage atherosclerosis. Qualitatively, there was a good agreement between the phantoms and the readers. Quantitatively, the low WSS regions directly distal to the stenosis were sensitive to small reconstruction errors.