Oscillatory Shear Index Research Articles

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.

Hemodynamics in the treatment of pseudoaneurysm caused by extreme constriction of aortic arch with coated stent.

To evaluate the changes in distal vascular morphology and hemodynamics in patients with extremely severe aortic coarctation (CoA) after covered palliative (CP) stent dilation with different surgical strategies. Perioperative computed tomography angiography and digital subtraction angiography were utilized to construct three aortic models with varying stenosis rates and one follow-up model in a patient with extremely severe CoA. The models included: an idealized non-stenosed model (A: 0%), a model post initial stent deployment (B: 28%), a model post balloon expansion (C: 39%), and a model 18 months after post-balloon expansion (D: 39%). Consistent boundary conditions were applied to all models, and hemodynamic simulation was conducted using the pure fluid method. The narrowest and distal diameter of the stent increased by 34.71% and 59.29%, respectively, from model B to C. Additionally, the distal diameter of the stent increased by -13.80% and +43.68% compared to the descending aorta diameter, respectively. Furthermore, the ellipticity of the maximum cross-section of the aneurysm region in model A to D continued to increase. The oscillatory shear index at the stenosis to the region of the aneurysm were found to be higher in Models A and B, and lower in Models C and D. At the moment of maximum flow velocity, the blood flow distribution in models A and B was more uniform in the widest section of the blood vessels at the distal end of the stenosis, whereas models C and D exhibited disturbed blood flow with more than 2 eddy currents. The time-averaged wall shear stress (TAWSS) decreased in the distal and basal aneurysms, while it significantly increased at the step position. The aneurysmal region exhibited an endothelial cell activation potential value lower than 0.4 Pa-1. In patients with extremely severe CoA, it is crucial to ensure that the expanded diameter at both ends of the CP stent does not exceed the native vascular diameter during deployment. Our simulation results demonstrate that overdilation leads to a decrease in the TAWSS above the injured vessel, creating an abnormal hemodynamic environment that may contribute to the development and enlargement of false aneurysms in the early postoperative period. ClinicalTrials.gov, (NCT02917980).

Numerical assessment of using various outlet boundary conditions on the hemodynamics of an idealized left coronary artery model.

Vascular diseases are greatly influenced by the hemodynamic parameters and the accuracy of determining these parameters depends on the use of correct boundary conditions. The present work carries out a two-way fluid-structure interaction (FSI) simulation to investigate the effects of outlet pressure boundary conditions on the hemodynamics through the left coronary artery bifurcation with moderate stenosis (50%) in the left anterior descending (LAD) branch. The Carreau viscosity model is employed to characterise the shear-thinning behaviour of blood. The results of the study reveal that the employment of zero pressure at the outlet boundaries significantly overestimates the values of hemodynamic variables like wall shear stress (WSS), and time-averaged wall shear stress (TAWSS) compared with human healthy and pulsatile pressure outlet conditions. However, the difference between these variables is marginally low for human healthy and pulsatile pressure outlets. The oscillatory shear index (OSI) remains the same across all scenarios, indicating independence from the outlet boundary condition. Furthermore, the magnitude of negative axial velocity and pressure drop across the plaque are found to be higher at the zero pressure outlet boundary condition.

Validation of classification system for isolated superior mesenteric artery dissections using image-based computational flow analysis

Background and ObjectiveThe isolated superior mesenteric artery dissection (ISMAD) is a rare but potentially fatal vascular disorder. Classifications for ISMAD were previously proposed based on morphometric features. However, the classification systems were not standardized and verified yet. This study conducted computational flow analysis to validate the latest classification system of ISMAD and aid clinical decision-making based on hemodynamic parameters. Methods62 patients with ISMAD were included and classified into different types according to false lumen structures (five types, Type I-V) and true lumen patency (two types, Type P and Type S) according to Qiu classification system. Computational fluid dynamics and three-dimensional structural analyses were conducted on the basis of computed tomography angiography datasets. Quantitative and qualitative functional analyses were performed via parameters of interest including volume flow of each minute, pressure drop, pressure gradient, the derivative parameters of wall shear stress such as time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and the relative residence time (RRT). Statistical analyses were conducted among different ISMAD types. ResultsTAWSS, OSI and RRT showed significant difference among different types when classified using false lumen structures. In detail, Type IV showed significantly higher TAWSS than other types (p = 0.007). OSI was obviously higher in Type II (p = 0.015). Type IV also presented the lowest RRT (p = 0.005). The pressure drop, pressure gradient, OSI and RRT showed higher value in Type S than that in Type P, demonstrating a statistical significance with p values of 0.017, 0.041, 0.001 and 0.012, respectively. While Type P had larger volume flow than Type S (p = 0.041). ConclusionsThe notable differences in hemodynamic features among different types demonstrated the feasibility of Qiu classification system. The evaluation based on hemodynamic simulation might also provide insights into risk identification and guide therapeutic decisions for ISMAD.

Can β-blockers prevent intracranial aneurysm rupture? - insights from Computational Fluid Dynamics analysis.

Hypertension is a risk factor for intracranial aneurysm rupture. We analyzed whether the intake of drugs from specific classes of anti-hypertensive medications affects hemodynamic parameters of intracranial aneurysm dome. We recorded medical history including medications and the in-hospital blood pressure values. We then obtained 3D reconstruction of each patients' aneurysm dome and the feeding artery. Using OpenFOAM software we performed Computational Fluid Dynamics analysis of blood flow through the modeled structures. Blood was modeled as Newtonian fluid, using the incompressible transient solver. As the inlet boundary condition we used the patient-specific Internal Carotid Artery blood velocity waves obtained with Doppler ultrasound. We calculated haemodynamic parameters of the aneurysm dome. All presented analyses are cross-sectional.We included 72 patients with a total of 91 unruptured intracranial aneurysms. The history of β-blocker intake significantly influenced hemodynamic parameters of aneurysm dome. The patients on β-blockers had significantly smaller aneurysm domes (5.09 ± 2.11 mm vs. 7.41 ± 5.89 mm; p = 0.03) and did not have aneurysms larger than 10 mm (0% vs 17.0%; p = 0.01). In the Computational Fluid Dynamics analysis, walls of aneurysms in patients who took β-blockers were characterized by lower Wall Shear Stress Gradient (1.67 ± 1.85 Pa vs. 4.3 ± 6.06 Pa; p = 0.03), Oscillatory Shear Index (0.03 ± 0.02 vs. 0.07 ± 0.10; p = 0.04) and Surface Vortex Fraction (16.2% ± 5.2% vs. 20.0% ± 6.8%; p<0.01). After controlling for covariates, we demonstrated difference of Surface Vortex Fraction (F[1, 48] = 4.36; p = 0.04) and Oscillatory Shear Index (F[1, 48] = 6.51; p = 0.01) between patients taking and not taking β-blockers, respectively. Intake of β-blockers might contribute to more favorable hemodynamics inside aneurysmal sac. Other antihypertensive medication classes were not associated with differences in intracranial aneurysm parameters.

In-silico analysis of hemodynamic indicators in idealized stented coronary arteries for varying stent indentation

In this work, we investigate the effects of stent indentation on hemodynamic indicators in stented coronary arteries. Our aim is to assess in-silico risk factors for in-stent restenosis (ISR) and thrombosis after stent implantation. The proposed model is applied to an idealized artery with Xience V stent for four indentation percentages and three mesh refinements. We analyze the patterns of hemodynamic indicators arising from different stent indentations and propose an analysis of time-averaged WSS (TAWSS), topological shear variation index (TSVI), oscillatory shear index (OSI), and relative residence time (RRT). We observe that higher indentations display higher frequency of critically low TAWSS, high TSVI, and non-physiological OSI and RRT. Furthermore, an appropriate mesh refinement is needed for accurate representation of hemodynamics in the stent vicinity. The results suggest that disturbed hemodynamics could play a role in the correlation between high indentation and ISR.

Quantification of morpho-hemodynamic changes in unruptured intracranial aneurysms with irregular pulsation during the cardiac cycle using 4D-CTA.

Previous studies predicting the rupture risk of intracranial aneurysms (IAs) have predominantly utilized static imaging data, overlooking the dynamic blood flow and biomechanical properties of the aneurysm wall. Irregular pulsation detected by 4D-CTA is a potential predictor of aneurysm rupture, albeit with uncertain clinical significance. This study aimed to analyze the changes in morpho-hemodynamic characteristics of IAs during the cardiac cycle to elucidate the dynamic changes and the associated hemodynamic mechanisms. A retrospective review was conducted on the 4D-CTA data of IA patients between January 2017 and September 2019. R-R intervals were segmented into 20-time phases, reconstructing 20 CT datasets to identify irregular pulsation and extract 3D aneurysm models. Computational fluid dynamics (CFD) simulations analyzed hemodynamic parameters such as oscillatory shear index (OSI) and relative residence time (RRT). Changes in morpho-hemodynamic characteristics were quantified in terms of the absolute change (parameter*) and relative change rate (parameter%). Rupture risk was assessed using the rupture resemblance model (RRS). Eleven UIAs from 10 patients were finally included, with five aneurysms showing irregular pulsation (45.45%). No significant differences in morpho-hemodynamic characteristics were observed between aneurysms with or without irregular pulsation. More remarkable changes in aneurysm size (size*: 0.59 ± 0.14 mm vs. 0.32 ± 0.12 mm, p = 0.010; size%: 10.49% ± 1.43% vs. 3.95% ± 1.79%, p < 0.001), volume (volume%: 13.72% vs. 6.39%, p = 0.009), OSI (OSI*: 0.02 ± 0.01 vs. 0.004 ± 0.005, p = 0.004; OSI%: 200% vs. 12.50%, p = 0.004) and RRT (RRT%: 97.14% vs. 43.95, p = 0.052) over the cardiac cycle were significantly linked to irregular pulsation. Aneurysms with irregular pulsation demonstrated a more unfavorable hemodynamic environment during the cardiac cycle, irrespective of the predicted rupture risk. Furthermore, irregular pulsation at the aneurysm dome exhibited higher hemodynamic instability than at the sidewall. Irregular pulsation may indicate hemodynamic instability within the aneurysm, leading to an increased rupture risk in the area where irregular pulsation occurs. This proof-of-concept study could enhance understanding of dynamic changes in UIAs during the cardiac cycle and the underlying hemodynamic mechanisms.

A comprehensive MRI-based computational model of blood flow in compliant aorta using radial basis function interpolation

BackgroundProperly understanding the origin and progression of the thoracic aortic aneurysm (TAA) can help prevent its growth and rupture. For a better understanding of this pathogenesis, the aortic blood flow has to be studied and interpreted in great detail. We can obtain detailed aortic blood flow information using magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) with a prescribed motion of the aortic wall.MethodsWe performed two different types of simulations—static (rigid wall) and dynamic (moving wall) for healthy control and a patient with a TAA. For the latter, we have developed a novel morphing approach based on the radial basis function (RBF) interpolation of the segmented 4D-flow MRI geometries at different time instants. Additionally, we have applied reconstructed 4D-flow MRI velocity profiles at the inlet with an automatic registration protocol.ResultsThe simulated RBF-based movement of the aorta matched well with the original 4D-flow MRI geometries. The wall movement was most dominant in the ascending aorta, accompanied by the highest variation of the blood flow patterns. The resulting data indicated significant differences between the dynamic and static simulations, with a relative difference for the patient of 7.47±14.18% in time-averaged wall shear stress and 15.97±43.32% in the oscillatory shear index (for the whole domain).ConclusionsIn conclusion, the RBF-based morphing approach proved to be numerically accurate and computationally efficient in capturing complex kinematics of the aorta, as validated by 4D-flow MRI. We recommend this approach for future use in MRI-based CFD simulations in broad population studies. Performing these would bring a better understanding of the onset and growth of TAA.

Understanding the role of the left atrial appendage on the flow in the atrium.

The exclusion/occlusion of the left atrial appendage (LAA) is a treatment option for atrial fibrillation (AF) patients who are at high risk of stroke and high risk of bleeding. As the role of the LAA is not well understood or explored, this study aims to assess its role on flow dynamics in the left atrium. Computational fluid dynamics (CFD) simulations were carried out for nine AF patients before and after LAA exclusion. The flow parameters investigated included the LA velocities, Time Averaged Wall Shear Stress (TAWSS), Oscillatory Shear Index (OSI), Relative Residence Time (RRT), and Pressure in the LA. This study shows that, on average, a decrease in TAWSS (1.82 ± 1.85 Pa to 1.27 ± 0.96 Pa, p < 0.05) and a slight increase in OSI (0.16 ± 0.10 to 0.17 ± 0.10, p < 0.05), RRT (1.87 ± 1.84 Pa-1 to 2.11 ± 1.78 Pa-1, p < 0.05), and pressure (-19.2 ± 6.8 mmHg to -15.3 ± 8.3 mmHg, p < 0.05) were observed in the LA after the exclusion of the LAA, with a decrease in low-magnitude velocities. The exclusion of the LAA seems to be associated with changes in LA flow dynamics. Further studies are needed to elucidate the clinical implications of these changes.

Non-invasive fractional flow reserve estimation in coronary arteries using angiographic images

Coronary artery disease is the leading global cause of mortality and Fractional Flow Reserve (FFR) is widely regarded as the gold standard for assessing coronary artery stenosis severity. However, due to the limitations of invasive FFR measurements, there is a pressing need for a highly accurate virtual FFR calculation framework. Additionally, it’s essential to consider local haemodynamic factors such as time-averaged wall shear stress (TAWSS), which play a critical role in advancement of atherosclerosis. This study introduces an innovative FFR computation method that involves creating five patient-specific geometries from two-dimensional coronary angiography images and conducting numerical simulations using computational fluid dynamics with a three-element Windkessel model boundary condition at the outlet to predict haemodynamic distribution. Furthermore, four distinct boundary condition methodologies are applied to each geometry for comprehensive analysis. Several haemodynamic features, including velocity, pressure, TAWSS, and oscillatory shear index are investigated and compared for each case. Results show that models with average boundary conditions can predict FFR values accurately and observed errors between invasive FFR and virtual FFR are found to be less than 5%.

Unsteady electroosmotic flow of Carreau–Newtonian fluids through a cylindrical tube

The present study aims to investigate the fully developed electroosmotic flow of Carreau–Newtonian fluids within a circular tube, considering two different boundary conditions of zeta potential. The unsteady behavior of Carreau–Newtonian fluids, is attributed to the pulsatile pressure-driven flow. Studying two distinct formulations, slip and no-slip zeta potentials, illustrates the relative movement of fluid particles and the charged surface. The proposed study’s framework is divided into two regions: a central region containing a non-Newtonian Carreau fluid exhibiting both shear-thinning and thickening behavior, and a peripheral region containing a Newtonian fluid with the effect of an electrical double layer at the solid wall of the cylindrical tube. The Poisson–Boltzmann equation under the assumption of Debye–Hückel approximation governs the potential distribution due to the electric double layer, facilitating the examination of electroosmotic flow through tube. Obtaining analytical solutions for the momentum equations in both regions becomes challenging due to the inclusion of pulsatile pressure-driven flow and a non-Newtonian Carreau fluid. Asymptotic series expansions are employed, utilizing small parameters such as the Weissenberg number (We≪1) and a pulsatile Reynolds number (α≪1), to derive velocity expressions for both regions. Wall shear stress fluctuations, as indicated by time-averaged wall shear stress and oscillatory shear index are assessed to gauge the temporal variation of wall shear stress from the dominant blood flow direction. By utilizing potential and hydrodynamic variables, an asymptotic formulation for the streaming potential is developed to ascertain the asymptotic expression of electrokinetic energy conversion efficiency and scrutinize the impacts of pertinent physical parameters on it. The proposed study prominently shows that the Debye–Hückel parameter, Weissenberg number, pulsatile Reynolds number, and time-dependent pressure gradient significantly influence the fluid velocity, flow rate, and flow resistance. The notable determination of the present study is that the velocity amplitude and electrokinetic energy conversion efficiency are enlarged in the slip-dependent case relative to the no-slip zeta potential. It is noteworthy that as inertial forces become more prominent compared to viscous forces, both time-average wall shear stress and oscillatory shear index experience a slight decrease, although the reduction in oscillatory shear index is minimal owing to the absence of obstruction within the tube wall. The NDSolve numerical scheme in Mathematica software is employed to calculate numerical solutions for the governing equations, which are subsequently verified against results obtained through asymptotic analysis. The findings of this study are expected to deepen our understanding of unsteady electroosmotic flows of Carreau fluid and contribute to the design and development of advanced microfluidic devices with enhanced performance for various applications.

Turbulent blood flow in a cerebral artery with an aneurysm

Unruptured intracranial aneurysms are common in the general population, and many uncertainties remain when predicting rupture risks and treatment outcomes. One of the cutting-edge tools used to investigate this condition is computational fluid dynamics (CFD). However, CFD is not yet mature enough to guide the clinical management of this disease. In addition, recent studies have reported significant flow instabilities when refined numerical methods are used. Questions remain as to how to properly simulate and evaluate this flow, and whether these instabilities are really turbulence. The purpose of the present study is to evaluate the impact of the simulation setup on the results and investigate the occurrence of turbulence in a cerebral artery with an aneurysm. For this purpose, direct numerical simulations were performed with up to 200 cardiac cycles and with data sampling rates of up to 100,000 times per cardiac cycle. Through phase-averaging or triple decomposition, the contributions of turbulence and of laminar pulsatile waves to the velocity, pressure and wall shear stress fluctuations were distinguished. For example, the commonly used oscillatory shear index was found to be closely related to the laminar waves introduced at the inlet, rather than turbulence. The turbulence energy cascade was evaluated through energy spectrum estimates, revealing that, despite the low flow rates and Reynolds number, the flow is turbulent near the aneurysm. Phase-averaging was shown to be an approach that can help researchers better understand this flow, although the results are highly dependent on simulation setup and post-processing choices.

AnXplore: a comprehensive fluid-structure interaction study of 101 intracranial aneurysms.

Advances in computational fluid dynamics continuously extend the comprehension of aneurysm growth and rupture, intending to assist physicians in devising effective treatment strategies. While most studies have first modelled intracranial aneurysm walls as fully rigid with a focus on understanding blood flow characteristics, some researchers further introduced Fluid-Structure Interaction (FSI) and reported notable haemodynamic alterations for a few aneurysm cases when considering wall compliance. In this work, we explore further this research direction by studying 101 intracranial sidewall aneurysms, emphasizing the differences between rigid and deformable-wall simulations. The proposed dataset along with simulation parameters are shared for the sake of reproducibility. A wide range of haemodynamic patterns has been statistically analyzed with a particular focus on the impact of the wall modelling choice. Notable deviations in flow characteristics and commonly employed risk indicators are reported, particularly with near-dome blood recirculations being significantly impacted by the pulsating dynamics of the walls. This leads to substantial fluctuations in the sac-averaged oscillatory shear index, ranging from -36% to +674% of the standard rigid-wall value. Going a step further, haemodynamics obtained when simulating a flow-diverter stent modelled in conjunction with FSI are showcased for the first time, revealing a 73% increase in systolic sac-average velocity for the compliant-wall setting compared to its rigid counterpart. This last finding demonstrates the decisive impact that FSI modelling can have in predicting treatment outcomes.

ASSESSMENT OF GROWTH AND RUPTURE OF CAROTID ATHEROSCLEROSIS AND PLAQUE WITH DIFFERENT DEGREES OF STENOSIS

In this study, the growth and rupture of plaque were investigated to predict the development of carotid atherosclerosis and plaque. Realistic 3D models of multi-component plaque and carotid artery were established in healthy, mild, moderate, and severe stenosis. Fluid–structure interaction (FSI) simulations were performed on these models, and the wall shear stress (WSS), time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and wall tensile stress (WTS) in a natural cardiac cycle were analyzed to assess the epidermal rupture, growth tendency, and internal rupture of plaques in a natural state. Results show that in severe stenosis, the WSS at the stenosis sites exceeds the threshold of rupture during almost the entire cardiac cycle, thus causing an epidermal injury, bringing about plaque detachment and thrombosis. TAWSS decreased and OSI and RRT increased in the bifurcation region and the downstream region of the plaque. The deposition is more likely to occur in these regions. A higher degree of stenosis will increase the OSI and RRT in the downstream region of the plaque, leading to continuous deterioration of the plaque, while the degree of stenosis has little effect on the upstream region of the plaque. The fibrous cap is the starting point for rupture within the plaque, and calcification increases the stress on the fibrous cap, thereby increasing the risk of rupture.

Quantitative evaluation using carotid ultrasonography-based high-frame-rate vector flow imaging in patients with low carotid stenosis.

To explore the role of quantitative evaluation using carotid ultrasonography (US)-based high-frame-rate vector flow (V Flow) imaging in patients with low carotid stenosis. This single-centre cross-sectional study consecutively recruited volunteers without carotid plaque and patients with low carotid stenosis from August 2022 to May 2023. Patients were divided into symptomatic and asymptomatic groups according to their head CT or MRI results within 8 weeks. All V Flow imaging examinations were performed using a Mindray Resona R9 US system. The wall shear stress (WSS) values, oscillatory shear index (OSI) values, and turbulence (Tur) indexes in the normal common carotid artery (CCA), normal carotid bifurcation (CB), and on the upstream and downstream surface of carotid plaque were measured. Pearson Chi-square test and Fisher exact test were used for counting data according to their type. For measurement data, independent sample t test and non-parametric rank sum test were used. The results proved that patients have higher WSS values and Tur indexes of CB than volunteers, and higher WSS values were detected on the surface of the plaques in symptomatic patients. What's more, the downstream side of the plaque was more vulnerable to plaque rupture than the upstream side due to more dynamic blood flow. Therefore, carotid US-based high-frame-rate V Flow imaging provides reliable mechanical biomarkers for assessing the haemodynamic change in patients with low stenosis. Our study may provide a new imaging tool for monitoring the progression of atherosclerosis and aiding the management of early atherosclerotic patients. Our study firstly investigated the difference of V Flow parameters on the surface of carotid plaques between symptomatic and asymptomatic patients with low carotid stenosis, which is expected to provide haemodynamic information and the mechanical basis for plaque rupture.

Ruptured Vertebral Artery Dissection Treated with Endovascular Intervention

A retrospective study of cases of endovascular treatment of dissection of the vertebral artery with subarachnoid hemorrhage was conducted. Data were 11 cases of vertebral artery dissecting aneurysm (VADA) among 291 consecutive subarachnoid hemorrhage patients who underwent clipping or endovascular treatment at Ota Memorial Hospital. Classified into 4 patterns based on the location of the dissection and posterior inferior cerebellar artery (PICA): pre-PICA, post-PICA, involved PICA, and non-PICA. And one of the cases had bilateral vertebral artery dissection, and computational fluid dynamics analysis was included in the study. Ruptured VADA occurred in 11 of the 291 patients (3.8%). Endovascular treatment was performed in 8 of these 11 patients. Postoperative diffusion-weighted imaging detected no high-intensity lesions and no postoperative ischemic complications or rebleeding occurred in any patient. In a case of bilateral VADA, computational fluid dynamics analysis of very low or high wall shear stress at the dissection, low aneurysm formation indicator, and high oscillatory shear index may be considered rupture factors. Treatment strategies for each branching pattern of PICA can prevent rupture and avoid ischemic complications. And prediction of the rupture side is important in patients with bilateral dissection to consider the appropriate treatment and timing.

Hemodynamic coupling between a primary atherosclerotic plaque and subsequent secondary lesions

Atherosclerosis is one of the most common diseases of the arterial tree, especially in the coronary arteries. Stenoses exceeding 50% area reduction are shown to alternate the downstream coronary flow, and hemodynamics will lead to further atherogenesis. Clinical evidence also confirms that vascular stenoses are not stagnant since they are usually associated with downstream lesions. In this study, it is hypothesized that the formation of secondary plaques, or aneurysms downstream of a primary stenosis, compensates for the abnormal ranges of hemodynamic forces caused by the primary stenosis. An experimental setup captured the hemodynamics of non-Newtonian blood flow in three-dimensional (3D) printed phantoms of coronary arteries with various sequences of lesions. Then, based on the collected data, in silico models of these lesions were simulated using computational fluid dynamics. For the proposed cases, time-averaged wall shear stress, velocity profile, oscillatory shear index, and relative residence time were extracted at the plaque side and the plaque front walls and compared to the reference model with only the primary plaque. The secondary plaque postulated the abnormal hemodynamic conditions to its downstream, which implies endothelial activation and onset of further pathologic events. However, the secondary aneurysm restored flow conditions to normal after its distal shoulder, preventing more damage to the endothelium. Examined angiograms of patients with developed atherosclerotic lesions unveiled that a sequence of plaques is formed over time, and most interestingly, the series stopped after the formation of an aneurysm.

A comparative study of Newtonian and non-Newtonian blood flow through Bi-Leaflet Mechanical Heart Valve

The present study examines flow through Bi-Leaflet Mechanical Heart Valves (BMHV) at physiological conditions considering both Newtonian and non-Newtonian fluid models for blood rheology. It is well known that the non-Newtonian effects of blood are pronounced in small diameter arteries. Most of the earlier works on Mechanical Heart Valves (MHV) have considered blood as a Newtonian fluid as the flow involves large-diameter artery such as the aorta. In this work, we have reported the predicted parameters, such as leaflet kinematics, vortex structures, wall shear stress, and blood damage index for both blood models. It is found that the leaflet attributes smaller asynchronous motion in the case of non-Newtonian Carreau fluid model with slightly reduced angular velocity compared to the Newtonian assumption. Predictions on the blood damage index suggest a 21% higher damage while using non-Newtonian model than Newtonian model, which may be attributed to higher levels of mechanical stress within the fluid. However, vortex structures, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) are found to be similar in predictions using both the fluid models. We have used an in-house sharp interface immersed boundary method with fluid–structure interaction to simulate the coupled action of moving valves and pulsatile blood flow. Our findings suggest that the general consensus of using Newtonian model in large arteries may not be appropriate for prediction of leaflet kinematics and blood damage index in Mechanical heart valves.

#806 Predicting arteriovenous fistula success using ferumoxytol-enhanced MRI (Fe-MRI) and computational flow dynamics (CFD)

Abstract Background and Aims Arteriovenous fistula (AVF) failure is common and is invariably associated with preexisting anatomic problems of the inflow and outflow circuit. Ferumoxytol, an iron oxide nanoparticle, provides a safe alternative to gadolinium contrast material for MR angiography in chronic kidney disease (CKD). Computational fluid dynamics (CFD) has been utilised to virtually demonstrate the haemodynamics within AVFs using high-quality vascular mapping obtained with ferumoxytol-enhanced MRI (Fe-MRI). Using Fe-MRI and CFD we aimed to assess AVF haemodynamics and identify flow patterns and geometries associated with AVF success. Method In this pilot study we used Fe-MRI datasets of CKD patients who had vascular mapping 6 weeks after AVF creation. Using CFD methods we performed vessel segmentation to construct three-dimensional (3D) patient-specific geometries from images obtained with Fe-MRI (Fig. 1). We then completed a detailed CFD analysis for each AVF geometry to establish the associations of haemodynamic (mean wall shear stress (WSS), oscillatory shear index (OSI), blood flow velocity, and areas of abnormal flow) and anatomical parameters (anastomosis angle, feeding artery diameter, feeding artery curvature, draining vein diameter, and draining vein curvature) with AVF outcomes. The primary outcome was AVF success defined as AVF use (assisted or unassisted) for dialysis for a minimum of 3 months. ROC analysis was used to assess anatomical predictors of AVF flows of ≥1000 ml/min. Results The study included 17 patients with CKD who had AVF creation between 2018 and 2020. From these, 13 AVFs were successfully used for dialysis and 4 failed to mature or were used for dialysis for &amp;lt;3 months before failure. Compared to failed fistulas, successful AVFs had higher mean WSS and OSI. Failed AVFs were characterised by different haemodynamic patterns including lower flow rates with less helical flows. On ROC analysis the three metrics with the highest area under the curve (AUC) values were the feeding artery curvature (0.82), feeding artery diameter (0.76), and draining vein diameter (0.74) with a combined AUC value of 0.83 (Fig. 2). Conclusion Using Fe-MRI and CFD in the post-surgical juxta-anastomotic environment higher mean WSS and OSI were predictive of AVF success. In addition, the feeding artery curvature, feeing artery and draining vein diameters are strongly correlated with high venous flow rates, one of the key factors of AVF success. Although vessel diameters are routinely measured as part of pre-operative vascular mapping, a greater emphasis should be given on mean WSS, OSI and feeding artery vessel curvature when selecting vessels for AVF creation.

Rupture point is associated with divergent hemodynamics in intracranial aneurysms.

Understanding the risk factors leading to intracranial aneurysm (IA) rupture have still not been fully clarified. They are vital for proper medical guidance of patients harboring unruptured IAs. Clarifying the hemodynamics associated with the point of rupture could help could provide useful information about some of the risk factors. Thus far, few studies have studied this issue with often diverging conclusions. We identified a point of rupture in patients operated for an IAs during surgery, using a combination of preoperative computed tomography (CT) and computed tomography angiography (CTA). Hemodynamic parameters were calculated both for the aneurysm sac as a whole and the point of rupture. In two cases, the results of CFD were compared with those of the experiment using particle image velocimetry (PIV). We were able to identify 6 aneurysms with a well-demarcated point of rupture. In four aneurysms, the rupture point was near the vortex with low wall shear stress (WSS) and high oscillatory shear index (OSI). In one case, the rupture point was in the flow jet with high WSS. In the last case, the rupture point was in the significant bleb and no specific hemodynamic parameters were found. The CFD results were verified in the PIV part of the study. Our study shows that different hemodynamic scenarios are associated with the site of IA rupture. The numerical simulations were confirmed by laboratory models. This study further supports the hypothesis that various pathological pathways may lead to aneurysm wall damage resulting in its rupture.