The impact of scaled boundary conditions on wall shear stress computations in atherosclerotic human coronary bifurcations

Main Coronary Flow Calculation With the Assistance of Physiological Side Branch Flow

OBJECTIVE

Fractional flow reserve (FFR) derived from coronary angiography, referred to as ICA-FFR, is a less-invasive alternative for invasive FFR measurement based on computational fluid dynamics. Blood flow into side branches influences the accuracy of ICA-FFR. However, properly compensating for side branch flow in ICA-FFR analysis is challenging. In this study, we proposed a physiological side branch flow model to comprehensively compensate side branch flow for ICA-FFR analysis with no need for reconstructing side branch geometry.

METHODOLOGY

the physiological side branch flow model employed a reduced-order model to calculate the pressure distribution in vessel segments. The main coronary artery (without side branches) was delineated and divided based on bifurcation nodes. The model compensates for flow to invisible side branches within each segment and flow to visible side branches at each bifurcation node. Lastly, ICA-FFR based on physiological side branch flow model (ICA-FFR) was calculated from a single angiographic view. Functional stenosis is defined by FFR 0.80.

RESULT

Our study involved 223 vessels from 172 patients. Using invasive FFR as a reference, the Pearson correlation coefficient of ICA-FFR was 0.93. ICA-FFR showed a high AUC (AUC = 0.96) and accuracy (91.9) in predicting functional stenosis.

CONCLUSION

The proposed model accurately compensates for flow to side branches without their geometry in ICA-FFR analysis. ICA-FFR analysis exhibits high feasibility and diagnostic performance in identifying functional stenosis.

A new understanding of coronary curvature and haemodynamic impact on the course of plaque onset and progression

The strong link between atherosclerosis and luminal biomechanical stresses is well established. Yet, this understanding has not translated into preventative coronary diagnostic imaging, particularly due to the under-explored role of coronary anatomy and haemodynamics in plaque onset, which we aim to address with this work. The left coronary trees of 20 non-stenosed (%diameter stenosis [%DS] = 0), 12 moderately (0 < %DS < 70) and 7 severely (%DS ≥ 70) stenosed cases were dissected into bifurcating and non-bifurcating segments for whole-tree and segment-specific comparisons, correlating nine three-dimensional coronary anatomical features, topological shear variation index (TSVI) and luminal areas subject to low time-average endothelial shear stress (%LowTAESS), high oscillatory shear index (%HighOSI) and high relative residence time (%HighRRT). We found that TSVI is the only metric consistently differing between non-stenosed and stenosed cases across the whole tree, bifurcating and non-bifurcating segments (p < 0.002, AUC = 0.876), whereas average curvature and %HighOSI differed only for the whole trees (p < 0.024) and non-bifurcating segments (p < 0.027), with AUC > 0.711. Coronary trees with moderate or severe stenoses differed only in %LowTAESS (p = 0.009) and %HighRRT (p = 0.012). This suggests TSVI, curvature and %HighOSI are potential factors driving plaque onset, with greater predictive performance than the previously recognized %LowTAESS and %HighRRT, which appears to play a role in plaque progression.

Integrating Computational and Biological Hemodynamic Approaches to Improve Modeling of Atherosclerotic Arteries

Atherosclerosis is the primary cause of cardiovascular disease, resulting in mortality, elevated healthcare costs, diminished productivity, and reduced quality of life for individuals and their communities. This is exacerbated by the limited understanding of its underlying causes and limitations in current therapeutic interventions, highlighting the need for sophisticated models of atherosclerosis. This review critically evaluates the computational and biological models of atherosclerosis, focusing on the study of hemodynamics in atherosclerotic coronary arteries. Computational models account for the geometrical complexities and hemodynamics of the blood vessels and stenoses, but they fail to capture the complex biological processes involved in atherosclerosis. Different in vitro and in vivo biological models can capture aspects of the biological complexity of healthy and stenosed vessels, but rarely mimic the human anatomy and physiological hemodynamics, and require significantly more time, cost, and resources. Therefore, emerging strategies are examined that integrate computational and biological models, and the potential of advances in imaging, biofabrication, and machine learning is explored in developing more effective models of atherosclerosis.

Evaluation of models of sequestration flow in coronary arteries-Physiology versus anatomy?

BACKGROUND

Myocardial ischaemia results from insufficient coronary blood flow. Computed virtual fractional flow reserve (vFFR) allows quantification of proportional flow loss without the need for invasive pressure-wire testing. In the current study, we describe a novel, conductivity model of side branch flow, referred to as 'leak'. This leak model is a function of taper and local pressure, the latter of which may change radically when focal disease is present. This builds upon previous techniques, which either ignore side branch flow, or rely purely on anatomical factors. This study aimed to describe a new, conductivity model of side branch flow and compare this with established anatomical models.

METHODS AND RESULTS

The novel technique was used to quantify vFFR, distal absolute flow (Qd) and microvascular resistance (CMVR) in 325 idealised 1D models of coronary arteries, modelled from invasive clinical data. Outputs were compared to an established anatomical model of flow. The conductivity model correlated and agreed with the reference model for vFFR (r = 0.895, p < 0.0001; ＋0.02, 95% CI 0.00 to ＋ 0.22), Qd (r = 0.959, p < 0.0001; -5.2 mL/min, 95% CI -52.2 to ＋13.0) and CMVR (r = 0.624, p < 0.0001; ＋50 Woods Units, 95% CI -325 to ＋2549).

CONCLUSION

Agreement between the two techniques was closest for vFFR, with greater proportional differences seen for Qd and CMVR. The conductivity function assumes vessel taper was optimised for the healthy state and that CMVR was not affected by local disease. The latter may be addressed with further refinement of the technique or inferred from complementary image data. The conductivity technique may represent a refinement of current techniques for modelling coronary side-branch flow. Further work is needed to validate the technique against invasive clinical data.

Turbulent blood dynamics in the left heart in the presence of mitral regurgitation: a computational study based on multi-series cine-MRI

In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution.

A spatiotemporal analysis of the left coronary artery biomechanics using fluid-structure interaction models

Biomechanics plays a critical role in coronary artery disease development. FSI simulation is commonly used to understand the hemodynamics and mechanical environment associated with atherosclerosis pathology. To provide a comprehensive characterization of patient-specific coronary biomechanics, an analysis of FSI simulation in the spatial and temporal domains was performed. In the current study, a three-dimensional FSI model of the LAD coronary artery was built based on a patient-specific geometry using COMSOL Multiphysics. The effect of myocardial bridging was simulated. Wall shear stress and its derivatives including time-averaged wall shear stress, wall shear stress gradient, and OSI were calculated across the cardiac cycle in multiple locations. Arterial wall strain (radial, circumferential, and longitudinal) and von Mises stress were calculated. To assess perfusion, vFFR was calculated. The results demonstrated the FSI model could identify regional and transient differences in biomechanical parameters within the coronary artery. The addition of myocardial bridging caused a notable change in von Mises stress and an increase in arterial strain during systole. The analysis performed in this manner takes greater advantage of the information provided in the space and time domains and can potentially assist clinical evaluation.

Modelling coronary flows: impact of differently measured inflow boundary conditions on vessel-specific computational hemodynamic profiles

BACKGROUND AND OBJECTIVES

The translation of hemodynamic quantities based on wall shear stress (WSS) or intravascular helical flow into clinical biomarkers of coronary atherosclerotic disease is still hampered by the assumptions/idealizations required by the computational fluid dynamics (CFD) simulations of the coronary hemodynamics. In the resulting budget of uncertainty, inflow boundary conditions (BCs) play a primary role. Accordingly, in this study we investigated the impact of the approach adopted for in vivo coronary artery blood flow rate assessment on personalized CFD simulations where blood flow rate is used as inflow BC.

METHODS

CFD simulations were carried out on coronary angiograms by applying personalized inflow BCs derived from four different techniques assessing in vivo surrogates of flow rate: continuous thermodilution, intravascular Doppler, frame count-based 3D contrast velocity, and diameter-based scaling law. The impact of inflow BCs on coronary hemodynamics was evaluated in terms of WSS- and helicity-based quantities.

RESULTS

As main findings, we report that: (i) coronary flow rate values may differ based on the applied flow derivation technique, as continuous thermodilution provided higher flow rate values than intravascular Doppler and diameter-based scaling law (p = 0.0014 and p = 0.0023, respectively); (ii) such intrasubject differences in flow rate values lead to different surface-averaged values of WSS magnitude and helical blood flow intensity (p<0.0020); (iii) luminal surface areas exposed to low WSS and helical flow topological features showed robustness to the flow rate values.

CONCLUSIONS

Although the absence of a clinically applicable gold standard approach prevents a general recommendation for one coronary blood flow rate derivation technique, our findings indicate that the inflow BC may impact computational hemodynamic results, suggesting that a standardization would be desirable to provide comparable results among personalized CFD simulations of the coronary hemodynamics.

Endothelial shear stress and vascular remodeling in bioresorbable scaffold and metallic stent

BACKGROUND AND AIMS

The impact of endothelial shear stress (ESS) on vessel remodeling in vessels implanted with bioresorbable scaffold (BRS) as compared to metallic drug-eluting stent (DES) remains elusive. The aim of this study was to determine whether the relationship between ESS and remodeling patterns differs in BRS from those seen in metallic DES at 3-year follow-up.

METHODS

In the ABSORB II randomized trial, lesions were investigated by serial coronary angiography and intravascular ultrasound (IVUS). Three-dimensional reconstructions of coronary arteries post-procedure and at 3 years were performed. ESS was quantified using non-Newtonian steady flow simulation. IVUS cross-sections in device segment were matched using identical landmarks.

RESULTS

Paired ESS calculations post-procedure and at 3 years were feasible in 57 lesions in 56 patients. Post-procedure, median ESS at frame level was higher in BRS than in DES, with marginal statistical significance (0.97 ± 0.48 vs. 0.75 ± 0.39 Pa, p = 0.063). In the BRS arm, vessel area and lumen area showed larger increases in the highest tercile of median ESS post-procedure as compared to the lowest tercile. In contrast, in DES, no significant relationship between median ESS post-procedure and remodeling was observed. In multivariate analysis, smaller vessel area, larger lumen area, higher plaque burden post-procedure, and higher median ESS post-procedure were independently associated with expansive remodeling in matched frames. Only in BRS, younger age was an additional significant predictor of expansive remodeling.

CONCLUSIONS

In a subset of lesions with large plaque burden, shear stress could be associated with expansive remodeling and late lumen enlargement in BRS, while ESS had no impact on vessel dimension in metallic DES.

The importance of side branches in modeling 3D hemodynamics from angiograms for patients with coronary artery disease

Genesis of atherosclerotic lesions in the human arterial system is critically influenced by the fluid mechanics. Applying computational fluid dynamic tools based on accurate coronary physiology derived from conventional biplane angiogram data may be useful in guiding percutaneous coronary interventions. The primary objective of this study is to build and validate a computational framework for accurate personalized 3-dimensional hemodynamic simulation across the complete coronary arterial tree and demonstrate the influence of side branches on coronary hemodynamics by comparing shear stress between coronary models with and without these included. The proposed novel computational framework based on biplane angiography enables significant arterial circulation analysis. This study shows that models that take into account flow through all side branches are required for precise computation of shear stress and pressure gradient whereas models that have only a subset of side branches are inadequate for biomechanical studies as they may overestimate volumetric outflow and shear stress. This study extends the ongoing computational efforts and demonstrates that models based on accurate coronary physiology can improve overall fidelity of biomechanical studies to compute hemodynamic risk-factors.

Impaired Flow-Mediated Dilation and Severity and Vulnerability of Culprit Plaque in Patients with Coronary Artery Disease

The association between endothelial function, evaluated using flow-mediated dilatation (FMD), and the severity of coronary artery disease remains to be elucidated.A total of 245 consecutive patients with stable angina were prospectively enrolled. FMD was evaluated in the brachial artery before percutaneous coronary intervention. Patients were divided into 2 groups according to the FMD value (lower FMD group [FMD < 2.0], n = 82; higher FMD group [FMD ≥ 2.0], n = 163). The severity of coronary artery disease was evaluated using findings of angiography and optical coherence tomography, and compared between the 2 groups.The prevalence of left main (LM) disease was significantly higher in the lower FMD group than in the higher FMD group (8.5% versus 2.5%, P = 0.046), although the prevalence of multivessel disease was comparable between the groups. Lower FMD was independently associated with a higher prevalence of LM disease (odds ratio, 3.89; 95% confidence interval, 1.12-15.5; P = 0.033). A general linear model with multiple variables revealed that the minimal lumen area (MLA) in the culprit lesion was significantly smaller in patients with lower FMD than in those with higher FMD (regression coefficient b, -0.249 mm²; 95% confidence interval, -0.479--0.018 mm²; P = 0.035). The prevalence ofvulnerable plaque characteristics was comparable between the 2 groups.Patients with lower FMD had a higher incidence of LM disease and a smaller MLA in the culprit lesion. FMD may be a useful, noninvasive indicator for identifying patients with severe coronary artery disease.

Coronary fractional flow reserve measurements of a stenosed side branch: a computational study investigating the influence of the bifurcation angle

Background

Coronary hemodynamics and physiology specific for bifurcation lesions was not well understood. To investigate the influence of the bifurcation angle on the intracoronary hemodynamics of side branch (SB) lesions computational fluid dynamics simulations were performed.

Methods

A parametric model representing a left anterior descending—first diagonal coronary bifurcation lesion was created according to the literature. Diameters obeyed fractal branching laws. Proximal and distal main branch (DMB) stenoses were both set at 60 %. We varied the distal bifurcation angles (40°, 55°, and 70°), the flow splits to the DMB and SB (55 %:45 %, 65 %:35 %, and 75 %:25 %), and the SB stenoses (40, 60, and 80 %), resulting in 27 simulations. Fractional flow reserve, defined as the ratio between the mean distal stenosis and mean aortic pressure during maximal hyperemia, was calculated for the DMB and SB (FFR_SB) for all simulations.

Results

The largest differences in FFR_SB comparing the largest and smallest bifurcation angles were 0.02 (in cases with 40 % SB stenosis, irrespective of the assumed flow split) and 0.05 (in cases with 60 % SB stenosis, flow split 55 %:45 %). When the SB stenosis was 80 %, the difference in FFR_SB between the largest and smallest bifurcation angle was 0.33 (flow split 55 %:45 %). By describing the ΔP_SB−Q_SB relationship using a quadratic curve for cases with 80 % SB stenosis, we found that the curve was steeper (i.e. higher flow resistance) when bifurcation angle increases (ΔP = 0.451*Q + 0.010*Q² and ΔP = 0.687*Q + 0.017*Q² for 40° and 70° bifurcation angle, respectively). Our analyses revealed complex hemodynamics in all cases with evident counter-rotating helical flow structures. Larger bifurcation angles resulted in more pronounced helical flow structures (i.e. higher helicity intensity), when 60 or 80 % SB stenoses were present. A good correlation (R² = 0.80) between the SB pressure drop and helicity intensity was also found.

Conclusions

Our analyses showed that, in bifurcation lesions with 60 % MB stenosis and 80 % SB stenosis, SB pressure drop is higher for larger bifurcation angles suggesting higher flow resistance (i.e. curves describing the ΔP_SB−Q_SB relationship being steeper). When the SB stenosis is mild (40 %) or moderate (60 %), SB resistance is minimally influenced by the bifurcation angle, with differences not being clinically meaningful. Our findings also highlighted the complex interplay between anatomy, pressure drops, and blood flow helicity in bifurcations.