Effect of common carotid artery inlet length on normal carotid bifurcation hemodynamics

Controversy exists regarding the suitability of fully developed versus measured inlet velocity profiles for image-based computational fluid dynamics (CFD) studies of carotid bifurcation hemodynamics. Here, we attempt to resolve this by investigating the impact of the reconstructed common carotid artery (CCA) inlet length on computed metrics of "disturbed" flow. Twelve normal carotid bifurcation geometries were reconstructed from contrast-enhanced angiograms acquired as part of the Vascular Aging--The Link That Bridges Age to Atherosclerosis study (VALIDATE). The right carotid artery lumen geometry was reconstructed from its brachiocephalic origin to well above the bifurcation, and the CCA was truncated objectively at locations one, three, five, and seven diameters proximal to where it flares into the bifurcation. Relative to the simulations carried out using the full CCA, models truncated at one CCA diameter strongly overestimated the amount of disturbed flow. Substantial improvement was offered by using three CCA diameters, with only minor further improvement using five CCA diameters. With seven CCA diameters, the amounts of disturbed flow agreed unambiguously with those predicted by the corresponding full-length models. Based on these findings, we recommend that image-based CFD models of the carotid bifurcation should incorporate at least three diameters of CCA length if fully developed velocity profiles are to be imposed at the inlet. The need for imposing measured inlet velocity profiles would seem to be relevant only for those cases where the CCA is severely truncated.     

Use of factor analysis to characterize arterial geometry and predict hemodynamic risk: application to the human carotid bifurcation

The detailed geometry of atherosclerosis-prone vascular segments may influence their susceptibility by mediating local hemodynamics. An appreciation of the role of specific geometric variables is complicated by the considerable correlation among the many parameters that can be used to describe arterial shape and size. Factor analysis is a useful tool for identifying the essential features of such an inter-related data set, as well as for predicting hemodynamic risk in terms of these features and for interpreting the role of specific geometric variables. Here, factor analysis is applied to a set of 14 geometric variables obtained from magnetic resonance images of 50 human carotid bifurcations. Two factors alone were capable of predicting 12 hemodynamic metrics related to shear and near-wall residence time with adjusted squared Pearson's correlation coefficient as high as 0.54 and P-values less than 0.0001. One factor measures cross-sectional expansion at the bifurcation; the other measures the colinearity of the common and internal carotid artery axes at the bifurcation. The factors explain the apparent lack of an effect of branch angle on hemodynamic risk. The relative risk among the 50 bifurcations, based on time-average wall shear stress, could be predicted with a sensitivity and specificity as high as 0.84. The predictability of the hemodynamic metrics and relative risk is only modestly sensitive to assumptions about flow rates and flow partitions in the bifurcation.     

Measurement of hemodynamics in human carotid artery using ultrasound and computational fluid dynamics.

The objective of the study was to investigate the feasibility of using computational fluid dynamic modeling (CFD) with noninvasive ultrasound measurements to determine time-variant three-dimensional (3D) carotid arterial hemodynamics in humans in vivo. The effects of hyperoxia and hypoxic hypercapnia on carotid artery local hemodynamics were examined by use of this approach. Six normotensive volunteers followed a double-blind randomized crossover design. Blood pressure, heart rate, and carotid blood flow were measured while subjects breathed normal air, a mixture of 5% CO(2) and 15% O(2) (hypoxic hypercapnia), and 100% O(2) (hyperoxia). Carotid artery geometry was reconstructed on the basis of B-mode ultrasound images by using purpose-built image processing software. Time-variant 3D carotid hemodynamics were estimated by using finite volume-based CFD. Systemic blood pressure was not significantly affected by hyperoxia or hypoxic hypercapnia, but heart rate decreased significantly with hyperoxia. There was an increase in diastolic flow velocity in the external carotid artery after hypoxic hypercapnia, but otherwise carotid blood flow velocities did not change significantly. Compared with normal air, hyperoxic conditions were associated with a decrease in the width of the region of flow separation in the external carotid artery. During hyperoxia, there was also an increase in the minimum and a decrease in maximum shear stress in the bifurcation and hence a reduction in cyclic variation in shear stress. Hypoxic hypercapnia was associated with a reduced duration of flow separation in the external carotid artery and an increase in the minimum shear stress without affecting the cyclic variation in shear stress. This study demonstrates the feasibility of using noninvasive ultrasound techniques in conjunction with CFD to describe time-variant 3D hemodynamics in the human carotid arterial bifurcation in vivo.     

Comparison of steady and pulsatile flow in a double branching arterial model

Data are presented to compare fluid flow parameters for steady flow with those for time-varying flow in a simplified two branch model which simulates the region of the abdominal aorta near the celiac and superior mesenteric branches of the dog. Measurements in the model included laser doppler anemometry velocity profiles during steady flow, sinusoidal flow with a superimposed mean flow (referred to as simple oscillatory flow) and arterial pulsatile flow. Shear rate measurements were made by an electrochemical technique during steady flow. Flow visualization studies were done during steady and pulsatile flow. Fluid flow effects in the simplified model during steady flow showed many similarities to the results from previous steady flow studies in a canine aortic cast. Shear rates in the region of the proximal (first, or celiac) branch were independent of flow rates in the distal (second, or mesenteric) branch, but the shear pattern within the proximal branch changed significantly as flow in the proximal branch increased. Shear rates on the proximal flow divider (leading edge into the distal branch) depended primarily on the flow rate to the proximal branch, but not on flow to the distal branch. At certain daughter branch flow ratios (approximately 2:1, proximal to distal), flow separation was promoted at the outer wall of the second branch, but flow separation did not occur in the first branch. In contrast to the canine aortic case results, flow separation was never detected on the distal (mesenteric) flow divider of the simplified model. This observation reflects the subtle effects of geometry on flow since the mesenteric flow divider in the canine cast protrudes into the main flow whereas the distal flow divider in the simplified model does not. There were distinct differences in the flow phenomena between steady, simple oscillatory and arterial pulsatile flow. Peak shear rates during pulsatile flow were as much as 10--100 times greater than steady flow shear rates at comparable mean flow rates. Particularly noteworthy for the pulsatile flow with a Womersley parameter of sixteen were very blunt velocity profiles throughout systole, and the absence of flow separation or reversal in those regions of the model that exhibited flow separation during steady flow. The shape of the waveform influences the nature of the flow during time-varying flows. Future studies of fluid dynamics in model systems must consider the pulsatile nature of the flow if a true interpretation of arterial flow phenomena is to be made.     

Hemodynamics analysis of a stenosed carotid bifurcation and its plaque-mitigating design

Considering transient two-dimensional laminar flow in a diseased carotid artery segment with realistic inlet and outflow conditions, detailed velocity profiles, pressure fields, wall shear stress distributions and coupled, localized plaque formations have been simulated. The type of outflow boundary condition influences to a certain degree the extent of plaque build-up, which in turn reduces "disturbed flow" phenomena such as flow separations, recirculation zones, and wavy flow patterns in the artery branches during portions of the pulse. Based on computer experiments varying key geometric factors, a plaque-mitigating design of a carotid artery bifurcation has been proposed. Elimination of the carotid bulb, a smaller bifurcation angle, lower area ratios, and smooth wall curvatures generated a design with favorable hemodynamics parameters, leading to reduced plaque build-up by factors of 10 and 2 in the internal carotid and in the external carotid, respectively.     

Flow in t-bifurcations: Effect of the sharpness of the flow divider

The local geometry of a bifurcation has been hypothesized to be a potential geometrical risk factor for the development of atherosclerosis. While flow division and branch area ratios clearly affect the flow field, the importance of the flow divider shape is not as clear. A fast spectral element computational fluid mechanics (CFD) solver was used to simulate flow through 90 ° T-bifurcations with three different flow divider shapes. Other factors, such as flow partition, area ratio, and bifurcation angle, were kept constant. A Reynolds number range of 15 to 350 was studied to bracket experimental results in the literature. The variation in the sharpness of the corners had a dramatic effect on both the flow field and wall shear stress distribution in the side branch, but little effect on the flow in the main tube. The magnitude of reverse velocities and wall shear stress in the side branch increased linearly over a physiological range of Reynolds number and corner shape. This paper verifies the accuracy and usefulness of spectral element CFD in studying three-dimensional hemodynamics.     

Targeted particle tracking in computational models of human carotid bifurcations

A significant and largely unsolved problem of computational fluid dynamics (CFD) simulation of flow in anatomically relevant geometries is that very few calculated pathlines pass through regions of complex flow. This in turn limits the ability of CFD-based simulations of imaging techniques (such as MRI) to correctly predict in vivo performance. In this work, I present two methods designed to overcome this filling problem, firstly, by releasing additional particles from areas of the flow inlet that lead directly to the complex flow region ("preferential seeding") and, secondly, by tracking particles both "downstream" and "upstream" from seed points within the complex flow region itself. I use the human carotid bifurcation as an example of complex blood flow that is of great clinical interest. Both idealized and healthy volunteer geometries are investigated. With uniform seeding in the inlet plane (in the common carotid artery (CCA)) of an idealized bifurcation geometry, approximately half the particles passed through the internal carotid artery (ICA) and half through the external carotid artery. However, of those particles entering the ICA, only 16% passed directly through the carotid bulb region. Preferential seeding from selected regions of the CCA was able to increase this figure to 47%. In the second method, seeding of particles within the carotid bulb region itself led to a very high proportion (97%) of pathlines running from CCA to ICA. Seeding of particles in the bulb plane of three healthy volunteer carotid bifurcation geometries led to much better filling of the bulb regions than by particles seeded at the inlet alone. In all cases, visualization of the origin and behavior of recirculating particles led to useful insights into the complex flow patterns. Both seeding methods produced significant improvements in filling the carotid bulb region with particle tracks compared with uniform seeding at the inlet and led to an improved understanding of the complex flow patterns. The methods described may be combined and are generally applicable to CFD studies of fluid and gas flow and are, therefore, of relevance in hemodynamics, respiratory mechanics, and medical imaging science.     

Numerical flow studies in human carotid artery bifurcations: basic discussion of the geometric factor in atherogenesis

In this study fluid dynamic variables are analysed numerically in different human carotid artery bifurcation models in order to clarify the geometric factor in carotid bifurcation atherogenesis. The geometric variations describe healthy human carotid bifurcation anatomy and concern the shape of the carotid sinus and the angle between the branches. The flow conditions remain unchanged. The governing Navier-Stokes equations describing incompressible, pulsatile, three-dimensional viscous flow are approximated using a pressure correction finite element procedure which has been developed for time-consuming, three-dimensional, time-dependent viscous flow problems. The study concentrates on flow velocity, on detailed analysis of flow separation and flow recirculation, and on wall shear stress distribution. The results show that the extension and the location of the recirculation zone in the sinus as well as the duration of separated flow during the pulse cycle are affected by the geometrical variations. In view of the significance of the reversed flow zones and of the accompanied low shear regions in atherogenesis the geometry-dependent flow separation characteristics in the sinus is of substantial interest.     

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.     

Numerical analysis of flow through a severely stenotic carotid artery bifurcation

The results of computational simulations may supplement MR and other in vivo diagnostic techniques to provide an accurate picture of the hemodynamics in particular vessels, which may help demonstrate the risks of embolism or plaque rupture posed by particular plaque deposits. In this study, a model based on an endarterectomy specimen of the plaque in a carotid bifurcation was examined. The flow conditions include steady flow at Reynolds numbers of 300, 600, and 900 as well as unsteady pulsatile flow. Both dynamic pressure and wall shear stress are very high, with shear values up to 70 N/m2, proximal to the stenosis throat in the internal carotid artery, and both vary significantly through the flow cycle. The wall shear stress gradient is also strong along the throat. Vortex shedding is observed downstream of the most severe occlusion. Two turbulence models, the Chien and Goldberg varieties of k-epsilon, are tested and evaluated for their relevance in this geometry. The Chien model better captures phenomena such as vortex shedding. The flow distal to stenosis is likely transitional, so a model that captures both laminar and turbulent behavior is needed.     

Effect of variations in pressure and flow on the geometry of isolated canine carotid arteries

A study is described in which the effects of hemodynamics on arterial geometry are investigated in vitro. A novel perfusion apparatus is employed to deliver pulsatile flow through excised canine carotid arteries under carefully controlled conditions. Data of perfused vessel diameter and arterial wall thickness are derived from the radial displacement of the pulsating vessel as measured using a scanning laser micrometer whose accuracy is determined to be 0.0125 mm (0.0005 in). The results of 30 perfusion experiments suggest that the hemodynamic variables of transmural pressure, pulse pressure and flow rate influence vessel size and radial strain. The physiologic implications of these findings are discussed.     

The effect of proximal artery flow on the hemodynamics at the distal anastomosis of a vascular bypass graft: computational study

The formation of distal anastomotic intimal hyperplasia (IH), one common mode of bypass graft failure, has been shown to occur in the areas of disturbed flow particular to this site. The nature of theflow in the segment of artery proximal to the distal anastomosis varies from case to case depending on the clinical situation presented. A partial stenosis of a bypassed arterial segment may allow residual prograde flow through the proximal artery entering the distal anastomosis of the graft. A complete stenosis may allow for zero flow in the proximal artery segment or retrograde flow due to the presence of small collateral vessels upstream. Although a number of investigations on the hemodynamics at the distal anastomosis of an end-to-side bypass graft have been conducted, there has not been a uniform treatment of the proximal artery flow condition. As a result, direct comparison of results from study to study may not be appropriate. The purpose of this work was to perform a three-dimensional computational investigation to study the effect of the proximal artery flow condition (i.e., prograde, zero, and retrograde flow) on the hemodynamics at the distal end-to-side anastomosis. We used the finite volume method to solve the full Navier-Stokes equations for steady flow through an idealized geometry of the distal anastomosis. We calculated the flow field and local wall shear stress (WSS) and WSS gradient (WSSG) everywhere in the domain. We also calculated the severity parameter (SP), a quantification of hemodynamic variation, at the anastomosis. Our model showed a marked difference in both the magnitude and spatial distribution of WSS and WSSG. For example, the maximum WSS magnitude on the floor of the artery proximal to the anastomosis for the prograde and zero flow cases is 1.8 and 3.9 dynes/cm2, respectively, while it is increased to 10.3 dynes/cm2 in the retrograde flow case. Similarly, the maximum value of WSSG magnitude on thefloor of the artery proximal to the anastomosis for the prograde flow case is 4.9 dynes/cm3, while it is increased to 13.6 and 24.2 dynes/cm3, respectively, in the zero and retrograde flow cases. The value of SP is highest for the retrograde flow case (13.7 dynes/cm3) and 8.1 and 12.1 percent lower than this for the prograde (12.6 dynes/cm3) and zero (12.0 dynes/cm3) flow cases, respectively. Our model results suggest that the flow condition in the proximal artery is an important determinant of the hemodynamics at the distal anastomosis of end-to-side vascular bypass grafts. Because hemodynamic forces affect the response of vascular endothelial cells, the flow situation in the proximal artery may affect IH formation and, therefore, long-term graft patency. Since surgeons have some control over the flow condition in the proximal artery, results from this study could help determine which flow condition is clinically optimal.     

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.     

Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements

The presence of atherosclerotic plaques has been shown to be closely related to the vessel geometry. Studies on postmortem human arteries and on the experimental animal show positive correlation between the presence of plaque thickness and low shear stress, departure of unidirectional flow and regions of flow separation and recirculation. Numerical simulations of arterial blood flow and direct blood flow velocity measurements by magnetic resonance imaging (MRI) are two approaches for the assessment of arterial blood flow patterns. In order to verify that both approaches give equivalent results magnetic resonance velocity data measured in a compliant anatomical carotid bifurcation model were compared to the results of numerical simulations performed for a corresponding computational vessel model. Cross sectional axial velocity profiles were calculated and measured for the midsinus and endsinus internal carotid artery. At both locations a skewed velocity profile with slow velocities at the outer vessel wall, medium velocities at the side walls and high velocities at the flow divider (inner) wall were observed. Qualitative comparison of the axial velocity patterns revealed no significant differences between simulations and in vitro measurements. Even quantitative differences such as for axial peak flow velocities were less than 10%. Secondary flow patterns revealed some minor differences concerning the form of the vortices but maximum circumferential velocities were in the same range for both methods.     

Local and global geometric influence on steady flow in distal anastomoses of peripheral bypass grafts

We consider the effect of geometrical configuration on the steady flow field of representative geometries from an in vivo anatomical data set of end-to-side distal anastomoses constructed as part of a peripheral bypass graft. Using a geometrical classification technique, we select the anastomoses of three representative patients according to the angle between the graft and proximal host vessels (GPA) and the planarity of the anastomotic configuration. The geometries considered include two surgically tunneled grafts with shallow GPAs which are relatively planar but have different lumen characteristics, one case exhibiting a local restriction at the perianastomotic graft and proximal host whilst the other case has a relatively uniform cross section. The third case is nonplanar and characterized by a wide GPA resulting from the graft being constructed superficially from an in situ vein. In all three models the same peripheral resistance was imposed at the computational outflows of the distal and proximal host vessels and this condition, combined with the effect of the anastomotic geometry, has been observed to reasonably reproduce the in vivo flow split. By analyzing the flow fields we demonstrate how the local and global geometric characteristics influences the distribution of wall shear stress and the steady transport of fluid particles. Specifically, in vessels that have a global geometric characteristic we observe that the wall shear stress depends on large scale geometrical factors, e.g., the curvature and planarity of blood vessels. In contrast, the wall shear stress distribution and local mixing is significantly influenced by morphology and location of restrictions, particular when there is a shallow GPA. A combination of local and global effects are also possible as demonstrated in our third study of an anastomosis with a larger GPA. These relatively simple observations highlight the need to distinguish between local and global geometric influences for a given reconstruction. We further present the geometrical evolution of the anastomoses over a series of follow-up studies and observe how the lumen progresses towards the faster bulk flow of the velocity in the original geometry. This mechanism is consistent with the luminal changes in recirculation regions that experience low wall shear stress. In the shallow GPA anastomoses the proximal part of the native host vessel occludes or stenoses earlier than in the case with wide GPA. A potential contribution to this behavior is suggested by the stronger mixing that characterizes anastomoses with large GPA.     

Carotid bifurcation position and branching angle in patients with atherosclerotic carotid disease

Carotid artery bifurcation (CB) is the preferred site for development of atherosclerosis (AS) in extracranial cerebral arteries; internal carotid artery stenosis is the most common cause of ischemic stroke. The frequent atherosclerotic disease of CB may best be explained by the hemodynamic influence of complex blood flow that results from the unique geometry of the bifurcation. Few papers analyze all possible geometric structural characteristics of this bifurcation. While performing many carotid endarterectomies, we noticed that a certain correlation between CB height in the neck and its angle existed, that a larger angle is accompanied with increased frequency of elongation and kinking and that CB shape influences distribution of atherosclerosis. The purpose of this paper is to quantify and evaluate these clinical observations. Radiogrametric analysis of 154 bi-plane orthogonal aortic arch arteriograms of patients with symptomatic atherosclerotic carotid artery disease was performed and a total of 289 CBs were analyzed. The CB height in relation to cervical spine segments was measured and real angles of each bifurcation were calculated. A positive linear correlation between CB height and angle exists: the CB angle increases /decreases 3.34 degrees for each third of the cervical vertebral body height or intervertebral space height. The CB is positioned a little higher on the left side. The proximal border of the atherosclerotic process is found at the level of intersection of the axes of the common carotid artery branches in 92.6% of examined CBs. In lower CBs (with smaller angles) the proximal border was located in the last segment of the common carotid artery, while in high bifurcations (wider angles) the proximal border of the AS process is more distally in the blood flow, in the beginning of the internal carotid artery, and the process was more extensive. High CBs are more suitable for eversion endarterectomy while normal and low CBs are more suitable for open (classic) carotid endarterectomy. The influence of the geometric risk factor demands further investigation.     

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.     

A numerical and experimental analysis of the flow field in a two-dimensional model of the human carotid artery bifurcation

Axial velocities were measured in an enlarged, two-dimensional, rigid model of the carotid artery bifurcation by means of a laser-Doppler anemometer, under both steady and unsteady flow conditions. Also a numerical model was developed, based on the finite element approximation of the Navier-Stokes and continuity equations. From this study it appeared that the numerically predicted velocities agree well with the experimentally obtained values. Besides, the bifurcation hardly influenced the upstream flow in the main branch (common carotid artery), high velocity gradients were observed at the divider walls of the daughter branches (internal and external carotid arteries) and large zones with reversed flow were present near the nondivider walls of these branches. For steady flow the maximal diameter of this zone at the entrance of the internal carotid artery (carotid sinus) was about 25% of the local diameter of this branch. For unsteady flow this zone was absent during the initial phase of flow acceleration and maximal at the end of flow deceleration with a maximal diameter of about 50% of the local diameter of the carotid sinus.     

An experimental analysis of the flow field in a three-dimensional model of the human carotid artery bifurcation

Steady flow measurements were carried out in a rigid three-dimensional model of the human carotid artery bifurcation at a Reynolds number of 640 and a flow division ratio of 50/50. Both axial and secondary velocities were measured with a laser-Doppler anemometer. In the bulb opposite to the flow divider a zone with negative axial velocities was found with a maximal diameter of about 60% of the local diameter of the branch and a cross-sectional extent of about 25% of the local cross-sectional area. In the bulb the maximum axial velocity shifted towards the divider wall and at the end of the bulb an axial velocity plateau arose near the non-divider wall. Halfway through the bulb, secondary flow showed a vortex through which fluid flowed towards the divider wall near the bifurcation plane and back towards the non-divider wall near the upper walls.