Phasic flow patterns at a hemodialysis venous anastomosis

A phase-by-phase analysis of local flow patterns at the venous anastomosis of an arteriovenous hemodialysis angioaccess loop graft (AVLG) was made. The study was carried out in an elastic, transparent Silastic in vitro flow model, which duplicates the detail geometry of the AVLG obtained from an animal model (30+ kg dogs with 12 weeks bilateral femoral AVLG implantation). The flow model was installed in a mock pulsatile flow loop system designed to simulate physiological conditions. Flow visualization was made in laser-illuminated flow fields using a high-speed cine camera. Analysis of the high-speed cine indicates there is a distinct separation region downstream of the anastomotic toe in the median plane and a stagnation region that oscillates along the opposite wall. During inward motion of the vessel wall, accumulation of particles in the separation region and the nearby stagnation region is observed. A large swirl appears in the distal vein during end-systolic period. A double-helical flow pattern occurs further down in the distal vein. Retrograde flow in the distal vein occurs in an "oscillating" manner following each cardiac cycle.     

Haemodynamics of angioaccess venous anastomoses

A large percentage of arteriovenous haemodialysis angioaccess loop grafts (AVLG) fail within the first year after surgery, the occlusive lesions being found predominantly at the venous anastomosis site. This paper presents a detailed flow dynamic study of the AVLG system using three elastic, transparent bench-top flow models, which were based on the geometry of silicone rubber casts obtained at different times from a chronic animal model. Each model thus represented a different stage of the lesion development. Flow visualization and laser Doppler anemometer surveys of the flow field confirmed that the hydrodynamic factors favour lesion development near the stagnation point opposite the anastomotic toe, where the momentum of the impinging jet stream, combined with the oscillating wall shear stress generated in the vicinity of the stagnation point, acts in both directions. The accumulation of tracer particles in the region of flow separation is believed to be a combined contribution from the hydraulic forces and the inward motion of the vessel wall. As these hydrodynamic factors are enhanced upon further development of the occlusive lesion, a vicious cycle may be formed.     

An in vitro comparative study of St. Jude Medical and Edwards-Duromedics bileaflet valves using laser anemometry

An in vitro comparative study of St. Jude (SJ) and Edwards-Duromedics (DM) Bileaflet valves was performed under steady and physiological pulsatile flow conditions in an axisymmetric chamber using Laser Doppler Anemometry (LDA). LDA measurements were conducted in two different orientations; in the first orientation, the LDA traverse was perpendicular and, in the second orientation, parallel to the tilt axis of the leaflets. The axial velocities were measured in both orientations at two different locations distal to the valves. The velocity profiles at peak systole show the presence of stronger vortex in the sinus region for flow past SJ valve in the first orientation compared to the DM valve. Velocity profile distal to the SJ valve in second orientation was relatively flat where as for the DM valve, a jet-like flow was present. The differences found in the velocity profiles between the two valves can be attributed to the differences in geometry with thicker leaflets, smaller angle of leaflets opening and the presence of the leaflet curvature for the DM valve. The results obtained in this study do not show any fluid dynamic advantages due to the curved leaflet geometry of the DM valve.     

The effect of angle and flow rate upon hemodynamics in distal vascular graft anastomoses: an in vitro model study.

A steady flow, in vitro model of distal arterial bypass graft junctions was used to examine the effects of junction angle and flow rate on the local velocity field. Three test sections were fabricated from Plexiglas tubing having anastomotic junction angles of either 30, 45, or 60 deg. Flow visualization revealed velocity profiles skewed toward the outer wall with a flow split around a clear stagnation point along the outer wall. Laser Doppler anemometry [LDA] measurements confirmed a distinct stagnation point at the outer wall and both reverse and forward shear were detected immediately upstream and downstream, respectively, of this site. Axial velocities and shear rates along the outer wall were higher than along the inner wall and occurred in the junction angle order: 45, 60, and 30 deg. This study clearly identified changes in wall shear which varied with the anastomotic angle and flow rate.     

Analysis of blood flow in an out-of-plane CABG model

Coronary artery bypass graft (CABG) is a routine surgical treatment for ischemic and infarcted myocardium. A large number of CABG fail postoperatively because of intimal hyperplasia within months or years. The cause of this failure is thought to be partly related to the flow patterns and shear stresses acting on the endothelial cells. An accurate representation of the flow field and associated wall shear stress (WSS) requires a detailed three-dimensional (3D) model of the CABG. The purpose of this study is to present a detailed analysis of blood flow in a 3D aorto/left CABG, bypassing the occluded left anterior descending coronary (LAD) artery. The analysis takes into account the influence of the out-of-plane geometry of the graft. The finite volume technique was employed to model the 3D blood flow pattern to determine the velocity and WSS distributions. This study presents the flow field distributions of the velocity and WSS at four instances of the cardiac cycle, two in systole and two in diastole. Our results reveal that the CABG geometry has a significant effect on the velocity distribution. The axial velocity profiles at different instances of the cardiac cycle exhibit strong skewing; significant secondary flow and vortex structures are seen in the in-plane velocity patterns. The maximum WSS on the bed of the occluded LAD artery opposite to the graft junction is 14 Pa in middiastole, whereas there is a significantly lower and more uniform distribution of WSS on the bed of the anastomosis. The present results indicate that nonplanarity of the blood vessel along with the inflow conditions has a substantial effect on the fluid mechanics of CABG that contribute to the patency of graft.     

Experimental study of physiological pulsatile flow past valve prostheses in a model of human aorta--I. Caged ball valves

Pulsatile flow development past a caged ball valve in a model human aorta was studied using laser Doppler anemometry. Velocity profiles measured in the ascending aorta and in the mid-arch region were strongly influenced by the geometry of the valve at the root of the aorta. Velocity profiles distal to the valve were asymmetric with jet-like flow in the peripheral region having larger velocity magnitudes towards the left lateral wall. In early diastole, a streamwise vortex motion was observed throughout the model aorta with fluid moving towards the downstream direction along the left lateral wall and reversed flow along the right lateral wall. With the caged ball valve at the root of the aorta, no reversed flow was observed along the inner wall of curvature in the mid-arch region.     

Quantitative comparison of CFD predicted and MRI measured velocity fields in a carotid bifurcation phantom

Steady flow of a blood mimicking fluid in a physiologically realistic model of the human carotid bifurcation was studied using both magnetic resonance imaging (MRI) and computational fluid dynamics (CFD) modelling techniques. Quantitative comparisons of the 3D velocity field in the bifurcation phantom were made between phase contrast MRI measurements and CFD predictions. The geometry for the CFD model was reconstructed from T(1) weighted MR imaging of the test phantom. It was found that the predicted velocity fields were in fair agreement with MR measured velocities. In both the internal and external carotid arteries, the agreement between CFD predictions and MRI measurements was better along the inner-outer wall axis with a correlation factor C>0.897 (average 0.939) where the velocity profiles were skewed, than along the anterior-posterior axis (average correlation factor 0.876) where the velocity profiles were in M-shape.     

An introduction to biofluid mechanics--basic models and applications

Cardiovascular disease is the primary cause of morbidity and mortality in the western world. Complex hemodynamics play a critical role in the development of atherosclerosis and the processes of aging, as well as many other disease processes. Biofluid mechanics play a major role in the cardiovascular system and it is important to understand the forces and movement of blood cells and whole blood as well as the interaction between blood cells and the vessel wall. Fundamental fluid mechanical, which are important for the understanding of the blood flow in the cardiovascular circulatory system of the human body aspects are presented. Measurement techniques for model studies such as LDA, ultrasound, and MRI studies will be discussed. Viscosity and flow behavior changes specifically the creation of vortices and flow disturbances can be used to show how medication can influence flow behavior. Experiments have shown that hemodynamics may have a strong influence on the creation of aneurysms and varicose veins. Other factors such as vessel wall structure are also important. In preliminary studies, it has been demonstrated that geometry and elasticity of vessel walls help determine flow behavior. High velocity fluctuations indicate flow disturbances that should be avoided. Health care practitioners must understand fluid dynamic factors such as flow rate ratio, pressure and velocity gradients, and flow behavior, velocity distribution, shear stress on the wall and on blood cells. These mechanical factors are largely responsible for the deposit of blood cells and lipids, a leading cause of atherosclerosis. The interaction between blood cells and of the cells with the vessel, leads to the formation of plaques and agglomerations. These deposits are found predominantly at arterial bends and bifurcations where blood flow is disturbed, where a secondary flow is created, and where flow separation regions are found. Experiments on hemodynamic effects in elastic silicon rubber models of the cardiovascular system with flow wire, stents, or patches for vessel surgery will be discussed. These studies can be important in improving diagnostics and therapeutic applications.     

Axial flow velocity patterns in a pulmonary artery model with varying degrees of valvular pulmonic stenosis: pulsatile in vitro studies

With the advent of noninvasive clinical techniques which can measure blood flow velocities (Doppler ultrasound), it is suggested that a fundamental knowledge of the axial flow velocity patterns in the pulmonary artery, and the changes caused by stenosis, may be used to support accurate diagnosis of valvular pulmonic stenosis. The present study was designed to characterize the axial flow velocity patterns in an in vitro model of a human adult pulmonary artery with varying degrees of valvular pulmonic stenosis. A two-dimensional laser Doppler anemometer (LDA) system was used to map the flow fields in the main (MPA), left (LPA), and right (RPA) branches of the pulmonary artery model. The study was conducted in the Georgia Tech. right heart pulse duplicator system. It was observed that the axial flow velocity patterns in the MPA and the LPA change dramatically with increasing degree of valvular stenosis. This indicates that the axial flow velocity patterns in these two branches are strongly influenced by the degree of valvular stenosis. The axial flow velocity patterns in the RPA, however, do not change much with varying degrees of valvular stenosis, indicating that the axial flow fields in the RPA are mainly influenced by the geometry of the bifurcation. It may be concluded therefore, that the changes in the axial flow velocity patterns in the MPA and LPA (rather than in the RPA) could be sensitive and reliable indicators of the severity of the defect.     

What measurements are necessary for a comprehensive evaluation of the peripheral arterial circulation?

Several methods are available to detect atherosclerotic lesions with a severe degree of stenosis (>70%), but the diagnosis of atherosclerotic lesions with no stenosis or with a minor degree of stenosis (<20%), is problematic. Hemodynamics associated with stenotic lesions are well described by the relationship of blood pressure and blood flow velocity, both as a function of time and localization (along the length and cross-section of the vessel). The use of this relationship in the clinic is difficult because no precise information is available about the geometry and branching of arteries, blood viscosity, and the velocity distribution over the cross-sectional area of the blood vessel. Besides, the invasiveness of the technique to measure arterial pressure as a function of time and localization does not allow routine application in patients. Because of these limitations, alternative methods have been developed. The degree and extensiveness of atherosclerotic disease can, for instance, be estimated from the changes in maximum blood flow velocity and in velocity profile, i.e., velocity distribution along the cross-section of the vessel. Moreover, the delay between simultaneously recorded arterial blood flow velocity tracings (pulse-wave velocity determination) is used to assess the elastic properties of the vessel. Changes in velocity profile occur at relatively slight degrees of arterial stenosis (around 20%), so that determination of these profiles along diseased arteries may contribute to the early diagnosis of atherosclerotic lesions. In man, transcutaneous information about the maximum and mean blood flow velocities over the cross-sectional area of the artery as an instantaneous function of time as well as the flow pattern can be obtained online with continuous wave Doppler flowmeters, at least when audio spectrum analysis is used as a processing technique. Velocity profiles can be determined with multichannel pulsed Doppler systems if the resolution of the system is adequate and a sufficient number of sample volumes can be obtained, limiting the interpolation between these samples. The on-line recording of velocity profiles can be facilitated by combining the pulsed Doppler device with either a velocity imaging system or a B-mode scan. In systems with a high resolution (sample distance 0.5 mm), one should be able to detect local disturbances in the velocity profile at the site of the lesion (due to local increases in shear stress) and proximal to the lesion (due to reflections), so that lesions with a minor degree of stenosis can be detected. In resistive systems (e.g., internal carotid arteries) in which the relationship between pressure and velocity changes during the cardiac cycle is relatively simple, the elasticity of the arterial wall can be determined by relating the relative diameter changes of the vessel, determined on-line with multichannel pulsed Doppler systems, to the instantaneous velocity pulse. Although the detection of atherosclerotic lesions at an early stage of the disease with sophisticated Doppler devices looks promising, further clinical evaluation is required.     

Blood flow and vessel mechanics in a physiologically realistic model of a human carotid arterial bifurcation

The pulsatile flow in an anatomically realistic compliant human carotid bifurcation was simulated numerically. Pressure and mass flow waveforms in the carotid arteries were obtained from an individual subject using non-invasive techniques. The geometry of the computational model was reconstructed from magnetic resonance angiograms. Maps of time-average wall shear stress, contours of velocity in the flow field as well as wall movement and tensile stress on the arterial wall are all presented. Inconsistent with previous findings from idealised geometry models, flow in the carotid sinus is dominated by a strong helical flow accompanied by a single secondary vortex motion. This type of flow is induced primarily by the asymmetry and curvature of the in vivo geometry. Flow simulations have been carried out under the rigid wall assumption and for the compliant wall, respectively. Comparison of the results demonstrates the quantitative influence of the vessel wall motion. Generally there is a reduction in the magnitude of wall shear stress, with its degree depending on location and phase of the cardiac cycle. The region of slow or reversed flow was greater, in both spatial and temporal terms in the compliant model, but the global characteristics of the flow and stress patterns remain unchanged. The analysis of mechanical stresses on the vessel surface shows a complicated stress field. Stress concentration occurs at both the anterior and posterior aspects of the proximal internal bulb. These are also regions of low wall shear stress. The comparison of computed and measured wall movement generally shows good agreement.     

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.     

A model investigation of the velocity and pressure spectra in vascular murmurs

A physical model consisting of an axisymmetrical jet in a rigid plexiglass pipe was used to study the flow and pressure fluctuations downstream from an aortic stenosis. The fluctuating velocity components, u and v, at several locations in the steady liquid jet were measured using a laser Doppler anemometer system. Simultaneous wall pressure fluctuations were monitored by an array of nine miniature pressure transducers wall mounted in the axial direction. This paper presents the detailed measurements of mean velocity profiles, turbulent intensity distributions and RMS pressure fluctuations. The energy spectra obtained for the pressure fluctuations and the u and v velocity components are compared. Contrary to earlier works, we found that the differences between peak frequencies of the pressure spectra and the characteristic frequencies of the velocity spectra vary with positions downstream from the nozzle. These differences are discussed in light of pseudosound generation by the eddy structures in the stenotic flow field.     

Non-Newtonian effects of blood flow on hemodynamics in distal vascular graft anastomoses

Non-Newtonian fluid flow in a stenosed coronary bypass is investigated numerically using the Carreau-Yasuda model for the shear thinning behavior of the blood. End-to-side coronary bypass anastomosis is considered in a simplified model geometry where the host coronary artery has a 75% severity stenosis. Different locations of the bypass graft to the stenosis and different flow rates in the graft and in the host artery are studied. Particular attention is given to the non-Newtonian effect of the blood on the primary and secondary flow patterns in the host coronary artery and the wall shear stress (WSS) distribution there. Interaction between the jet flow from the stenosed artery and the flow from the graft is simulated by solving the three-dimensional Navier-Stokes equation coupled with the non-Newtonian constitutive model. Results for the non-Newtonian flow, the Newtonian flow and the rescaled Newtonian flow are presented. Significant differences in axial velocity profiles, secondary flow streamlines and WSS between the non-Newtonian and Newtonian fluid flows are revealed. However, reasonable agreement between the non-Newtonian and the rescaled Newtonian flows is found. Results from this study support the view that the residual flow in a partially occluded coronary artery interacts with flow in the bypass graft and may have significant hemodynamic effects in the host vessel downstream of the graft. Non-Newtonian property of the blood alters the flow pattern and WSS distribution and is an important factor to be considered in simulating hemodynamic effects of blood flow in arterial bypass grafts.     

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

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study 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 cell hydrodynamics and their effects on E. coli biofilm formation under different nutrient conditions and turbulent flow

Biofilm formation is a major factor in the growth and spread of both desirable and undesirable bacteria as well as in fouling and corrosion. In order to simulate biofilm formation in industrial settings a flow cell system coupled to a recirculating tank was used to study the effect of a high (550 mg glucose l?1) and a low (150 mg glucose l?1) nutrient concentration on the relative growth of planktonic and attached biofilm cells of Escherichia coli JM109(DE3). Biofilms were obtained under turbulent flow (a Reynolds number of 6000) and the hydrodynamic conditions of the flow cell were simulated by using computational fluid dynamics. Under these conditions, the flow cell was subjected to wall shear stresses of 0.6 Pa and an average flow velocity of 0.4 m s?1 was reached. The system was validated by studying flow development on the flow cell and the applicability of chemostat model assumptions. Full development of the flow was assessed by analysis of velocity profiles and by monitoring the maximum and average wall shear stresses. The validity of the chemostat model assumptions was performed through residence time analysis and identification of biofilm forming areas. These latter results were obtained through wall shear stress analysis of the system and also by assessment of the free energy of interaction between E. coli and the surfaces. The results show that when the system was fed with a high nutrient concentration, planktonic cell growth was favored. Additionally, the results confirm that biofilms adapt their architecture in order to cope with the hydrodynamic conditions and nutrient availability. These results suggest that until a certain thickness was reached nutrient availability dictated biofilm architecture but when that critical thickness was exceeded mechanical resistance to shear stress (ie biofilm cohesion) became more important.     

Steady flow development past valve prostheses in a model human aorta—I. Centrally occluding valves

In this paper, laser-Doppler anemometry measurement of steady flow development in a model human aorta has been reported. Studies were made with uniform entry flow at the root of the aorta and our measurements showed the establishment of a pair of Dean vortices in the mid-arch region. Subsequently, the nature of flow development past centrally occluding caged ball valves in the model aorta was investigated. Our studies showed that in the ascending aorta, an asymmetric velocity profile is obtained with larger velocity gradients towards the inner wall of tertiary curvature (anatomically the left lateral wall) with centrally occluding valves. The peripheral flow past these valves prevented the development of Dean vortices in the mid-arch region. The caged ball valves at the root of the aorta had no discernible effect on the velocity profiles in the brachio-cephalic artery.     

A two-dimensional numerical analysis of unsteady flow in the carotid artery bifurcation. A comparison with three-dimensional in-vitro measurements and the influence of minor stenoses

In the present study a two-dimensional finite element model for incompressible Newtonian flow is applicated to the modelling of carotid artery flow. In earlier studies, the numerical model was validated experimentally for several flow configurations. In general the pulsatile flow is characterized by reversed flow regions at the non-divider side walls of both the internal and external carotid arteries. The unsteadiness of the flow is associated with rather complex spatial and temporal velocity distributions and leads to temporal variations of the location and length of the reversed flow regions. As a consequence, pronounced spatial and temporal variations in the wall shear stresses are found. At the non-divider side walls, wall shear stresses are relatively low and exhibits an oscillatory behaviour in space and time. At the divider side walls, wall shear stresses are relatively high and approximately follow the flow rate distribution in time. The aim of this study is not only to present two-dimensional calculations but also to compare the calculated two-dimensional velocity profiles with those from three-dimensional experiments. It is observed that in the common carotid artery and in the proximal parts of the internal and external carotid arteries, the two-dimensional numerical model provides valuable information with respect to the three-dimensional configuration. In the more distal parts of especially the internal carotid artery, deviations are found between the two-dimensional numerical and three-dimensional experimental model. These deviations can mainly be attributed to the neglect of the secondary velocity distribution in the two-dimensional model. In the two-dimensional numerical model the influence of a minor stenosis in the internal carotid artery is hardly distinguishable from a minor geometrical variation without stenosis. Full three-dimensional analyses of the influence of minor stenoses are needed to prove numerically whether in-vivo measurements of the axial velocity distribution are useful in the detection of minor stenoses.     

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.     

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.