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.     

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.     

Non-Newtonian blood flow in human right coronary arteries: transient simulations

This study looks at pulsatile blood flow through four different right coronary arteries, which have been reconstructed from biplane angiograms. A non-Newtonian blood model (the Generalised Power Law), as well as the usual Newtonian model of blood viscosity, is used to study the wall shear stress in each of these arteries over the entire cardiac cycle. The difference between Newtonian and non-Newtonian blood models is also studied over the whole cardiac cycle using the recently generalised global non-Newtonian importance factor. In addition, the flow is studied by considering paths of massless particles introduced into the flow field. The study shows that, when studying the wall shear stress distribution for transient blood flow in arteries, the use of a Newtonian blood model is a reasonably good approximation. However, to study the flow within the artery in greater detail, a non-Newtonian model is more appropriate.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Non-Newtonian blood flow in human right coronary arteries: steady state simulations

This study looks at blood flow through four different right coronary arteries, which have been reconstructed from bi-plane angiograms. Five non-Newtonian blood models, as well as the usual Newtonian model of blood viscosity, are used to study the wall shear stress in each of these arteries at a particular point in the cardiac cycle. It was found that in the case of steady flow in a given artery, the pattern of wall shear stress is consistent across all models. The magnitude of wall shear stress, however, is influenced by the model used and correlates with graphs of shear stress versus strain for each model. For mid-range velocities of around 0.2 m s(-1) the models are virtually indistinguishable. Local and global non-Newtonian importance factors are introduced, in an attempt to quantify the types of flows where non-Newtonian behaviour is significant. It is concluded that, while the Newtonian model of blood viscosity is a good approximation in regions of mid-range to high shear, it is advisable to use the Generalised Power Law model (which tends to the Newtonian model in those shear ranges in any case) in order to achieve better approximation of wall shear stress at low shear.     

Numerical simulation of non-Newtonian blood flow dynamics in human thoracic aorta

Three non-Newtonian blood viscosity models plus the Newtonian one are analysed for a patient-specific thoracic aorta anatomical model under steady-state flow conditions via wall shear stress (WSS) distribution, non-Newtonian importance factors, blood viscosity and shear rate. All blood viscosity models yield a consistent WSS distribution pattern. The WSS magnitude, however, is influenced by the model used. WSS is found to be the lowest in the vicinity of the three arch branches and along the distal walls of the branches themselves. In this region, the local non-Newtonian importance factor and the blood viscosity are elevated, and the shear rate is low. The present study revealed that the Newtonian assumption is a good approximation at mid-and-high flow velocities, as the greater the blood flow, the higher the shear rate near the arterial wall. Furthermore, the capabilities of the applied non-Newtonian models appeared at low-flow velocities. It is concluded that, while the non-Newtonian power-law model approximates the blood viscosity and WSS calculations in a more satisfactory way than the other non-Newtonian models at low shear rates, a cautious approach is given in the use of this blood viscosity model. Finally, some preliminary transient results are presented.     

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.     

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 simulation of LDL mass transfer in a common carotid artery under pulsatile flows

Symmetrical 30-60% stenosis in carotid artery with a semi-permeable wall under steady/unsteady flows for Newtonian/non-Newtonian fluids is investigated numerically. The results show that the unsteadiness of blood flow, blood pressure rise and LDL component size increase the luminal concentration, LC, of the surface. The maximum LC occurring immediately after the separation point and the non-Newtonian fluid predicts higher LDL accumulation. LC decreased as the recirculation length is increased and reaches maximum at 40% stenosis. This process is used to estimate the time-dependent growth of the arterial wall.     

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

Axial dispersion in packed bed reactors involving viscoinelastic and viscoelastic non-Newtonian fluids

Axial dispersion is an important parameter in the performance of packed bed reactors. A lot of fluids exhibit non-Newtonian behaviour but the effect of rheological parameters on axial dispersion is not available in literature. The effect of rheology on axial dispersion has been analysed for viscoinelastic and viscoelastic non-Newtonian fluids. Aqueous solutions of carboxymethyl cellulose and polyacrylamide have been chosen to represent viscoinelastic and viscoelastic liquid-phases. Axial dispersion has been measured in terms of Bo_L number. The single parameter axial dispersion model has been applied to analyse RTD response curve. The Bo_L numbers were observed to increase with increase in liquid flow rate and consistency index ‘K’ for viscoinelastic as well as viscoelastic fluids. Bodenstein correlation for Newtonian fluids proposed has been modified to account for the effect of fluid rheology. Further, Weissenberg number is introduced to quantify the effect of viscoelasticity.