The effect of inlet waveforms on computational hemodynamics of patient-specific intracranial aneurysms

Due to the lack of patient-specific inlet flow waveform measurements, most computational fluid dynamics (CFD) simulations of intracranial aneurysms usually employ waveforms that are not patient-specific as inlet boundary conditions for the computational model. The current study examined how this assumption affects the predicted hemodynamics in patient-specific aneurysm geometries. We examined wall shear stress (WSS) and oscillatory shear index (OSI), the two most widely studied hemodynamic quantities that have been shown to predict aneurysm rupture, as well as maximal WSS (MWSS), energy loss (EL) and pressure loss coefficient (PLc). Sixteen pulsatile CFD simulations were carried out on four typical saccular aneurysms using 4 different waveforms and an identical inflow rate as inlet boundary conditions. Our results demonstrated that under the same mean inflow rate, different waveforms produced almost identical WSS distributions and WSS magnitudes, similar OSI distributions but drastically different OSI magnitudes. The OSI magnitude is correlated with the pulsatility index of the waveform. Furthermore, there is a linear relationship between aneurysm-averaged OSI values calculated from one waveform and those calculated from another waveform. In addition, different waveforms produced similar MWSS, EL and PLc in each aneurysm. In conclusion, inlet waveform has minimal effects on WSS, OSI distribution, MWSS, EL and PLc and a strong effect on OSI magnitude, but aneurysm-averaged OSI from different waveforms has a strong linear correlation with each other across different aneurysms, indicating that for the same aneurysm cohort, different waveforms can consistently stratify (rank) OSI of aneurysms.     

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.     

Inlet conditions for image-based CFD models of the carotid bifurcation: is it reasonable to assume fully developed flow?

BACKGROUND: Computational fluid dynamics tools are useful for their ability to model patient specific data relevant to the genesis and progression of atherosclerosis, but unavailable to measurement tools. The sensitivity of the physiologically relevant parameters of wall shear stress (WSS) and the oscillatory shear index (OSI) to secondary flow in the inlet velocity profiles was investigated in three realistic models of the carotid bifurcation. METHOD OF APPROACH: Secondary flow profiles were generated using sufficiently long entrance lengths, to which curvature and helical pitch were added. The differences observed were contextualized with respect to effect of the uncertainty of the models' geometry on the same parameters. RESULTS: The effects of secondary velocities in the inlet profile on WSS and OSI break down within a few diameters of the inlet. Overall, the effect of secondary inlet flow on these models was on average more than 3.5 times smaller than the effect of geometric variability, with 13% and 48% WSS variability induced by inlet secondary flow and geometric differences, respectively. CONCLUSIONS: The degree of variation is demonstrated to be within the range of the other computational assumptions, and we conclude that given a sufficient entrance length of realistic geometry, simplification to fully developed axial (i.e., Womersley) flow may be made without penalty. Thus, given a choice between measuring three components of inlet velocity or a greater geometric extent, we recommend effort be given to more accurate and detailed geometric reconstructions, as being of primary influence on physiologically significant indicators.     

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.     

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.     

Wall shear stresses remain elevated in mature arteriovenous fistulas: a case study

Maintaining vascular access (VA) patency continues to be the greatest challenge for dialysis patients. VA dysfunction, primarily due to venous neointimal hyperplasia development and stenotic lesion formation, is mainly attributed to complex hemodynamics within the arteriovenous fistula (AVF). The effect of VA creation and the subsequent geometrical remodeling on the hemodynamics and shear forces within a mature patient-specific AVF is investigated. A 3D reconstructed geometry of a healthy vein and a fully mature patient-specific AVF was developed from a series of 2D magnetic resonance image scans. A previously validated thresholding technique for region segmentation and lumen cross section contour creation was conducted in MIMICS 10.01, allowing for the creation of a 3D reconstructed geometry. The healthy vein and AVF computational models were built, subdivided, and meshed in GAMBIT 2.3. The computational fluid dynamic (CFD) code FLUENT 6.3.2 (Fluent Inc., Lebanon, NH) was employed as the finite volume solver to determine the hemodynamics and shear forces within the healthy vein and patient-specific AVF. Geometrical alterations were evaluated and a CFD analysis was conducted. Substantial geometrical remodeling was observed, following VA creation with an increase in cross-sectional area, out of plane curvature (maximum angle of curvature in AVF=30?deg), and angle of blood flow entry. The mean flow velocity entering the vein of the AVF is dramatically increased. These factors result in complex three-dimensional hemodynamics within VA junction (VAJ) and efferent vein of the AVF. Complex flow patterns were observed and the maximum and mean wall shear stress (WSS) magnitudes are significantly elevated. Flow reversal was found within the VAJ and efferent vein. Extensive geometrical remodeling during AVF maturation does not restore physiological hemodynamics to the VAJ and venous conduit of the AVF, and high WSS and WSS gradients, and flow reversal persist. It is theorized that the vessel remodelling and the continued non-physiological hemodynamics within the AVF compound to result in stenotic lesion development.     

Approximating hemodynamics of cerebral aneurysms with steady flow simulations

Computational fluid dynamics (CFD) simulations can be employed to gain a better understanding of hemodynamics in cerebral aneurysms and improve diagnosis and treatment. However, introduction of CFD techniques into clinical practice would require faster simulation times. The aim of this study was to evaluate the use of computationally inexpensive steady flow simulations to approximate the aneurysm's wall shear stress (WSS) field. Two experiments were conducted. Experiment 1 compared for two cases the time-averaged (TA), peak systole (PS) and end diastole (ED) WSS field between steady and pulsatile flow simulations. The flow rate waveform imposed at the inlet was varied to account for variations in heart rate, pulsatility index, and TA flow rate. Consistently across all flow rate waveforms, steady flow simulations accurately approximated the TA, but not the PS and ED, WSS field. Following up on experiment 1, experiment 2 tested the result for the TA WSS field in a larger population of 20 cases covering a wide range of aneurysm volumes and shapes. Steady flow simulations approximated the space-averaged WSS with a mean error of 4.3%. WSS fields were locally compared by calculating the absolute error per node of the surface mesh. The coefficient of variation of the root-mean-square error over these nodes was on average 7.1%. In conclusion, steady flow simulations can accurately approximate the TA WSS field of an aneurysm. The fast computation time of 6 min per simulation (on 64 processors) could help facilitate the introduction of CFD into clinical practice.     

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 and wall shear stress distribution in human internal carotid arteries: the geometric effect on the risk for stenoses

It has been widely observed that atherosclerotic stenosis occurs at sites with complex hemodynamics, such as arteries with high curvature or bifurcations. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present study, 3D sinusoidally pulsatile blood flow through the models of internal carotid artery (ICA) with different geometries was investigated with computational simulation. Three preferred sites of stenoses were found along the carotid siphon with low and highly oscillatory WSS. The risk for stenoses at these sites was scaled with the values of time-averaged WSS and oscillating shear index (OSI). The local risk for stenoses at every preferred site of stenoses was found different between 3 types of ICA, indicating that the geometry of the blood vessel plays significant roles in the atherogenesis. Specifically, the large curvature and planarity of the vessel were found to increase the risk for stenoses, because they tend to lower WSS and elevate OSI. Therefore, the geometric study makes it possible to estimate the stenosis location in the ICA siphon as long as the shape of ICA was measured.     

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.     

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.     

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.     

Influence of model boundary conditions on blood flow patterns in a patient specific stenotic right coronary artery

Background

In literature, the effect of the inflow boundary condition was investigated by examining the impact of the waveform and the shape of the spatial profile of the inlet velocity on the cardiac hemodynamics. However, not much work has been reported on comparing the effect of the different combinations of the inlet/outlet boundary conditions on the quantification of the pressure field and flow distribution patterns in stenotic right coronary arteries.

Method

Non-Newtonian models were used to simulate blood flow in a patient-specific stenotic right coronary artery and investigate the influence of different boundary conditions on the phasic variation and the spatial distribution patterns of blood flow. The 3D geometry of a diseased artery segment was reconstructed from a series of IVUS slices. Five different combinations of the inlet and the outlet boundary conditions were tested and compared.

Results

The temporal distribution patterns and the magnitudes of the velocity, the wall shear stress (WSS), the pressure, the pressure drop (PD), and the spatial gradient of wall pressure (WPG) were different when boundary conditions were imposed using different pressure/velocity combinations at inlet/outlet. The maximum velocity magnitude in a cardiac cycle at the center of the inlet from models with imposed inlet pressure conditions was about 29% lower than that from models using fully developed inlet velocity data. Due to the fact that models with imposed pressure conditions led to blunt velocity profile, the maximum wall shear stress at inlet in a cardiac cycle from models with imposed inlet pressure conditions was about 29% higher than that from models with imposed inlet velocity boundary conditions. When the inlet boundary was imposed by a velocity waveform, the models with different outlet boundary conditions resulted in different temporal distribution patterns and magnitudes of the phasic variation of pressure. On the other hand, the type of different boundary conditions imposed at the inlet and the outlet did not have significant effect on the spatial distribution patterns of the PD, the WPG and the WSS on the lumen surface, regarding the locations of the maximum and the minimum of each quantity.

Conclusions

The observations from this study indicated that the ways how pressure and velocity boundary conditions are imposed in computational models have considerable impact on flow velocity and shear stress predictions. Accuracy of in vivo measurements of blood pressure and velocity is of great importance for reliable model predictions.

     

The influence of inlet velocity profile on predicted flow in type B aortic dissection

In order for computational fluid dynamics to provide quantitative parameters to aid in the clinical assessment of type B aortic dissection, the results must accurately mimic the hemodynamic environment within the aorta. The choice of inlet velocity profile (IVP) therefore is crucial; however, idealised profiles are often adopted, and the effect of IVP on hemodynamics in a dissected aorta is unclear. This study examined two scenarios with respect to the influence of IVP—using (a) patient-specific data in the form of a three-directional (3D), through-plane (TP) or flat IVP; and (b) non-patient-specific flow waveform. The results obtained from nine simulations using patient-specific data showed that all forms of IVP were able to reproduce global flow patterns as observed with 4D flow magnetic resonance imaging. Differences in maximum velocity and time-averaged wall shear stress near the primary entry tear were up to 3% and 6%, respectively, while pressure differences across the true and false lumen differed by up to 6%. More notable variations were found in regions of low wall shear stress when the primary entry tear was close to the left subclavian artery. The results obtained with non-patient-specific waveforms were markedly different. Throughout the aorta, a 25% reduction in stroke volume resulted in up to 28% and 35% reduction in velocity and wall shear stress, respectively, while the shape of flow waveform had a profound influence on the predicted pressure. The results of this study suggest that 3D, TP and flat IVPs all yield reasonably similar velocity and time-averaged wall shear stress results, but TP IVPs should be used where possible for better prediction of pressure. In the absence of patient-specific velocity data, effort should be made to acquire patient’s stroke volume and adjust the applied IVP accordingly.

     

Tetrahedral vs. polyhedral mesh size evaluation on flow velocity and wall shear stress for cerebral hemodynamic simulation

Haemodynamic factors, in particular wall shear stresses (WSSs) may have significant impact on growth and rupture of cerebral aneurysms. Without a means to measure WSS reliably in vivo, computational fluid dynamic (CFD) simulations are frequently employed to visualise and quantify blood flow from patient-specific computational models. With increasing interest in integrating these CFD simulations into pretreatment planning, a better understanding of the validity of the calculations in respect to computation parameters such as volume element type, mesh size and mesh composition is needed. In this study, CFD results for the two most common aneurysm types (saccular and terminal) are compared for polyhedral- vs. tetrahedral-based meshes and discussed regarding future clinical applications. For this purpose, a set of models were constructed for each aneurysm with spatially varying surface and volume mesh configurations (mesh size range: 5119-258, 481 volume elements). WSS distribution on the model wall and point-based velocity measurements were compared for each configuration model. Our results indicate a benefit of polyhedral meshes in respect to convergence speed and more homogeneous WSS patterns. Computational variations of WSS values and blood velocities are between 0.84 and 6.3% from the most simple mesh (tetrahedral elements only) and the most advanced mesh design investigated (polyhedral mesh with boundary layer).     

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.     

Circumferential vascular deformation after stent implantation alters wall shear stress evaluated with time-dependent 3D computational fluid dynamics models.

The success of vascular stents in the restoration of blood flow is limited by restenosis. Recent data generated from computational fluid dynamics (CFD) models suggest that stent geometry may cause local alterations in wall shear stress (WSS) that have been associated with neointimal hyperplasia and subsequent restenosis. However, previous CFD studies have ignored histological evidence of vascular straightening between circumferential stent struts. We tested the hypothesis that consideration of stent-induced vascular deformation may more accurately predict alterations in indexes of WSS that may subsequently account for histological findings after stenting. We further tested the hypothesis that the severity of these alterations in WSS varies with the degree of vascular deformation after implantation. Steady-state and time-dependent simulations of three-dimensional CFD arteries based on canine coronary artery measurements of diameter and blood flow were conducted, and WSS and WSS gradients were calculated. Circumferential straightening introduced areas of high WSS between stent struts that were absent in stented vessels of circular cross section. The area of vessel exposed to low WSS was dependent on the degree of circumferential vascular deformation and axial location within the stent. Stents with four vs. eight struts increased the intrastrut area of low WSS in vessels, regardless of cross-sectional geometry. Elevated WSS gradients were also observed between struts in vessels with polygonal cross sections. The results obtained using three-dimensional CFD models suggest that changes in vascular geometry after stent implantation are important determinants of WSS distributions that may be associated with subsequent neointimal hyperplasia.     

Errors in power-law estimations of inflow rates for intracranial aneurysm CFD

Patient-specific inflow rates are rarely available for computational fluid dynamics (CFD) studies of intracranial aneurysms. Instead, inflow rates are often estimated from parent artery diameters via power laws, i.e. Q ∝ Dⁿ, reflecting adaptation of conduit arteries to demanded flow. The present study aimed to validate the accuracy of these power laws. Internal carotid artery (ICA) flow rates were measured from 25 ICA aneurysm patients via 2D phase contrast MRI. ICA diameters, derived from 3D segmentation of rotational angiograms, were used to estimate inflow rates via power laws from the aneurysm CFD literature assuming the same inlet wall shear stress (WSS) (n = 3), velocity (n = 2) or flow rate (n = 0) for all cases. To illustrate the potential impact of errors in flow rate estimates, pulsatile CFD was carried out for four cases having large errors for at least one power law. Flow rates estimated by n = 3 and n = 0 power laws had significant (p < 0.01) mean biases of −22% to +32%, respectively, but with individual errors ranging from −78% to +120%. The n = 2 power law had no significant bias, but had non-negligible individual errors of −58% to +71%. CFD showed similarly large errors for time-averaged sac WSS; however, these were reduced after normalizing by parent artery WSS. High frequency WSS fluctuations, evident in 2/4 aneurysms, were also sensitive to inflow rate errors. Care should therefore be exercised in the interpretation of aneurysm CFD studies that rely on power law estimates of inflow rates, especially if absolute (vs. normalized) WSS, or WSS instabilities, are of interest.     

MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models

Pulsatile flow was studied in physiologically realistic models of a normal and a moderately stenosed (30% diameter reduction) human carotid bifurcation. Time-resolved velocity measurements were made using magnetic resonance imaging, from which wall shear stress (WSS) vectors were calculated. Velocity measurements in the inflow and outflow regions were also used as boundary conditions for a computational fluid dynamics (CFD) model. Experimental flow patterns and derived WSS vectors were compared qualitatively with the corresponding CFD predictions. In the stenosed phantom, flow in the bulb region of the "internal carotid artery" was concentrated along the outer wall, with a region of low and recirculating flow near the inner wall. In the normal phantom, the converse was found, with a low flow region near the outer wall of the bulb. Time-averaged WSS and oscillatory shear index were also markedly different for the two phantoms.