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.     

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.     

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.     

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.     

Multiphase hemodynamic simulation of pulsatile flow in a coronary artery

A multiphase transient non-Newtonian three-dimensional (3-D) computational fluid dynamics (CFD) simulation has been performed for pulsatile hemodynamics in an idealized curved section of a human coronary artery. We present the first prediction, to the authors' knowledge, of particulate buildup on the inside curvature using the multiphase theory of dense suspension hemodynamics. In this study, the particulates are red blood cells (RBCs). The location of RBC buildup on the inside curvature correlates with lower wall shear stress (WSS) relative to the outside curvature. These predictions provide insight into how blood-borne particulates interact with artery walls and hence, have relevance for understanding atherogenesis since clinical observations show that atherosclerotic plaques generally form on the inside curvatures of arteries. The buildup of RBCs on the inside curvature is driven by the secondary flow and higher residence times. The higher viscosity in the central portion of the curved vessel tends to block their flow, causing them to migrate preferentially through the boundary layer. The reason for this is the nearly neutrally buoyant nature of the dense two-phase hemodynamic flow. The two-phase non-Newtonian viscosity model predicts greater shear thinning than the single-phase non-Newtonian model. Consequently, the secondary flow induced in the curvature is weaker. The waveforms for computed hemodynamic parameters, such as hematocrit, WSS, and viscosity, follow the prescribed inlet velocity waveforms. The lower oscillatory WSS produced on the inside curvature has implications for understanding thickening of the intimal layer.     

Fluid flow structure in arterial bypass anastomosis

The fluid flow through a stenosed artery and its bypass graft in an anastomosis can substantially influence the outcome of bypass surgery. To help improve our understanding of this and related issues, the steady Navier-Stokes flows are computed in an idealized arterial bypass system with partially occluded host artery. Both the residual flow issued from the stenosis--which is potentially important at an earlier stage after grafting--and the complex flow structure induced by the bypass graft are investigated. Seven geometric models, including symmetric and asymmetric stenoses in the host artery, and two major aspects of the bypass system, namely, the effects of area reduction and stenosis asymmetry, are considered. By analyzing the flow characteristics in these configurations, it is found that (1) substantial area reduction leads to flow recirculation in both upstream and downstream of the stenosis and in the host artery near the toe, while diminishes the recirculation zone in the bypass graft near the bifurcation junction, (2) the asymmetry and position of the stenosis can affect the location and size of these recirculation zones, and (3) the curvature of the bypass graft can modify the fluid flow structure in the entire bypass system.     

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.     

Hemodynamic Analysis in an Idealized Artery Tree: Differences in Wall Shear Stress between Newtonian and Non-Newtonian Blood Models

Development of many conditions and disorders, such as atherosclerosis and stroke, are dependent upon hemodynamic forces. To accurately predict and prevent these conditions and disorders hemodynamic forces must be properly mapped. Here we compare a shear-rate dependent fluid (SDF) constitutive model, based on the works by Yasuda et al in 1981, against a Newtonian model of blood. We verify our stabilized finite element numerical method with the benchmark lid-driven cavity flow problem. Numerical simulations show that the Newtonian model gives similar velocity profiles in the 2-dimensional cavity given different height and width dimensions, given the same Reynolds number. Conversely, the SDF model gave dissimilar velocity profiles, differing from the Newtonian velocity profiles by up to 25% in velocity magnitudes. This difference can affect estimation in platelet distribution within blood vessels or magnetic nanoparticle delivery. Wall shear stress (WSS) is an important quantity involved in vascular remodeling through integrin and adhesion molecule mechanotransduction. The SDF model gave a 7.3-fold greater WSS than the Newtonian model at the top of the 3-dimensional cavity. The SDF model gave a 37.7-fold greater WSS than the Newtonian model at artery walls located immediately after bifurcations in the idealized femoral artery tree. The pressure drop across arteries reveals arterial sections highly resistive to flow which correlates with stenosis formation. Numerical simulations give the pressure drop across the idealized femoral artery tree with the SDF model which is approximately 2.3-fold higher than with the Newtonian model. In atherosclerotic lesion models, the SDF model gives over 1 Pa higher WSS than the Newtonian model, a difference correlated with over twice as many adherent monocytes to endothelial cells from the Newtonian model compared to the SDF model.     

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.     

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.     

Numerical modeling of simulated blood flow in idealized composite arterial coronary grafts: steady state simulations

This paper presents a comparative study of simulated blood flow in different configurations of simplified composite arterial coronary grafts (CACGs). Even though the composite arterial grafting is increasingly used in cardiac surgery, it is still questionable whether or not the blood flow in such grafts can adequately meet the demands of the native myocardial circulation. A computational fluid dynamics (CFD) model was developed to conduct computer-based studies of simulated blood flow in four different geometric configurations of CACGs, corresponding to routinely used networks in cardiac surgery coronary grafts (T, Y, Pi and sequential). The flow was assumed three-dimensional, laminar and steady and the fluid as Newtonian, while the vessel walls were considered as inelastic and impermeable. It was concluded that local haemodynamics, practically described by velocity, pressure drop, wall shear stress (WSS) and flow rates, may be strongly influenced by the local geometry, especially at the anastomotic sites. The computations were made at mean flow rates of 37.5, 75 and 150ml/min. The side-branch outflow rates, computed for each bypass graft, showed noticeable differences. The results, which were found both qualitatively and quantitatively consistent with other studies, indicate that the Pi-graft exhibits significantly less uniform distribution of outflow rates than the other geometric configurations. Moreover, prominent variations in WSS and velocity distribution among the assessed CACGs were predicted, showing remarkable flow interactions among the arterial branches. The lowest shear stress regions were found on the lateral walls of bifurcations, which are predominantly susceptible to the occurrence of coronary artery disease (CAD). In contrast, the highest WSS were observed at the turn of the arterial branches.     

Hemodynamic Based Surgical Decision on Sequential Graft and Y-Type Graft in Coronary Artery Bypass Grafting

Purpose: Sequential graft and Y-type graft are two different surgical procedures in coronary artery bypass grafting (CABG). The hemodynamic environment of them are different, that may cause different short-term surgical result and long-term patency. In this study, the short-term and long-term result of sequential and Y-type graft was discussed by comparing the hemodynamics of them. Materials and Methods: Two postoperative 3-dimensional (3D) models were built by applying different graft on a patient-specific 3D model with serious stenosis. Then zero-dimensional (0D)/3D coupled simulation was carried out by coupling the postoperative 3D models with a 0D lumped parameter model of the cardiovascular system. Results: The flow rate of native coronary arteries and grafts are all calculated and illustrated in this paper. No significant difference of the native coronary arteries flow and graft flow exists between two surgical procedures. The wall shear stress (WSS) and streamline were also depicted. The graft WSS of sequential graft is 19.1% higher than Y-type graft. While flow separation appears at the bifurcation of Y-type graft. Conclusion: The short-term outcomes of sequential graft and Y-type graft are almost the same. But it can be found from the hemodynamics factors that the longterm patency of the sequential graft is better.     

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.     

Three-dimensional numerical simulations of flow through a stenosed coronary bypass

This work analyzes the flow patterns at the anastomosis of a stenosed coronary bypass. Three-dimensional numerical simulations are performed using a finite elements method. We consider a geometrical model of the host coronary artery with and without a 75% severity stenosis for three different locations from the anastomosis. The flow features - velocity profiles, secondary motions and wall shear stresses - are compared for different configurations of the flow rate and of the distance of the anastomosis from the site of occlusion (called distance of grafting). The combination of the junction flow effects - counter rotating vortices - with the stenosis effects - confined jet flow - is particularly important when the distance of grafting is short. Given that the residual flow issued from the pathologic stenosis being non-negligible after two weeks grafting, models without stenosis cannot predict the evolution of the wall shear stress in the vicinity of the anastomosis.     

Flow and wall shear stress in end-to-side and side-to-side anastomosis of venous coronary artery bypass grafts

Purpose  

Coronary artery bypass graft (CABG) surgery represents the standard treatment of advanced coronary artery disease. Two major types of anastomosis exist to connect the graft to the coronary artery, i.e., by using an end-to-side or a side-to-side anastomosis. There is still controversy because of the differences in the patency rates of the two types of anastomosis. The purpose of this paper is to non-invasively quantify hemodynamic parameters, such as mass flow and wall shear stress (WSS), in end-to-side and side-to-side anastomoses of patients with CABG using computational fluid dynamics (CFD).     

Analysis of blood flow in an out-of-plane CABG model

Coronary artery bypass graft (CABG) is a routine surgical treatment for ischemic and infarcted myocardium. A large number of CABG fail postoperatively because of intimal hyperplasia within months or years. The cause of this failure is thought to be partly related to the flow patterns and shear stresses acting on the endothelial cells. An accurate representation of the flow field and associated wall shear stress (WSS) requires a detailed three-dimensional (3D) model of the CABG. The purpose of this study is to present a detailed analysis of blood flow in a 3D aorto/left CABG, bypassing the occluded left anterior descending coronary (LAD) artery. The analysis takes into account the influence of the out-of-plane geometry of the graft. The finite volume technique was employed to model the 3D blood flow pattern to determine the velocity and WSS distributions. This study presents the flow field distributions of the velocity and WSS at four instances of the cardiac cycle, two in systole and two in diastole. Our results reveal that the CABG geometry has a significant effect on the velocity distribution. The axial velocity profiles at different instances of the cardiac cycle exhibit strong skewing; significant secondary flow and vortex structures are seen in the in-plane velocity patterns. The maximum WSS on the bed of the occluded LAD artery opposite to the graft junction is 14 Pa in middiastole, whereas there is a significantly lower and more uniform distribution of WSS on the bed of the anastomosis. The present results indicate that nonplanarity of the blood vessel along with the inflow conditions has a substantial effect on the fluid mechanics of CABG that contribute to the patency of graft.     

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 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 simulations of pulsatile flow in an end-to-side anastomosis model

A potential interaction between the local hemodynamics and the artery wall response has been suggested for vascular graft failure by intimal hyperplasia (IH). Among the various hemodynamic factors, wall shear stress has been implicated as the primary factor responsible for the development of IH. In order to explore the role of hemodynamics in the formation of IH in end-to-side anastomosis, computational fluid dynamics is employed. To validate the numerical simulations, comparisons with existing experimental data are performed for both steady and pulsatile flows. Generally, good agreement is observed with the velocity profiles whereas some discrepancies are found in wall shear stress (WSS) distributions. Using the same end-to-side anastomosis geometry, numerical simulations are extended using a femoral artery waveform to identify the possible role of unsteady hemodynamics. In the current simulations, Carreau-Yasuda model is used to account for the non-Newtonian nature of blood. Computations indicated a disturbed flow field at the artery-graft junction leading to locally elevated shear stresses on the vascular wall. Furthermore, the shear stress distribution followed the same behavior with oscillating magnitude over the entire flow cycle. Thus, distal IH observed in end-to-side artery-graft models may be caused by the fluctuations in WSS's along the wall.