Numerical analysis of non-Newtonian blood flow and wall shear stress in realistic single, double and triple aorto-coronary bypasses

A hemodynamic analysis of energy loss in abdominal aortic aneurysm using three-dimension idealized model

Objective: The aim of this study is to perform specific hemodynamic simulations of idealized abdominal aortic aneurysm (AAA) models with different diameters, curvatures and eccentricities and evaluate the risk of thrombosis and aneurysm rupture. Methods: Nine idealized AAA models with different diameters (3 cm or 5 cm), curvatures (0° or 30°) and eccentricities (centered on or tangent to the aorta), as well as a normal model, were constructed using commercial software (Solidworks; Dassault Systemes S.A, Suresnes, France). Hemodynamic simulations were conducted with the same time-varying volumetric flow rate extracted from the literature and 3-element Windkessel model (3 EWM) boundary conditions were applied at the aortic outlet. Several hemodynamic parameters such as time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), endothelial cell activation potential (ECAP) and energy loss (EL) were obtained to evaluate the risk of thrombosis and aneurysm rupture under different conditions. Results: Simulation results showed that the proportion of low TAWSS region and high OSI region increases with the rising of aneurysm diameter, whereas decreases in the curvature and eccentric models of the corresponding diameters, with the 5 cm normal model having the largest low TAWSS region (68.5%) and high OSI region (40%). Similar to the results of TAWSS and OSI, the high ECAP and high RRT areas were largest in the 5 cm normal model, with the highest wall-averaged value (RRT: 5.18 s, ECAP: 4.36 Pa-1). Differently, the increase of aneurysm diameter, curvature, and eccentricity all lead to the increase of mean flow EL and turbulent EL, such that the highest mean flow EL (0.82 W · 10-3) and turbulent EL (1.72 W · 10-3) were observed in the eccentric 5 cm model with the bending angle of 30°. Conclusion: Collectively, increases in aneurysm diameter, curvature, and eccentricity all raise mean flow EL and turbulent flow EL, which may aggravate the damage and disturbance of flow in aneurysm. In addition, it can be inferred by conventional parameters (TAWSS, OSI, RRT and ECAP) that the increase of aneurysm diameter may raise the risk of thrombosis, whereas the curvature and eccentricity appeared to have a protective effect against thrombosis.

A hemodynamic model of artery bypass graft considering microcirculation function

BACKGROUND

The incidence of arterial stenosis is increasing year by year. In order to better diagnose and treat arterial stenosis, numerical simulation technology has become a popular method.

OBJECTIVE

A novel model is constructed to investigate the influence of microcirculation on the hemodynamics of artery bypass graft.

METHODS

In this paper, a severely narrow artery bypass graft model is considered. The geometric shape includes a narrow artery tube and a bypass graft of the same diameter with a 45° suture angle. The fluid-structure interaction model is considered by finite element numerical calculation, and the flow is simulated with microcirculation as the outlet boundary condition. The changes of blood flow velocity, pressure and wall shear stress are analyzed.

RESULTS

The results show that blood almost entirely flows into the graft tube and there is no recirculation area at the anastomosis.

CONCLUSION

The artery bypass graft model considering microcirculation function could simulate the physiological characteristics of blood flow more reasonably, and it provide helps for clinicians to diagnose and treat arterial stenosis.

Multiphase Flow Hemodynamic Evaluation of Vertebral Artery Stenosis Lesions and Plaque Stability

BACKGROUND

Atherosclerosis is one of the main causes of vertebral artery stenosis, which reduces blood supply to the posterior circulation, resulting in cerebral infarction or death.

OBJECTIVE

To investigate stenosis rates and locations on the development of vertebral artery plaques.

METHODS

Stenosis models with varying degrees and positions of stenosis were established. The stenosis area was comprehensively analyzed using multiphase flow numerical simulation. Wall shear stress (WSS), blood flow velocity, and red blood cell (RBC) volume fraction were calculated.

RESULTS

Blood flow velocity in 30-70% stenosis of each segment tended to increase significantly higher than normal. Downstream of 50% stenosis exhibited turbulent flow; downstream of 70% displayed reflux. Severe stenosis increases the WSS and distribution area. The mixed area of high and low WSS appeared downstream of the stenosis. The RBC volume fraction at the stenosis increased (maximum value: 0.487 at 70% stenosis in the V4), which was 1.08 times the normal volume fraction. Turbulent and backflow regions exhibited complex RBC volume fraction distributions.

CONCLUSION

Flow velocity, WSS, and RBC volume fraction at the stenosis increase with stenosis severity, increasing plaque shedding. Narrow downstream spoiler and reflux areas possess low WSS and high erythrocyte volume fractions, accelerating plaque growth.

Effects of different bypass surgeries on LSA revascularization in patients with left subclavian occlusion

INTRO

Left subclavian artery bypass surgery is mainly carried out for patients with severe left subclavian occlusion. This paper aimed to evaluate the hemodynamic effects of different surgical bypass modes on left subclavian artery revascularization.

METHODS

Three-dimensional models of the aorta were reconstructed from CTA images of a patient with left subclavian artery occlusion, a patient with type B aortic dissection with left subclavian artery coverage during thoracic endovascular aortic repair, and a healthy 74 year-old man, resulting in six modes for each person: healthy left subclavian artery mode, left subclavian artery occlusion mode and four bypass modes. Hemodynamic parameters, including flow field, flow distribution, pressure gradient, and wall shear stress, were calculated using computational fluid dynamics.

RESULTS

After left subclavian artery bypass surgery, distal left subclavian artery blood flow resulting from left common carotid artery to distal left subclavian artery bypass was 100% of that in the healthy mode, while the other modes yielded flows at least 91%. Moreover, reversed flow only completely disappeared with left common carotid artery to distal left subclavian artery bypass, whereas reverse flow was observed in the other three modes in early systole.

CONCLUSION

Left common carotid artery to distal left subclavian artery bypass can effectively reduce reverse blood flow in the left vertebral artery, and it is a feasible, effective, and safe option for left subclavian artery revascularization in patients with left subclavian occlusion.

A Hemodynamic Analysis of the Thrombosis Within Occluded Coronary Arterial Fistulas With Terminal Aneurysms Using a Blood Stasis Model

Objective: The aim of this study is to numerically evaluate thrombosis risk within occluded coronary arterial fistulas (CAF) with terminal aneurysms, and provide guidance in choosing occlusion positions, with clinical observations as reference. Method: Four patients with CAF were studied, with different occlusion positions in actual treatments. Hemodynamics simulations were conducted, with blood residue predicted using the blood stasis model. Three types of models (untreated model, aneurysm-reserved model and aneurysm-removed model) were studeid for each patient. Four metrics, i.e., proportion of high oscillatory shear index (OSI), area of high OSI, old blood volume fraction (OBVF)) and old blood volume (OBV) was obtained to distinguish the thrombosis risk of different treatments (proximal or distal occlusion), comparing with the follow-up CTA. Results: For all the postopertive models, the high OBVF, high OSI(>0.3) and low time-averaged wall shear stress (TAWSS) regions were mainly at the distal fistula, indicating these regions were prone to thrombosis. The regions where blood residue remains are roughly regions of high OSI, corresponding well with clinical observations. In contrast, TAWSS failed to distinguish the difference in thrombosis risk. Absolute values (area of high OSI, OBV) can better reflect the degree of thrombosis risk between treatment types compared with percentage values (proportion of high OSI, OBVF). By comparing with the actual clinical treatments and observations, the OBV is superior to the area of high OSI in determining treatment type. Conclusion: The OBV, a volumetric parameter for blood stasis, can better account for the CAF thrombosis and reflect the degree of blood stasis compared with OSI or TAWSS, is a more appropriate metric for thrombosis in the fistula. Together with morphological parameters, the OBV could guide clinicians to formulate more appropriate surgical plans, which is of great significance for the preoperative evaluation and treatment prognosis of CAF patients.

Effect of non-Newtonian fluid rheology on an arterial bypass graft: A numerical investigation guided by constructal design

In post-operative scenarios of arterial graft surgeries to bypass coronary artery stenosis, fluid dynamics plays a crucial role. Problems such as intimal hyperplasia have been related to fluid dynamics and wall shear stresses near the graft junction. This study focused on the question of the use of Newtonian and non-Newtonian models to represent blood in this type of problem in order to capture important flow features, as well as an analysis of the performance of geometry from the view of Constructive Theory. The objective of this study was to investigate the effects rheology on the steady-state flow and on the performance of a system consisting of an idealized version of a partially obstructed coronary artery and bypass graft. The Constructal Design Method was employed with two degrees of freedom: the ratio between bypass and artery diameters and the junction angle at the bypass inlet. The flow problem was solved numerically using the Finite Volume Method with blood modeled employing the Carreau equation for viscosity. The Computational Fluid Dynamics model associated with the Sparse Grid method generated eighteen response surfaces, each representing a severe stenosis degree of 75% for specific combinations of rheological parameters, dimensionless viscosity ratio, Carreau number and flow index at two distinct Reynolds numbers of 150 and 250. There was a considerable dependence of the pressure drop on rheological parameters. For the two Reynolds numbers studied, the Newtonian case presented the lowest value of the dimensionless pressure drop, suggesting that the choice of applying Newtonian blood may underestimate the value of pressure drop in the system by about 12.4% (Re =150) and 7.8% (Re = 250). Even so, results demonstrated that non-Newtonian rheological parameters did not influence either the shape of the response surfaces or the optimum bypass geometry, which consisted of a diameter ratio of 1 and junction angle of 30°. However, the viscosity ratio and the flow index had the greatest impact on pressure drop, recirculation zones and wall shear stress. Rheological parameters also affected the recirculation zones downstream of stenosis, where intimal hyperplasia is more prevalent. Newtonian and most non-Newtonian results had similar wall shear stresses, except for the non-Newtonian case with high viscosity ratio. In the view of Constructal Design, the geometry of best performance was independent of the rheological model. However, rheology played an important role on pressure drop and flow dynamics, allowing the prediction of recirculation zones that were not captured by a Newtonian model.

On the relevance of boundary conditions and viscosity models in blood flow simulations in patient-specific aorto-coronary bypass models

Physiologically realistic results are the aim of every blood flow simulation. This is not different in aorto-coronary bypasses where the properties of the coronary circulation may significantly affect the relevance of the performed simulations. By considering three patient-specific bypass geometries, the present article focuses on two aspects of the coronary blood flow - its phasic flow pattern and its behaviour affected by blood rheology. For the phasic flow property, a multiscale modelling approach is chosen as a means to assess the ability of five different types of coronary boundary conditions (mean arterial pressure, Windkessel model and three lumped parameter models) to attain realistic coronary haemodynamics. From the analysed variants of boundary conditions, the best option in terms of physiological characteristics and its potential for use in patient-based simulations, is utilised to account for the effect of shear-dependent viscosity on the resulting haemodynamics and wall shear stress stimulation. Aside from the Newtonian model, the blood rheology is approximated by two non-Newtonian models in order to determine whether the choice of a viscosity model is important in simulations involving coronary circulation. A comprehensive analysis of obtained results demonstrated notable superiority of all lumped parameter models, especially in comparison to the constant outlet pressure, which regardless of bypass type gave overestimated and physiologically misleading results. In terms of rheology, it was noted that blood in undamaged coronary arteries behaves as a Newtonian fluid, whereas in vessels with atypical lumen geometry, such as that of anastomosis or stenosis, its shear-thinning behaviour should not be ignored.

Numerical study of hemodynamics in a complete coronary bypass with venous and arterial grafts and different degrees of stenosis

Cardiovascular diseases are among the leading causes of death in the world. The coronary blockage is one of most common types of these diseases that in the majority of cases has been treated by bypass surgery. In the bypass surgery, a graft is implemented to alter the blocked coronary and allow the blood supply process. The hemodynamic characteristics of the bypass strongly depend on the geometry and mechanical properties of the graft. In the present study, the fluid-structure interaction (FSI) analysis is conducted to investigate the bypass performance for a thoracic artery as well as a saphenous vein graft. Blood flow introduces a pressure on the walls of the graft which behaves as a hyperelastic material. A complete coronary bypass with stenosis degrees of 70% and 100% is modeled. To consider the nonlinear stress-strain behavior of the grafts, a five parameter Mooney-Rivlin hyperplastic model is implemented for the structural analysis and blood is assumed to behave as a Newtonian fluid. The simulations are performed for a structured grid to solve the governing equations using finite element method (FEM). The results show that wall shear stress (WSS) for saphenous vein is larger than that of thoracic artery while the total deformation of the thoracic artery is larger compared to the saphenous vein. Also, for the venous grafts or lower stenosis degree, the oscillatory shear index (OSI) is higher at both left and right anastomoses meaning that venous grafts as well as lower degree of stenosis are more critical in terms of restenosis.

Computational fluid dynamic study of multiple sequential coronary artery bypass anastomoses in a native coronary stenosis model

BACKGROUND

The objective of this study was to evaluate the hemodynamic characteristics of multiple sequential coronary artery bypass grafting using a computational fluid dynamics study.

METHODS

First anastomosis was configured into parallel and diamond anastomoses, and the second anastomosis was set as end-side anastomosis. The anastomosis incision lengths were fixed at 2 mm. Various combinations of the degree of first and second stenoses were studied. The diameter of both the native and graft vessels was set at 2 mm. The inlet boundary condition was set by a sample of the transient time flow measurement, which was measured intraoperatively.

RESULTS

Both swirl and stagnation were observed at the outlets of the stenosis and the anastomosis sites. When the severity of the second stenosis was larger than that of the first, the flow at the outlet of the second stenosis was more unstable. Higher wall shear stress and larger oscillatory shear index regions were observed when the severe stenosis was bypassed by the first anastomosis, especially with diamond anastomoses. Less energy loss and higher energy efficiency were present when the vessel with more severe stenosis was bypassed as the second anastomosis. Energy loss was lower and energy efficiency was higher with parallel anastomosis than diamond anastomosis when the severity of the two stenoses was the same.

CONCLUSIONS

It is ideal to bypass the less severe stenosis vessel first with a parallel anastomosis method when employing multiple sequential bypass grafting. This improves hemodynamic stability and energy efficiency, according to a computational fluid dynamics model.

Computational fluid dynamic study of different incision length of coronary artery bypass grafting in a native coronary stenosis model

Background

The objective of this study was to evaluate hemodynamic patterns in end-side coronary artery bypass grafting with different anastomosis length by computational fluid dynamic study in the native coronary stenosis model.

Methods

The fluid dynamic computations were carried out using ANSYS CFX. Incision length was set to be 2, 4, 6, 8, 10 mm. The angle between the two blood vessels corresponded to the length of the incision. Native vessels were set to be 90% stenosis. The radius of both native and graft vessels was set to be 2 mm. The inlet boundary condition was set by the sample of the transient time flow which was measured intraoperatively.

Results

The energy efficiency was higher and energy loss was lower when the anastomosis length was longer until 8 mm. However, energy efficiency was lowest and energy loss was highest in the 10-mm model. In the 10-mm incision model, the streamline showed the scanty bypass flow in the bottom. Vortex showed that only 10-mm model showed the vortex just distal to the stenosis in the native inlet, and more vortex in native outlet than other length models. The oscillatory shear index (OSI) was higher in the outlet top in all models. And only 10-mm model showed high oscillatory index just distal to the stenosis.

Conclusions

In the end-side anastomosis, an anastomosis length of 8 mm was the ideal length with less flow complexity, low OSI, and less energy loss and high energy efficiency in the native 90% stenosis model.

Investigating the impact of non-Newtonian blood models within a heart pump

A detailed computational fluid dynamics (CFD) study of transient, turbulent blood flow through a positive displacement left ventricular assist device is performed. Two common models for non-Newtonian blood flow are compared to the Newtonian model to investigate their impact on predicted levels of shear rate and wall shear stress. Given that both parameters are directly relevant to the evaluation of risk from thrombus and haemolysis, there is a need to assess the sensitivity to modelling non-Newtonian flow effects within a pulsatile turbulent flow, in order to identify levels of uncertainly in CFD. To capture the effects of turbulence, the elliptic blending Reynolds stress model is used in the present study, on account of superior performance of second moment closure schemes previously identified by the present authors. The CFD configuration includes two cyclically rotating valves and a moving pusher plate to periodically vary the chamber volume. An overset mesh algorithm is used for each instance of mesh motion, and a zero gap technique was employed to ensure full valve closure. The left ventricular assist device was operated at a pumping rate of 86 BPM (beats per minute) and a systolic duration of 40% of the pumping cycle, in line with existing experimental data to which comparisons are made. The sensitivity of the variable viscosity models is investigated in terms of mean flow field, levels of turbulence and global shear rate, and a non-dimensional index is used to directly evaluate the impact of non-Newtonian effects. The clinical relevance of the results is reported along with a discussion of modelling uncertainties, observing that the turbulent kinetic energy is generally predicted to be higher in non-Newtonian flow than that observed in Newtonian flow. Copyright © 2016 John Wiley & Sons, Ltd.

The numerical analysis of non-Newtonian blood flow in human patient-specific left ventricle

Recently, various non-invasive tools such as the magnetic resonance image (MRI), ultrasound imaging (USI), computed tomography (CT), and the computational fluid dynamics (CFD) have been widely utilized to enhance our current understanding of the physiological parameters that affect the initiation and the progression of the cardiovascular diseases (CVDs) associated with heart failure (HF). In particular, the hemodynamics of left ventricle (LV) has attracted the attention of the researchers due to its significant role in the heart functionality. In this study, CFD owing its capability of predicting detailed flow field was adopted to model the blood flow in images-based patient-specific LV over cardiac cycle. In most published studies, the blood is modeled as Newtonian that is not entirely accurate as the blood viscosity varies with the shear rate in non-linear manner. In this paper, we studied the effect of Newtonian assumption on the degree of accuracy of intraventricular hemodynamics. In doing so, various non-Newtonian models and Newtonian model are used in the analysis of the intraventricular flow and the viscosity of the blood. Initially, we used the cardiac MRI images to reconstruct the time-resolved geometry of the patient-specific LV. After the unstructured mesh generation, the simulations were conducted in the CFD commercial solver FLUENT to analyze the intraventricular hemodynamic parameters. The findings indicate that the Newtonian assumption cannot adequately simulate the flow dynamic within the LV over the cardiac cycle, which can be attributed to the pulsatile and recirculation nature of the flow and the low blood shear rate.