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.     

Three-dimensional numerical analysis of pulsatile flow and wall shear stress in the carotid artery bifurcation.

To analyse the pulsatile flow field and the mechanical stresses in a three-dimensional carotid artery bifurcation model, computer simulation is applied. The approximation of the Navier-Stokes equations uses a pressure correction finite element method. Numerical results are presented for axial and secondary flow velocity and wall shear stresses with special emphasis on the fluid dynamics in the carotid sinus. This region is of major interest because it is affected preferentially by lesions. Detailed local flow studies as carried out here should lead to a further insight into the mechanisms of atherogenesis. The flow conditions used in the study were chosen according to Ku et al. (Arteriosclerosis 5, 293-302, 1985). The results of this numerical analysis agree in the essential features with their experimental results.     

A numerical calculation of flow in a curved tube model of the left main coronary artery

The flow pattern in the left main coronary artery has been calculated using an idealized geometry and by numerically solving the full Navier-Stokes equations for a Newtonian fluid. Two different forms for the entrance velocity profile were used, one a time-varying, flat profile and the other a time-varying, less flat velocity profile. The results obtained demonstrate the presence of secondary motions for conditions simulating flow in the left main coronary artery, with maximum secondary flow velocities being on the order of three to four percent of the maximum axial velocity. This secondary flow phenomenon has an important influence on the wall shear stress distribution, in spite of the fact that there is virtually no alteration in the axial velocity profile. The maximum ratio of the outer wall shear stress to that on the inner wall is 1.4 at a Reynolds number of Re = 270, and it increases with increasing Reynolds number, reaching a value of 1.7 at Re = 810. Although there are significant differences in the results in the immediate vicinity of the inlet for the two different forms of the entrance velocity profile used, this difference does not persist far into the tube. Independent of the choice of the entrance velocity profile, it appears that there will be significant secondary flow effects on the wall shear stress.     

Three-dimensional steady flow through a bifurcation

Steady flow of an incompressible, Newtonian fluid through a symmetric bifurcated rigid channel was numerically analyzed by solving the three-dimensional Navier-Stokes equations. The upstream Reynolds number ranged from 100 to 1500. The bifurcation was symmetrical with a branch angle of 60 deg and the area ratio of the daughter to the mother vessel was 2.0. The numerical procedure utilized a coordinate transformation and a control volume approach to discretize the equations to finite difference form and incorporated the SIMPLE algorithm in performing the calculation. The predicted velocity pattern was in qualitative agreement with experimental measurements available in the literature. The results also showed the effect of secondary flow which can not be predicted using previous two-dimensional simulations. A region of reversed flow was observed near the outer wall of the branch except for the case of the lowest Reynolds number. Particle trajectory was examined and it was found that no fluid particles remained within the recirculation zone. The shear stress was calculated on both the inner and the outer wall of the branch. The largest wall shear stress, located in the vicinity of the apex of the branch, was of the same order of magnitude as the level that can cause damage to the vessel wall as reported in a recent study.     

Oscillatory flow and gas transport through a symmetrical bifurcation

Axial gas transport due to the interaction between radial mixing and radially nonuniform axial velocities is responsible for gas transport in thick airways during High-frequency oscillatory ventilation (HFO). Because the airways can be characterized by a bifurcating tube network, the secondary flow in the curved portion of a bifurcating tube contributes to cross-stream mixing. In this study the oscillatory flow and concentration fields through a single symmetrical airway bifurcating tube model were numerically analyzed by solving three-dimensional Navier-Stokes and mass concentration equations with the SIMPLER algorithm. The simulation conditions were for a Womersley number, alpha = 9.1 and Reynolds numbers in the parent tube between 200 and 1000, corresponding to Dn2/alpha 4 in the curved portion between 2 and 80, where Dn is Dean number. For comparison with the results from the bifurcating tube, we calculated the velocity and concentration fields for fully developed oscillatory flow through a curved tube with a curvature rate of 1/10, which is identical to the curved portion of the bifurcating tube. For Dn2/alpha 4 < or = 10 in the curved portion of the bifurcating tube, the flow divider and area changes dominate the axial gas transport, because the effective diffusivity is greater than in either a straight or curved tube, in spite of low secondary velocities. However, for Dn2/alpha 4 > or = 20, the gas transport characteristics in a bifurcation are similar to a curved tube because of the significant effect of secondary flow.     

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

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.     

Comparison of axisymmetric and three-dimensional models for gas uptake in a single bifurcation during steady expiration

Reactive gas uptake is predicted and compared in a single bifurcation at steady expiratory flow in terms of Sherwood number using an axisymmetric single-path model (ASPM) and a three-dimensional computational fluid dynamics model (CFDM). ASPM is validated in a two-generation geometry by comparing the average gas-phase mass transfer coefficients with the experimental values. ASPM predicted mass transfer coefficients within 20% of the experimental values. The flow and concentration variables in the ASPM were solved using Galerkin finite element method and in the CFDM using commercial finite element software FIDAP. The simulations were performed for reactive gas flowing at Reynolds numbers ranging from 60 to 350 in both symmetric bifurcation for three bifurcation angles, 30 deg, 70 deg, and 90 deg, and in an asymmetric bifurcation. The numerical models compared with each other qualitatively but quantitatively they were within 0.4-8% due to nonfully developed flow in the parent branch predicted by the CFDM. The radially averaged concentration variation along the axial location matched qualitatively between the CFDM and ASPM but quantitatively they were within 32% due to differences in the flow field. ASPM predictions compared well with the CFDM predictions for an asymmetric bifurcation. These results validate the simplified ASPM and the complex CFDM. ASPM predicts higher Sherwood number with a flat velocity inlet profile compared to a parabolic inlet velocity profile. Sherwood number increases with the inlet average velocity, wall mass transfer coefficient, and bifurcation angle since the boundary layer grows slower in the parent and daughter branches.     

Measurement and calculations of laminar flow in a ninety degree bifurcation

Measurements and numericaL calculations of laminar flow in a plane 90 degrees bifurcation are presented. The corresponding two-dimensional steady flow Navier-Stokes equations solved by a finite-difference procedure employing pressure and velocity as dependent variables. The influence of Reynolds number and mass flow ratio on the velocity field, streamlines, local shear stress and pressure drop are quantified and shown to be substantial. The circulation patterns and shear stresses are examined in view of available data regarding the formation of atherotic plaques in the human circulatory system. The calculated velocity profiles are compared with measurements obtained with laser Doppler anemometry and the agreement is shown to be satisfactory. Calculations outside the range of measurements which are of value to biomechanics are also presented.     

Computational study of the effect of geometric and flow parameters on the steady flow field at the rabbit aorto-celiac bifurcation

Arterial hemodynamic forces may play a role in the localization of early atherosclerotic lesions. We have been developing numerical techniques based on overset or "Chimera" type formulations to solve the Navier-Stokes equations in complex geometries simulating arterial bifurcations. This paper presents three-dimensional steady flow computations in a model of the rabbit aorto-celiac bifurcation. The computational methods were validated by comparing the numerical results to previously-obtained flow visualization data. Once validated, the numerical algorithms were used to investigate the sensitivity of the computed flow field and resulting wall shear stress distribution to various geometric and hemodynamic parameters. The results demonstrated that a decrease in the extent of aortic taper downstream of the celiac artery induced looping fluid motion along the lateral walls of the aorta and shifted the peak wall shear stress from downstream of the celiac artery to upstream. Increasing the flow Reynolds number led to a sharp increase in spatial gradients of wall shear stress. The flow field was highly sensitive to the flow division ratio, i.e., the fraction of total flow rate that enters the celiac artery, with larger values of this ratio leading to the occurrence of flow separation along the dorsal wall of the aorta. Finally, skewness of the inlet velocity profile had a profound impact on the wall shear stress distribution near the celiac artery. While not physiological due to the assumption of steady flow, these results provide valuable insight into the fluid physics at geometries simulating arterial bifurcations.     

Numerical study on the effect of secondary flow in the human aorta on local shear stresses in abdominal aortic branches

Flow in the aortic arch is characterized primarily by the presence of a strong secondary flow superimposed over the axial flow, skewed axial velocity profiles and diastolic flow reversals. A significant amount of helical flow has also been observed in the descending aorta of humans and in models. In this study a computational model of the abdominal aorta complete with two sets of outflow arteries was adapted for three-dimensional steady flow simulations. The flow through the model was predicted using the Navier-Stokes equations to study the effect that a rotational component of flow has on the general flow dynamics in this vascular segment. The helical velocity profile introduced at the inlet was developed from magnetic resonance velocity mappings taken from a plane transaxial to the aortic arch. Results showed that flow division ratios increased in the first set of branches and decreased in the second set with the addition of rotational flow. Shear stress varied in magnitude with the addition of rotational flow, but the shear stress distribution did not change. No regions of flow separation were observed in the iliac arteries for either case. Helical flow may have a stabilizing effect on the flow patterns in branches in general, as evidenced by the decreased difference in shear stress between the inner and outer walls in the iliac arteries.     

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.     

Effect of geometric variations on pressure loss for a model bifurcation of the human lung airway

Characteristics of pressure loss (ΔP) in human lung airways were numerically investigated using a realistic model bifurcation. Flow equations were numerically solved for the steady inspiratory condition with the tube length, the branching angle and flow velocity being varied over a wide range. In general, the ΔP coefficient K showed a power-law dependence on Reynolds number (Re) and length-to-diameter ratio with a different exponent for Re≥100 than for Re<100. The effect of different branching angles on pressure loss was very weak in the smooth-branching airways.     

Secondary velocity fields in the conducting airways of the human lung

An understanding of flow and dispersion in the human respiratory airways is necessary to assess the toxicological impact of inhaled particulate matter as well as to optimize the design of inhalable pharmaceutical aerosols and their delivery systems. Secondary flows affect dispersion in the lung by mixing solute in the lumen cross section. The goal of this study is to measure and interpret these secondary velocity fields using in vitro lung models. Particle image velocimetry experiments were conducted in a three-generational, anatomically accurate model of the conducting region of the lung. Inspiration and expiration flows were examined under steady and oscillatory flow conditions. Results illustrate secondary flow fields as a function of flow direction, Reynolds number, axial location with respect to the bifurcation junction, generation, branch, phase in the oscillatory cycle, and Womersley number. Critical Dean number for the formation of secondary vortices in the airways, as well as the strength and development length of secondary flow, is characterized. The normalized secondary velocity magnitude was similar on inspiration and expiration and did not vary appreciably with generation or branch. Oscillatory flow fields were not significantly different from corresponding steady flow fields up to a Womersley number of 1 and no instabilities related to shear were detected on flow reversal. These observations were qualitatively interpreted with respect to the simple streaming, augmented dispersion, and steady streaming convective dispersion mechanisms.     

Turbulence modeling in three-dimensional stenosed arterial bifurcations

Under normal healthy conditions, blood flow in the carotid artery bifurcation is laminar. However, in the presence of a stenosis, the flow can become turbulent at the higher Reynolds numbers during systole. There is growing consensus that the transitional k-omega model is the best suited Reynolds averaged turbulence model for such flows. Further confirmation of this opinion is presented here by a comparison with the RNG k-epsilon model for the flow through a straight, nonbifurcating tube. Unlike similar validation studies elsewhere, no assumptions are made about the inlet profile since the full length of the experimental tube is simulated. Additionally, variations in the inflow turbulence quantities are shown to have no noticeable affect on downstream turbulence intensity, turbulent viscosity, or velocity in the k-epsilon model, whereas the velocity profiles in the transitional k-omega model show some differences due to large variations in the downstream turbulence quantities. Following this validation study, the transitional k-omega model is applied in a three-dimensional parametrically defined computer model of the carotid artery bifurcation in which the sinus bulb is manipulated to produce mild, moderate, and severe stenosis. The parametric geometry definition facilitates a powerful means for investigating the effect of local shape variation while keeping the global shape fixed. While turbulence levels are generally low in all cases considered, the mild stenosis model produces higher levels of turbulent viscosity and this is linked to relatively high values of turbulent kinetic energy and low values of the specific dissipation rate. The severe stenosis model displays stronger recirculation in the flow field with higher values of vorticity, helicity, and negative wall shear stress. The mild and moderate stenosis configurations produce similar lower levels of vorticity and helicity.     

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.     

The influence of the non-Newtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model.

Laser Doppler anemometry experiments and finite element simulations of steady flow in a three dimensional model of the carotid bifurcation were performed to investigate the influence of non-Newtonian properties of blood on the velocity distribution. The axial velocity distribution was measured for two fluids: a non-Newtonian blood analog fluid and a Newtonian reference fluid. Striking differences between the measured flow fields were found. The axial velocity field of the non-Newtonian fluid was flattened, had lower velocity gradients at the divider wall, and higher velocity gradients at the non-divider wall. The flow separation, as found with the Newtonian fluid, was absent. In the computations, the shear thinning behavior of the analog blood fluid was incorporated through the Carreau-Yasuda model. The viscoelastic properties of the fluid were not included. A comparison between the experimental and numerical results showed good agreement, both for the Newtonian and the non-Newtonian fluid. Since only shear thinning was included, this seems to be the dominant non-Newtonian property of the blood analog fluid under steady flow conditions.