Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

The effect of inlet waveforms on computational hemodynamics of patient-specific intracranial aneurysms

Due to the lack of patient-specific inlet flow waveform measurements, most computational fluid dynamics (CFD) simulations of intracranial aneurysms usually employ waveforms that are not patient-specific as inlet boundary conditions for the computational model. The current study examined how this assumption affects the predicted hemodynamics in patient-specific aneurysm geometries. We examined wall shear stress (WSS) and oscillatory shear index (OSI), the two most widely studied hemodynamic quantities that have been shown to predict aneurysm rupture, as well as maximal WSS (MWSS), energy loss (EL) and pressure loss coefficient (PLc). Sixteen pulsatile CFD simulations were carried out on four typical saccular aneurysms using 4 different waveforms and an identical inflow rate as inlet boundary conditions. Our results demonstrated that under the same mean inflow rate, different waveforms produced almost identical WSS distributions and WSS magnitudes, similar OSI distributions but drastically different OSI magnitudes. The OSI magnitude is correlated with the pulsatility index of the waveform. Furthermore, there is a linear relationship between aneurysm-averaged OSI values calculated from one waveform and those calculated from another waveform. In addition, different waveforms produced similar MWSS, EL and PLc in each aneurysm. In conclusion, inlet waveform has minimal effects on WSS, OSI distribution, MWSS, EL and PLc and a strong effect on OSI magnitude, but aneurysm-averaged OSI from different waveforms has a strong linear correlation with each other across different aneurysms, indicating that for the same aneurysm cohort, different waveforms can consistently stratify (rank) OSI of aneurysms.     

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.     

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.     

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.     

Comparison between a generalized Newtonian model and a network-type multiscale model for hemodynamic behavior in the aortic arch: Validation with 4D MRI data for a case study

Blood is a complex fluid in which the presence of the various constituents leads to significant changes in its rheological properties. Thus, an appropriate non-Newtonian model is advisable; and we choose a Modified version of the rheological model of Phan-Thien and Tanner (MPTT). The different parameters of this model, derived from the rheology of polymers, allow characterization of the non-Newtonian nature of blood, taking into account the behavior of red blood cells in plasma. Using the MPTT model that we implemented in the open access software OpenFOAM, numerical simulations have been performed on blood flow in the thoracic aorta for a healthy patient. We started from a patient-specific model which was constructed from medical images. Exiting flow boundary conditions have been developped, based on a 3-element Windkessel model to approximate physiological conditions. The parameters of the Windkessel model were calibrated with in vivo measurements of flow rate and pressure. The influence of the selected viscosity of red blood cells on the flow and wall shear stress (WSS) was investigated. Results obtained from this model were compared to those of the Newtonian model, and to those of a generalized Newtonian model, as well as to in vivo dynamic data from 4D MRI during a cardiac cycle. Upon evaluating the results, the MPTT model shows better agreement with the MRI data during the systolic and diastolic phases than the Newtonian or generalized Newtonian model, which confirms our interest in using a complex viscoelastic model.     

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.     

Inlet conditions for image-based CFD models of the carotid bifurcation: is it reasonable to assume fully developed flow?

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.     

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.     

Analysis of non-Newtonian effects within an aorta-iliac bifurcation region

The geometry of the arteries at or near arterial bifurcation influences the blood flow field, which is an important factor affecting arteriogenesis. The blood can act sometimes as a non-Newtonian fluid. However, many studies have argued that for large and medium arteries, the blood flow can be considered to be Newtonian. In this work a comprehensive investigation of non-Newtonian effects on the blood fluid dynamic behavior in an aorta-iliac bifurcation is presented. The aorta-iliac geometry is reconstructed with references to the values reported in Shah et al. (1978); the 3D geometrical model consists of three filleted cylinders of different diameters. Governing equations with the appropriate boundary conditions are solved with a finite-element code. Different rheological models are used for the blood flow through the lumen and detailed comparisons are presented for the aorta-iliac bifurcation. Results are presented in terms of the velocity profiles in the bifurcation zone and Wall Shear Stress (WSS) for different sides of the bifurcation both for male and female geometries, showing that the Newtonian fluid assumption can be made without any particular loss in terms of accuracy with respect to the other more complex rheological models.     

Realistic non-Newtonian viscosity modelling highlights hemodynamic differences between intracranial aneurysms with and without surface blebs

Most computational fluid dynamic (CFD) simulations of aneurysm hemodynamics assume constant (Newtonian) viscosity, even though blood demonstrates shear-thinning (non-Newtonian) behavior. We sought to evaluate the effect of this simplifying assumption on hemodynamic forces within cerebral aneurysms, especially in regions of low wall shear stress, which are associated with rupture. CFD analysis was performed for both viscosity models using 3D rotational angiography volumes obtained for 26 sidewall aneurysms (12 with blebs, 12 ruptured), and parametric models incorporating blebs at different locations (inflow/outflow zone). Mean and lowest 5% values of time averaged wall shear stress (TAWSS) computed over the dome were compared using Wilcoxon rank-sum test. Newtonian modeling not only resulted in higher aneurysmal TAWSS, specifically in areas of low flow and blebs, but also showed no difference between aneurysms with or without blebs. In contrast, for non-Newtonian analysis, bleb-bearing aneurysms showed significantly lower 5% TAWSS compared to those without (p=0.005), despite no significant difference in mean dome TAWSS (p=0.32). Non-Newtonian modeling also accentuated the differences in dome TAWSS between ruptured and unruptured aneurysms (p<0.001). Parametric models further confirmed that realistic non-Newtonian viscosity resulted in lower bleb TAWSS and higher focal viscosity, especially when located in the outflow zone. The results show that adopting shear-thinning non-Newtonian blood viscosity in CFD simulations of intracranial aneurysms uncovered hemodynamic differences induced by bleb presence on aneurysmal surfaces, and significantly improved discriminant statistics used in risk stratification. These findings underline the possible implications of using a realistic model of blood viscosity in predictive computational hemodynamics.     

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.     

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 common carotid artery inlet length on normal carotid bifurcation hemodynamics

Controversy exists regarding the suitability of fully developed versus measured inlet velocity profiles for image-based computational fluid dynamics (CFD) studies of carotid bifurcation hemodynamics. Here, we attempt to resolve this by investigating the impact of the reconstructed common carotid artery (CCA) inlet length on computed metrics of "disturbed" flow. Twelve normal carotid bifurcation geometries were reconstructed from contrast-enhanced angiograms acquired as part of the Vascular Aging--The Link That Bridges Age to Atherosclerosis study (VALIDATE). The right carotid artery lumen geometry was reconstructed from its brachiocephalic origin to well above the bifurcation, and the CCA was truncated objectively at locations one, three, five, and seven diameters proximal to where it flares into the bifurcation. Relative to the simulations carried out using the full CCA, models truncated at one CCA diameter strongly overestimated the amount of disturbed flow. Substantial improvement was offered by using three CCA diameters, with only minor further improvement using five CCA diameters. With seven CCA diameters, the amounts of disturbed flow agreed unambiguously with those predicted by the corresponding full-length models. Based on these findings, we recommend that image-based CFD models of the carotid bifurcation should incorporate at least three diameters of CCA length if fully developed velocity profiles are to be imposed at the inlet. The need for imposing measured inlet velocity profiles would seem to be relevant only for those cases where the CCA is severely truncated.     

Tetrahedral vs. polyhedral mesh size evaluation on flow velocity and wall shear stress for cerebral hemodynamic simulation

Haemodynamic factors, in particular wall shear stresses (WSSs) may have significant impact on growth and rupture of cerebral aneurysms. Without a means to measure WSS reliably in vivo, computational fluid dynamic (CFD) simulations are frequently employed to visualise and quantify blood flow from patient-specific computational models. With increasing interest in integrating these CFD simulations into pretreatment planning, a better understanding of the validity of the calculations in respect to computation parameters such as volume element type, mesh size and mesh composition is needed. In this study, CFD results for the two most common aneurysm types (saccular and terminal) are compared for polyhedral- vs. tetrahedral-based meshes and discussed regarding future clinical applications. For this purpose, a set of models were constructed for each aneurysm with spatially varying surface and volume mesh configurations (mesh size range: 5119-258, 481 volume elements). WSS distribution on the model wall and point-based velocity measurements were compared for each configuration model. Our results indicate a benefit of polyhedral meshes in respect to convergence speed and more homogeneous WSS patterns. Computational variations of WSS values and blood velocities are between 0.84 and 6.3% from the most simple mesh (tetrahedral elements only) and the most advanced mesh design investigated (polyhedral mesh with boundary layer).     

Wall shear stress at the initiation site of cerebral aneurysms

Hemodynamics are believed to play an important role in the initiation of cerebral aneurysms. In particular, studies have focused on wall shear stress (WSS), which is a key regulator of vascular biology and pathology. In line with the observation that aneurysms predominantly occur at regions of high WSS, such as bifurcation apices or outer walls of vascular bends, correlations have been found between the aneurysm initiation site and high WSS. The aim of our study was to analyze the WSS field at an aneurysm initiation site that was neither a bifurcation apex nor the outer wall of a vascular bend. Ten cases with aneurysms on the A1 segment of the anterior cerebral artery were analyzed and compared with ten controls. Aneurysms were virtually removed from the vascular models of the cases to mimic the pre-aneurysm geometry. Computational fluid dynamics (CFD) simulations were created to assess the magnitude, gradient, multidirectionality, and pulsatility of the WSS. To aid the inter-subject comparison of hemodynamic variables, we mapped the branch surfaces onto a two-dimensional parametric space. This approach made it possible to view the whole branch at once for qualitative evaluation. It also allowed us to empirically define a patch for quantitative analysis, which was consistent among subjects and encapsulated the aneurysm initiation sites in our dataset. To test the sensitivity of our results, CFD simulations were repeated with a second independent observer virtually removing the aneurysms and with a 20 % higher flow rate at the inlet. We found that branches harboring aneurysms were characterized by high WSS and high WSS gradients. Among all assessed variables, the aneurysm initiation site most consistently coincided with peaks of temporal variation in the WSS magnitude.     

Abdominal aortic hemodynamics in young healthy adults at rest and during lower limb exercise: quantification using image-based computer modeling

Localization of atherosclerotic lesions in the abdominal aorta has been previously correlated to areas of adverse hemodynamic conditions, such as flow recirculation, low mean wall shear stress, and high temporal oscillations in shear. Along with its many systemic benefits, exercise is also proposed to have local benefits in the vasculature via the alteration of these regional flow patterns. In this work, subject-specific models of the human abdominal aorta were constructed from magnetic resonance angiograms of five young, healthy subjects, and computer simulations were performed under resting and exercise (50% increase in resting heart rate) pulsatile flow conditions. Velocity fields and spatial variations in mean wall shear stress (WSS) and oscillatory shear index (OSI) are presented. When averaged over all subjects, WSS increased from 4.8 +/- 0.6 to 31.6 +/- 5.7 dyn/cm2 and OSI decreased from 0.22 +/- 0.03 to 0.03 +/- 0.02 in the infrarenal aorta between rest and exercise. WSS significantly increased, whereas OSI decreased between rest and exercise at the supraceliac, infrarenal, and suprabifurcation levels, and significant differences in WSS were found between anterior and posterior sections. These results support the hypothesis that exercise provides localized benefits to the cardiovascular system through acute mechanical stimuli that trigger longer-term biological processes leading to protection against the development or progression of atherosclerosis.