Uncertainty quantification of simulated biomechanical stimuli in coronary artery bypass grafts

Personalized and uncertainty-aware coronary hemodynamics simulations: From Bayesian estimation to improved multi-fidelity uncertainty quantification

Background

Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline.

Objective

We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data.

Methods

We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs.

Results

Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance.

Conclusions

The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

Quantification of total uncertainty in the physics-informed reconstruction of CVSim-6 physiology

When predicting physical phenomena through simulation, quantification of the total uncertainty due to multiple sources is as crucial as making sure the underlying numerical model is accurate. Possible sources include irreducible aleatoric uncertainty due to noise in the data, epistemic uncertainty induced by insufficient data or inadequate parameterization, and model-form uncertainty related to the use of misspecified model equations. In addition, recently proposed approaches provide flexible ways to combine information from data with full or partial satisfaction of equations that typically encode physical principles. Physics-based regularization interacts in nontrivial ways with aleatoric, epistemic and model-form uncertainty and their combination, and a better understanding of this interaction is needed to improve the predictive performance of physics-informed digital twins that operate under real conditions. To better understand this interaction, with a specific focus on biological and physiological models, this study investigates the decomposition of total uncertainty in the estimation of states and parameters of a differential system simulated with MC X-TFC, a new physics-informed approach for uncertainty quantification based on random projections and Monte-Carlo sampling. After an introductory comparison between approaches for physics-informed estimation, MC X-TFC is applied to a six-compartment stiff ODE system, the CVSim-6 model, developed in the context of human physiology. The system is first analyzed by progressively removing data while estimating an increasing number of parameters, and subsequently by investigating total uncertainty under model-form misspecification of non-linear resistance in the pulmonary compartment. In particular, we focus on the interaction between the formulation of the discrepancy term and quantification of model-form uncertainty, and show how additional physics can help in the estimation process. The method demonstrates robustness and efficiency in estimating unknown states and parameters, even with limited, sparse, and noisy data. It also offers great flexibility in integrating data with physics for improved estimation, even in cases of model misspecification.

Credibility assessment of in silico clinical trials for medical devices

In silico clinical trials (ISCTs) are an emerging method in modeling and simulation where medical interventions are evaluated using computational models of patients. ISCTs have the potential to provide cost-effective, time-efficient, and ethically favorable alternatives for evaluating the safety and effectiveness of medical devices. However, ensuring the credibility of ISCT results is a significant challenge. This paper aims to identify unique considerations for assessing the credibility of ISCTs and proposes an ISCT credibility assessment workflow based on recently published model assessment frameworks. First, we review various ISCTs described in the literature, carefully selected to showcase the range of methodological options available. These studies cover a wide variety of devices, reasons for conducting ISCTs, patient model generation approaches including subject-specific versus 'synthetic' virtual patients, complexity levels of devices and patient models, incorporation of clinician or clinical outcome models, and methods for integrating ISCT results with real-world clinical trials. We next discuss how verification, validation, and uncertainty quantification apply to ISCTs, considering the range of ISCT approaches identified. Based on our analysis, we then present a hierarchical workflow for assessing ISCT credibility, using a general credibility assessment framework recently published by the FDA's Center for Devices and Radiological Health. Overall, this work aims to promote standardization in ISCTs and contribute to the wider adoption and acceptance of ISCTs as a reliable tool for evaluating medical devices.

Integrated modeling and simulation of recruitment of myocardial perfusion and oxygen delivery in exercise

While exercise-mediated vasoregulation in the myocardium is understood to be governed by autonomic, myogenic, and metabolic-mediated mechanisms, we do not yet understand the spatial heterogeneity of vasodilation or its effects on microvascular flow patterns and oxygen delivery. This study uses a simulation and modeling approach to explore the mechanisms underlying the recruitment of myocardial perfusion and oxygen delivery in exercise. The simulation approach integrates model components representing: whole-body cardiovascular hemodynamics, cardiac mechanics and myocardial work; myocardial perfusion; and myocardial oxygen transport. Integrating these systems together, model simulations reveal: (1.) To match expected flow and transmural flow ratios at increasing levels of exercise, a greater degree of vasodilation must occur in the subendocardium compared to the subepicardium. (2.) Oxygen extraction and venous oxygenation are predicted to substantially decrease with increasing exercise level preferentially in the subendocardium, suggesting that an oxygen-dependent error signal driving metabolic mediated recruitment of flow would be operative only in the subendocardium. (3.) Under baseline physiological conditions approximately 4% of the oxygen delivered to the subendocardium may be supplied via retrograde flow from coronary veins.

Influence of material parameter variability on the predicted coronary artery biomechanical environment via uncertainty quantification

Central to the clinical adoption of patient-specific modeling strategies is demonstrating that simulation results are reliable and safe. Indeed, simulation frameworks must be robust to uncertainty in model input(s), and levels of confidence should accompany results. In this study, we applied a coupled uncertainty quantification-finite element (FE) framework to understand the impact of uncertainty in vascular material properties on variability in predicted stresses. Univariate probability distributions were fit to material parameters derived from layer-specific mechanical behavior testing of human coronary tissue. Parameters were assumed to be probabilistically independent, allowing for efficient parameter ensemble sampling. In an idealized coronary artery geometry, a forward FE model for each parameter ensemble was created to predict tissue stresses under physiologic loading. An emulator was constructed within the UncertainSCI software using polynomial chaos techniques, and statistics and sensitivities were directly computed. Results demonstrated that material parameter uncertainty propagates to variability in predicted stresses across the vessel wall, with the largest dispersions in stress within the adventitial layer. Variability in stress was most sensitive to uncertainties in the anisotropic component of the strain energy function. Moreover, unary and binary interactions within the adventitial layer were the main contributors to stress variance, and the leading factor in stress variability was uncertainty in the stress-like material parameter that describes the contribution of the embedded fibers to the overall artery stiffness. Results from a patient-specific coronary model confirmed many of these findings. Collectively, these data highlight the impact of material property variation on uncertainty in predicted artery stresses and present a pipeline to explore and characterize forward model uncertainty in computational biomechanics.

Influence of Material Parameter Variability on the Predicted Coronary Artery Biomechanical Environment via Uncertainty Quantification

Central to the clinical adoption of patient-specific modeling strategies is demonstrating that simulation results are reliable and safe. Indeed, simulation frameworks must be robust to uncertainty in model input(s), and levels of confidence should accompany results. In this study, we applied a coupled uncertainty quantification-finite element (FE) framework to understand the impact of uncertainty in vascular material properties on variability in predicted stresses. Univariate probability distributions were fit to material parameters derived from layer-specific mechanical behavior testing of human coronary tissue. Parameters were assumed to be probabilistically independent, allowing for efficient parameter ensemble sampling. In an idealized coronary artery geometry, a forward FE model for each parameter ensemble was created to predict tissue stresses under physiologic loading. An emulator was constructed within the UncertainSCI software using polynomial chaos techniques, and statistics and sensitivities were directly computed. Results demonstrated that material parameter uncertainty propagates to variability in predicted stresses across the vessel wall, with the largest dispersions in stress within the adventitial layer. Variability in stress was most sensitive to uncertainties in the anisotropic component of the strain energy function. Moreover, unary and binary interactions within the adventitial layer were the main contributors to stress variance, and the leading factor in stress variability was uncertainty in the stress-like material parameter that describes the contribution of the embedded fibers to the overall artery stiffness. Results from a patient-specific coronary model confirmed many of these findings. Collectively, these data highlight the impact of material property variation on uncertainty in predicted artery stresses and present a pipeline to explore and characterize forward model uncertainty in computational biomechanics.

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.

Reducing Geometric Uncertainty in Computational Hemodynamics by Deep Learning-Assisted Parallel-Chain MCMC

Computational hemodynamic modeling has been widely used in cardiovascular research and healthcare. However, the reliability of model predictions is largely dependent on the uncertainties of modeling parameters and boundary conditions, which should be carefully quantified and further reduced with available measurements. In this work, we focus on propagating and reducing the uncertainty of vascular geometries within a Bayesian framework. A novel deep learning (DL)-assisted parallel Markov chain Monte Carlo (MCMC) method is presented to enable efficient Bayesian posterior sampling and geometric uncertainty reduction. A DL model is built to approximate the geometry-to-hemodynamic map, which is trained actively using online data collected from parallel MCMC chains and utilized for early rejection of unlikely proposals to facilitate convergence with less expensive full-order model evaluations. Numerical studies on two-dimensional aortic flows are conducted to demonstrate the effectiveness and merit of the proposed method.

A new noninvasive and patient-specific hemodynamic index for the severity of renal stenosis and outcome of interventional treatment

Renal arterial stenosis (RAS) often causes renovascular hypertension, which may result in kidney failure and life-threatening consequences. Direct assessment of the hemodynamic severity of RAS has yet to be addressed. In this work, we present a computational concept to derive a new, noninvasive, and patient-specific index to assess the hemodynamic severity of RAS and predict the potential benefit to the patient from a stenting therapy. The hemodynamic index is derived from a functional relation between the translesional pressure indicator (TPI) and lumen volume reduction (S) through a parametric deterioration of the RAS. Our in-house computational platform, InVascular, for image-based computational hemodynamics is used to compute the TPI at given S. InVascular integrates unified computational modeling for both image processing and computational hemodynamics with graphic processing unit parallel computing technology. The TPI-S curve reveals a pair of thresholds of S indicating mild or severe RAS. The TPI at S = 0 represents the pressure improvement following a successful stenting therapy. Six patient cases with a total of 6 aortic and 12 renal arteries are studied. The computed blood pressure waveforms have good agreements with the in vivo measured ones and the systolic pressure is statistical equivalence to the in-vivo measurements with p < .001. Uncertainty quantification provides the reliability of the computed pressure through the corresponding 95% confidence interval. The severity assessments of RAS in four cases are consistent with the medical practice. The preliminary results inspire a more sophisticated investigation for real medical insights of the new index. This computational concept can be applied to other arterial stenoses such as iliac stenosis. Such a noninvasive and patient-specific hemodynamic index has the potential to aid in the clinical decision-making of interventional treatment with reduced medical cost and patient risks.

Surrogate-based Bayesian calibration of biomechanical models with isotropic material behavior

This work introduces a computational methodology to calibrate material models in biomechanical applications under uncertainty. We adopt a Bayesian approach, which estimates the probability distributions of hyperelastic material parameters, based on force-strain measurements. We approximate the parametric biomechanical model by combining a reduced order representation of the force response with a Polynomial Chaos expansion. The surrogate model allows to employ sampling-intensive Markov chain Monte Carlo methods and provides an efficient way to estimate (generalized) Sobol coefficients. We use a Sobol sensitivity analysis to assess the influence of material parameters and present an iterative procedure to quantify the accuracy of the surrogate model as additional uncertainty during Bayesian updating. The methodology is illustrated with three cases, tensile experiments on heat-induced whey protein gel, indentation experiments for oocytes and a manufactured example. Real experimental data are used for the calibration.

Geometric Uncertainty in Patient-Specific Cardiovascular Modeling with Convolutional Dropout Networks

We propose a novel approach to generate samples from the conditional distribution of patient-specific cardiovascular models given a clinically aquired image volume. A convolutional neural network architecture with dropout layers is first trained for vessel lumen segmentation using a regression approach, to enable Bayesian estimation of vessel lumen surfaces. This network is then integrated into a path-planning patient-specific modeling pipeline to generate families of cardiovascular models. We demonstrate our approach by quantifying the effect of geometric uncertainty on the hemodynamics for three patient-specific anatomies, an aorto-iliac bifurcation, an abdominal aortic aneurysm and a sub-model of the left coronary arteries. A key innovation introduced in the proposed approach is the ability to learn geometric uncertainty directly from training data. The results show how geometric uncertainty produces coefficients of variation comparable to or larger than other sources of uncertainty for wall shear stress and velocity magnitude, but has limited impact on pressure. Specifically, this is true for anatomies characterized by small vessel sizes, and for local vessel lesions seen infrequently during network training.

Uncertainty quantification implementations in human hemodynamic flows

BACKGROUND AND OBJECTIVE

Human hemodynamic modeling is usually influenced by uncertainties occurring from a considerable unavailability of information linked to the boundary conditions and the physical properties used in the numerical models. Calculating the effect of these uncertainties on the numerical findings along the cardiovascular system is a demanding process due to the complexity of the morphology of the body and the area dynamics. To cope with all these difficulties, Uncertainty Quantification (UQ) methods seem to be an ideal tool.

RESULTS

This study focuses on analyzing and summarizing some of the recent research efforts and directions of implementing UQ in human hemodynamic flows by analyzing 139 research papers. Initially, the suitability of applying this approach is analyzed and demonstrated. Then, an overview of the most significant research work in various fields of biomedical hemodynamic engineering is presented. Finally, it is attempted to identify any possible forthcoming directions for research and methodological progress of UQ in biomedical sciences.

CONCLUSION

This review concludes that by finding the best statistical methods and parameters to represent the propagated uncertainties, while achieving a good interpretation of the interaction between input-output, is crucial for implementing UQ in biomedical sciences.

Multiscale modeling meets machine learning: What can we learn?

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learning to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

Computational Evaluation of Venous Graft Geometries in Coronary Artery Bypass Surgery

Cardiothoracic surgeons are faced with a choice of different revascularization techniques and diameters for saphenous vein grafts (SVG) in coronary artery bypass graft surgery . Using computational simulations, we virtually investigate the effect of SVG geometry on hemodynamics of both venous grafts and the target coronary arteries. We generated patient-specific 3-dimensional anatomic models of coronary artery bypass graft surgery patients and quantified mechanical stimuli. We performed virtual surgery on 3 patient-specific models by modifying the geometry vein grafts to reflect single, Y, and sequential surgical configurations with SVG diameters ranging from 2 mm to 5 mm. Our study demonstrates that the coronary artery runoffs are relatively insensitive to the choice of SVG revascularization geometry. We observe a 10% increase in runoff when the SVG diameter is changed from 2 mm to 5 mm. The wall shear stress of SVG increases dramatically when the diameter drops, following an inverse power scaling with diameter. For a fixed diameter, the average wall shear stress on the vein graft varies in ascending order as single, Y, and sequential graft in the patient cohort. The runoff to the target coronary arteries changes marginally due to the choice of graft configuration or diameter. The shear stress on the vein graft depends on both flow rate and diameter and follows an inverse power scaling consistent with a Poiseuille flow assumption. Given the similarity in runoff with different surgical configurations, choices of SVG geometries can be informed by propensity for graft failure using shear stress evaluations.

Global Sensitivity Analysis for Patient-Specific Aortic Simulations: The Role of Geometry, Boundary Condition and Large Eddy Simulation Modeling Parameters

Numerical simulations for computational hemodynamics in clinical settings require a combination of many ingredients, mathematical models, solvers and patient-specific data. The sensitivity of the solutions to these factors may be critical, particularly when we have a partial or noisy knowledge of data. Uncertainty quantification is crucial to assess the reliability of the results. We present here an extensive sensitivity analysis in aortic flow simulations, to quantify the dependence of clinically relevant quantities to the patient-specific geometry and the inflow boundary conditions. Geometry and inflow conditions are generally believed to have a major impact on numerical simulations. We resort to a global sensitivity analysis, (i.e., not restricted to a linearization around a working point), based on polynomial chaos expansion (PCE) and the associated Sobol' indices. We regard the geometry and the inflow conditions as the realization of a parametric stochastic process. To construct a physically consistent stochastic process for the geometry, we use a set of longitudinal-in-time images of a patient with an abdominal aortic aneurysm (AAA) to parametrize geometrical variations. Aortic flow is highly disturbed during systole. This leads to high computational costs, even amplified in a sensitivity analysis -when many simulations are needed. To mitigate this, we consider here a large Eddy simulation (LES) model. Our model depends in particular on a user-defined parameter called filter radius. We borrowed the tools of the global sensitivity analysis to assess the sensitivity of the solution to this parameter too. The targeted quantities of interest (QoI) include: the total kinetic energy (TKE), the time-average wall shear stress (TAWSS), and the oscillatory shear index (OSI). The results show that these indexes are mostly sensitive to the geometry. Also, we find that the sensitivity may be different during different instants of the heartbeat and in different regions of the domain of interest. This analysis helps to assess the reliability of in silico tools for clinical applications.

Exploring the potential of transfer learning for metamodels of heterogeneous material deformation

From the nano-scale to the macro-scale, biological tissue is spatially heterogeneous. Even when tissue behavior is well understood, the exact subject specific spatial distribution of material properties is often unknown. And, when developing computational models of biological tissue, it is usually prohibitively computationally expensive to simulate every plausible spatial distribution of material properties for each problem of interest. Therefore, one of the major challenges in developing accurate computational models of biological tissue is capturing the potential effects of this spatial heterogeneity. Recently, machine learning based metamodels have gained popularity as a computationally tractable way to overcome this problem because they can make predictions based on a limited number of direct simulation runs. These metamodels are promising, but they often still require a high number of direct simulations to achieve an acceptable performance. Here we show that transfer learning, a strategy where knowledge gained while solving one problem is transferred to solving a different but related problem, can help overcome this limitation. Critically, transfer learning can be used to leverage both low-fidelity simulation data and simulation data that is the outcome of solving a different but related mechanical problem. In this paper, we extend Mechanical MNIST, our open source benchmark dataset of heterogeneous material undergoing large deformation, to include a selection of low-fidelity simulation results that require ≈ 2 - 4 orders of magnitude less CPU time to run. Then, we show that transferring the knowledge stored in metamodels trained on these low-fidelity simulation results can vastly improve the performance of metamodels used to predict the results of high-fidelity simulations. In the most dramatic examples, metamodels trained on 100 high fidelity simulations but pre-trained on 60,000 low-fidelity simulations achieves nearly the same test error as metamodels trained on 60,000 high-fidelity simulations (1 - 1.5% mean absolute percent error). In addition, we show that transfer learning is an effective method for leveraging data from different load cases, and for leveraging low-fidelity two-dimensional simulations to predict the outcomes of high-fidelity three-dimensional simulations. Looking forward, we anticipate that transfer learning will enable us to better capture the influence of tissue spatial heterogeneity on the mechanical behavior of biological materials across multiple different domains.

Validation of Wall Shear Stress Assessment in Non-invasive Coronary CTA versus Invasive Imaging: A Patient-Specific Computational Study

Endothelial shear stress (ESS) identifies coronary plaques at high risk for progression and/or rupture leading to a future acute coronary syndrome. In this study an optimized methodology was developed to derive ESS, pressure drop and oscillatory shear index using computational fluid dynamics (CFD) in 3D models of coronary arteries derived from non-invasive coronary computed tomography angiography (CTA). These CTA-based ESS calculations were compared to the ESS calculations using the gold standard with fusion of invasive imaging and CTA. In 14 patients paired patient-specific CFD models based on invasive and non-invasive imaging of the left anterior descending (LAD) coronary arteries were created. Ten patients were used to optimize the methodology, and four patients to test this methodology. Time-averaged ESS (TAESS) was calculated for both coronary models applying patient-specific physiological data available at the time of imaging. For data analysis, each 3D reconstructed coronary artery was divided into 2 mm segments and each segment was subdivided into 8 arcs (45°).TAESS and other hemodynamic parameters were averaged per segment as well as per arc. Furthermore, the paired segment- and arc-averaged TAESS were categorized into patient-specific tertiles (low, medium and high). In the ten LADs, used for optimization of the methodology, we found high correlations between invasively-derived and non-invasively-derived TAESS averaged over segments (n = 263, r = 0.86) as well as arcs (n = 2104, r = 0.85, p < 0.001). The correlation was also strong in the four testing-patients with r = 0.95 (n = 117 segments, p = 0.001) and r = 0.93 (n = 936 arcs, p = 0.001).There was an overall high concordance of 78% of the three TAESS categories comparing both methodologies using the segment- and 76% for the arc-averages in the first ten patients. This concordance was lower in the four testing patients (64 and 64% in segment- and arc-averaged TAESS). Although the correlation and concordance were high for both patient groups, the absolute TAESS values averaged per segment and arc were overestimated using non-invasive vs. invasive imaging [testing patients: TAESS segment: 30.1(17.1-83.8) vs. 15.8(8.8-63.4) and TAESS arc: 29.4(16.2-74.7) vs 15.0(8.9-57.4) p < 0.001]. We showed that our methodology can accurately assess the TAESS distribution non-invasively from CTA and demonstrated a good correlation with TAESS calculated using IVUS/OCT 3D reconstructed models.

High precision invasive FFR, low-cost invasive iFR, or non-invasive CFR?: optimum assessment of coronary artery stenosis based on the patient-specific computational models

The objective of this paper is to apply computational fluid dynamic (CFD) as a complementary tool for clinical tests to not only predict the present and future status of left coronary artery stenosis but also to evaluate some clinical hypotheses. In order to assess the present status of the coronary artery stenosis severity, and thereby selecting the most appropriate type of treatment for each patient, fractional flow reserve (FFR), instantaneous wave free-ratio (iFR), and coronary flow reserve (CFR) are calculated. To examine FFR, iFR, and CFR results, the effect of geometric features of stenoses, including diameter reduction (%), lesion length (LL), and minimum lumen diameter (MLD), is studied on them. It is observed that FFR is a more conservative index than iFR and CFR to assess the severity of coronary stenosis. In addition, it is seen that FFR, iFR, and CFR decrease by increasing LL and decreasing MLD. Therefore, the morphological indices, LL/MLD and LL/MLD̂4, with the calculated conservative cut-off values equal to 5.5 and 3.6, are considered. Next, some controversial clinical hypotheses about the assessment of the severity of coronary stenosis are evaluated numerically. These include the examination of FFR, iFR, and CFR accuracies, investigating the effect of coronary hyperemia on iFR, as well as the reliability of the hybrid iFR-FFR decision-making strategy. The presented numerical model can also be used as a predictive tool to identify the atherosuseptible sites of arteries by calculating the time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

The effects of clinically-derived parametric data uncertainty in patient-specific coronary simulations with deformable walls

Cardiovascular simulations are increasingly used for noninvasive diagnosis of cardiovascular disease, to guide treatment decisions, and in the design of medical devices. Quantitative assessment of the variability of simulation outputs due to input uncertainty is a key step toward further integration of cardiovascular simulations in the clinical workflow. In this study, we present uncertainty quantification in computational models of the coronary circulation to investigate the effect of uncertain parameters, including coronary pressure waveform, intramyocardial pressure, morphometry exponent, and the vascular wall Young's modulus. We employ a left coronary artery model with deformable vessel walls, simulated via an Arbitrary-Lagrangian-Eulerian framework for fluid-structure interaction, with a prescribed inlet pressure and open-loop lumped parameter network outlet boundary conditions. Stochastic modeling of the uncertain inputs is determined from intra-coronary catheterization data or gathered from the literature. Uncertainty propagation is performed using several approaches including Monte Carlo, Quasi Monte Carlo sampling, stochastic collocation, and multi-wavelet stochastic expansion. Variabilities in the quantities of interest, including branch pressure, flow, wall shear stress, and wall deformation are assessed. We find that uncertainty in inlet pressures and intramyocardial pressures significantly affect all resulting QoIs, while uncertainty in elastic modulus only affects the mechanical response of the vascular wall. Variability in the morphometry exponent used to distribute the total downstream vascular resistance to the single outlets, has little effect on coronary hemodynamics or wall mechanics. Finally, we compare convergence behaviors of statistics of QoIs using several uncertainty propagation methods on three model benchmark problems and the left coronary simulations. From the simulation results, we conclude that the multi-wavelet stochastic expansion shows superior accuracy and performance against Quasi Monte Carlo and stochastic collocation methods.

Multilevel and multifidelity uncertainty quantification for cardiovascular hemodynamics

Standard approaches for uncertainty quantification in cardiovascular modeling pose challenges due to the large number of uncertain inputs and the significant computational cost of realistic three-dimensional simulations. We propose an efficient uncertainty quantification framework utilizing a multilevel multifidelity Monte Carlo (MLMF) estimator to improve the accuracy of hemodynamic quantities of interest while maintaining reasonable computational cost. This is achieved by leveraging three cardiovascular model fidelities, each with varying spatial resolution to rigorously quantify the variability in hemodynamic outputs. We employ two low-fidelity models (zero- and one-dimensional) to construct several different estimators. Our goal is to investigate and compare the efficiency of estimators built from combinations of these two low-fidelity model alternatives and our high-fidelity three-dimensional models. We demonstrate this framework on healthy and diseased models of aortic and coronary anatomy, including uncertainties in material property and boundary condition parameters. Our goal is to demonstrate that for this application it is possible to accelerate the convergence of the estimators by utilizing a MLMF paradigm. Therefore, we compare our approach to single fidelity Monte Carlo estimators and to a multilevel Monte Carlo approach based only on three-dimensional simulations, but leveraging multiple spatial resolutions. We demonstrate significant, on the order of 10 to 100 times, reduction in total computational cost with the MLMF estimators. We also examine the differing properties of the MLMF estimators in healthy versus diseased models, as well as global versus local quantities of interest. As expected, global quantities such as outlet pressure and flow show larger reductions than local quantities, such as those relating to wall shear stress, as the latter rely more heavily on the highest fidelity model evaluations. Similarly, healthy models show larger reductions than diseased models. In all cases, our workflow coupling Dakota's MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.

Effect of Wall Elasticity on Hemodynamics and Wall Shear Stress in Patient-Specific Simulations in the Coronary Arteries

Wall shear stress (WSS) has been shown to be associated with myocardial infarction (MI) and progression of atherosclerosis. Wall elasticity is an important feature of hemodynamic modeling affecting WSS calculations. The objective of this study was to investigate the role of wall elasticity on WSS, and justify use of either rigid or elastic models in future studies. Digital anatomic models of the aorta and coronaries were created based on coronary computed tomography angiography (CCTA) in four patients. Hemodynamics was computed in rigid and elastic models using a finite element flow solver. WSS in five timepoints in the cardiac cycle and time averaged wall shear stress (TAWSS) were compared between the models at each 3 mm subsegment and 4 arcs in cross sections along the centerlines of coronaries. In the left main (LM), proximal left anterior descending (LAD), left circumflex (LCX), and proximal right coronary artery (RCA) of the elastic model, the mean percent radial increase 5.95 ± 1.25, 4.02 ± 0.97, 4.08 ± 0.94, and 4.84 ± 1.05%, respectively. WSS at each timepoint in the cardiac cycle had slightly different values; however, when averaged over the cardiac cycle, there were negligible differences between the models. In both the subsegments (n = 704) and subarc analysis, TAWSS in the two models were highly correlated (r = 0.99). In investigation on the effect of coronary wall elasticity on WSS in CCTA-based models, the results of this study show no significant differences in TAWSS justifying using rigid wall models for future larger studies.

The effect of anastomotic angle and diameter ratio on flow field in the distal end-to-side anastomosis

Flow fields in the distal end-to-side anastomosis of coronary artery bypass graft are associated with intimal hyperplasia and bypass failure. This work aims to demonstrate the effect of anastomotic angle and diameter ratio on flow field of coronary artery bypass graft. The flow fields inside polydimethylsiloxane models of coronary artery bypass graft with various anastomotic angles (α = 30°, 45°, 60° and 75°) and different diameter ratios (Φ = 0.78 and 1.11) are investigated using particle image velocimetry and computational fluid dynamics method under pulsatile flow condition. The results show that the anastomotic angle is positively correlated with the number and area of the recirculation zone, and the flow field disturbance at the anastomosis will develop in the same direction. Compared with that of Φ = 0.78, when Φ = 1.11, the flow fields at the anastomosis are relatively smoother with less turbulence.