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.     

Pulsatile flow of non-Newtonian fluid in distensible models of human arteries

In addition to biochemical factors, hydromechanical influences are responsible for atherogenesis and deposits of blood platelets at bends and bifurcations of human arteries. Hence the flow patterns were simulated in a true-to-scale three-dimensional bifurcation of a human renal artery model and of an arterial femoralis with Newtonian and non-Newtonian blood like fluid. Investigations were made with steady and pulsatile flow. The velocity profiles (at physiological Re-numbers) were measured after the bifurcations with a laser-Doppler-anemometer. In previous works Newtonian fluids were used to investigate the flow in bends and bifurcations of rigid and elastic simplified models. In this paper, emphasis is placed on the difference between rigid and elastic models and also Newtonian and non Newtonian flow behavior. Differences between Newtonian and non Newtonian fluids may especially be expected to occur after branches where the flow has local strong convective elements such as in reverse zones and flow separation points.     

Flow investigations in a model of a three-dimensional human artery with Newtonian and non-Newtonian fluids. Part I

Together with biochemical factors, fluid mechanical factors play a role in atherogenesis and the deposition of blood platelets at bends and bifurcations in human arteries. Hence, flow patterns were investigated in a simplified 3-dimensional model of a human renal artery bifurcation using Newtonian (aqueous glycerol) and non-Newtonian (aqueous solution of polyacrylamide) fluids. Studies were carried out in steady as well as pulsatile flow at inflow Reynolds numbers of 498 and 951 with flow rate ratios main tube V1: right branch V4: left branch V3 of 1: 0.25: 0.25 and 1: 0.18: 0.18 respectively. The velocity distribution proximal and distal to the bifurcations was measured using a laser-Doppler anemometer. In steady flow, zones of flow separation and reverse flow were observed distal to the bifurcations. In pulsatile flow using non-Newtonian fluids, there was a significant enlargement of these zones. Differences between the Newtonian and non-Newtonian fluids occurred especially distal to the bifurcations. Shear stresses along all measuring positions were computed from the velocity gradients.     

Modelling and simulation of the behaviour of a biofluid in a microchannel biochip separator

This paper reports an investigation into the flow behaviour of a biofluid in a microchannel systems through conceptual analysis and modelling. The application is the design of a microfluidic chip developed for the separation of plasma from blood. The effect of key design features of the microchannels on the flow behaviour of the biofluid is explored. These include geometric features such as the constriction, bending channel, bifurcation and the channel length ratio between the main and side channels. The performance of each design is discussed in terms of separation efficiency of the red blood cells with respect to the rest of the medium. Particular phenomena such as the Fahraeus and Fahraeus–Lindqvist effects, the Zweifach–Fung bifurcation law and the cell-free layer are discussed. In this paper, the fluid is modelled as a single-phase flow assuming either Newtonian or Non-Newtonian behaviour to investigate the effect of the fluid viscosity on both flow and separation efficiency. For a flow rate-controlled Newtonian flow system, it is found that viscosity and outlet pressure have little effect on the velocity distribution through each of the microchannels. For a diluted fluid where the flow in the whole channel system is modelled with a uniform viscosity, less plasma is separated from blood than observed in the non-Newtonian case. This results in an increase in the flow rate ratio between the main and side channels. A comparison of Newtonian and non-Newtonian flows shows that both flows tend to behave identically with an increase in the shear strain rate.     

Modelling and simulation of the behaviour of a biofluid in a microchannel biochip separator

This paper reports an investigation into the flow behaviour of a biofluid in a microchannel systems through conceptual analysis and modelling. The application is the design of a microfluidic chip developed for the separation of plasma from blood. The effect of key design features of the microchannels on the flow behaviour of the biofluid is explored. These include geometric features such as the constriction, bending channel, bifurcation and the channel length ratio between the main and side channels. The performance of each design is discussed in terms of separation efficiency of the red blood cells with respect to the rest of the medium. Particular phenomena such as the Fahraeus and Fahraeus-Lindqvist effects, the Zweifach-Fung bifurcation law and the cell-free layer are discussed. In this paper, the fluid is modelled as a single-phase flow assuming either Newtonian or Non-Newtonian behaviour to investigate the effect of the fluid viscosity on both flow and separation efficiency. For a flow rate-controlled Newtonian flow system, it is found that viscosity and outlet pressure have little effect on the velocity distribution through each of the microchannels. For a diluted fluid where the flow in the whole channel system is modelled with a uniform viscosity, less plasma is separated from blood than observed in the non-Newtonian case. This results in an increase in the flow rate ratio between the main and side channels. A comparison of Newtonian and non-Newtonian flows shows that both flows tend to behave identically with an increase in the shear strain rate.     

Flush mounted hot film anemometer measurement of wall shear stress distal to a tri-leaflet valve for Newtonian and non-Newtonian blood analog fluids

Wall shear stress has been measured by flush-mounted hot film anemometry distal to an Ionescu-Shiley tri-leaflet valve under pulsatile flow conditions. Both Newtonian (aqueous glycerol) and non-Newtonian (aqueous polyacrylamide) blood analog fluids were investigated. Significant differences in the axial distribution of wall shear stress between the two fluids are apparent in flows having nearly identical Reynolds numbers. The Newtonian fluid exhibits a (peak) wall shear rate which is maximized near the valve seat (30 mm) and then decays to a fully developed flow value (by 106 mm). In contrast, the shear rate of the non-Newtonian fluid at 30 mm is less than half that of the Newtonian fluid and at 106 mm is more than twice that of the Newtonian fluid. It is suggested that non-Newtonian rheology influences valve flow patterns either through alterations in valve opening associated with low shear separation zones behind valve leaflets, or because of variations in the rate of jet spreading. More detailed studies are required to clarify the mechanisms. The Newtonian wall shear stresses for this valve are low. The highest value observed anywhere in the aortic chamber was 2.85 N/m2 at a peak Reynolds number of 3694.     

Development of a general method for designing microvascular networks using distribution of wall shear stress

In the present study, theoretical formulations for calculation of optimal bifurcation angle and relationship between the diameters of mother and daughter vessels using the power law model for non-Newtonian fluids are developed. The method is based on the distribution of wall shear stress in the mother and daughter vessels. Also, the effect of distribution of wall shear stress on the minimization of energy loss and flow resistance is considered. It is shown that constant wall shear stress in the mother and daughter vessels provides the minimum flow resistance and energy loss of biological flows. Moreover, the effects of different wall shear stresses in the mother and daughter branches, different lengths of daughter branches in the asymmetric bifurcations and non-Newtonian effect of biological fluid flows on the bifurcation angle and the relationship between the diameters of mother and daughter branches are considered. Using numerical simulations for non-Newtonian models such as power law and Carreau models, the effects of optimal bifurcation angle on the pressure drop and flow resistance of blood flow in the symmetric bifurcation are investigated. Numerical simulations show that optimal bifurcation angle decreases the pressure drop and flow resistance especially for bifurcations at large Reynolds number.     

Pulsatile non-Newtonian flow characteristics in a three-dimensional human carotid bifurcation model.

Numerical analysis of flow phenomena and wall shear stresses in the human carotid artery bifurcation has been carried out using a three-dimensional geometrical model. The primary aim of this study is the detailed discussion of non-Newtonian flow velocity and wall shear stress during the pulse cycle. A comparison of non-Newtonian and Newtonian results is also presented. The applied non-Newtonian behavior of blood is based on measured dynamic viscosity. In the foreground of discussion are the flow characteristics in the carotid sinus. The investigation shows complex flow patterns especially in the carotid sinus where flow separation occurs at the outer wall throughout the systolic deceleration phase. The changing sign of the velocity near the outer sinus wall results in oscillating shear stress during the pulse cycle. At the outer wall of the sinus at maximum diameter level the shear stress ranges from -1.92 N/m2 to 1.22 N/m2 with a time-averaged value of 0.04 N/m2. At the inner wall of the sinus at maximum diameter level the shear stress range is from 1.16 N/m2 to 4.18 N/m2 with a mean of 1.97 N/m2. The comparison of non-Newtonian and Newtonian results indicates unchanged flow phenomena and rather minor differences in the basic flow characteristics.     

An LDA and flow visualization study of pulsatile flow in an aortic bifurcation model

The structure of pulsatile flow in a rigid aortic bifurcation model was studied by means of a flow visualization technique and three-dimensional laser-Doppler anemometry. The model was made of glass, having the same shape as that of the average human aortic bifurcation. It was installed into a mock circulatory loop that generated physiological pulsatile flow. Flow separation was observed during accelerated and decelerated flow periods. Double helical flow existed inside the flow separation in the early accelerated flow period. In the decelerated flow period, disturbed flow appeared behind the separation zone. Flow was strongly disturbed during the back flow period, and then was gradually stabilized in the forward flow period. The flow separation and the disturbances released from the flow separation zone greatly influenced near-wall velocities along the lateral wall. The wave form of the near-wall velocity in the flow separation zone was much different from that observed in the aortic portion and behind the separation zone; for example, the magnitude of the negative peak velocity in the direction of the tube axis was larger than that of the positive one, and mean velocity in a cycle was very low. This abnormal phasic change of the near-wall velocity may be associated with atherogenesis. The three-dimensional velocity measurement is very useful for the detailed analysis of near-wall velocity patterns.     

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.     

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

The pulsatile flow of non-Newtonian fluid in a bifurcation model with a non-planar daughter branch is investigated numerically by using the Carreau-Yasuda model to take into account the shear thinning behavior of the analog blood fluid. The objective of this study is to deal with the influence of the non-Newtonian property of fluid and of out-of-plane curvature in the non-planar daughter vessel on wall shear stress (WSS), oscillatory shear index (OSI), and flow phenomena during the pulse cycle. The non-Newtonian property in the daughter vessels induces a flattened axial velocity profile due to its shear thinning behavior. The non-planarity deflects flow from the inner wall of the vessel to the outer wall and changes the distribution of WSS along the vessel, in particular in systole phase. Downstream of the bifurcation, the velocity profiles are shifted toward the flow divider, and low WSS and high shear stress temporal oscillations characterized by OSI occur on the outer wall region of the daughter vessels close to the bifurcation. Secondary motions become stronger with the addition of the out-of-plane curvature induced by the bending of the vessel, and the secondary flow patterns swirl along the non-planar daughter vessel. A significant difference between the non-Newtonian and the Newtonian pulsatile flow is revealed during the pulse cycle; however, reasonable agreement between the non-Newtonian and the rescaled Newtonian flow is found. Calculated results for the pulsatile flow support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Unsteady flow of Casson fluids through the eccentric circular tube.

As a model of the flow through the false lumen of the dissecting aortic aneurysm, the starting and stopping flows of the Casson fluid through the eccentric circular tube are studied numerically. The calculations are performed by using the one-direction fully developed unsteady flow model. The calculation results show that the non-Newtonian property of the Casson fluid has the following two effects. One is the decrease in the flow rate through the false lumen, the other is the effect to strengthen the non-uniformity of the velocity distribution in the false lumen. Moreover, the calculations based on a triphasic waveform measured by McDonald are performed and the development of the dissecting aortic aneurysm is discussed on the basis of the calculated wall shear stress.     

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.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

The influence of the non-Newtonian properties of blood on the flow in large arteries: unsteady flow in a 90 degrees curved tube.

A numerical and experimental investigation of unsteady entry flow in a 90 degrees curved tube is presented to study the impact of the non-Newtonian properties of blood on the velocity distribution. The time-dependent flow rate for the Newtonian and the non-Newtonian blood analog fluid were identical. For the numerical computation, a Carreau-Yasuda model was employed to accommodate the shear thinning behavior of the Xanthan gum solution. The viscoelastic properties were not taken into account. The experimental results indicate that significant differences between the Newtonian and non-Newtonian fluid are present. The numerical results for both the Newtonian and the non-Newtonian fluid agree well with the experimental results. Since viscoelasticity was not included in the numerical code, shear thinning behavior of the blood analog fluid seems to be the dominant non-Newtonian property, even under unsteady flow conditions. Finally, a comparison between the non-Newtonian fluid model and a Newtonian fluid at a rescaled Reynolds number is presented. The rescaled Reynolds number, based on a characteristic rather than the high-shear rate viscosity of the Xanthan gum solution, was about three times as low as the original Reynolds number. Comparison reveals that the character of flow of the non-Newtonian fluid is simulated quite well by using the appropriate Reynolds number.     

Linear and nonlinear analyses of pulsatile blood flow in a cylindrical tube

A non-Newtonian shear-thinning constitutive relation is proposed to study pulsatile flow of whole blood in a cylindrical tube. The constitutive relation, which satisfies the principle of material frame indifference, is derived from viscometric data obtained from whole blood over a range of hematocrits. Assuming axisymmetric flow in a rigid cylindrical tube of constant diameter, a second-order, nonlinear partial differential equation governing the axial velocity component is obtained. Imposing a periodic pressure gradient, the governing equation was solved numerically using finite difference methods over a range of Stokes values and hematocrits. For a forcing frequency of 1 Hz, results are presented over tube diameters ranging between 0.1 and 2 cm and over hematocrits ranging between 10 and 80%. For a given hematocrit, velocity profiles predicted for the non-Newtonian model under sinusoidal forcing reveal attenuated volume flow rate and enhanced vorticity transport over the tube cross-section relative to a Newtonian fluid having a viscosity corresponding to the high shear-rate limit. For moderate to high Stokes numbers, consistent with flow in large arteries, our results revealed a viscosity distribution that was nearly time invariant. An analytic solution was obtained for a fluid having arbitrarily prescribed radially varying, temporally invariant viscosity and density distributions under arbitrary periodic pressure forcing. Close agreement was observed between our numerical and analytical results when the imposed viscosity distribution was chosen to approximate the time-averaged viscosity distribution predicted by the shear-thinning non-Newtonian model. For St > or approximately= 100, the disparity between our results and those of a Newtonian fluid of constant viscosity grows with a decreasing ratio of the DC to AC components of the pressure-gradient amplitude below 50%. In particular, for any purely oscillatory pressure-gradient (vanishing DC component), the Womersley solution is a particularly poor predictor of the amplitude and phase of wall shear rate for over half of the flow cycle. Under such circumstances, the analytical models presented here provide a simple and accurate means of estimating instantaneous wall shear rate, knowing only the pressure gradient and hematocrit.     

The effect of compliance on wall shear in casts of a human aortic bifurcation

Rigid and compliant casts of a human aortic bifurcation were subjected to physiologically realistic pulsatile fluid flows. At a number of sites near the wall in the approximate median plane of the bifurcation of these models, fluid velocity was measured with a laser Doppler velocimeter, and wall motion (in the case of the compliant cast) was determined with a Reticon linescan camera. The velocity and wall motion data were combined to estimate the instantaneous shear rates at the cast wall. Analysis showed that at the outer walls the cast compliance reduced shear rates, while at the walls of the flow divider the shear rate was increased.     

Pulsatile non-newtonian blood flow in three-dimensional carotid bifurcation models: a numerical study of flow phenomena under different bifurcation angles

Flow and stress patterns in human carotid artery bifurcation models, which differ in the bifurcation angle, are analysed numerically under physiologically relevant flow conditions. The governing Navier-Stokes equations describing pulsatile, three-dimensional flow of an incompressible non-Newtonian fluid are approximated using a pressure correction finite element method, which has been developed recently. The non-Newtonian behaviour of blood is modelled using Casson's relation, based on measured dynamic viscosity. The study concentrates on flow and stress characteristics in the carotid sinus. The results show that the complex flow in the sinus is affected by the angle variation. The magnitude of reversed flow, the extension of the recirculation zone in the outer sinus region and the duration of flow separation during the pulse cycle as well as the resulting wall shear stress are clearly different in the small angle and in the large angle bifurcation. The haemodynamic phenomena, which are important in atherogenesis, are more pronounced in the large angle bifurcation.