Curvature effect on hemodynamic conditions at the inner bend of the carotid siphon and its relation to aneurysm formation

Although high-impact hemodynamic forces are thought to lead to cerebral aneurysmal change, little is known about the aneurysm formation on the inner aspect of vascular bends such as the internal carotid artery (ICA) siphon where wall shear stress (WSS) is expected to be low. This study evaluates the effect of vessel curvature and hemodynamics on aneurysm formation along the inner carotid siphon. Catheter 3D-rotational angiographic volumes of 35 ICA (10 aneurysms, 25 controls) were evaluated in 3D for radius of curvature and peak curvature of the siphon bend, followed by univariate statistical analysis. Computational fluid dynamic (CFD) simulations were performed on patient-derived models after aneurysm removal and on synthetic variants of increasing curvature. Peak focal siphon curvature was significantly higher in aneurysm bearing ICAs (0.36±0.045 vs. 0.30±0.048 mm⁻¹, p=0.003), with no difference in global radius of curvature (p=0.36). In CFD simulations, increasing parametric curvature tightness (from 5 to 3 mm radius) resulted in dramatic increase of WSS and WSS gradient magnitude (WSSG) on the inner wall of the bend. In patient-derived data, the location of aneurysms coincided with regions of low WSS (<4 Pa) flanked by high WSS and WSSG peaks. WSS peaks correlated with the aneurysm neck. In contrast, control siphon bends displayed low, almost constant, WSS and WSSG profiles with little spatial variation. High bend curvature induces dynamically fluctuating high proximal WSS and WSSG followed by regions of flow stasis and recirculation, leading to local conditions known to induce destructive vessel wall remodeling and aneurysmal initiation.     

Carotid geometry effects on blood flow and on risk for vascular disease

It has been widely observed that atherosclerotic diseases occur at sites with complex hemodynamics, such as artery bifurcations, junctions, and regions of high curvature. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present work, 3D pulsatile blood flow through a model of the carotid artery bifurcation was simulated using a finite volume numerical method. The goal was to quantify the risk of atherogenesis associated with different carotid artery geometries. A risk scale based on the average WSS on the sinus wall of the internal carotid artery was proposed-a scale that can be used to quantify the effect of the carotid geometry on the relative risk for developing vascular disease. It was found that the bifurcation angle and the out-of-plane angle of the internal carotid artery affect the formation of low stress regions on the carotid walls. The main conclusions are: (a) larger internal carotid artery angles (theta(IC)) generally increase the frequency and the area of blood recirculation and lower the WSS on the sinus wall, hence increasing the risk of plaque build-up; (b) off-plane angles were found to lower the WSS on the sinus for geometries with theta(IC)25 degrees . Larger off-plane angles generally increase the danger of plague build-up; (c) for theta(IC) < 25 degrees , the off-plane angle does not have an obvious effect on the hemodynamic WSS; (d) symmetric bifurcations were found to increase the WSS on the sinus wall and ease the risk of vascular disease.     

Numerical investigation of the non-Newtonian pulsatile blood flow in a bifurcation model with a non-planar branch

The pulsatile flow of non-Newtonian fluid in a bifurcation model with a non-planar daughter branch is investigated numerically by using the Carreau-Yasuda model to take into account the shear thinning behavior of the analog blood fluid. The objective of this study is to deal with the influence of the non-Newtonian property of fluid and of out-of-plane curvature in the non-planar daughter vessel on wall shear stress (WSS), oscillatory shear index (OSI), and flow phenomena during the pulse cycle. The non-Newtonian property in the daughter vessels induces a flattened axial velocity profile due to its shear thinning behavior. The non-planarity deflects flow from the inner wall of the vessel to the outer wall and changes the distribution of WSS along the vessel, in particular in systole phase. Downstream of the bifurcation, the velocity profiles are shifted toward the flow divider, and low WSS and high shear stress temporal oscillations characterized by OSI occur on the outer wall region of the daughter vessels close to the bifurcation. Secondary motions become stronger with the addition of the out-of-plane curvature induced by the bending of the vessel, and the secondary flow patterns swirl along the non-planar daughter vessel. A significant difference between the non-Newtonian and the Newtonian pulsatile flow is revealed during the pulse cycle; however, reasonable agreement between the non-Newtonian and the rescaled Newtonian flow is found. Calculated results for the pulsatile flow support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Simulation of particle-hemodynamics in a partially occluded artery segment with implications to the initiation of microemboli and secondary stenoses.

Computational results of laminar incompressible blood-particle flow analyses in an axisymmetric artery segment with a smooth local area constriction of 75 percent have been presented. The flow input waveform was sinusoidal with a nonzero average. The non-Newtonian behavior of blood was simulated with a modified Quemada model, platelet concentrations were calculated with a drift-flux model, and monocyte trajectories were described and compared for both Newtonian and Quemada rheologies. Indicators of "disturbed flow" included the time-averaged wall shear stress (WSS), the oscillatory shear index (OSI), and the wall shear stress gradient (WSSG). Implications of the vortical flow patterns behind the primary stenosis to the formation of microemboli and downstream stenoses are as follows. Elevated platelet concentrations due to accumulation in recirculation zones mixed with thrombin and ADP complexes assumed to be released upstream in high wall shear stress regions, could form microemboli, which are convected downstream. Distinct near-wall vortices causing a local increase in the WSSG and OSI as well as blood-particle entrainment with possible wall deposition, indicate sites susceptible to the onset of an additional stenosis proximal to the initial geometric disturbance.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

MRI-based quantification of outflow boundary conditions for computational fluid dynamics of stenosed human carotid arteries

Accurate assessment of wall shear stress (WSS) is vital for studies on the pathogenesis of atherosclerosis. WSS distributions can be obtained by computational fluid dynamics (CFD) using patient-specific geometries and flow measurements. If patient-specific flow measurements are unavailable, in- and outflow have to be estimated, for instance by using Murray’s Law. It is currently unknown to what extent this law holds for carotid bifurcations, especially in cases where stenoses are involved. We performed flow measurements in the carotid bifurcation using phase-contrast MRI in patients with varying degrees of stenosis. An empirical relation between outflow and degree of area stenosis was determined and the outflow measurements were compared to estimations based on Murray’s Law. Furthermore, the influence of outflow conditions on the WSS distribution was studied.For bifurcations with an area stenosis smaller than 65%, the outflow ratio of the internal carotid artery (ICA) to the common carotid artery (CCA) was 0.62±0.12 while the outflow ratio of the external carotid artery (ECA) was 0.35±0.13. If the area stenosis was larger than 65%, the flow to the ICA decreased linearly to zero at 100% area stenosis. The empirical relation fitted the flow data well (R²=0.69), whereas Murray’s Law overestimated the flow to the ICA substantially for larger stenosis, resulting in an overestimation of the WSS. If patient-specific flow measurements of the carotid bifurcation are unavailable, estimation of the outflow ratio by the presented empirical relation will result in a good approximation of calculated WSS using CFD.     

Can temporal fluctuation in spatial wall shear stress gradient initiate a cerebral aneurysm? A proposed novel hemodynamic index,the gradient oscillatory number (GON)

We propose a new hemodynamic index for the initiation of a cerebral aneurysm, defined by the temporal fluctuations of tension/compression forces acting on endothelial cells. We employed a patient-specific geometry of a human internal carotid artery (ICA) with an aneurysm, and reconstructed the geometry of the ICA before aneurysm formation by artificially removing the aneurysm. We calculated the proposed hemodynamic index and five other hemodynamic indices (wall shear stress (WSS) at peak systole, time-averaged WSS, time-averaged spatial WSS gradient, oscillatory shear index (OSI), and potential aneurysm formation indicator (AFI)) for the geometry before aneurysm formation using a computational fluid dynamics technique. By comparing the distribution of each index at the location of aneurysm formation, we discussed the validity of each. The results showed that only the proposed hemodynamic index had a significant correlation with the location of aneurysm formation. Our findings suggest that the proposed index may be useful as a hemodynamic index for the initiation of cerebral aneurysms.     

Mis-sizing of stent promotes intimal hyperplasia: impact of endothelial shear and intramural stress

Stent can cause flow disturbances on the endothelium and compliance mismatch and increased stress on the vessel wall. These effects can cause low wall shear stress (WSS), high wall shear stress gradient (WSSG), oscillatory shear index (OSI), and circumferential wall stress (CWS), which may promote neointimal hyperplasia (IH). The hypothesis is that stent-induced abnormal fluid and solid mechanics contribute to IH. To vary the range of WSS, WSSG, OSI, and CWS, we intentionally mismatched the size of stents to that of the vessel lumen. Stents were implanted in coronary arteries of 10 swine. Intravascular ultrasound (IVUS) was used to size the coronary arteries and stents. After 4 wk of stent implantation, IVUS was performed again to determine the extent of IH. In conjunction, computational models of actual stents, the artery, and non-Newtonian blood were created in a computer simulation to yield the distribution of WSS, WSSG, OSI, and CWS in the stented vessel wall. An inverse relation (R(2) = 0.59, P < 0.005) between WSS and IH was found based on a linear regression analysis. Linear relations between WSSG, OSI, and IH were observed (R(2) = 0.48 and 0.50, respectively, P < 0.005). A linear relation (R(2) = 0.58, P < 0.005) between CWS and IH was also found. More statistically significant linear relations between the ratio of CWS to WSS (CWS/WSS), the products CWS × WSSG and CWS × OSI, and IH were observed (R(2) = 0.67, 0.54, and 0.56, respectively, P < 0.005), suggesting that both fluid and solid mechanics influence the extent of IH. Stents create endothelial flow disturbances and intramural wall stress concentrations, which correlate with the extent of IH formation, and these effects were exaggerated with mismatch of stent/vessel size. These findings reveal the importance of reliable vessel and stent sizing to improve the mechanics on the vessel wall and minimize IH.     

Wall shear stress and near-wall flows in the stenosed femoral artery

Stenotic artery hemodynamics are often characertised by metrics including oscillatory shear index (OSI) and residence time (RT). This analysis was conducted to clarify the link between the near-wall flow behaviour and these resultant flow metrics. A computational simulation was conducted of a stenosed femoral artery, with an idealised representative geometry and a physiologically realistic inlet profile. The overall flow behaviour was characterised through consideration of the axial flow, which was non-dimensionalised against mean flow velocity. The OSI and RT metrics, which are a useful indicator of likely atherosclerotic sites, were explained through a discussion of the WSS values at different time points, the velocity behaviour and velocity profiles, with a particular focus on the near-wall behaviour which influences wall shear stress and the transient evolution of the wall shear stress. While, the stenosis throat experiences high values of wall shear stress, the smooth flow through this contracted region results in low variation in wall shear stress vectors and limited opportunity for any particle stasis. However, regions were noted distal and proximal (though to a lesser extent), where the change in recirculation zones over the cycle created highly elevated regions of both OSI and RT.     

Hemodynamics and wall mechanics in human carotid bifurcation and its consequences for atherogenesis: investigation of inter-individual variation

Finite element simulations of fluid-solid interactions were used to investigate inter-individual variations in flow dynamics and wall mechanics at the carotid artery bifurcation, and its effects on atherogenesis, in three healthy humans (normal volunteers: NV1, NV2, NV4). Subject-specific calculations were based on MR images of structural anatomy and ultrasound measurements of flow at domain boundaries. For all subjects, the largest contiguous region of low wall shear stress (WSS) occurred at the carotid bulb, WSS was high (6-10 Pa) at the apex, and a small localized region of WSS > 10 Pa occurred close to the inner wall of the external carotid artery (ECA). NV2 and NV4 had a "spot" of low WSS distal to the bifurcation at the inner wall of the ECA. Low WSS patches in the common carotid artery (CCA) were contiguous with the carotid bulb low WSS region in NV1 and NV2, but not in NV4. In all three subjects, areas of high oscillatory shear index (OSI) were confined to regions of low WSS. Only NV4 exhibited high levels of OSI on the external adjoining wall of the ECA and CCA. For all subjects, the maximum wall shear stress temporal gradient (WSSTG) was highest at the flow divider (reaching 1,000 Pa/s), exceeding 300 Pa/s at the walls connecting the ECA and CCA, but remaining below 250 Pa/s outside of the ECA. In all subjects, (maximum principle) cyclic strain (CS) was greatest at the apex (NV1: 14%; NV2: 11%; NV4: 6%), and a second high CS region occurred at the ECA-CCA adjoining wall (NV1: 11%, NV2: 9%, NV4: 5%). Wall deformability was included in one simulation (NV2) to verify that it had little influence on the parameters studied. Location and magnitude of low WSS were similar, except for the apex (differences of up to 25%). Wall distensibility also influenced OSI, doubling it in most of the CCA, separating the single high OSI region of the carotid bulb into two smaller regions, and shrinking the ECA internal and external walls' high OSI regions. These observations provide further evidence that significant intra-subject variability exists in those factors thought to impact atherosclerosis.     

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73 Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.     

The effects of stenosis severity on the hemodynamic parameters-assessment of the correlation between stress phase angle and wall shear stress

To study the effects of increase in the degree of stenosis severity and subsequent complexity of hemodynamic patterns on hemodynamic parameters, experimental investigations and numerical simulations were performed. The correlations between the large negative Stress Phase Angle (SPA), the low mean Wall Shear Stress (WSS) and high Oscillatory Shear Index (OSI) were investigated at the distal shoulder and post-stenotic regions as the outcomes of elevated stenosis severity. Models included non-Newtonian fluid flow in stenotic arteries with 30-80% symmetrical stenoses. To study the interactions between pulsatile WSS and pulsatile wall circumferential stress (WCS) acting on endothelial cells, SPA as the phase difference between WSS and WCS waves was used. Moreover, the distribution of SPA on the lumen axis was compared to the distributions of the mean WSS and OSI that have been regarded until now as the determinants of atherosclerosis-prone regions. Results indicate that an increase in stenosis severity, not only affects the mean WSS, mean WCS and pulse amplitudes, but also influences the phase difference between them. The SPA is large negative on the distal shoulder and post-stenotic areas where atherosclerotic plaque develops. The increasing stenosis severity and the subsequent increasing complexity of hemodynamic patterns affect the correlation between any of the low mean WSS and high OSI with large negative SPA, such that it not only leads to create and develop some regions where the correlation between any of the low mean WSS and high OSI with large negative SPA is well but also leads to create and develop other regions where such correlations fail.     

Direct numerical simulation of transitional flow in a stenosed carotid bifurcation

The blood flow dynamics of a stenosed, subject-specific, carotid bifurcation were numerically simulated using the spectral element method. Pulsatile inlet conditions were based on in vivo color Doppler ultrasound measurements of blood velocity. The results demonstrated the transitional or weakly turbulent state of the blood flow, which featured rapid velocity and pressure fluctuations in the post-stenotic region of the internal carotid artery (ICA) during systole and laminar flow during diastole. High-frequency vortex shedding was greatest downstream of the stenosis during the deceleration phase of systole. Velocity fluctuations had a frequency within the audible range of 100-300Hz. Instantaneous wall shear stress (WSS) within the stenosis was relatively high during systole ( approximately 25-45Pa) compared to that in a healthy carotid. In addition, high spatial gradients of WSS were present due to flow separation on the inner wall. Oscillatory flow reversal and low pressure were observed distal to the stenosis in the ICA. This study predicts the complex flow field, the turbulence levels and the distribution of the biomechanical stresses present in vivo within a stenosed carotid artery.     

Cerebral perfusion changes in patients with internal carotid artery stenoses after surgical revascularization

Changes of cerebral perfusion and the condition of collateral blood supply in patients with internal carotid artery stenoses may have a prognostic value for effective blood flow restoration after revascularization of the internal carotid arteries (ICAs). To determine the patterns of cerebral perfusion changes in patients with ICA stenoses before and after surgical treatment, a clinical CT perfusion study of 41 patients with moderate to severe ICA stenoses was performed. Perfusion CT (PCT) had been conducted in 17 patients with moderate ICA stenoses (50–69%) and in 24 patients with severe ICA stenoses (70–99%) 3 times: before intervention (balloon angioplasty with stenting or carotid endarterectomy), on the 3rd to 7th day, and within 1 to 3 months after surgery. Scanning was performed at the level of the basal ganglia and semioval centers. In patients without ICA stenosis (control group of 39 individuals), PCT was conducted once. We found that surgical recanalization of the ICA leads to normalization of the cerebral blood flow parameters in the perfusion area of the middle cranial artery, as evidenced by a decrease of MTT and CBV and an increase of CBF to values comparable to those in the control group. However, blood flow restoration in the anterior and posterior watershed areas, which are known to be mostly affected under chronic hypoperfusion conditions, was observed only in patients with a complete circle of Willis (CW) and moderate ICA stenosis. Therefore, severe stenosis (>70%) and the incomplete CW are the prognostic factors for inadequate blood flow restoration after revascularization in patients with ICA stenoses.     

Influence of head position on carotid hemodynamics in young adults

Studies in adults have shown marked changes in geometry and relative positions of the carotid arteries when rotating the head. The aim of this study was to quantify the change in geometry and analyze its effect on carotid hemodynamics as a result of head rotation. The right carotid arteries of nine young adult subjects were investigated in supine position with straight and left turned head positions, respectively. The three-dimensional (3D) carotid geometry was reconstructed by using 3D ultrasound (3D US), and the carotid hemodynamics were calculated by combining 3D US with computational fluid dynamics. It was observed that cross-sectional areas and shapes did not change markedly with head rotation, but carotid vessel center lines altered with planarification of the common carotid artery as a main feature (P < 0.05). Measured common carotid flow rates changed significantly at the individual level when the head was turned, but on the average, the change in mean common carotid flow rate was relatively small (0.37 +/- 1.11 ml/s). The effect of the altered center lines and flow rates on the atherogenic nature of the carotid bifurcation was evaluated by using calculated hemodynamic wall parameters, such as wall shear stress (WSS) and oscillatory shear index (OSI). It was found that WSS and OSI patterns changed significantly with head rotation, but the variations were very subject dependent and could not have been predicted without assessing the altered geometry and flow of the carotid bifurcation for individual cases. This study suggests that there is a need for standardization of the choice of head position in the 3D US scan protocol, and that carotid stents and emboli diverters should be studied in different head positions.     

Flow patterns in three-dimensional porcine epicardial coronary arterial tree

The branching pattern of epicardial coronary arteries is clearly three-dimensional, with correspondingly complex flow patterns. The objective of the present study was to perform a detailed hemodynamic analysis using a three-dimensional finite element method in a left anterior descending (LAD) epicardial arterial tree, including main trunk and primary branches, based on computed tomography scans. The inlet LAD flow velocity was measured in an anesthetized pig, and the outlet pressure boundary condition was estimated based on scaling laws. The spatial and temporal wall shear stress (WSS), gradient of WSS (WSSG), and oscillatory shear index (OSI) were calculated and used to identify regions of flow disturbances in the vicinity of primary bifurcations. We found that low WSS and high OSI coincide with disturbed flows (stagnated, secondary, and reversed flows) opposite to the flow divider and lateral to the junction orifice of the main trunk and primary branches. High time-averaged WSSG occurs in regions of bifurcations, with the flow divider having maximum values. Low WSS and high OSI were found to be related through a power law relationship. Furthermore, zones of low time-averaged WSS and high OSI amplified for larger diameter ratio and high inlet flow rate. Hence, different focal atherosclerotic-prone regions may be explained by different physical mechanism associated with certain critical levels of low WSS, high OSI, and high WSSG, which are strongly affected by the diameter ratio. The implications of the flow patterns for atherogenesis are enumerated.     

Computational study of pulsatile blood flow in prototype vessel geometries of coronary segments

The spatial and temporal distributions of wall shear stress (WSS) in prototype vessel geometries of coronary segments are investigated via numerical simulation, and the potential association with vascular disease and specifically atherosclerosis and plaque rupture is discussed. In particular, simulation results of WSS spatio-temporal distributions are presented for pulsatile, non-Newtonian blood flow conditions for: (a) curved pipes with different curvatures, and (b) bifurcating pipes with different branching angles and flow division. The effects of non-Newtonian flow on WSS (compared to Newtonian flow) are found to be small at Reynolds numbers representative of blood flow in coronary arteries. Specific preferential sites of average low WSS (and likely atherogenesis) were found at the outer regions of the bifurcating branches just after the bifurcation, and at the outer-entry and inner-exit flow regions of the curved vessel segment. The drop in WSS was more dramatic at the bifurcating vessel sites (less than 5% of the pre-bifurcation value). These sites were also near rapid gradients of WSS changes in space and time – a fact that increases the risk of rupture of plaque likely to develop at these sites. The time variation of the WSS spatial distributions was very rapid around the start and end of the systolic phase of the cardiac cycle, when strong fluctuations of intravascular pressure were also observed. These rapid and strong changes of WSS and pressure coincide temporally with the greatest flexion and mechanical stresses induced in the vessel wall by myocardial motion (ventricular contraction). The combination of these factors may increase the risk of plaque rupture and thrombus formation at these sites.     

On the relative importance of rheology for image-based CFD models of the carotid bifurcation

BACKGROUND: Patient-specific computational fluid dynamics (CFD) models derived from medical images often require simplifying assumptions to render the simulations conceptually or computationally tractable. In this study, we investigated the sensitivity of image-based CFD models of the carotid bifurcation to assumptions regarding the blood rheology. METHOD OF APPROACH: CFD simulations of three different patient-specific models were carried out assuming: a reference high-shear Newtonian viscosity, two different non-Newtonian (shear-thinning) rheology models, and Newtonian viscosities based on characteristic shear rates or, equivalently, assumed hematocrits. Sensitivity of wall shear stress (WSS) and oscillatory shear index (OSI) were contextualized with respect to the reproducibility of the reconstructed geometry, and to assumptions regarding the inlet boundary conditions. RESULTS: Sensitivity of WSS to the various rheological assumptions was roughly 1.0 dyn/cm(2) or 8%, nearly seven times less than that due to geometric uncertainty (6.7 dyn/cm(2) or 47%), and on the order of that due to inlet boundary condition assumptions. Similar trends were observed regarding OSI sensitivity. Rescaling the Newtonian viscosity based on time-averaged inlet shear rate served to approximate reasonably, if overestimate slightly, non-Newtonian behavior. CONCLUSIONS: For image-based CFD simulations of the normal carotid bifurcation, the assumption of constant viscosity at a nominal hematocrit is reasonable in light of currently available levels of geometric precision, thus serving to obviate the need to acquire patient-specific rheological data.     

Can sites prone to flow induced vascular complications in a-v fistulas be assessed using computational fluid dynamics?

Arterio-venous fistulas (shunts between arteries and veins) are the preferred vascular access for hemodialysis. Despite their superior patency, compared with synthetic tubes and grafts, functional problems and inadequate flow rates are the common complications. Local flow conditions, in particular low and oscillating wall shear stresses (WSS), are central to vascular problems and a robust framework for analyzing flow conditions in vascular structures could provide an understanding of the mechanisms leading to vascular complications, such as stenoses, aneurisms, and thromboses.We hypothesize that a validated computational fluid dynamics (CFD) framework can be used to identify critical fistula configurations with elevated risk of complications. Therefore, the aim of the present study was to develop a CFD framework for analyzing fluid flow in complex vascular structures, such as arterio-venous fistulas validated by comparisons of in vitro volume flows with CFD results and flow fields from ultrasound scans with CFD simulations.Volume flows measured in vitro and CFD data differed quantitatively. However, good relative correlations exist between the data using logarithmic scales. Qualitatively, visual comparisons between ultrasound and CFD images showed good agreement between the two methods. In addition, WSS levels and the oscillatory shear index (OSI) were calculated and visualized on the model surface. The method was successfully validated and the method is deemed suitable for more thorough investigations into the field of vascular complications in a-v fistulas.     

Numerical investigation of physiologically realistic pulsatile flow through arterial stenosis

Numerical simulations of pulsatile blood flow in straight tube stenosis models were performed to investigate the poststenotic flow phenomena. In this study, three axisymmetrical and three asymmetrical stenosis models with area reduction of 25%, 50% and 75% were constructed. A measured human common carotid artery blood flow waveform was used as the upstream flow condition which has a mean Reynold's number of 300. All calculations were performed with high spatial and temporal resolutions. Flow features such as velocity profiles, flow separation zone (FSZ), and wall shear stress (WSS) distributions in the poststenotic region for all models are presented. The results have demonstrated that the formation and development of FSZs in the poststenotic region are very complex, especially in the flow deceleration phase. In axisymmetric stenoses the poststenotic flow is more sensitive to changes in the degree of stenosis than in asymmetric models. For severe stenoses, the stenosis influence length is shorter in asymmetrical models than in axisymmetrical cases. WSS oscillations (between positive and negative values) have been observed at various downstream locations in some models. The amplitude of the oscillation depends strongly on the axial location and the degree of stenosis.