Low Time-averaged Wall Shear Stress Research Articles

Modelling lower-limb peripheral arterial disease using clinically available datasets: impact of inflow boundary conditions on hemodynamic indices for restenosis prediction

Background and ObjectivesThe integration of hemodynamic markers as risk factors in restenosis prediction models for lower-limb peripheral arteries is hindered by fragmented clinical datasets. Computed tomography (CT) scans enable vessel geometry reconstruction and can be obtained at different times than the Doppler ultrasound (DUS) images, which provide information on blood flow velocity. Computational fluid dynamics (CFD) simulations allow the computation of near-wall hemodynamic indices, whose accuracy depends on the prescribed inlet boundary condition (BC), derived from the DUS images. This study aims to: (i) investigate the impact of different DUS-derived velocity waveforms on CFD results; (ii) test whether the same vessel areas, subjected to altered hemodynamics, can be detected independently of the applied inlet BC; (iii) suggest suitable DUS images to obtain reliable CFD results. MethodsCFD simulations were conducted on three patients treated with bypass surgery, using patient-specific DUS-derived inlet BCs recorded at either the same or different time points than the CT scan. The impact of the chosen inflow condition on bypass hemodynamics was assessed in terms of wall shear stress (WSS)-derived quantities. Patient-specific critical thresholds for the hemodynamic indices were applied to identify critical luminal areas and compare the results with a reference obtained with a DUS image acquired in close temporal proximity to the CT scan. ResultsThe main findings indicate that: (i) DUS-derived inlet velocity waveforms acquired at different time points than the CT scan led to statistically significantly different CFD results (p<0.001); (ii) the same luminal surface areas, exposed to low time-averaged WSS, could be identified independently of the applied inlet BCs; (iii) similar outcomes were observed for the other hemodynamic indices if the prescribed inlet velocity waveform had the same shape and comparable systolic acceleration time to the one recorded in close temporal proximity to the CT scan. ConclusionsDespite a lack of standardised data collection for diseased lower-limb peripheral arteries, an accurate estimation of luminal areas subjected to altered near-wall hemodynamics is possible independently of the applied inlet BC. This holds if the applied inlet waveform shares some characteristics – derivable from the DUS report – as one matching the acquisition time of the CT scan.

Assessing the impact of tear direction in coronary artery dissection on thrombosis development: A hemodynamic computational study

ObjectiveIatrogenic coronary artery dissection is a complication of coronary intimal injury and dissection due to improper catheter manipulation. The impact of tear direction on the prognosis of coronary artery dissection (CAD) remains unclear. This study examines the hemodynamic effects of different tear directions (transverse and longitudinal) of CAD and evaluates the risk of thrombosis, rupture and further dilatation of CAD. MethodsTwo types of CAD models (Type I: transverse tear, Type II: longitudinal tear) were reconstructed from the aorto-coronary CTA dataset of 8 healthy cases. Four WSS-based indicators were analyzed, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and cross flow index (CFI). A thrombus growth model was also introduced to predict the trend of thrombus growth in CAD with two different tear directions. ResultsFor most of the WSS-based indicators, including TAWSS, RRT, and CFI, no statistically significant differences were observed across the CAD models with varying tear directions, except for OSI, where a significant difference was noted (p < 0.05). Meanwhile, in terms of thrombus growth, the thrombus growing at the tear of the Type I (transverse tear) CAD model extended into the true lumen earlier than that of the Type II (longitudinal tear) model. ConclusionsNumerical simulations suggest that: (1) The CAD with transverse tear have a high risk of further tearing of the dissection at the distal end of the tear. (2) The CAD with longitudinal tear create a hemodynamic environment characterized by low TAWSS and high OSI in the false lumen, which may additionally increase the risk of vessel wall injury. (3) The CAD with transverse tear may have a higher risk of thrombosis and coronary obstruction and myocardial ischemia in the early phase of the dissection.

Spatiotemporal analysis of the effects of exercise on the hemodynamics of the aorta in hypertensive rats using fluid-structure interaction simulation

Abstract Background and Objective Hemodynamic changes that lead to increased blood pressure represent the main drivers of organ damage in hypertension. Prolonged increases to blood pressure can lead to vascular remodeling, which also affects vascular hemodynamics during the pathogenesis of hypertension. Exercise is beneficial for relieving hypertension, however the mechanistic link between exercise training and how it influences hemodynamics in the context of hypertension is not well understood. Methods n exercise model was developed using spontaneously hypertensive rats (SHR) subject to a 12-week treadmill training regime. The heart rates and blood pressures of rats were measured using the tail cuff method, while micro-computed tomography (CT) scanning was used to develop three-dimensional structures of rat aorta, and ultrasound was used to detect rat aortic blood flow and changes to vessel wall structures. Computational fluid dynamics (CFD) and fluid-structure interaction (FSI) models were used to simulate and measure hemodynamic parameters of the rat aortic vessels. In parallel, Masson staining was performed on fixed samples of blood vessels to investigate collagen volume fraction. Hypertensive rats in the sedentary and long-term exercise training groups were subjected to a single bout exercise training, and their aortic hemodynamic parameters were analyzed before, 5 min, 24 h, and 72 h after the single bout exercise. Results Of the two models, in comparison to actual ultrasonic measurement values recorded, we found that numerical simulation results from the FSI model could more accurately model blood flow in the ascending aorta of hypertensive rats, compared to the CFD model. Moreover, longterm exercise training improved local hemodynamic parameters of blood vessels, and led to improvements in adverse hemodynamic features documented, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT). Longterm exercise training of SHR also improved local vascular collagen deposition in the aorta, while improvements in vascular remodeling were also correlated with favorable hemodynamic parameters. Compared with sedentary SHR, signals for low TAWSS regions of the aortic arch in SHR on the long-term exercise regime shifted to the position of the ascending aorta after a single bout of exercise. Conclusions This study demonstrates that FSI is informative to study the spatiotemporal effects of long-term exercise training on hemodynamic changes within the aortas of hypertensive rats, and that long-term exercise is beneficial through its effects to modulate vascular hemodynamics in hypertension.

Surgical Modulation of Pulmonary Artery Shear Stress: A Patient-Specific CFD Analysis of the Norwood Procedure.

This study created 3D CFD models of the Norwood procedure for hypoplastic left heart syndrome (HLHS) using standard angiography and echocardiogram data to investigate the impact of shunt characteristics on pulmonary artery (PA) hemodynamics. Leveraging routine clinical data offers advantages such as availability and cost-effectiveness without subjecting patients to additional invasive procedures. Patient-specific geometries of the intrathoracic arteries of two Norwood patients were generated from biplane cineangiograms. "Virtual surgery" was then performed to simulate the hemodynamics of alternative PA shunt configurations, including shunt type (modified Blalock-Thomas-Taussig shunt (mBTTS) vs. right ventricle-to-pulmonary artery shunt (RVPAS)), shunt diameter, and pulmonary artery anastomosis angle. Left-right pulmonary flow differential, Qp/Qs, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) were evaluated. There was strong agreement between clinically measured data and CFD model output throughout the patient-specific models. Geometries with a RVPAS tended toward more balanced left-right pulmonary flow, lower Qp/Qs, and greater TAWSS and OSI than models with a mBTTS. For both shunt types, larger shunts resulted in a higher Qp/Qs and higher TAWSS, with minimal effect on OSI. Low TAWSS areas correlated with regions of low flow and changing the PA-shunt anastomosis angle to face toward low TAWSS regions increased TAWSS. Excellent correlation between clinically measured and CFD model data shows that 3D CFD models of HLHS Norwood can be developed using standard angiography and echocardiographic data. The CFD analysis also revealed consistent changes in PA TAWSS, flow differential, and OSI as a function of shunt characteristics.

Effects of preoperative aortic tortuosity and postoperative hypertension on patient-specific hemodynamics of abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) is a serious dilated vascular disease. The risk factors of aneurysm rupture and postoperative blood pressure are the major clinical concerns. The purpose of this work is to analyze the hemodynamic difference between preoperative and postoperative AAA using patient-specific boundary conditions and explore the rupture risk of AAA with different torsion and to simulate the blood flow of different degrees of hypertension. A three-element Windkessel model is utilized as the outlet boundary condition. 20-sim software and self-developed user-defined functions are used to calibrate the parameters. To analyze the influence of aortic tortuosity on hemodynamics, five AAA geometries with different torsion degrees are modified and obtained. We also grade the hypertension and explore the effect of hypertension after operation. The analysis results are consistent with the patient-specific situation. The results show that stent implantation for aneurysms reduces the pressure at the lesion site and increase the speed and wall shear stress. Vessels with high torsion have lower time-averaged wall shear stress (TAWSS), higher oscillatory shear stress index (OSI) and relative retention time, and more disordered blood flow. After stent implantation, with the increase in blood pressure, TAWSS and OSI continue to increase, and the blood flow speed is also faster. AAA with high torsion is more likely to rupture. In clinical practice, attention should be paid to the degree of torsion of the lesion site. Furthermore, blood pressure should be monitored and controlled in time to prevent postoperative complications.

Left main coronary artery morphological phenotypes and its hemodynamic properties.

Atherosclerosis may be linked to morphological defects that lead to variances in coronary artery hemodynamics. Few objective strategies exit at present for generalizing morphological phenotypes of coronary arteries in terms of hemodynamics. We used unsupervised clustering (UC) to classify the morphology of the left main coronary artery (LM) and looked at how hemodynamic distribution differed between phenotypes. In this study, 76 LMs were obtained from 76 patients. After LMs were reconstructed with coronary computed tomography angiography, centerlines were used to extract the geometric characteristics. Unsupervised clustering was carried out using these characteristics to identify distinct morphological phenotypes of LMs. The time-averaged wall shear stress (TAWSS) for each phenotype was investigated by means of computational fluid dynamics (CFD) analysis of the left coronary artery. We identified four clusters (i.e., four phenotypes): Cluster 1 had a shorter stem and thinner branches (n = 26); Cluster 2 had a larger bifurcation angle (n = 10); Cluster 3 had an ostium at an angulation to the coronary sinus and a more curved stem, and thick branches (n = 10); and Cluster 4 had an ostium at an angulation to the coronary sinus and a flatter stem (n = 14). TAWSS features varied widely across phenotypes. Nodes with low TAWSS (L-TAWSS) were typically found around the branching points of the left anterior descending artery (LAD), particularly in Cluster 2. Our findings demonstrated that UC is a powerful technique for morphologically classifying LMs. Different LM phenotypes exhibited distinct hemodynamic characteristics in certain regions. This morphological clustering method could aid in identifying people at high risk for developing coronary atherosclerosis, hence facilitating early intervention.

A hemodynamic analysis of energy loss in abdominal aortic aneurysm using three-dimension idealized model.

Objective: The aim of this study is to perform specific hemodynamic simulations of idealized abdominal aortic aneurysm (AAA) models with different diameters, curvatures and eccentricities and evaluate the risk of thrombosis and aneurysm rupture. Methods: Nine idealized AAA models with different diameters (3cm or 5cm), curvatures (0° or 30°) and eccentricities (centered on or tangent to the aorta), as well as a normal model, were constructed using commercial software (Solidworks; Dassault Systemes S.A, Suresnes, France). Hemodynamic simulations were conducted with the same time-varying volumetric flow rate extracted from the literature and 3-element Windkessel model (3 EWM) boundary conditions were applied at the aortic outlet. Several hemodynamic parameters such as time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), endothelial cell activation potential (ECAP) and energy loss (EL) were obtained to evaluate the risk of thrombosis and aneurysm rupture under different conditions. Results: Simulation results showed that the proportion of low TAWSS region and high OSI region increases with the rising of aneurysm diameter, whereas decreases in the curvature and eccentric models of the corresponding diameters, with the 5cm normal model having the largest low TAWSS region (68.5%) and high OSI region (40%). Similar to the results of TAWSS and OSI, the high ECAP and high RRT areas were largest in the 5cm normal model, with the highest wall-averaged value (RRT: 5.18 s, ECAP: 4.36 Pa-1). Differently, the increase of aneurysm diameter, curvature, and eccentricity all lead to the increase of mean flow EL and turbulent EL, such that the highest mean flow EL (0.82W · 10-3) and turbulent EL (1.72W · 10-3) were observed in the eccentric 5cm model with the bending angle of 30°. Conclusion: Collectively, increases in aneurysm diameter, curvature, and eccentricity all raise mean flow EL and turbulent flow EL, which may aggravate the damage and disturbance of flow in aneurysm. In addition, it can be inferred by conventional parameters (TAWSS, OSI, RRT and ECAP) that the increase of aneurysm diameter may raise the risk of thrombosis, whereas the curvature and eccentricity appeared to have a protective effect against thrombosis.

Computational analysis of patient-specific pulsatile blood flow: The influence of non-Newtonian models on wall shear stress assessment

Blood is a sophisticated biological fluid with components like erythrocytes that give it non-Newtonian behavior. Hemodynamic factors such as velocity magnitude, pressure, and wall shear stress descriptors are the most important factors in the development of atherosclerosis. The wall shear stress descriptors are regulated not only by flow geometry but also by blood rheological properties. In the current study, we carried out a numerical analysis of the non-Newtonian pulsatile blood flow while taking into account a patient-specific geometry and transient boundary conditions. Non-Newtonian blood flow is modeled using the four non-Newtonian models: the power-law model, the Carreau model, the Casson model, and the Quemada model, and compared with the Newtonian model. Streamline analysis vividly illustrates velocity patterns, revealing the presence of recirculation zones near sinus regions. The study suggests the significance of selecting appropriate viscosity models for accurate assessments, particularly in regions with low time-average wall shear stress values, such as those associated with atherosclerotic plaques. The differences in the time-averaged wall shear stress between the four non-Newtonian models were found to be the highest in the Quemada model. The study concluded that the non-Newtonian model is required when the focus is on the low-time-averaged wall shear stress area.

Hemodynamic influence of mild stenosis morphology in different coronary arteries: a computational fluid dynamic modelling study.

Introduction: Mild stenosis [degree of stenosis (DS) < 50%] is commonly labeled as nonobstructive lesion. Some lesions remain stable for several years, while others precipitate acute coronary syndromes (ACS) rapidly. The causes of ACS and the factors leading to diverse clinical outcomes remain unclear. Method: This study aimed to investigate the hemodynamic influence of mild stenosis morphologies in different coronary arteries. The stenoses were modeled with different morphologies based on a healthy individual data. Computational fluid dynamics analysis was used to obtain hemodynamic characteristics, including flow waveforms, fractional flow reserve (FFR), flow streamlines, time-average wall shear stress (TAWSS), and oscillatory shear index (OSI). Results: Numerical simulation indicated significant hemodynamic differences among different DS and locations. In the 20%-30% range, significant large, low-velocity vortexes resulted in low TAWSS (<4dyne/cm2) around stenoses. In the 30%-50% range, high flow velocity due to lumen area reduction resulted in high TAWSS (>40dyne/cm2), rapidly expanding the high TAWSS area (averagely increased by 0.46cm2) in left main artery and left anterior descending artery (LAD), where high OSI areas remained extensive (>0.19 cm2). Discussion: While mild stenosis does not pose any immediate ischemic risk due to a FFR > 0.95, 20%-50% stenosis requires attention and further subdivision based on location is essential. Rapid progression is a danger for lesions with 20%-30% DS near the stenoses and in the proximal LAD, while lesions with 30%-50% DS can cause plaque injury and rupture. These findings support clinical practice in early assessment, monitoring, and preventive treatment.

Effects of hematocrit levels on flow structures and stress levels in the healthy and diseased carotid arteries

Flow structures and wall shear stress (WSS) distribution in human arteries depend strongly on blood rheology. Blood viscosity undergoes significant alterations due to varying hematocrit (Hct) levels. Hematocrit levels can change from the normal level (45%) due to numerous disorders such as leukemia, lymphoma, and congenital heart disease. However, the effects of Hct levels on the flow features and WSS levels in large stenotic arteries are not well reported in literature. The present study computationally assesses the variation in flow features in healthy and stenosed carotid arteries due to changes in blood Hct levels. We have reported variations in axial and secondary flow characteristics, vortex structures, and WSS-based parameters for different Hct levels. Specifically, we consider high (65%) and low (25%) Hct levels in addition to a normal Hct level. Our results reveal that both high and low Hct levels have adverse effects on the hemodynamic features inside a healthy carotid artery model. The low Hct level displays a much more widespread low time-averaged WSS (TAWSS) area, while the high Hct shows an increased oscillatory shear index (OSI) zone at the carotid sinus wall. In stenotic arteries, low Hct enhances the spread of the low TAWSS zone in the post-stenotic carotid sinus, whereas high Hct increases the risk of plaque rupture. The high Hct level also increases the spread of high OSI zone at the inner wall of internal carotid artery for the moderately stenosed artery.

Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstructs flow and increases pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the small pulmonary arteries due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we applied a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) approach to investigate the impact of CTEPH on pulmonary artery stiffening, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) in extralobar (main, right, and left) pulmonary arteries and intralobar (distal to the extralobar) arteries. Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar and intralobar arteries. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI, quantified by the novel parameter φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document}, which is related to thrombogenicity. Our analysis reveals a strong positive correlation between increases in mean pulmonary artery pressure (mPAP) and φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document} from baseline to CTEPH in individual subjects, which supports the suggestion that increased φ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varphi$$\\end{document} drives disease severity. This subject-specific experimental–computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.

Hemodynamics analysis and computational fluid dynamics simulations to identify potential location of rupture in abdominal aortic aneurysm

Abstract Background Abdominal Aortic Aneurysm (AAA) rupture is the most severe consequence of aneurysmal disease, with an overall death rate of 80%. Early diagnosis and appropriate follow-up and management of AAA will significantly reduce the risk of rupture and mortality. Rupture risk assessment of AAA has a critical role for the early warning. Purpose To Predict of the course of an AAA and determine the rupture risk and location from a biomechanical point of view. AAA hemodynamics were investigated to clarify the relation between the hemodynamic parameters and the AAA rupture. Methods Computed tomography (CT) images of one healthy aorta and 23 patients’ aorta with aneurysms were analyzed. WSS-related hemodynamic parameters of time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell activation potential (ECAP), and relative residence time (RRT) are investigated in patient-specific AAA geometries by performing computational fluid dynamics (CFD) simulations (Figure 1). Patient-specific AAA geometries are identified using medical imaging techniques and realistic boundary conditions are determined using Doppler ultrasonography. Results AAA rupture locations predicted by CFD simulations were in good agreement with the actual rupture locations of retrospective data. The rupture locations had relatively low TAWSS, high OSI, high ECAP, and high RTT levels. High risk of rupture is expected in regions with a local minimum of TAWSS, and a local maximum of OSI, ECAP, and RRT. Conclusions Examining the spatial distribution of the hemodynamic parameters was a powerful approach for rupture risk assessment of AAA. Predicting the rupture location using a single hemodynamic parameter is not feasible, and the accuracy of AAA rupture risk assessment can be improved by using multiple WSS-related parameters by identifying the common intersection regions of low TAWSS, high OSI, high ECAP, and high RRT.Figure 1

Abstract 008: Morphomechanical Analysis of Saccular Aneurysms

Introduction An irregular morphology is considered a high‐risk feature of unruptured intracranial aneurysms. Some of these aneurysms may have blebs where there is a high likelihood of rupture. The process of aneurysm rupture is dynamic and involves factors such as aneurysm geometry, flow dynamics, and aneurysmal wall enhancement. In our study, we aimed to conduct a comprehensive analysis of brain aneurysms, specifically focusing on irregular and regular aneurysms. This analysis included the examination of aneurysmal wall enhancement (AWE), computational fluid dynamics (CFD), and finite element analysis (FEA). Methods The study focused exclusively on unruptured aneurysms and employed high‐resolution vessel wall imaging using 3T MRI. Aneurysms were classified as irregular if they displayed a bleb or multiple lobes. A 3D segmentation of the aneurysm and vessel was generated using 3D Slicer software. Orthogonal spokes were extended from the segmentation into the vessel wall, enabling the creation of intensity (SI) heatmaps. Aneurysmal wall enhancement (AWE) was calculated based on the post‐contrast SI. Furthermore, wall tension (WT) was determined using FEA, and time‐averaged wall shear stress (TAWSS) was computed through CFD utilizing the 3D segmentation. To identify regions of notable activity within the aneurysm, an analysis was conducted by examining values exceeding the 50th percentile threshold. This approach facilitated the inclusion of areas demonstrating heightened activity within the aneurysm. The study comprised two primary analyses: first, a comparison between regular and irregular aneurysms, and second, a secondary analysis that focused specifically on aneurysms with blebs. During this analysis, the bleb, body, and neck regions were isolated for comparative evaluation among different compartments within the aneurysms. Results The study included a total of forty‐seven unruptured saccular aneurysms, with 38% (18/47) classified as having irregular morphology. Comparing aneurysms with irregular morphology to those with regular morphology, it was observed that irregular aneurysms exhibited higher WT (p=0.05) and lower TAWSS (p=0.04). Among the irregular aneurysms, 61% (11/18) were found to have blebs. Within aneurysms containing blebs, significant differences were observed in terms of WT and AWE between the bleb, neck, and aneurysm body (p&lt;0.001 and p=0.03, respectively). Specifically, the blebs exhibited the highest AWE within the aneurysm, followed by the aneurysm body and the neck. The neck exhibited the highest value of WT, followed by the aneurysm body and the bleb. Finally, the neck also demonstrated the highest value of TAWSS, followed by the aneurysm body and the bleb. Conclusion Irregular aneurysms had a higher WT and lower TAWSS. Furthermore, the morphomechanical characteristics of irregular aneurysms varied between the aneurysm neck and body, indicating distinct properties in these regions. Unstable regions of the aneurysm, such as blebs, exhibited elevated AWE along with low TAWSS and WT.

Numerical study of hemodynamic changes in the Circle of Willis after stenosis of the internal carotid artery

Background and objectivesIn clinical practice a large number of patients with ischemic stroke have internal carotid artery (ICA) stenosis accompanied by Circle of Willis (CoW) stenosis. In the presence of carotid artery stenosis, CoW atherosclerosis may cause cerebral blood flow decompensation and may promote the development of ischemic stroke. The reason for the concomitant stenosis at both sites is unknown. This study investigated the hemodynamic effects of ICA stenosis on the CoW. MethodsWe developed a three-dimensional/zero-dimensional (3D/0D) closed-loop geometric multiscale model of the cerebral artery to quantify the hemodynamic indicators, including time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI). Mild (<50 %), moderate (50–69 %) and severe (>69 %) ICA stenoses were established based on 3D models of cerebral arteries in two volunteers. Geometric multiscale computing models were numerically evaluated to obtain local hemodynamic changes in the CoW in order to assess the risk of stenosis in the CoW. ResultsModel calculations showed that for all 3D models the A1 segment of the anterior cerebral artery (ACA) or the posterior communicating artery (PCA) within the CoW exhibited a hemodynamic environment with high OSI (>0.2) and low TAWSS (<1 Pa) when the ICA had a moderate stenosis. While in the case of mild and severe stenosis in ICA, there is no such phenomenon. The proportion of the surface area possessing high OSI and low TAWSS in the A1 segment of the ACA or in the PCA was mostly greater than 60 %, which might potentially cause the formation and development of atherosclerosis in CoW and finally lead to CoW stenosis. ConclusionsTherefore, although moderate carotid artery stenosis may not cause ischemic stroke, it may cause hemodynamic changes in the CoW, which in turn may promote CoW stenosis and cause CoW decompensation. In clinical treatment attention should be paid not only to stenosis of the carotid arteries but also to changes in the hemodynamic environment within the CoW, in order to prevent the adverse effects of CoW stenosis.

Numerical analysis of hemodynamic parameters in stenosed arteries under pulsatile flow conditions

This research studies the changes in flow patterns and hemodynamic parameters of diverse shapes and sizes of stenosis. Six different shapes and sizes of stenosis are constructed to investigate the variations in hemodynamics as the morphology changes. Changes in shape (trapezoidal and bell-shaped) and sizes of stenosis change the stresses on the walls and their flow patterns. TAWSS and OSI results specify that trapezoidal stenosis exerts greater stress than bell-shaped stenosis. Also, as the length of the trapezoidal stenosis increases, the TAWSS increases, whereas the trend is the opposite for bell-shaped stenosis. Later, this paper also studies different degrees of stenosis extracted from real images. Changes in velocity flow patterns, wall shear stress (WSS), Time-averaged wall shear stress (TAWSS) and Oscillatory shear index (OSI) have been studied for these images. Results illustrate that the peak velocity rises drastically as the stenosis percentage increases. Negative velocity is seen close to the artery's walls, indicating flow separation. This flow separation region is seen throughout the cycle except in the accelerating flow region. An increase in stenosis also increases WSS and TAWSS drastically. Negative WSS is seen downstream of stenosis, indicating flow recirculation. Such negative WSS in the blood vessels also promotes endothelial dysfunction. OSI values greater than 0.2 are seen near the stenosis region, indicating atherosclerosis growth. Regions of high OSI and low TAWSS are also identified, indicating probable regions of plaque development.

Assessment of Complex Flow Patterns in Patients With Carotid Webs, Patients With Carotid Atherosclerosis, and Healthy Subjects Using 4D Flow MRI.

Carotid webs (CaWs) are fibromuscular projections in the internal carotid artery (ICA) that cause mild luminal narrowing (<50%), but may be causative in up to one-third of seemingly cryptogenic strokes. Understanding hemodynamic alterations caused by CaWs is imperative to assessing stroke risk. Time-Average Wall Shear Stress (TAWSS) and Oscillatory Shear Index (OSI) are hemodynamic parameters linked to vascular dysfunction and thrombosis. To test the hypothesis: "CaWs are associated with lower TAWSS and higher OSI than mild atherosclerosis or healthy carotid bifurcation." Prospective study. A total of 35 subjects (N = 14 bifurcations with CaW, 11F, age: 49 ± 10, 10 mild atherosclerosis 6F, age: 72 ± 9, 11 healthy 9F, age: 42 ± 13). 4D flow/STAR-MATCH/3D TOF/3T MRI, CTA. 4D Flow velocity data were analyzed in two ways: 1) 3D ROI in the ICA bulbar segment (complex flow patterns are expected) was used to quantify the regions with low TAWSS and high OSI. 2) 2D planes were placed perpendicular to the centerline of the carotid bifurcation for detailed analysis of TAWSS and OSI. Independent-samples Kruskal-Wallis-H test with 0.05 used for statistical significance. The percent surface area where low TAWSS was present in the ICA bulb was 12.3 ± 8.0% (95% CI: 7.6-16.9) in CaW subjects, 1.6 ± 1.9% (95% CI: 0.2-2.9) in atherosclerosis, and 8.5 ± 7.7% (95% CI: 3.6-13.4) in healthy subjects, all differences were statistically significant (ƞ2 = 0.3 [95% CI: 0.05-0.5], P-value CaW vs. healthy = 0.2). OSI had similar values in the CCA between groups (ƞ2 = 0.07 [95% CI: 0.0-0.2], P-value = 0.5), but OSI was significantly higher downstream of the bifurcation in CaW subjects compared to atherosclerosis and normal subjects. OSI returned to similar values between groups 1.5 diameters distal to the bifurcation (ƞ2 = 0.03 [95% CI: 0.0-0.2], P-value = 0.7). Lower TAWSS and higher OSI are present in the ICA bulb in patients with CaW when compared to patients with atherosclerotic or healthy subjects. 2 TECHNICAL EFFICACY: Stage 2.

Haemodynamic changes in visceral hybrid repairs of type III and type V thoracoabdominal aortic aneurysms

The visceral hybrid procedure combining retrograde visceral bypass grafting and completion endovascular stent grafting is a feasible alternative to conventional open surgical or wholly endovascular repairs of thoracoabdominal aneurysms (TAAA). However, the wide variability in visceral hybrid configurations means that a priori prediction of surgical outcome based on haemodynamic flow profiles such as velocity pattern and wall shear stress post repair remain challenging. We sought to appraise the clinical relevance of computational fluid dynamics (CFD) analyses in the setting of visceral hybrid TAAA repairs. Two patients, one with a type III and the other with a type V TAAA, underwent successful elective and emergency visceral hybrid repairs, respectively. Flow patterns and haemodynamic parameters were analysed using reconstructed pre- and post-operative CT scans. Both type III and type V TAAAs showed highly disturbed flow patterns with varying helicity values preoperatively within their respective aneurysms. Low time-averaged wall shear stress (TAWSS) and high endothelial cell action potential (ECAP) and relative residence time (RRT) associated with thrombogenic susceptibility was observed in the posterior aspect of both TAAAs preoperatively. Despite differing bypass configurations in the elective and emergency repairs, both treatment options appear to improve haemodynamic performance compared to preoperative study. However, we observed reduced TAWSS in the right iliac artery (portending a theoretical risk of future graft and possibly limb thrombosis), after the elective type III visceral hybrid repair, but not the emergency type V repair. We surmise that this difference may be attributed to the higher neo-bifurcation of the aortic stent graft in the type III as compared to the type V repair. Our results demonstrate that CFD can be used in complicated visceral hybrid repair to yield potentially actionable predictive insights with implications on surveillance and enhanced post-operative management, even in patients with complicated geometrical bypass configurations.

Abstract P1048: Lattice-Boltzmann Method Achieves Feasible Hemodynamic Analysis In Carotid Atherosclerosis With Low Computational Cost

Introduction: Atherosclerosis is a common cardiovascular pathophysiology that can lead to plaque-related complications, including stroke and myocardial infarction. Non-invasive hemodynamic analysis can provide valuable insights into the plaque's rupture vulnerability and growth progression. Computational fluid dynamics is widely studied to understand hemodynamic pathophysiology further; however, traditional computational fluid dynamics methods often require daily-based analysis time to visualize hemodynamic parameters. This study investigated the feasibility of using the lattice Boltzmann method (LBM) for hemodynamic analysis of atherosclerosis lesions. Methods: From ten carotid artery stenosis patients, 20 carotid artery models were generated and analyzed using LBM. Hemodynamic parameters, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT), were compared between the stenotic lesion and contralateral control sides. The patient-specific boundary conditions were applied. Their cardiac cycle was subdivided into three periods: systolic acceleration, systolic deceleration, and diastolic, and analyzed to investigate their relationship with the hemodynamic parameters. Results: LBM effectively captured the hemodynamic parameters associated with carotid atherosclerosis with an hourly-based computational cost (2.765 [2.435-3.233] hours). The analysis visualized high TAWSS and low RRT distribution at stenotic areas during the systolic deceleration period and low TAWSS and high RRT distribution at the post-stenotic distal lesion during the diastolic period. These findings were reproducible among the ten patients. Conclusion: LBM-based hemodynamic analysis has the potential to serve as a non-invasive diagnostic tool that can help understand atherosclerosis pathophysiology with low computational cost.

Hemodynamic parameters impact the stability of distal stent graft-induced new entry

Stent graft-induced new entry tear (SINE) is a serious complication in aortic dissection patients caused by the stent-graft itself after thoracic endovascular aortic repair (TEVAR). The stability of SINE is a key indicator for the need and timing of reinterventions. This study aimed to understand the role of hemodynamics in SINE stability by means of computational fluid dynamics (CFD) analysis based on patient-specific anatomical information. Four patients treated with TEVAR who developed a distal SINE (dSINE) were included; two patients had a stable dSINE and two patients experienced expansion of the dSINE upon follow-up examinations. CFD simulations were performed on geometries reconstructed from computed tomography scans acquired upon early detection of dSINE in these patients. Computational results showed that stable dSINEs presented larger regions with low time-averaged wall shear stress (TAWSS) and high relative residence time (RRT), and partial thrombosis was observed at subsequent follow-ups. Furthermore, significant systolic antegrade flow was observed in the unstable dSINE which also had a larger retrograde flow fraction (RFF) on the SINE plane. In conclusion, this pilot study suggested that high RRT and low TAWSS may indicate stable dSINE by promoting thrombosis, whereas larger RFF and antegrade flows inside dSINE might be associated with its expansion.

Abnormal flow pattern of low wall shear stress and high oscillatory shear index in spontaneous vertebral artery dissection with vertebral artery hypoplasia.

Spontaneous vertebral artery dissection (sVAD) might tend to develop in vertebral artery hypoplasia (VAH) with hemodynamic dysfunction and it is crucial to assess hemodynamics in sVAD with VAH to investigate this hypothesis. This retrospective study aimed to quantify hemodynamic parameters in patients with sVAD with VAH. Patients who had suffered ischemic stroke due to an sVAD of VAH were enrolled in this retrospective study. The geometries of 14 patients (28 vessels) were reconstructed using Mimics and Geomagic Studio software from CT angiography (CTA). ANSYS ICEM and ANSYS FLUENT were utilized for mesh generation, set boundary conditions, solve governing equations, and perform numerical simulations. Slices were obtained at the upstream area, dissection or midstream area and downstream area of each VA. The blood flow patterns were visualized through instantaneous streamline and pressure at peak systole and late diastole. The hemodynamic parameters included pressure, velocity, time-averaged blood flow, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell action potential (ECAP), relative residence time (RRT) and time-averaged nitric oxide production rate (TARNO). Significant focal increased velocity was present in the dissection area of steno-occlusive sVAD with VAH compared to other nondissected areas (0.910 m/s vs. 0.449 vs. 0.566, p < 0.001), while focal slow flow velocity was observed in the dissection area of aneurysmal dilatative sVAD with VAH according to velocity streamlines. Steno-occlusive sVAD with VAH arteries had a lower time-averaged blood flow (0.499 cm3/s vs. 2.268, p < 0.001), lower TAWSS (1.115 Pa vs. 2.437, p = 0.001), higher OSI (0.248 vs. 0.173, p = 0.006), higher ECAP (0.328 Pa-1 vs. 0.094, p = 0.002), higher RRT (3.519 Pa-1 vs. 1.044, p = 0.001) and deceased TARNO (104.014 nM/s vs. 158.195, p < 0.001) than the contralateral VAs. Steno-occlusive sVAD with VAH patients had abnormal blood flow patterns of focal increased velocity, low time-averaged blood flow, low TAWSS, high OSI, high ECAP, high RRT and decreased TARNO. These results provide a good basis for further investigation of sVAD hemodynamics and support the applicability of the CFD method in testing the hemodynamic hypothesis of sVAD. More detailed hemodynamic conditions with different stages of sVAD are warranted in the future.