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.     

Blood flow dynamics in saccular aneurysm models of the basilar artery

Blood flow dynamics under physiologically realistic pulsatile conditions plays an important role in the growth, rupture, and surgical treatment of intracranial aneurysms. The temporal and spatial variations of wall pressure and wall shear stress in the aneurysm are hypothesized to be correlated with its continuous expansion and eventual rupture. In addition, the assessment of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils. This paper describes the flow dynamics in two representative models of a terminal aneurysm of the basilar artery under Newtonian and non-Newtonian fluid assumptions, and compares their hemodynamics with that of a healthy basilar artery. Virtual aneurysm models are investigated numerically, with geometric features defined by beta = 0 deg and beta = 23.2 deg, where beta is the tilt angle of the aneurysm dome with respect to the basilar artery. The intra-aneurysmal pulsatile flow shows complex ring vortex structures for beta = 0 deg and single recirculation regions for beta = 23.2 deg during both systole and diastole. The pressure and shear stress on the aneurysm wall exhibit large temporal and spatial variations for both models. When compared to a non-Newtonian fluid, the symmetric aneurysm model (beta = 0 deg) exhibits a more unstable Newtonian flow dynamics, although with a lower peak wall shear stress than the asymmetric model (beta = 23.2 deg). The non-Newtonian fluid assumption yields more stable flows than a Newtonian fluid, for the same inlet flow rate. Both fluid modeling assumptions, however, lead to asymmetric oscillatory flows inside the aneurysm dome.     

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.     

On the effect of parent-aneurysm angle on flow patterns in basilar tip aneurysms: towards a surrogate geometric marker of intra-aneurismal hemodynamics

As a part of previous computational fluid dynamic (CFD) validation studies, particle image velocimetry (PIV) of two anatomically realistic basilar artery tip aneurysm models revealed two distinct types of flow (one of which has yet to be reported in the literature), characterized by the location and strength of the intra-aneurismal vortex. We hypothesized that these distinct "hemodynamic phenotypes" could be anticipated by a simple geometric parameter: the angle of the aneurysm bulb relative to the parent artery. An idealized basilar tip aneurysm model was constructed to allow independent control of this angle, and CFD simulations were carried out for angles ranging from 2 degrees to 30 degrees , these extremes corresponding to the angles measured from the two anatomically realistic models. The gross hemodynamics predicted by the idealized model for 2 degrees and 30 degrees were consistent with those seen in the corresponding anatomically realistic models. For the idealized model, the flow type switched at an angle between 8 degrees and 12 degrees . Sensitivity studies suggested that, near these angles, the hemodynamic phenotype was sensitive to inflow momentum. Outside this range, however, the parent-bulb angle appeared to be a robust predictor of hemodynamic phenotype. Our findings suggest that blood flow dynamics in basilar artery tip aneurysms fall into one of the two broad phenotypes, each subject to distinct hemodynamic forces. That the general features of these flow types may be anticipated by a relatively simple-to-measure geometric parameter could help ease the introduction of hemodynamic information into routine clinical decision-making.     

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.     

Hemodynamic stress in lateral saccular aneurysms

The flow velocities in glass and silastic lateral aneurysm models were quantitatively measured with the non-invasive laser Doppler method. The influences of the elasticity of the wall, the pulse wave and the properties of the perfusion medium on the intra-aneurysmal circulation were investigated. As shown previously, the inflow into the aneurysm arose from the downstream lip and was directed toward the center of the fundus. Backflow to the parent vessel took place along the walls of the fundus. With non-pulsatile perfusion, flow velocities in the center of the standardized aneurysms varied between 0.4 and 2% of the maximum velocity in the parent vessel. With pulsatile perfusion, flow velocities in the center of the fundus ranged between 8 and 13% of the flow velocity in the axis of the parent vessel. Flow velocities in the aneurysms were slower with a polymer suspension with blood-like properties compared to a glycerol/water solution. Flow velocity measurements near the aneurysmal wall allowed the estimation of the shear stresses at critical locations. The maximum shear stresses at the downstream lip of the aneurysm were in the range of the stresses measured at the flow divider of an arterial bifurcation. The present results suggest that in human saccular aneurysms intra-aneurysmal flow and shear stress on the wall are directly related to the pulsatility of perfusion, i.e. the systolic/diastolic pressure difference and that the tendency to spontaneous thrombosis depends on the viscoelastic properties of the blood, namely the hematocrit.     

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.     

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.     

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.     

Morphological and Hemodynamic Discriminators for Rupture Status in Posterior Communicating Artery Aneurysms

Background and Purpose

The conflicting findings of previous morphological and hemodynamic studies on intracranial aneurysm rupture may be caused by the relatively small sample sizes and the variation in location of the patient-specific aneurysm models. We aimed to determine the discriminators for aneurysm rupture status by focusing on only posterior communicating artery (PCoA) aneurysms.

Materials and Methods

In 129 PCoA aneurysms (85 ruptured, 44 unruptured), clinical, morphological and hemodynamic characteristics were compared between the ruptured and unruptured cases. Multivariate logistic regression analysis was performed to determine the discriminators for rupture status of PCoA aneurysms.

Results

While univariate analyses showed that the size of aneurysm dome, aspect ratio (AR), size ratio (SR), dome-to-neck ratio (DN), inflow angle (IA), normalized wall shear stress (NWSS) and percentage of low wall shear stress area (LSA) were significantly associated with PCoA aneurysm rupture status. With multivariate analyses, significance was only retained for higher IA (OR = 1.539, p < 0.001) and LSA (OR = 1.393, p = 0.041).

Conclusions

Hemodynamics and morphology were related to rupture status of intracranial aneurysms. Higher IA and LSA were identified as discriminators for rupture status of PCoA aneurysms.     

Effects of stent porosity on hemodynamics in a sidewall aneurysm model

Computation and experiment have been complementarily performed to study the fluid flow inside a stented lateral aneurysm anchored on the straight parent vessel. The implicit solver was based on the time-dependent incompressible Navier-Stokes equations of laminar flow. Solutions were generated by a cell-center finite-volume method that used second-order upwind and second-order center flux difference splitting for the convection and diffusion term, respectively. The second-order Crank-Nicolson method was used in the time integration term. Experimental techniques used were flow visualization (FV) and particle tracking velocimetry (PTV). Experimentally, the straight afferent vessel had an inner diameter 10mm. The diameters of the aneurysmal orifice, neck, and fundus were 14, 10, and 15 mm, respectively, and the distance between the orifice and dome measured 20mm. A 30 mm long helix-shaped stent was tested. Four stent porosities of 100%, 70%, 50%, and 25% were examined. Volume-flow rate waveform of a cerebral artery was considered with a maximum Reynolds number of 250 and Womersley number of 3.9. Results are presented in terms of the pulsatile main and secondary flow velocity vector fields, the volume inflow rates into the aneurysm, and the wall shear stress (WSS) and wall pressure at the aneurysm dome. Some comparisons of computed results with the present FV and PTV results and with the data available from the literature are also made. The maximum flow velocity inside the aneurysm ostium and the WSS in the dome region at the peak flow can, respectively, be suppressed to less than 5% of the parent vessel bulk velocity (or 20% of the unstented case) and 8% of the unstented case if the stent porosity is smaller than 40% (about the porosity of the two-layer stents). In general, the three-layer stents seem not as effective as the two-layer stents in reducing the magnitude of aneurysm inflow rate and WSS.     

犬囊状动脉瘤模型的制作

本文报道用“静脉囊镶嵌技术·制成犬的囊状动脉瘤模型。18个模型(6个单侧型,6个分叉型,6个末梢型)造型后2周经IA DSA检查。本模型在分型、血流动力学改变方面与人类囊状脑动脉瘤类似。不同类型的动脉瘤模型既有相同的血流动力学特征,又有各自的特点,这与动脉瘤与载瘤动脉的角度有关。我们认为该模型可应用于研究动脉瘤的血流动力学与血管内栓塞治疗。     

Onset and nature of flow-induced vibrations in cerebral aneurysms via fluid–structure interaction simulations

Clinical, experimental, and recent computational studies have demonstrated the presence of wall vibrations in cerebral aneurysms, thought to be induced by blood flow instability. These vibrations could induce irregular, high-rate deformation of the aneurysm wall, and potentially disrupt regular cell behavior and promote deleterious wall remodeling. In order to elucidate, for the first time, the onset and nature of such flow-induced vibrations, in this study we imposed a linearly increasing flow rate on high-fidelity fluid–structure interaction models of three anatomically realistic aneurysm geometries. Prominent narrow-band vibrations in the range of 100–500 Hz were found in two out of the three aneurysm geometries tested, while the case that did not exhibit flow instability did not vibrate. Aneurysm vibrations consisted mostly of fundamental modes of the entire aneurysm sac, with the vibrations exhibiting more frequency content at higher frequencies than the flow instabilities driving those vibrations. The largest vibrations occurred in the case which exhibited strongly banded fluid frequency content, and the vibration amplitude was highest when the strongest fluid frequency band was an integer multiple of one of the natural frequencies of the aneurysm sac. Lower levels of vibration occurred in the case which exhibited turbulent-like flow with no distinct frequency bands. The current study provides a plausible mechanistic explanation for the high-frequency sounds observed in cerebral aneurysms, and suggests that narrow-band (vortex-shedding type) flow might stimulate the wall more, or at least at lower flow rates, than broad-band, turbulent-like flow.

     

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.     

Numerical simulation of hemodynamics in stented internal carotid aneurysm based on patient-specific model

There is still a considerable lack of quantitative information concerning the effects of stent structures on blood flow in an aneurismal cavity. In this paper, five virtual stents with different structures and wire cross-sections were designed for incorporation into the same patient-specific aneurysm model. Computational fluid dynamics simulations were performed so as to study how these five types of stents modified hemodynamic parameters. Numerical results demonstrated that the mean flow rate in the aneurismal cavity decreased the most in the model that used a stent with a rectangular wire cross-section, and that the wall shear stresses at the dome and neck of the aneurysm decreased more in models that used a stent with a circular wire cross-section or a spiral stent with a rectangular wire cross-section compared to other models. In addition, the wall pressure on the aneurysm increased slightly after implantation of the stent in all five models. This result differs from that previously published, and may help guide the design and assist clinicians in selecting an appropriate stent for treating cerebral aneurysms.     

Flow instability detected in ruptured versus unruptured cerebral aneurysms at the internal carotid artery

Flow instability has emerged as a new hemodynamic metric hypothesized to have potential value in assessing the rupture risk of cerebral aneurysms. However, diverse findings have been reported in the literature. In the present study, high-resolution hemodynamic simulations were performed retrospectively on 35 aneurysms (10 ruptured & 25 unruptured) located at the internal carotid artery (ICA). Simulated hemodynamic parameters were statistically compared between the ruptured and unruptured aneurysms, with emphasis on examining the correlation of flow instability with the status of aneurysm rupture. Pronounced flow instability was detected in 20% (2 out of 10) of the ruptured aneurysms, whereas in 44% (11 out of 25) of the unruptured aneurysms. Statistically, the flow instability metric (quantified by the temporally and spatially averaged fluctuating kinetic energy over the aneurysm sac) did not differ significantly between the ruptured and unruptured aneurysms. In contrast, low wall shear stress area (LSA) and pressure loss coefficient (PLC) exhibited significant correlations with the status of aneurysm rupture. In conclusion, the present study suggests that the presence of flow instability may not correlate closely with the status of aneurysm rupture, at least for ICA aneurysms. On the other hand, the retrospective nature of the study and the small sample size may have to some extent compromised the reliability of the conclusion, and therefore large-scale prospective studies would be needed to further address the issue.     

Pulsatile flow effects on the hemodynamics of intracranial aneurysms

High-resolution numerical simulations are carried out to systematically investigate the effect of the incoming flow waveform on the hemodynamics and wall shear stress patterns of an anatomic sidewall intracranial aneurysm model. Various wave forms are constructed by appropriately scaling a typical human waveform such that the waveform maximum and time-averaged Reynolds numbers, the Womersley number (α), and the pulsatility index (PI) are systematically varied within the human physiologic range. We show that the waveform PI is the key parameter that governs the vortex dynamics across the aneurysm neck and the flow patterns within the dome. At low PI, the flow in the dome is similar to a driven cavity flow and is characterized by a quasi-stationary shear layer that delineates the parent artery flow from the recirculating flow within the dome. At high PI, on the other hand, the flow is dominated by vortex ring formation, transport across the neck, and impingement and breakdown at the distal wall of the aneurysm dome. We further show that the spatial and temporal characteristics of the wall shear stress field on the aneurysm dome are strongly correlated with the vortex dynamics across the neck. We finally argue that the ratio between the characteristic time scale of transport by the mean flow across the neck and the time scale of vortex ring formation can be used to predict for a given sidewall aneurysm model the critical value of the waveform PI for which the hemodynamics will transition from the cavity mode to the vortex ring mode.     

Numerical and experimental studies on pulsatile flow in aneurysms arising laterally from a curved parent vessel at various angles

Both numerical and experimental studies have been performed to characterize the fluid flow inside the lateral aneurysms arising from the curved parent vessels at various angles gamma. The implicit solver was based on the time-dependent Navier-Stokes equations of incompressible laminar flow. Solutions were generated by a cell-center finite-volume method that used second order upwind and second order center flux difference splitting for the convection and diffusion term, respectively. The second order Crank-Nicolson method was used in the time integration term while the SIMPLEC algorithm was adopted to handle the pressure-velocity coupling. Complementarily, the particle tracking velocimetry (PTV) was used to measure the velocity fields. The conditions selected were to simulate an internal carotid artery with a diameter of 5 mm by similarity rules. The values of gamma explored were 0 degrees, 45 degrees, 90 degrees, and 135 degrees. Pulsatile flow with Wormersley number 3.9 and Reynolds numbers varying from 350 to 850 was considered. The computed results are firstly verified by the PTV measured ones. Discussion of the results is in terms of pulsatile main and secondary velocity vector fields, inflow rates into the aneurysm, and the distributions of wall shear stress and static pressure. It is found that among the angles examined gamma=45( composite function) is the riskiest angle from a fluid dynamics point of view and the aneurysmal dome is at risk.