On the Relative Importance of Rheology for Image-Based CFD Models of the Carotid Bifurcation

2006 ◽  
Vol 129 (2) ◽  
pp. 273-278 ◽  
Author(s):  
Sang-Wook Lee ◽  
David A. Steinman

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.0dyn∕cm2 or 8%, nearly seven times less than that due to geometric uncertainty (6.7dyn∕cm2 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.

2006 ◽  
Vol 128 (3) ◽  
pp. 371-379 ◽  
Author(s):  
Keri R. Moyle ◽  
Luca Antiga ◽  
David A. Steinman

Background: Computational fluid dynamics tools are useful for their ability to model patient specific data relevant to the genesis and progression of atherosclerosis, but unavailable to measurement tools. The sensitivity of the physiologically relevant parameters of wall shear stress (WSS) and the oscillatory shear index (OSI) to secondary flow in the inlet velocity profiles was investigated in three realistic models of the carotid bifurcation. Method of Approach: Secondary flow profiles were generated using sufficiently long entrance lengths, to which curvature and helical pitch were added. The differences observed were contextualized with respect to effect of the uncertainty of the models’ geometry on the same parameters. Results: The effects of secondary velocities in the inlet profile on WSS and OSI break down within a few diameters of the inlet. Overall, the effect of secondary inlet flow on these models was on average more than 3.5 times smaller than the effect of geometric variability, with 13% and 48% WSS variability induced by inlet secondary flow and geometric differences, respectively. Conclusions: The degree of variation is demonstrated to be within the range of the other computational assumptions, and we conclude that given a sufficient entrance length of realistic geometry, simplification to fully developed axial (i.e., Womersley) flow may be made without penalty. Thus, given a choice between measuring three components of inlet velocity or a greater geometric extent, we recommend effort be given to more accurate and detailed geometric reconstructions, as being of primary influence on physiologically significant indicators.


2021 ◽  
Author(s):  
Alejandro Gonzalo ◽  
Manuel Garcia-Villalba ◽  
Lorenzo Rossini ◽  
Eduardo Duran ◽  
David Vigneault ◽  
...  

Atrial fibrillation (AF) is the most common arrhythmia, affecting ~35M people worldwide. The irregular beating of the left atrial (LA) caused by AF impacts the LA hemodynamics increasing the risk of thrombosis and ischemic stroke. Most LA thrombi appear in its appendage (LAA), a narrow sac of varied morphology where blood is prone to stagnate. In the LAA, the combination of slow blood flow and low shear rates (<100 [1/s]) promotes the formation of red blood cell aggregations called rouleaux. Blood experiences a non-Newtonian behavior when rouleaux formed that has not been considered in previous CFD analysis of the LA. We model the anatomy and motion of the LA from 4D-CT images and solve the blood flow inside the LA geometry with our CFD in-house code, which models Non-Newtonian rheology with the shear-hematocrit-dependent Carreau-Yasuda equation. We cover a wide range of non-Newtonian effects considering a small and a large hematocrit, including an additional constitutive relation to account for themrouleaux formation time, and we compare our results with Newtonian simulations. Blood rheology influence in LAA hemostasis is studied in 6 patient-specific anatomies. Two subjects had an LAA thrombus (digitally removed before running the simulations), another had a history of TIAs, and the remaining three had normal atrial function. In our simulations, the shear rate remains below 50 [1/s] in the LAA for all non-Newtonian models considered. This triggers an increase of viscosity that alters the flow behavior in that site, which exhibits different flow patterns than Newtonian simulations. These hemodynamic changes translate into differences in the LAA hemostasis, calculated with the residence time.


Author(s):  
D. Keith Walters ◽  
Greg W. Burgreen ◽  
Robert L. Hester ◽  
David S. Thompson ◽  
David M. Lavallee ◽  
...  

Computational fluid dynamics (CFD) simulations were performed for unsteady periodic breathing conditions, using large-scale models of the human lung airway. The computational domain included fully coupled representations of the orotracheal region and large conducting zone up to generation four (G4) obtained from patient-specific CT data, and the small conducting zone (to G16) obtained from a stochastically generated airway tree with statistically realistic geometrical characteristics. A reduced-order geometry was used, in which several airway branches in each generation were truncated, and only select flow paths were retained to G16. The inlet and outlet flow boundaries corresponded to the oronasal opening (superior), the inlet/outlet planes in terminal bronchioles (distal), and the unresolved airway boundaries arising from the truncation procedure (intermediate). The cyclic flow was specified according to the predicted ventilation patterns for a healthy adult male at three different activity levels, supplied by the whole-body modeling software HumMod. The CFD simulations were performed using Ansys FLUENT. The mass flow distribution at the distal boundaries was prescribed using a previously documented methodology, in which the percentage of the total flow for each boundary was first determined from a steady-state simulation with an applied flow rate equal to the average during the inhalation phase of the breathing cycle. The distal pressure boundary conditions for the steady-state simulation were set using a stochastic coupling procedure to ensure physiologically realistic flow conditions. The results show that: 1) physiologically realistic flow is obtained in the model, in terms of cyclic mass conservation and approximately uniform pressure distribution in the distal airways; 2) the predicted alveolar pressure is in good agreement with previously documented values; and 3) the use of reduced-order geometry modeling allows accurate and efficient simulation of large-scale breathing lung flow, provided care is taken to use a physiologically realistic geometry and to properly address the unsteady boundary conditions.


2018 ◽  
Vol 34 (6) ◽  
pp. e2972 ◽  
Author(s):  
S.I.S. Pinto ◽  
J.B.L.M. Campos ◽  
E. Azevedo ◽  
C.F. Castro ◽  
L.C. Sousa

Author(s):  
Megan Cummins ◽  
Jenn S. Rossmann

The hemodynamics and fluid mechanical forces in blood vessels have long been implicated in the deposition and growth of atherosclerotic plaque. Detailed information about the hemodynamics in vessels affected by significant plaque deposits can provide insight into the mechanisms and likelihood of plaque weakening and rupture. In the current study, the governing equations are solved in their finite volume formulation in several patient-specific geometries. Recirculation zones, vortex shedding, and secondary flows are captured. The forces on vessel walls are shown to correlate with unstable plaque deposits. The results of these simulations suggest morphological features that may usefully supplement percent stenosis as a predictor of plaque vulnerability.


Blood ◽  
2017 ◽  
Vol 130 (Suppl_1) ◽  
pp. 964-964
Author(s):  
Erdem Kucukal ◽  
Jane A. Little ◽  
Umut A. Gurkan

Abstract The pathophysiology of sickle cell disease (SCD) involves altered biophysical properties of red blood cells (RBCs) and increased cellular adhesion, which can synergistically trigger recurrent and painful vaso-occlusive events in the microcirculatory network. RBC adhesion to the endothelial wall is heterogeneous and may initiate such occlusions by disrupting the local flow thus activating platelets and promoting subsequent cell-cell interactions. Moreover, these episodic events take place within a wide range of dynamically changing shear rates at the microscale. In order to better understand the role of shear rate on this process, we quantified shear-dependent RBC adhesion to endothelial proteins fibronectin (FN) and laminin (LN) utilizing a microfluidic system that can simulate physiologically relevant shear gradients of microcirculatory blood flow at a single flow rate. Whole blood samples were collected from 20 patients (10 males and 10 females) with homozygous SCD (HbSS). Samples were perfused through FN and LN immobilized shear-gradient microchannels (Fig. 1A) in which the shear rate continuously changes along flow direction. Computational simulations characterized the flow dynamics near the adherent RBCs (Fig. 1B). Based on the numerical results, a rectangular "field of interest (FOI)", along which the shear rate dropped approximately three-fold, was chosen for quantification of shear-dependent RBC adhesion. We observed changes in RBC adhesion to LN and FN in the shear gradient flow. Figure 1C and 1D show typical adhesion curves of surface adherent RBCs for an individual SCD sample within the FOI. To assess patient specific shear-dependent adhesion, we defined a parameter, "shear dependent adhesion rate (SDAR)", which is the slope of the adhesion curves based on normalized RBC adhesion numbers. A higher SDAR value was indicative of marked numbers of adherent RBCs that detach at higher shear rates whereas the effect of shear rate on RBC detachment was less for a lower SDAR. We observed an inverse relationship between SDAR and number of persistently adherent RBCs at high shear rates. Shear-dependent RBC adhesion to LN was heterogeneous among SCD patients. Patients with higher WBC counts constituted the low SDAR population with a threshold SDAR value of 60 (Fig. 1E, p=0.005, ANOVA). WBCs from patients with higher SDARs (and fewer persistently adhered cells) were all within the normal range. Patients in the low SDAR group also had significantly elevated absolute neutrophil counts (Fig. 1F, p=0.006, ANOVA), and ferritin levels (Fig. 1G, p=0.007, ANOVA). The mean ferritin level of those with low SDAR was nearly ten times greater than normal (mean= [3272.3 ± 791.9] μg/L vs. [784.5±219.6] μg/L). No white blood cell (WBC) adhesion was observed in the experiments. Here, we report a novel shear dependent adhesion ratio of sickle RBCs utilizing LN and FN functionalized microchannels. The approach presented here enabled us to create a shear gradient throughout the channel which may simulate the physiological flow conditions in the post-capillary venules. We further analyzed shear-dependent RBC adhesion in a patient specific manner and identified patient groups with low and high SDAR. The findings also suggested a link between lower shear dependent sickle RBC adhesion to LN and patient clinical phenotypes including inflammation and iron overload. Acknowledgments: This work was supported by grant #2013126 from the Doris Duke Charitable Foundation, National Heart Lung and Blood Institute R01HL133574, and National Science Foundation CAREER Award 1552782. Figure 1: Shear-dependent sickle RBC adhesion in microscale flow. (A) Macroscopic image of the shear-gradient microchannel with the arrow indicating flow direction. (B) Velocity and shear rate contours on a 2D plane above the bottom surface. The dashed rectangular area indicates the field of interest (FOI) where the experimental data were obtained. (C, D) Typical distribution of adherent deformable and non-deformable RBCs in LN and FN functionalized microchannels with the shear gradient. Dashed lines represent the adhesion curves and the corresponding equations were used to quantify shear dependent adhesion data. Shear-dependent RBC adhesion was lower (nSDAR&lt;60) in patients with elevated white blood cell counts (E), absolute neutrophil counts (F), and serum ferritin levels (G). The dashed rectangles indicate the normal clinical values. Figure 1 Figure 1. Disclosures Little: Hemex Health: Equity Ownership. Gurkan: Hemex Health: Employment, Equity Ownership.


Author(s):  
Mathias Vermeulen ◽  
Cedric Van Holsbeke ◽  
Tom Claessens ◽  
Jan De Backer ◽  
Peter Van Ransbeeck ◽  
...  

An experimental and numerical platform was developed to investigate the fluidodynamics in human airways. A pre operative patient specific geometry was used to create an identical experimental and numerical model. The experimental results obtained from Particle Image Velocimetry (PIV) measurements were compared to Computational Fluid Dynamics (CFD) simulations under stationary and pulsatile flow regimes. Together these results constitute the first step in predicting the clinical outcome of patients after lung surgeries such as Lung Volume Reduction.


Agronomy ◽  
2019 ◽  
Vol 9 (11) ◽  
pp. 769 ◽  
Author(s):  
Jorge Badules ◽  
Mariano Vidal ◽  
Antonio Boné ◽  
Emilio Gil ◽  
F. Javier García-Ramos

A computational fluid dynamics (CFD) model of the fluid velocities generated by the agitation system of an air-assisted sprayer was developed and validated by practical experiments in a laboratory. The model was developed considering different settings of the agitation system: Three water levels in the tank (1000, 2000, and 3000 L); two different numbers of active nozzles (2 or 4); and three working pressures of the agitation circuit (8, 10, or 12 bar). Actual measurements of the fluid velocity into the tank were taken using an acoustic Doppler velocimeter (ADV). CFD simulations made it possible to estimate fluid velocities at 38% of the measuring points with relative errors of less than 30%. Additionally, the CFD models have allowed the correct prediction of the general behavior of the fluid in the tank considering mean velocities depending on the setting parameters of the agitation system (water level in the tank, hydraulic circuit pressure, and number of active nozzles).


2016 ◽  
Vol 11 (4) ◽  
pp. 721-727
Author(s):  
Johnny Ong King Ngu ◽  
Darrien Yau Seng Mah ◽  
Charles Hin Joo Bong

In this paper the flow characteristics of stormwater are analyzed as it travels from a roof gutter down-pipe and the turbulent flow generated on entering an individual lot on-site stormwater detention (OSD) unit beneath a residential carport. Comparison was made between a full-scale model and computational fluid dynamic (CFD) simulations to determine the flow characteristics. These modular tanks with multi-unit chambers can capture the roof run-off from a 15-minute, 10-year return period storm. The results from the physical and CFD models matched well, suggesting that turbulent flow occurs when stormwater is directed to an individual lot stormwater detention tank. However, turbulence in the OSD was concentrated around the inlet, after which the pattern changed from turbulent to laminar flow. This work implies that the use of modular underground storage tanks is practical for managing stormwater from a roof.


Sign in / Sign up

Export Citation Format

Share Document