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.     

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.     

Shape optimization in unsteady blood flow: a numerical study of non-Newtonian effects

This paper presents a numerical study of non-Newtonian effects on the solution of shape optimization problems involving unsteady pulsatile blood flow. We consider an idealized two dimensional arterial graft geometry. Our computations are based on the Navier-Stokes equations generalized to non-Newtonian fluid, with the modified Cross model employed to account for the shear-thinning behavior of blood. Using a gradient-based optimization algorithm, we compare the optimal shapes obtained using both the Newtonian and generalized Newtonian constitutive equations. Depending on the shear rate prevalent in the domain, substantial differences in the flow as well as in the computed optimal shape are observed when the Newtonian constitutive equation is replaced by the modified Cross model. By varying a geometric parameter in our test case, we investigate the influence of the shear rate on the solution.     

Shape optimization in unsteady blood flow: A numerical study of non-Newtonian effects

This paper presents a numerical study of non-Newtonian effects on the solution of shape optimization problems involving unsteady pulsatile blood flow. We consider an idealized two dimensional arterial graft geometry. Our computations are based on the Navier–Stokes equations generalized to non-Newtonian fluid, with the modified Cross model employed to account for the shear-thinning behavior of blood. Using a gradient-based optimization algorithm, we compare the optimal shapes obtained using both the Newtonian and generalized Newtonian constitutive equations. Depending on the shear rate prevalent in the domain, substantial differences in the flow as well as in the computed optimal shape are observed when the Newtonian constitutive equation is replaced by the modified Cross model. By varying a geometric parameter in our test case, we investigate the influence of the shear rate on the solution.     

Numerical study of the impact of non-Newtonian blood behavior on flow over a two-dimensional backward facing step

Endothelial cell (EC) responsiveness to shear stress is essential for vasoregulation and plays a role in atherogenesis. Although blood is a non-Newtonian fluid, EC flow studies in vitro are typically performed using Newtonian fluids. The goal of the present study was to determine the impact of non-Newtonian behavior on the flow field within a model flow chamber capable of producing flow disturbance and whose dimensions permit Reynolds and Womersley numbers comparable to those present in vivo. We performed two-dimensional computational fluid dynamic simulations of steady and pulsatile laminar flow of Newtonian and non-Newtonian fluids over a backward facing step. In the non-Newtonian simulations, the fluid was modeled as a shear-thinning Carreau fluid. Steady flow results demonstrate that for Re in the range 50-400, the flow recirculation zone downstream of the step is 22-63% larger for the Newtonian fluid than for the non-Newtonian fluid, while spatial gradients of shear stress are larger for the non-Newtonian fluid. In pulsatile flow, the temporal gradients of shear stress within the flow recirculation zone are significantly larger for the Newtonian fluid than for the non-Newtonian fluid. These findings raise the possibility that in regions of flow disturbance, EC mechanotransduction pathways stimulated by Newtonian and non-Newtonian fluids may be different.     

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.     

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.     

Comparative study of Newtonian and non-Newtonian simulations of drug transport in a model drug-eluting stent

To elucidate the difference between Newtonian and shear thinning non-Newtonian assumptions of blood in the analysis of DES drug delivery, we numerically simulated the local flow pattern and the concentration distribution of the drug at the lumen-tissue interface for a structurally simplified DES deployed in a curved segment of an artery under pulsatile blood flow conditions. The numerical results showed that when compared with the Newtonian model, the Carreau (shear thinning) model could lead to some differences in the luminal surface drug concentration in certain areas along the outer wall of the curved vessel. In most areas of the vessel, however, there were no significant differences between the 2 models. Particularly, no significant difference between the two models was found in terms of the area-averaged luminal surface drug concentration. Therefore, we believe that the shear thinning property of blood may play little roles in DES drug delivery. Nevertheless, before we draw the conclusion that Newtonian assumption of blood can be used to replace its non-Newtonian one for the numerical simulation of drug transport in the DES implanted coronary artery, other more complex mechanical properties of blood such as its thixotropic behavior should be tested.     

Modeling pulsatile flow in aortic aneurysms: effect of non-Newtonian properties of blood

Pulsatile flow in an axisymmetric rigid-walled model of an abdominal aorta aneurysm was analyzed numerically for various aneurysm dilations using physiologically realistic resting waveform at time-averaged Reynolds number of 300 and peak Reynolds number of 1607. Discretization of the governing equations was achieved using a finite element scheme based on the Galerkin method of weighted residuals. Comparisons with previously published work on the basis of special cases were performed and found to be in excellent agreement. Our findings indicate that the velocity fields are significantly affected by non-Newtonian properties in pathologically altered configurations. Non-Newtonian fluid shear stress is found to be greater than Newtonian fluid shear stress during peak systole. Further, the maximum shear stress is found to occur near the distal end of AAA during peak systole. The impact of non-Newtonian blood flow characteristics on pressure compared to Newtonian model is found insignificant under resting conditions. Viscous and inertial forces associated with blood flow are responsible for the changes in the wall that result in thrombus deposition and dilation while rupture of AAA is more likely determined by much larger mechanical stresses imposed by pulsatile pressure on the wall of AAA.     

Wall shear stress in backward-facing step flow of a red blood cell suspension.

An experimental investigation of the wall shear stress distribution downstream of a backward-facing step is carried out. The wall shear stress distribution was determined by measuring the deformation of a gel layer, attached to the wall downstream of the step. Speckle pattern interferometry was applied to measure the deformation of the gel layer. The measured deformation, combined with the properties of the gel layer, served as an input for a finite element solid mechanics computation to determine the stress distribution in the gel layer. The wall shear stress, required to generate the measured deformation of the gel layer, was determined from these computations. A Newtonian buffer solution and a non-Newtonian red blood cell suspension were used as measuring fluids. The deformation of the gel layer was determined for a Newtonian buffer solution to evaluate the method and to obtain the properties of the gel layer. Subsequently, the wall shear stress distribution for the non-Newtonian red blood cell suspension was determined for three different flow rates. The inelastic non-Newtonian Carreau-Yasuda model served as constitutive model for the red blood cell suspension. Using this model, the velocity and wall shear stress distribution were computed by means of a finite element fluid mechanics computation. From the comparison between the numerical and the experimental results, it can be concluded that wall shear stresses, induced by the red blood cell suspension, can be modeled accurately by employing a Carreau-Yasuda model.     

Effect of non-Newtonian and pulsatile blood flow on mass transport in the human aorta

To investigate the effects of both non-Newtonian behavior and the pulsation of blood flow on the distributions of luminal surface LDL concentration and oxygen flux along the wall of the human aorta, we numerically compared a non-Newtonian model with the Newtonian one under both steady flow and in vivo pulsatile flow conditions using a human aorta model constructed from MRI images. The results showed that under steady flow conditions, although the shear thinning non-Newtonian nature of blood could elevate wall shear stress (WSS) in most regions of the aorta, especially areas with low WSS, it had little effect on luminal surface LDL concentration (c(w)) in most regions of the aorta. Nevertheless, it could significantly enhance c(w) in areas with high luminal surface LDL concentration through the shear dependent diffusivity of LDLs. For oxygen transport, the shear thinning non-Newtonian nature of blood could slightly reduce oxygen flux in most regions of the aorta, but this effect became much more apparent in areas with already low oxygen flux. The pulsation of blood flow could significantly reduce c(w) and enhance oxygen flux in these disturbed places. In most other regions of the aorta, the oxygen flux was also significantly higher than that for the steady flow simulation. In conclusion, the shear shining non-Newtonian nature of blood has little effect on LDL and oxygen transport in most regions of the aorta, but in the atherogenic-prone areas where luminal surface LDL concentration is high and oxygen flux is low, its effect is apparent. Similar is for the effect of pulsatile flow on the transport of LDLs. But, the pulsation of blood flow can apparently affect oxygen flux in the aorta, especially in areas with low oxygen flux.     

Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows

Effects of the non-Newtonian viscosity of blood on a flow in a coronary arterial casting of man were studied numerically using a finite element method. Various constitutive models were examined to model the non-Newtonian viscosity of blood and their model constants were summarized. A method to incorporate the non-Newtonian viscosity of blood was introduced so that the viscosity could be calculated locally. The pressure drop, wall shear stress and velocity profiles for the case of blood viscosity were compared for the case of Newtonian viscosity (0.0345 poise). The effect of the non-Newtonian viscosity of blood on the overall pressure drop across the arterial casting was found to be significant at a flow of the Reynolds number of 100 or less. Also in the region of flow separation or recirculation, the non-Newtonian viscosity of blood yields larger wall shear stress than the Newtonian case. The origin of the non-Newtonian viscosity of blood was discussed in relation to the viscoelasticity and yield stress of blood.     

The flow of power law fluids in elastic networks and porous media

The flow of power law fluids, which include shear thinning and shear thickening as well as Newtonian as a special case, in networks of interconnected elastic tubes is investigated using a residual-based pore scale network modeling method with the employment of newly derived formulae. Two relations describing the mechanical interaction between the local pressure and local cross-sectional area in distensible tubes of elastic nature are considered in the derivation of these formulae. The model can be used to describe shear dependent flows of mainly viscous nature. The behavior of the proposed model is vindicated by several tests in a number of special and limiting cases where the results can be verified quantitatively or qualitatively. The model, which is the first of its kind, incorporates more than one major nonlinearity corresponding to the fluid rheology and conduit mechanical properties, that is non-Newtonian effects and tube distensibility. The formulation, implementation, and performance indicate that the model enjoys certain advantages over the existing models such as being exact within the restricting assumptions on which the model is based, easy implementation, low computational costs, reliability, and smooth convergence. The proposed model can, therefore, be used as an alternative to the existing Newtonian distensible models; moreover, it stretches the capabilities of the existing modeling approaches to reach non-Newtonian rheologies.     

Computational fluid dynamics in abdominal aorta bifurcation: non-Newtonian versus Newtonian blood flow in a real case study

Hemodynamic in abdominal aorta bifurcation was investigated in a real case using computational fluid dynamics. A Newtonian and non-Newtonian (Walburn-Schneck) viscosity models were compared. The geometrical model was obtained by 3D reconstruction from CT-scan and hemodynamic parameters obtained by laser-Doppler. Blood was assumed incompressible fluid, laminar flow in transient regime and rigid vessel wall. Finite volume-based was used to study the velocity, pressure, wall shear stress (WSS) and viscosity throughout cardiac cycle. Results obtained with Walburn-Schneck’s model, during systole, present lower viscosity due to shear thinning behavior. Furthermore, there is a significant difference between the results obtained by the two models for a specific patient. During the systole, differences are more pronounced and are preferably located in the tortuous regions of the artery. Throughout the cardiac cycle, the WSS amplitude between the systole and diastole is greater for the Walburn-Schneck’s model than for the Newtonian model. However, the average viscosity along the artery is always greater for the non-Newtonian model, except in the systolic peak. The hemodynamic model is crucial to validate results obtained with CFD and to explore clinical potential.     

An experimental study of Newtonian and non-Newtonian flow dynamics in a ventricular assist device

The fluid dynamic behavior of a Newtonian water/glycerol solution, a non-Newtonian polymer (separan) solution, and bovine blood were compared in the Penn State Electrical Ventricular Assist Device (EVAD). Pulsed doppler ultrasound velocimetry was used to measure velocities in the near wall region (0.95-2.7 mm) along the perimeter of the pump. Mean velocity, turbulence intensity, local and convective acceleration, and shear rate were calculated from the PDU velocity measurements. Flow visualization provided qualitative information about the general flow patterns in the EVAD. Results indicate that water/glycerol does not accurately model the flow characteristics of bovine blood in the EVAD. The non-Newtonian separan solution produced results closer to those of the bovine blood than did the water/glycerol solution. Near wall velocity magnitudes for the separan were similar to those of the bovine blood, but the profile shapes differed for portions of the pump cycle. All three fluids exhibited periods of stagnation. Bovine blood results indicated the presence of a desired rotational washout pattern at midsystole, while results with the other fluids did not show this feature.     

Accurate prediction of wall shear stress in a stented artery: newtonian versus non-newtonian models

A significant amount of evidence linking wall shear stress to neointimal hyperplasia has been reported in the literature. As a result, numerical and experimental models have been created to study the influence of stent design on wall shear stress. Traditionally, blood has been assumed to behave as a Newtonian fluid, but recently that assumption has been challenged. The use of a linear model; however, can reduce computational cost, and allow the use of Newtonian fluids (e.g., glycerine and water) instead of a blood analog fluid in an experimental setup. Therefore, it is of interest whether a linear model can be used to accurately predict the wall shear stress caused by a non-Newtonian fluid such as blood within a stented arterial segment. The present work compares the resulting wall shear stress obtained using two linear and one nonlinear model under the same flow waveform. All numerical models are fully three-dimensional, transient, and incorporate a realistic stent geometry. It is shown that traditional linear models (based on blood's lowest viscosity limit, 3.5 Pa s) underestimate the wall shear stress within a stented arterial segment, which can lead to an overestimation of the risk of restenosis. The second linear model, which uses a characteristic viscosity (based on an average strain rate, 4.7 Pa s), results in higher wall shear stress levels, but which are still substantially below those of the nonlinear model. It is therefore shown that nonlinear models result in more accurate predictions of wall shear stress within a stented arterial segment.     

Numerical simulation of blood pulsatile flow in a stenosed carotid artery using different rheological models

Symmetrical 30-60% stenosis in a common carotid artery under unsteady flow condition for Newtonian and six non-Newtonian viscosity models are investigated numerically. Results show power-law model produces higher deviations, in terms of velocity and wall shear stress in comparison with other models while generalized power-law and modified-Casson models are more prone to Newtonian state. Comparing separation length of recirculation region at different critical points of cardiac cycle confirms the necessity of considering blood flow in unsteady mode. Increasing stenosis intensity causes flow patterns more disturbed downstream of the stenosis and WSS appear to develop remarkably at the stenosis throat.     

Flow of a blood analogue fluid in a compliant abdominal aortic aneurysm model: Experimental modelling

The aim of this work is to develop a unique in vitro set-up in order to analyse the influence of the shear thinning fluid-properties on the flow dynamics within the bulge of an abdominal aortic aneurysm (AAA). From an experimental point of view, the goals are to elaborate an analogue shear thinning fluid mimicking the macroscopic blood behaviour, to characterise its rheology at low shear rates and to propose an experimental device able to manage such an analogue fluid without altering its feature while reproducing physiological flow rate and pressure, through compliant AAA. Once these experimental prerequisites achieved, the results obtained in the present work show that the flow dynamics is highly dependent on the fluid rheology. The main results point out that the propagation of the vortex ring, generated in the AAA bulge, is slower for shear thinning fluids inducing a smaller travelled distance by the vortex ring so that it never impacts the anterior wall in the distal region, in opposition to Newtonian fluids. Moreover, scalar shear rate values are globally lower for shear thinning fluids inducing higher maximum stress values than those for the Newtonian fluids. Consequently, this work highlights that a Newtonian fluid model is finally inadequate to obtain a reliable prediction of the flow dynamics within AAA.     

Non-Newtonian Blood Flow in Left Coronary Arteries with Varying Stenosis: A Comparative Study

This paper presents Computational fluid dynamic (CFD) analysis of blood flow in three different 3-D models of left coronary artery (LCA). A comparative study of flow parameters (pressure distribution, velocity distribution and wall shear stress) in each of the models is done for a non-Newtonian (Carreau) as well as the Newtonian nature of blood viscosity over a complete cardiac cycle. The difference between these two types of behavior of blood is studied for both transient and steady states of flow. Additionally, flow parameters are compared for steady and transient boundary conditions considering blood as non-Newtonian fluid. The study shows that the highest wall shear stress (WSS), velocity and pressure are found in artery having stenosis in all the three branches of LCA. The use of Newtonian blood model is a good approximation for steady as well as transient blood flow boundary conditions if shear rate is above 100 s-1. However, the assumption of steady blood flow results in underestimating the values of flow parameters such as wall shear stress, pressure and velocity.