Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were 'treated' with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were ‘treated’ with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Flow dynamics in Models of Intracranial Terminal Aneurysms

Flow dynamics play an important role in the pathogenesis and treatment of intracranial aneurysms. The evaluation of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils, and the temporal and spatial variations of wall shear stress in the aneurysm are correlated with its growth and rupture. This numerical investigation describes the hemodynamic in two models of terminal aneurysm of the basilar artery. Aneurysm models with a aspect ratio of 1.0 and 1.67 were studied. Each model was subject to physiological representative waveform of inflow for a mean Reynolds number of 560. The effects of symmetric and asymmetric outflow conditions in the branches were studied.     

Modeling pulsatile flow in aortic aneurysms: effect of non-Newtonian properties of blood

Pulsatile flow in an axisymmetric rigid-walled model of an abdominal aorta aneurysm was analyzed numerically for various aneurysm dilations using physiologically realistic resting waveform at time-averaged Reynolds number of 300 and peak Reynolds number of 1607. Discretization of the governing equations was achieved using a finite element scheme based on the Galerkin method of weighted residuals. Comparisons with previously published work on the basis of special cases were performed and found to be in excellent agreement. Our findings indicate that the velocity fields are significantly affected by non-Newtonian properties in pathologically altered configurations. Non-Newtonian fluid shear stress is found to be greater than Newtonian fluid shear stress during peak systole. Further, the maximum shear stress is found to occur near the distal end of AAA during peak systole. The impact of non-Newtonian blood flow characteristics on pressure compared to Newtonian model is found insignificant under resting conditions. Viscous and inertial forces associated with blood flow are responsible for the changes in the wall that result in thrombus deposition and dilation while rupture of AAA is more likely determined by much larger mechanical stresses imposed by pulsatile pressure on the wall of AAA.     

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.     

Non-Newtonian blood flow in human right coronary arteries: transient simulations

This study looks at pulsatile blood flow through four different right coronary arteries, which have been reconstructed from biplane angiograms. A non-Newtonian blood model (the Generalised Power Law), as well as the usual Newtonian model of blood viscosity, is used to study the wall shear stress in each of these arteries over the entire cardiac cycle. The difference between Newtonian and non-Newtonian blood models is also studied over the whole cardiac cycle using the recently generalised global non-Newtonian importance factor. In addition, the flow is studied by considering paths of massless particles introduced into the flow field. The study shows that, when studying the wall shear stress distribution for transient blood flow in arteries, the use of a Newtonian blood model is a reasonably good approximation. However, to study the flow within the artery in greater detail, a non-Newtonian model is more appropriate.     

Exploring high frequency temporal fluctuations in the terminal aneurysm of the basilar bifurcation

Cerebral aneurysms are a common cause of death and disability. Of all the cardiovascular diseases, aneurysms are perhaps the most strongly linked with the local fluid mechanic environment. Aside from early in vivo clinical work that hinted at the possibility of high-frequency intra-aneurysmal velocity oscillations, flow in cerebral aneurysms is most often assumed to be laminar. This work investigates, through the use of numerical simulations, the potential for disturbed flow to exist in the terminal aneurysm of the basilar bifurcation. The nature of the disturbed flow is explored using a series of four idealized basilar tip models, and the results supported by four patient specific terminal basilar tip aneurysms. All four idealized models demonstrated instability in the inflow jet through high frequency fluctuations in the velocity and the pressure at approximately 120?Hz. The instability arises through a breakdown of the inflow jet, which begins to oscillate upon entering the aneurysm. The wall shear stress undergoes similar high-frequency oscillations in both magnitude and direction. The neck and dome regions of the aneurysm present 180 deg changes in the direction of the wall shear stress, due to the formation of small recirculation zones near the shear layer of the jet (at the frequency of the inflow jet oscillation) and the oscillation of the impingement zone on the dome of the aneurysm, respectively. Similar results were observed in the patient-specific models, which showed high frequency fluctuations at approximately 112 Hz in two of the four models and oscillations in the magnitude and direction of the wall shear stress. These results demonstrate that there is potential for disturbed laminar unsteady flow in the terminal aneurysm of the basilar bifurcation. The instabilities appear similar to the first instability mode of a free round jet.     

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patient's average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.     

Numerical study of the impact of non-Newtonian blood behavior on flow over a two-dimensional backward facing step

Endothelial cell (EC) responsiveness to shear stress is essential for vasoregulation and plays a role in atherogenesis. Although blood is a non-Newtonian fluid, EC flow studies in vitro are typically performed using Newtonian fluids. The goal of the present study was to determine the impact of non-Newtonian behavior on the flow field within a model flow chamber capable of producing flow disturbance and whose dimensions permit Reynolds and Womersley numbers comparable to those present in vivo. We performed two-dimensional computational fluid dynamic simulations of steady and pulsatile laminar flow of Newtonian and non-Newtonian fluids over a backward facing step. In the non-Newtonian simulations, the fluid was modeled as a shear-thinning Carreau fluid. Steady flow results demonstrate that for Re in the range 50-400, the flow recirculation zone downstream of the step is 22-63% larger for the Newtonian fluid than for the non-Newtonian fluid, while spatial gradients of shear stress are larger for the non-Newtonian fluid. In pulsatile flow, the temporal gradients of shear stress within the flow recirculation zone are significantly larger for the Newtonian fluid than for the non-Newtonian fluid. These findings raise the possibility that in regions of flow disturbance, EC mechanotransduction pathways stimulated by Newtonian and non-Newtonian fluids may be different.     

Flush mounted hot film anemometer measurement of wall shear stress distal to a tri-leaflet valve for Newtonian and non-Newtonian blood analog fluids

Wall shear stress has been measured by flush-mounted hot film anemometry distal to an Ionescu-Shiley tri-leaflet valve under pulsatile flow conditions. Both Newtonian (aqueous glycerol) and non-Newtonian (aqueous polyacrylamide) blood analog fluids were investigated. Significant differences in the axial distribution of wall shear stress between the two fluids are apparent in flows having nearly identical Reynolds numbers. The Newtonian fluid exhibits a (peak) wall shear rate which is maximized near the valve seat (30 mm) and then decays to a fully developed flow value (by 106 mm). In contrast, the shear rate of the non-Newtonian fluid at 30 mm is less than half that of the Newtonian fluid and at 106 mm is more than twice that of the Newtonian fluid. It is suggested that non-Newtonian rheology influences valve flow patterns either through alterations in valve opening associated with low shear separation zones behind valve leaflets, or because of variations in the rate of jet spreading. More detailed studies are required to clarify the mechanisms. The Newtonian wall shear stresses for this valve are low. The highest value observed anywhere in the aortic chamber was 2.85 N/m2 at a peak Reynolds number of 3694.     

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.     

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.     

Structure of pulsatile flow in a model of elastic cerebral aneurysm

This study examines the effect of aneurysmal wall elasticity on the structure of flow within an elastic aneurysm during pulsatile flow. We visualized flow structure in a model of an elastic saccular aneurysm located at the bifurcation of the anterior cerebral artery and extending to the anterior communicating artery, and measured changes in the diameter of the aneurysm wall during pulsatile flow using particle imaging velocimetry (PIV). We similarly measured these features during steady flow by PIV and found that dilation of the aneurysmal wall absorbed the dynamic energy within the aneurysm. Accordingly, aneurysm wall elasticity functions as a biocompatible reaction that relieves wall shear stress acting on the vascular wall during pulsatile flow, and should thus inhibit the development and rupture of an aneurysm.     

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.     

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.     

Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm

In experiments turbulence has previously been shown to occur in intracranial aneurysms. The effects of turbulence induced oscillatory wall stresses could be of great importance in understanding aneurysm rupture. To investigate the effects of turbulence on blood flow in an intracranial aneurysm, we performed a high resolution computational fluid dynamics (CFD) simulation in a patient specific middle cerebral artery (MCA) aneurysm using a realistic, pulsatile inflow velocity. The flow showed transition to turbulence just after peak systole, before relaminarization occurred during diastole. The turbulent structures greatly affected both the frequency of change of wall shear stress (WSS) direction and WSS magnitude, which reached a maximum value of 41.5Pa. The recorded frequencies were predominantly in the range of 1-500Hz. The current study confirms, through properly resolved CFD simulations that turbulence can occur in intracranial aneurysms.     

Can aspect ratio be used to categorize intra-aneurysmal hemodynamics?-A study of elastase induced aneurysms in rabbit

Clinical studies suggest that aneurysm aspect ratio (AR) is an important indicator of rupture likelihood. The importance of AR is hypothesized to arise from its influence on intra-aneurysmal hemodynamics. It has been conjectured that slower flow in high AR sacs leads to a cascade of biological activities that weaken the aneurysm wall (Ujiie et al.,1999). However, the connection between AR, hemodynamics and wall weakening has never been proven. Animal models of saccular aneurysms provide a venue for evaluating this conjecture. The focus of this work was to evaluate whether a commonly used elastase induced aneurysm model in rabbits is suitable for a study of this kind from a hemodynamic perspective. In particular, to assess whether hemodynamic factors in low and high AR sacs are statistically different. To achieve this objective, saccular aneurysms were created in 51 rabbits and pulsatile computational fluid dynamics (CFD) studies were performed using rabbit specific inflows. Distinct hemodynamics were found in the low AR (AR<1.8, n=25), and high AR (AR>2.2, n=18) models. A single, stable recirculation zone was present in all low AR aneurysms, whereas a second, transient recirculation zone was also found in the superior aspect of the aneurysm dome for all high AR cases. Aneurysms with AR between 1.8 and 2.2 displayed transitional flow patterns. Differences in values and distributions of hemodynamic parameters were found between low and high AR cases including time averaged wall shear stress, oscillatory shear index, relative residence time and non-dimensional inflow rate. This work lays the foundation for future studies of the dependence of growth and remodeling on AR in the rabbit model and provides a motivation for further studies of the coupling between AR and hemodynamics in human aneurysms.     

Non-Newtonian effects of blood flow on hemodynamics in distal vascular graft anastomoses

Non-Newtonian fluid flow in a stenosed coronary bypass is investigated numerically using the Carreau-Yasuda model for the shear thinning behavior of the blood. End-to-side coronary bypass anastomosis is considered in a simplified model geometry where the host coronary artery has a 75% severity stenosis. Different locations of the bypass graft to the stenosis and different flow rates in the graft and in the host artery are studied. Particular attention is given to the non-Newtonian effect of the blood on the primary and secondary flow patterns in the host coronary artery and the wall shear stress (WSS) distribution there. Interaction between the jet flow from the stenosed artery and the flow from the graft is simulated by solving the three-dimensional Navier-Stokes equation coupled with the non-Newtonian constitutive model. Results for the non-Newtonian flow, the Newtonian flow and the rescaled Newtonian flow are presented. Significant differences in axial velocity profiles, secondary flow streamlines and WSS between the non-Newtonian and Newtonian fluid flows are revealed. However, reasonable agreement between the non-Newtonian and the rescaled Newtonian flows is found. Results from this study support the view that the residual flow in a partially occluded coronary artery interacts with flow in the bypass graft and may have significant hemodynamic effects in the host vessel downstream of the graft. Non-Newtonian property of the blood alters the flow pattern and WSS distribution and is an important factor to be considered in simulating hemodynamic effects of blood flow in arterial bypass grafts.     

Numerical simulation of saccular aneurysm hemodynamics: influence of morphology on rupture risk

The governing equations for pulsatile fluid flow were solved in their finite volume formulation in order to simulate blood flow in a variety of three-dimensional aneurysm geometries. The influence of geometric factors on flow patterns and fluid mechanical forces was studied with the goal of identifying the risk of aneurysm rupture. Aneurysm morphology was characterized by quantitative shape indices reflecting the three dimensionality of the vasculature derived from clinical studies. Recirculation zones and secondary flows were observed in aneurysms and arteries. Regions of extreme and alternating shear stress were observed and identified as sites for potential aneurysm rupture. The ellipticity of an aneurysm was observed to be strongly correlated with wall shear stress at the aneurysm fundus, while its non-sphericity, volume, and degree of undulation were more weakly correlated.     

Comparative study of Newtonian and non-Newtonian simulations of drug transport in a model drug-eluting stent

To elucidate the difference between Newtonian and shear thinning non-Newtonian assumptions of blood in the analysis of DES drug delivery, we numerically simulated the local flow pattern and the concentration distribution of the drug at the lumen-tissue interface for a structurally simplified DES deployed in a curved segment of an artery under pulsatile blood flow conditions. The numerical results showed that when compared with the Newtonian model, the Carreau (shear thinning) model could lead to some differences in the luminal surface drug concentration in certain areas along the outer wall of the curved vessel. In most areas of the vessel, however, there were no significant differences between the 2 models. Particularly, no significant difference between the two models was found in terms of the area-averaged luminal surface drug concentration. Therefore, we believe that the shear thinning property of blood may play little roles in DES drug delivery. Nevertheless, before we draw the conclusion that Newtonian assumption of blood can be used to replace its non-Newtonian one for the numerical simulation of drug transport in the DES implanted coronary artery, other more complex mechanical properties of blood such as its thixotropic behavior should be tested.