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.     

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.     

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

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.     

Influence of non-Newtonian properties of blood on the wall shear stress in human atherosclerotic right coronary arteries

The objective of this work is to investigate the effect of non-Newtonian properties of blood on the wall shear stress (WSS) in atherosclerotic coronary arteries using both Newtonian and non-Newtonian models. Numerical simulations were performed to examine how the spatial and temporal WSS distributions are influenced by the stenosis size, blood viscosity, and flow rate. The computational results demonstrated that blood viscosity properties had considerable effect on the magnitude of the WSS, especially where disturbed flow was observed. The WSS distribution is highly non-uniform both temporally and spatially, especially in the stenotic region. The maximum WSS occurred at the proximal side of the stenosis, near the outer wall in the curved artery with no stenosis. The lumen area near the inner wall distal to the stenosis region experienced a lower WSS during the entire cardiac cycle. Among the factors of stenosis size, blood viscosity, and flow rate, the size of the stenosis has the most significant effect on the spatial and temporal WSS distributions qualitatively and quantitatively.     

A study of non-Newtonian aspects of blood flow through stenosed arteries and its applications in arterial diseases

Blood flow through a stenosed artery has been investigated in this paper. Blood has been represented by a non-Newtonian fluid obeying Herschel-Bulkley equation. This model has been used to study the influence of the fluid behaviour index n, shear-dependent nonlinear viscosity K and the yield stress tau H in blood flow through stenosed arteries. The variation of the wall shear stress and the flow resistance with n, K and tau H has been shown graphically. It is observed that the wall shear stress and the flow resistance increase in Herschel-Bulkley fluid in comparison with corresponding Newtonian fluid. It is of interest to note that, in the present model, the thickness of the plug core varies with the axial distance z in the stenotic region. Finally, some biological implications of the present model for some arterial diseases have been briefly discussed.     

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.     

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.     

Realistic non-Newtonian viscosity modelling highlights hemodynamic differences between intracranial aneurysms with and without surface blebs

Most computational fluid dynamic (CFD) simulations of aneurysm hemodynamics assume constant (Newtonian) viscosity, even though blood demonstrates shear-thinning (non-Newtonian) behavior. We sought to evaluate the effect of this simplifying assumption on hemodynamic forces within cerebral aneurysms, especially in regions of low wall shear stress, which are associated with rupture. CFD analysis was performed for both viscosity models using 3D rotational angiography volumes obtained for 26 sidewall aneurysms (12 with blebs, 12 ruptured), and parametric models incorporating blebs at different locations (inflow/outflow zone). Mean and lowest 5% values of time averaged wall shear stress (TAWSS) computed over the dome were compared using Wilcoxon rank-sum test. Newtonian modeling not only resulted in higher aneurysmal TAWSS, specifically in areas of low flow and blebs, but also showed no difference between aneurysms with or without blebs. In contrast, for non-Newtonian analysis, bleb-bearing aneurysms showed significantly lower 5% TAWSS compared to those without (p=0.005), despite no significant difference in mean dome TAWSS (p=0.32). Non-Newtonian modeling also accentuated the differences in dome TAWSS between ruptured and unruptured aneurysms (p<0.001). Parametric models further confirmed that realistic non-Newtonian viscosity resulted in lower bleb TAWSS and higher focal viscosity, especially when located in the outflow zone. The results show that adopting shear-thinning non-Newtonian blood viscosity in CFD simulations of intracranial aneurysms uncovered hemodynamic differences induced by bleb presence on aneurysmal surfaces, and significantly improved discriminant statistics used in risk stratification. These findings underline the possible implications of using a realistic model of blood viscosity in predictive computational hemodynamics.     

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.     

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.     

Local hemodynamic analysis after coronary stent implantation based on Euler-Lagrange method

Coronary stents are deployed to treat the coronary artery disease (CAD) by reopening stenotic regions in arteries to restore blood flow, but the risk of the in-stent restenosis (ISR) is high after stent implantation. One of the reasons is that stent implantation induces changes in local hemodynamic environment, so it is of vital importance to study the blood flow in stented arteries. Based on regarding the red blood cell (RBC) as a rigid solid particle and regarding the blood (including RBCs and plasma) as particle suspensions, a non-Newtonian particle suspensions model is proposed to simulate the realistic blood flow in this work. It considers the blood’s flow pattern and non-Newtonian characteristic, the blood cell-cell interactions, and the additional effects owing to the bi-concave shape and rotation of the RBC. Then, it is compared with other four common hemodynamic models (Newtonian single-phase flow model, Newtonian Eulerian two-phase flow model, non-Newtonian single-phase flow model, non-Newtonian Eulerian two-phase flow model), and the comparison results indicate that the models with the non-Newtonian characteristic are more suitable to describe the realistic blood flow. Afterwards, based on the non-Newtonian particle suspensions model, the local hemodynamic environment in stented arteries is investigated. The result shows that the stent strut protrusion into the flow stream would be likely to produce the flow stagnation zone. And the stent implantation can make the pressure gradient distribution uneven. Besides, the wall shear stress (WSS) of the region adjacent to every stent strut is lower than 0.5 Pa, and along the flow direction, the low-WSS zone near the strut behind is larger than that near the front strut. What’s more, in the regions near the struts in the proximal of the stent, the RBC particle stagnation zone is easy to be formed, and the erosion and deposition of RBCs are prone to occur. These hemodynamic analyses illustrate that the risk of ISR is high in the regions adjacent to the struts in the proximal and the distal ends of the stent when compared with struts in other positions of the stent. So the research can provide a suggestion on the stent design, which indicates that the strut structure in these positions of a stent should be optimized further.

     

Non-Newtonian flow patterns associated with an arterial stenosis.

A non-Newtonian constitutive equation for blood has been introduced in this paper. Using this equation, blood flow attributes such as velocity profiles, flowrate, pressure gradient, and wall shear stress in both straight and stenotic (constricted) tubes have been examined. Results showed that compared with Newtonian flow at the same flowrate, the non-Newtonian normally features larger pressure gradient, higher wall shear stress, and different velocity profile, especially in stenotic tube. In addition, the non-Newtonian stenotic flow appears to be more stable than Newtonian flow.     

Effect of hematocrit on wall shear rate in oscillatory flow: do the elastic properties of blood play a role?

Porcine blood was used to examine the relationship between hematocrit levels and wall shear rate patterns in straight and curved artery models under fixed oscillatory flow conditions characteristic of larger arteries. It is demonstrated that porcine blood models both the viscous and elastic components of the 2 Hz complex viscosity of human blood quite accurately over a broad range of shear rates (1-1000 s-1) and hematocrits (20%-80%). For a fixed oscillatory flow waveform (Poiseuille peak shear rate = 168 s-1; mean shear rate 84 s-1), increases in hematocrit produced a decrease in the peak wall shear rate in both the straight and curved artery models and a corresponding decrease in wall shear rate reversal on the inside wall of the curved artery model. The same trends were also observed for oscillatory flows of aqueous glycerin solutions of increasing viscosity in the range of viscosity of the blood samples tested. Aqueous glycerin solutions produced wall shear rate waveforms of the same magnitude and shape as the porcine blood. This indicates that variations in the shear rate, and therefore the shear stress, were caused primarily by changes in the viscous and not the elastic properties of blood. The results suggest that simple Newtonian fluids may be sufficient for in vitro determination of the first order effects to be expected of human blood flow in large vessels having complex geometries and shear rates in or above the range of the present study.     

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.     

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.