Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm

In experiments turbulence has previously been shown to occur in intracranial aneurysms. The effects of turbulence induced oscillatory wall stresses could be of great importance in understanding aneurysm rupture. To investigate the effects of turbulence on blood flow in an intracranial aneurysm, we performed a high resolution computational fluid dynamics (CFD) simulation in a patient specific middle cerebral artery (MCA) aneurysm using a realistic, pulsatile inflow velocity. The flow showed transition to turbulence just after peak systole, before relaminarization occurred during diastole. The turbulent structures greatly affected both the frequency of change of wall shear stress (WSS) direction and WSS magnitude, which reached a maximum value of 41.5Pa. The recorded frequencies were predominantly in the range of 1-500Hz. The current study confirms, through properly resolved CFD simulations that turbulence can occur in intracranial aneurysms.     

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.     

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.     

Wall shear stress at the initiation site of cerebral aneurysms

Hemodynamics are believed to play an important role in the initiation of cerebral aneurysms. In particular, studies have focused on wall shear stress (WSS), which is a key regulator of vascular biology and pathology. In line with the observation that aneurysms predominantly occur at regions of high WSS, such as bifurcation apices or outer walls of vascular bends, correlations have been found between the aneurysm initiation site and high WSS. The aim of our study was to analyze the WSS field at an aneurysm initiation site that was neither a bifurcation apex nor the outer wall of a vascular bend. Ten cases with aneurysms on the A1 segment of the anterior cerebral artery were analyzed and compared with ten controls. Aneurysms were virtually removed from the vascular models of the cases to mimic the pre-aneurysm geometry. Computational fluid dynamics (CFD) simulations were created to assess the magnitude, gradient, multidirectionality, and pulsatility of the WSS. To aid the inter-subject comparison of hemodynamic variables, we mapped the branch surfaces onto a two-dimensional parametric space. This approach made it possible to view the whole branch at once for qualitative evaluation. It also allowed us to empirically define a patch for quantitative analysis, which was consistent among subjects and encapsulated the aneurysm initiation sites in our dataset. To test the sensitivity of our results, CFD simulations were repeated with a second independent observer virtually removing the aneurysms and with a 20 % higher flow rate at the inlet. We found that branches harboring aneurysms were characterized by high WSS and high WSS gradients. Among all assessed variables, the aneurysm initiation site most consistently coincided with peaks of temporal variation in the WSS magnitude.     

Requirements for mesh resolution in 3D computational hemodynamics

Computational techniques are widely used for studying large artery hemodynamics. Current trends favor analyzing flow in more anatomically realistic arteries. A significant obstacle to such analyses is generation of computational meshes that accurately resolve both the complex geometry and the physiologically relevant flow features. Here we examine, for a single arterial geometry, how velocity and wall shear stress patterns depend on mesh characteristics. A well-validated Navier-Stokes solver was used to simulate flow in an anatomically realistic human right coronary artery (RCA) using unstructured high-order tetrahedral finite element meshes. Velocities, wall shear stresses (WSS), and wall shear stress gradients were computed on a conventional "high-resolution" mesh series (60,000 to 160,000 velocity nodes) generated with a commercial meshing package. Similar calculations were then performed in a series of meshes generated through an adaptive mesh refinement (AMR) methodology. Mesh-independent velocity fields were not very difficult to obtain for both the conventional and adaptive mesh series. However, wall shear stress fields, and, in particular, wall shear stress gradient fields, were much more difficult to accurately resolve. The conventional (nonadaptive) mesh series did not show a consistent trend towards mesh-independence of WSS results. For the adaptive series, it required approximately 190,000 velocity nodes to reach an r.m.s. error in normalized WSS of less than 10 percent. Achieving mesh-independence in computed WSS fields requires a surprisingly large number of nodes, and is best approached through a systematic solution-adaptive mesh refinement technique. Calculations of WSS, and particularly WSS gradients, show appreciable errors even on meshes that appear to produce mesh-independent velocity fields.     

Influence of non-Newtonian properties of blood on the wall shear stress in human atherosclerotic right coronary arteries

The objective of this work is to investigate the effect of non-Newtonian properties of blood on the wall shear stress (WSS) in atherosclerotic coronary arteries using both Newtonian and non-Newtonian models. Numerical simulations were performed to examine how the spatial and temporal WSS distributions are influenced by the stenosis size, blood viscosity, and flow rate. The computational results demonstrated that blood viscosity properties had considerable effect on the magnitude of the WSS, especially where disturbed flow was observed. The WSS distribution is highly non-uniform both temporally and spatially, especially in the stenotic region. The maximum WSS occurred at the proximal side of the stenosis, near the outer wall in the curved artery with no stenosis. The lumen area near the inner wall distal to the stenosis region experienced a lower WSS during the entire cardiac cycle. Among the factors of stenosis size, blood viscosity, and flow rate, the size of the stenosis has the most significant effect on the spatial and temporal WSS distributions qualitatively and quantitatively.     

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.     

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.     

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.     

Approaches for treatment of aortic arch aneurysm,a numerical study

Aortic arch aneurysm is a complex pathology which requires coverage of one or more aortic arch vessels. In this study we explore the hemodynamic behavior of the aortic arch in aneurysmatic and treated cases with three currently available treatment approaches: Surgery Graft, hybrid Stent-Graft and chimney Stent Graft. The analysis included four models of the time-dependent fluid domains of aneurysmatic arch and of the surgery, hybrid and chimney endovascular techniques. Dimensions of the models are based on typical anatomy, and boundary conditions are based on typical physiological flow.The simulations used computational fluid dynamics (CFD) methods to delineate the time-dependent flow dynamics in the four geometric models.Results of velocity vectors, flow patterns, blood pressure and wall shear stress distributions are presented.The results delineate disturbed and recirculating flow in the aortic arch aneurysm accompanied with low wall shear stress and velocities, compared to a uniformly directed flow and nominal wall shear stress (WSS) in the model of Surgery graft. Out of the two endograft procedures, the hybrid procedure clearly exhibits better hemodynamic performances over the chimney model, with lower WSS, lower pressure drop and less disturbed and vortical flow regions. Although the chimney procedure requires less manufacturing time and cost, it is associated with higher risk rates, and therefore, it is recommended only for emergency cases. This study may shed light on the hemodynamic factors for these complications and provide insight into ways to improve the procedure.     

Optimization of cardiovascular stent design using computational fluid dynamics

Coronary stent design affects the spatial distribution of wall shear stress (WSS), which can influence the progression of endothelialization, neointimal hyperplasia, and restenosis. Previous computational fluid dynamics (CFD) studies have only examined a small number of possible geometries to identify stent designs that reduce alterations in near-wall hemodynamics. Based on a previously described framework for optimizing cardiovascular geometries, we developed a methodology that couples CFD and three-dimensional shape-optimization for use in stent design. The optimization procedure was fully-automated, such that solid model construction, anisotropic mesh generation, CFD simulation, and WSS quantification did not require user intervention. We applied the method to determine the optimal number of circumferentially repeating stent cells (N(C)) for slotted-tube stents with various diameters and intrastrut areas. Optimal stent designs were defined as those minimizing the area of low intrastrut time-averaged WSS. Interestingly, we determined that the optimal value of N(C) was dependent on the intrastrut angle with respect to the primary flow direction. Further investigation indicated that stent designs with an intrastrut angle of approximately 40 deg minimized the area of low time-averaged WSS regardless of vessel size or intrastrut area. Future application of this optimization method to commercially available stent designs may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes.     

Hemodynamic Changes Caused by Flow Diverters in Rabbit Aneurysm Models: Comparison of Virtual and Realistic FD Deployments Based on Micro-CT Reconstruction

Adjusting hemodynamics via flow diverter (FD) implantation is emerging as a novel method of treating cerebral aneurysms. However, most previous FD-related hemodynamic studies were based on virtual FD deployment, which may produce different hemodynamic outcomes than realistic (in vivo) FD deployment. We compared hemodynamics between virtual FD and realistic FD deployments in rabbit aneurysm models using computational fluid dynamics (CFD) simulations. FDs were implanted for aneurysms in 14 rabbits. Vascular models based on rabbit-specific angiograms were reconstructed for CFD studies. Real FD configurations were reconstructed based on micro-CT scans after sacrifice, while virtual FD configurations were constructed with SolidWorks software. Hemodynamic parameters before and after FD deployment were analyzed. According to the metal coverage (MC) of implanted FDs calculated based on micro-CT reconstruction, 14 rabbits were divided into two groups (A, MC >35%; B, MC <35%). Normalized mean wall shear stress (WSS), relative residence time (RRT), inflow velocity, and inflow volume in Group A were significantly different (P<0.05) from virtual FD deployment, but pressure was not (P>0.05). The normalized mean WSS in Group A after realistic FD implantation was significantly lower than that of Group B. All parameters in Group B exhibited no significant difference between realistic and virtual FDs. This study confirmed MC-correlated differences in hemodynamic parameters between realistic and virtual FD deployment.     

Patient-specific computational haemodynamics: generation of structured and conformal hexahedral meshes from triangulated surfaces of vascular bifurcations

Measuring the blood flow is still limited by current imaging technologies and is generally overcome using computational fluid dynamics (CFD) which, because of the complex geometry of blood vessels, has widely relied on tetrahedral meshes. Hexahedral meshes offer more accurate results with lower-density meshes and faster computation as compared to tetrahedral meshes, but their use is limited by the far more complex mesh generation. We present a robust methodology for conformal and structured hexahedral mesh generation - applicable to complex arterial geometries as bifurcating vessels - starting from triangulated surfaces. Cutting planes are used to slice the lumen surface and to construct longitudinal Bezier splines. Afterwards, an isoparametric transformation is used to map a parametrically defined quadrilateral surface mesh into the vessel volume, resulting in stacks of sections which can then be used for sweeping. Being robust and open source based, this methodology may improve the current standard in patient-specific mesh generation and enhance the reliability of CFD to patient-specific haemodynamics.     

Patient-specific computational haemodynamics: generation of structured and conformal hexahedral meshes from triangulated surfaces of vascular bifurcations

Measuring the blood flow is still limited by current imaging technologies and is generally overcome using computational fluid dynamics (CFD) which, because of the complex geometry of blood vessels, has widely relied on tetrahedral meshes. Hexahedral meshes offer more accurate results with lower-density meshes and faster computation as compared to tetrahedral meshes, but their use is limited by the far more complex mesh generation. We present a robust methodology for conformal and structured hexahedral mesh generation – applicable to complex arterial geometries as bifurcating vessels – starting from triangulated surfaces. Cutting planes are used to slice the lumen surface and to construct longitudinal Bezier splines. Afterwards, an isoparametric transformation is used to map a parametrically defined quadrilateral surface mesh into the vessel volume, resulting in stacks of sections which can then be used for sweeping. Being robust and open source based, this methodology may improve the current standard in patient-specific mesh generation and enhance the reliability of CFD to patient-specific haemodynamics.     

Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations

Aortic dissecting aneurysm is one of the most catastrophic cardiovascular emergencies that carries high mortality. It was pointed out from clinical observations that the aneurysm development is likely to be related to the hemodynamics condition of the dissected aorta. In order to gain more insight on the formation and progression of dissecting aneurysm, hemodynamic parameters including flow pattern, velocity distribution, aortic wall pressure and shear stress, which are difficult to measure in vivo, are evaluated using numerical simulations. Pulsatile blood flow in patient-specific dissecting aneurismal aortas before and after the formation of lumenal aneurysm (pre-aneurysm and post-aneurysm) is investigated by computational fluid dynamics (CFD) simulations. Realistic time-dependent boundary conditions are prescribed at various arteries of the complete aorta models. This study suggests the helical development of false lumen around true lumen may be related to the helical nature of hemodynamic flow in aorta. Narrowing of the aorta is responsible for the massive recirculation in the poststenosis region in the lumenal aneurysm development. High pressure difference of 0.21 kPa between true and false lumens in the pre-aneurismal aorta infers the possible lumenal aneurysm site in the descending aorta. It is also found that relatively high time-averaged wall shear stress (in the range of 4-8 kPa) may be associated with tear initiation and propagation. CFD modeling assists in medical planning by providing blood flow patterns, wall pressure and wall shear stress. This helps to understand various phenomena in the development of dissecting aneurysm.     

Influence of surface model extraction parameter on computational fluid dynamics modeling of cerebral aneurysms

Threshold image intensity for reconstructing patient-specific vascular models is generally determined subjectively. We assessed the effects of threshold image intensity differences on computational fluid dynamics (CFD) using a simple method of threshold determination. This study included 11 consecutive patients with internal carotid artery aneurysms collected retrospectively between April 2009 and March 2010. In 3-dimensional rotational angiography image data, we set a line probe across the coronal cross-section of the parent internal carotid artery, and calculated a profile curve of the image intensity along this line. We employed the threshold coefficient (C(thre)) value in this profile curve, in order to determine the threshold image intensity objectively. We assessed the effects of C(thre) value differences on vascular model configuration and the wall shear stress (WSS) distribution of the aneurysm. The threshold image intensity increased as the C(thre) value increased. The frequency of manual editing increased as the C(thre) value decreased, while disconnection of the posterior communicating artery occurred more frequently as the C(thre) value increased. The volume of the vascular model decreased and WSS increased according to the C(thre) value increase. The pattern of WSS distribution changed remarkably in one case. Threshold image intensity differences can produce profound effects on CFD. Our results suggest the uniform setting of C(thre) value is important for objective CFD.     

Towards the improved quantification of in vivo abnormal wall shear stresses in BAV-affected patients from 4D-flow imaging: Benchmarking and application to real data

Bicuspid aortic valve (BAV), i.e. the fusion of two aortic valve cusps, is the most frequent congenital cardiac malformation. Its progression is often characterized by accelerated leaflet calcification and aortic wall dilation. These processes are likely enhanced by altered biomechanical stimuli, including fluid-dynamic wall shear stresses (WSS) acting on both the aortic wall and the aortic valve. Several studies have proposed the exploitation of 4D-flow magnetic resonance imaging sequences to characterize abnormal in vivo WSS in BAV-affected patients, to support prognosis and timing of intervention. However, current methods fail to quantify WSS peak values.On this basis, we developed two new methods for the improved quantification of in vivo WSS acting on the aortic wall based on 4D-flow data.We tested both methods separately and in combination on synthetic datasets obtained by two computational fluid-dynamics (CFD) models of the aorta with healthy and bicuspid aortic valve. Tests highlighted the need for data spatial resolution at least comparable to current clinical guidelines, the low sensitivity of the methods to data noise, and their capability, when used jointly, to compute more realistic peak WSS values as compared to state-of-the-art methods.The integrated application of the two methods on the real 4D-flow data from a preliminary cohort of three healthy volunteers and three BAV-affected patients confirmed these indications. In particular, quantified WSS peak values were one order of magnitude higher than those reported in previous 4D-flow studies, and much closer to those computed by highly time- and space-resolved CFD simulations.     

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.     

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 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.