The effects of non-Newtonian blood modeling and pulsatility on hemodynamics in the food and drug administration’s benchmark nozzle model

BACKGROUND: Computational fluid dynamics (CFD) is an important tool for predicting cardiovascular device performance. The FDA developed a benchmark nozzle model in which experimental and CFD data were compared, however, the studies were limited by steady flows and Newtonian models. OBJECTIVE: Newtonian and non-Newtonian blood models will be compared under steady and pulsatile flows to evaluate their influence on hemodynamics in the FDA nozzle. METHODS: CFD simulations were validated against the FDA data for steady flow with a Newtonian model. Further simulations were performed using Newtonian and non-Newtonian models under both steady and pulsatile flows. RESULTS: CFD results were within the experimental standard deviations at nearly all locations and Reynolds numbers. The model differences were most evident at Re = 500, in the recirculation regions, and during diastole. The non-Newtonian model predicted blunter upstream velocity profiles, higher velocities in the throat, and differences in the recirculation flow patterns. The non-Newtonian model also predicted a greater pressure drop at Re = 500 with minimal differences observed at higher Reynolds numbers. CONCLUSIONS: An improved modeling framework and validation procedure were used to further investigate hemodynamics in geometries relevant to cardiovascular devices and found that accounting for blood’s non-Newtonian and pulsatile behavior can lead to large differences in predictions in hemodynamic parameters.

Comparison of Non-Newtonian Models of One-Dimensional Hemodynamics

The paper is devoted to the comparison of different one-dimensional models of blood flow. In such models, the non-Newtonian property of blood is considered. It is demonstrated that for the large arteries, the small parameter is observed in the models, and the perturbation method can be used for the analytical solution. In the paper, the simplified nonlinear problem for the semi-infinite vessel with constant properties is solved analytically, and the solutions for different models are compared. The effects of the flattening of the velocity profile and hematocrit value on the deviation from the Newtonian model are investigated.

The Influence of Non-Newtonian Model on Properties of Blood Flow Through a Left Coronary Artery with Presence of Different Double Stenosis

Cardiovascular diseases were the main cause for loosing lives in the last decades due to the restricted blood flow states in the blood vessels areas. Numerical investigations have been conducted as the aim of this work to examine the blood flow, and wall shear stresses adjacent to the mono stenosis up to different degrees involved in the main, side and distal main branches as well as observe the pulsatile flow of blood in the left coronary artery through various percentage of stenosis. Both the Carreau non-Newtonian rheological model and the Newtonian model were utilized to model the blood fluid and wall shear stresses of left coronary artery, in a row, all the calculated data were validated with the previously published papers. It was found that the blood flow inside areas of the artery lie within the range of non-Newtonian rheological effects can be present, verifying the need to treat blood as non-Newtonian fluid; especially, with the case of 90% blockage.

Non-Newtonian Blood Modeling in Intracranial Aneurysm Hemodynamics: Impact On the WSS and OSI Metrics for Ruptured and Unruptured Cases

When simulating blood flow in intracranial aneurysms, the Newtonian model seems to be ubiquitous. However, analyzing the results from the few studies on this subject, the doubt remains on whether it is necessary to use non-Newtonian models in wall shear stress (WSS) simulations of cerebral vascular flows. Another open question related to this topic is whether different rheology models would influence the flow parameters for ruptured and unruptured cases, especially because ruptured aneurysms normally have morphological features that could trigger non-Newtonian phenomena in the blood flow due to low shear rates. The objective of this study is to investigate such flows. By using Computational Fluid Dynamics (CFD) in an open-source framework, we simulated an equal number of ruptured and unruptured patient-specific aneurysms to assess the influence of the blood modeling on the main hemodynamic variables associated with aneurysm formation, growth, and rupture. Results for wall shear stress and oscillatory shear index and their metrics were obtained using Casson and Carreau-Yasuda non-Newtonian models and were compared with those obtained using the Newtonian model. We found that the wall shear stress at peak systole is overestimated by more than 50% by using the non-Newtonian models, but its metrics based on time and surface averaged values remain unaffected. On the other hand, the surface-averaged oscillatory shear index (OSI) is underestimated by more than 40% by the non-Newtonian models. In addition, all differences were consistent among all aneurysms cases irrespective of their rupture status.

Newtonian and Non-Newtonian Blood Rheology Inside a Model of Stenosis

In order to imitate the atherosclerosis artery disease and determine the key issues, Computational Fluid Dynamics (CFD) is able to play a leading rule in the analysis of flow physics within the clogged arteries, in particular the stenosis artery. The problem of blood flow blockage through the blood vessel has been investigated numerically within a stenosis artery. In this work, a CFD technique was employed to solve the three-dimensional, steady, laminar and non-Newtonian Carreau model blood flow through a stenosis artery using Star-CCM+ software. The shape of stenosis that has been selected is a trapezoidal with two cases (70% and 90% blockage). Shear rate, streamlines, vorticity and importance factor are examined to assess the influence of non-Newtonian model through the test section, the Carreau model was compared with Newtonian model. The clinical significance of the shear rate is reported for the examined cases, observing that the levels of non-Newtonian model are predicted to be higher in the 90% blockage than that observed within the 70%; the same finding as related with the axial velocity and vorticity. The levels of re-circulation areas and vorticity are showed to be enlarged in the Carreau model compared with the case of Newtonian.

A Modular, Non-Newtonian, Model, Library Framework (DebrisLib) for Post-Wildfire Flood Risk Management

Wildfires increase flow and sediment load through removal of vegetation, alteration of soils, decreasing infiltration, and production of ash commonly generating a wide variety of geophysical flows (i.e., hyperconcentrated flows, mudflows, debris flows, etc.). Numerical modellers have developed a variety of Non-Newtonian algorithms to simulate each of these processes, and therefore, it can be difficult to understand the assumptions and limitations in any given model or replicate work. This diversity in the processes and approach to non-Newtonian simulations makes a modular computation library approach advantageous. A computational library consolidates the algorithms for each process and discriminates between these processes and algorithms with quantitative non-dimensional thresholds. This work presents a flexible numerical library framework (DebrisLib) to simulate large-scale, post-wildfire, non-Newtonian geophysical flows using both kinematic wave and shallow-water models. DebrisLib is derived from a variety of non-Newtonian closure approaches that predict a range of non-Newtonian flow conditions. It is a modular code designed to operate with any Newtonian, shallow-water parent code architecture. This paper presents the non-Newtonian model framework and demonstrates its effectiveness by calling it from two very different modelling frameworks developed by the U.S. Army Corp of Engineers (USACE), specifically, within the one-dimensional and two-dimensional Hydrologic Engineering Centre River Analysis System (HEC-RAS) and two-dimensional Adaptive Hydraulics (AdH) numerical models. The development and linkage-architecture were verified and validated using two non-Newtonian flume experiments selected to represent a range of non-Newtonian flow conditions (i.e., hyperconcentrated flow, mudflow, debris flow) commonly associated with post-wildfire flooding.

P2428Comparison of Newtonian and non-Newtonian rheology in calculation of endothelial shear stress

Background Although blood is a non-Newtonian fluid, most clinical computational fluid dynamic (CFD) studies assume blood to be a Newtonian fluid with constant viscosity. At higher blood flow rates in larger arteries, the two models should present similar results, and the Newtonian assumption can be considered acceptable. However, whether the Newtonian assumption is valid in patient-specific coronary arteries under pulsatile flow has not been evaluated. Purpose To compare CFD results using Newtonian and non-Newtonian models of blood in order to determine whether the Newtonian assumption can be considered valid in patient-specific coronary arteries. Methods Coronary arteries of 16 patients were reconstructed from fusion of angiography and intracoronary optical coherence tomography imaging. Pulsatile CFD simulations using Newtonian and non-Newtonian models were performed to calculate endothelial shear stress (ESS). The absolute and percent difference in time-averaged and instantaneous ESS values (calculated as non-Newtonian minus Newtonian) were compared on a point-to-point basis. The percent area of the vessel exposed to proatherogenic ESS values (considered <1 Pa) in each model was also calculated. Results The Newtonian and non-Newtonian models produce similar qualitative distributions of ESS. However, quantitative comparison shows that compared to the Newtonian results, the non-Newtonian model estimates significantly higher time-averaged ESS (2.04±0.63Pa versus 1.59±0.54Pa, 95% CI 0.39–0.49, p<0.001) throughout the cardiac cycle. This results in significantly greater estimate of area exposed to ESS <1Pa in the Newtonian model (50.43±14.16% versus 37.20±13.57%, 95% CI 11.28–15.18, p<0.001). Instantaneous ESS plotted through the cardiac cycle indicates the greatest divergence in ESS values occurs at the transition between end-systole and early diastole, at approximately 0.35 seconds (FIGURE). Conclusions Despite similar qualitative ESS distributions, Newtonian and non-Newtonian simulations provide significantly different quantitative ESS values. This suggests that in patient-specific simulations of coronary blood flow, the non-Newtonian model may increase accuracy of ESS measurements. We hypothesize that using a non-Newtonian model may improve the diagnostic accuracy of abnormal ESS to predict clinically significant progression of atherosclerosis, however further study is necessary.