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.     

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.     

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.     

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

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.     

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.     

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.     

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.     

On the relative importance of rheology for image-based CFD models of the carotid bifurcation

BACKGROUND: Patient-specific computational fluid dynamics (CFD) models derived from medical images often require simplifying assumptions to render the simulations conceptually or computationally tractable. In this study, we investigated the sensitivity of image-based CFD models of the carotid bifurcation to assumptions regarding the blood rheology. METHOD OF APPROACH: CFD simulations of three different patient-specific models were carried out assuming: a reference high-shear Newtonian viscosity, two different non-Newtonian (shear-thinning) rheology models, and Newtonian viscosities based on characteristic shear rates or, equivalently, assumed hematocrits. Sensitivity of wall shear stress (WSS) and oscillatory shear index (OSI) were contextualized with respect to the reproducibility of the reconstructed geometry, and to assumptions regarding the inlet boundary conditions. RESULTS: Sensitivity of WSS to the various rheological assumptions was roughly 1.0 dyn/cm(2) or 8%, nearly seven times less than that due to geometric uncertainty (6.7 dyn/cm(2) or 47%), and on the order of that due to inlet boundary condition assumptions. Similar trends were observed regarding OSI sensitivity. Rescaling the Newtonian viscosity based on time-averaged inlet shear rate served to approximate reasonably, if overestimate slightly, non-Newtonian behavior. CONCLUSIONS: For image-based CFD simulations of the normal carotid bifurcation, the assumption of constant viscosity at a nominal hematocrit is reasonable in light of currently available levels of geometric precision, thus serving to obviate the need to acquire patient-specific rheological data.     

Computational study of pulsatile blood flow in prototype vessel geometries of coronary segments

The spatial and temporal distributions of wall shear stress (WSS) in prototype vessel geometries of coronary segments are investigated via numerical simulation, and the potential association with vascular disease and specifically atherosclerosis and plaque rupture is discussed. In particular, simulation results of WSS spatio-temporal distributions are presented for pulsatile, non-Newtonian blood flow conditions for: (a) curved pipes with different curvatures, and (b) bifurcating pipes with different branching angles and flow division. The effects of non-Newtonian flow on WSS (compared to Newtonian flow) are found to be small at Reynolds numbers representative of blood flow in coronary arteries. Specific preferential sites of average low WSS (and likely atherogenesis) were found at the outer regions of the bifurcating branches just after the bifurcation, and at the outer-entry and inner-exit flow regions of the curved vessel segment. The drop in WSS was more dramatic at the bifurcating vessel sites (less than 5% of the pre-bifurcation value). These sites were also near rapid gradients of WSS changes in space and time – a fact that increases the risk of rupture of plaque likely to develop at these sites. The time variation of the WSS spatial distributions was very rapid around the start and end of the systolic phase of the cardiac cycle, when strong fluctuations of intravascular pressure were also observed. These rapid and strong changes of WSS and pressure coincide temporally with the greatest flexion and mechanical stresses induced in the vessel wall by myocardial motion (ventricular contraction). The combination of these factors may increase the risk of plaque rupture and thrombus formation at these sites.     

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.     

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.     

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.     

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.     

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.