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.     

PIV-measured versus CFD-predicted flow dynamics in anatomically realistic cerebral aneurysm models

Computational fluid dynamics (CFD) modeling of nominally patient-specific cerebral aneurysms is increasingly being used as a research tool to further understand the development, prognosis, and treatment of brain aneurysms. We have previously developed virtual angiography to indirectly validate CFD-predicted gross flow dynamics against the routinely acquired digital subtraction angiograms. Toward a more direct validation, here we compare detailed, CFD-predicted velocity fields against those measured using particle imaging velocimetry (PIV). Two anatomically realistic flow-through phantoms, one a giant internal carotid artery (ICA) aneurysm and the other a basilar artery (BA) tip aneurysm, were constructed of a clear silicone elastomer. The phantoms were placed within a computer-controlled flow loop, programed with representative flow rate waveforms. PIV images were collected on several anterior-posterior (AP) and lateral (LAT) planes. CFD simulations were then carried out using a well-validated, in-house solver, based on micro-CT reconstructions of the geometries of the flow-through phantoms and inlet/outlet boundary conditions derived from flow rates measured during the PIV experiments. PIV and CFD results from the central AP plane of the ICA aneurysm showed a large stable vortex throughout the cardiac cycle. Complex vortex dynamics, captured by PIV and CFD, persisted throughout the cardiac cycle on the central LAT plane. Velocity vector fields showed good overall agreement. For the BA, aneurysm agreement was more compelling, with both PIV and CFD similarly resolving the dynamics of counter-rotating vortices on both AP and LAT planes. Despite the imposition of periodic flow boundary conditions for the CFD simulations, cycle-to-cycle fluctuations were evident in the BA aneurysm simulations, which agreed well, in terms of both amplitudes and spatial distributions, with cycle-to-cycle fluctuations measured by PIV in the same geometry. The overall good agreement between PIV and CFD suggests that CFD can reliably predict the details of the intra-aneurysmal flow dynamics observed in anatomically realistic in vitro models. Nevertheless, given the various modeling assumptions, this does not prove that they are mimicking the actual in vivo hemodynamics, and so validations against in vivo data are encouraged whenever possible.     

Hemodynamic study of overlapping bare-metal stents intervention to aortic aneurysm

To investigate the hemodynamic performance of overlapping bare-metal stents intervention treatment to thoracic aortic aneurysms (TAA), three simplified TAA models, representing, no stent, with a single stent and 2 overlapped stents deployed in the aneurismal sac, were studied and compared in terms of flow velocity, wall shear stress (WSS) and pressure distributions by means of computational fluid dynamics. The results showed that overlapping stents intervention induced a flow field of slow velocity near the aneurismal wall. Single stent deployment in the sac reduced the jet-like flow formed prior to the proximal neck of the aneurysm, which impinged on the internal wall of the aneurysm. This jet-like flow vanished completely in the overlapping double stents case. Overlapping stents intervention led to an evident decrease in WSS; meanwhile, the pressure acting on the wall of the aneurysm was reduced slightly and presented more uniform distribution. The results therefore indicated that overlapping stents intervention may effectively isolate the thoracic aortic aneurysm, protecting it from rupture. In conclusion, overlapping bare-metal stents may serve a purpose similar to that of the multilayer aneurysm repair system (MARS) manufactured by Cardiatis SA (Isnes, Belgium).     

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.     

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.     

Stability of pulsatile blood flow at the ostium of cerebral aneurysms

The strength and direction of blood flow into and within a cerebral aneurysm are important issues in developing effective interventional strategies to stabilize the aneurysm. We tested the hypothesis that there are significant major hemodynamic features that are common to many aneurysm flows of the type studied here. This was investigated by performing computational fluid dynamic simulations of flow near 7 cerebral aneurysms using geometrical data obtained from clinical CT scans. Our numerical simulations of flow across the ostium plane of an aneurysm show that in many cases there is relatively stable flow structure that is maintained over the phase of the pulsatile flow cycle. The two main features of this flow are (1) quasi-permanent regions of flow influx and efflux across the ostium plane exist, separated by a “virtual boundary”, and (2) a helical vortex flow pattern within the aneurismal sac with swirl in two orthogonal cross-sectional planes. These numerical observations are consistent with in vitro experimental data from ultrasound color-Doppler velocimetry and other numerical and experimental studies. The observed flow patterns are found to occur in different types of aneurysms (bifurcation and sidewall), and can persist even after flow parameters are perturbed beyond the normal range of physiological flow conditions. These results suggest that in many cases, major aspects of the behavior of aneurismal hemodynamics for important classes of aneurysms can be learned from an analysis of steady, non-pulsatile flow, which is simpler and faster to simulate than time-dependent, pulsatile flow. An understanding of this fluid dynamical behavior may also prove useful in the design of stents, coils, and various other endovascular flow diverting devices.     

Influence of stent configuration on cerebral aneurysm fluid dynamics

Embolic coiling is the most popular endovascular treatment available for cerebral aneurysms. Nevertheless, the embolic coiling of wide-neck aneurysms is challenging and, in many cases, ineffective. Use of highly porous stents to support coiling of wide-neck aneurysms has become a common procedure in recent years. Several studies have also demonstrated that high porosity stents alone can significantly alter aneurysmal hemodynamics, but differences among different stent configurations have not been fully characterized. As a result, it is usually unclear which stent configuration is optimal for treatment. In this paper, we present a flow study that elucidates the influence of stent configuration on cerebral aneurysm fluid dynamics in an idealized wide-neck basilar tip aneurysm model. Aneurysmal fluid dynamics for three different stent configurations (half-Y, Y and, cross-bar) were first quantified using particle image velocimetry and then compared. Computational fluid dynamics (CFD) simulations were also conducted for selected stent configurations to facilitate validation and provide more detailed characterizations of the fluid dynamics promoted by different stent configurations. In vitro results showed that the Y stent configuration reduced cross-neck flow most significantly, while the cross-bar configuration reduced velocity magnitudes within the aneurysmal sac most significantly. The half-Y configuration led to increased velocity magnitudes within the aneurysmal sac at high parent-vessel flow rates. Experimental results were in strong agreement with CFD simulations. Simulated results indicated that differences in fluid dynamic performance among the different stent configurations can be attributed primarily to protruding struts within the bifurcation region.     

Posterior Cerebral Artery Angle and the Rupture of Basilar Tip Aneurysms

Since the initial publication of the International Study of Unruptured Intracranial Aneurysms (ISUIA), management of unruptured intracranial aneurysms has been mainly based on the size of the aneurysm. The contribution of morphological characteristics to treatment decisions of unruptured aneurysms has not been well studied in a systematic and location specific manner. We present a large sample of basilar artery tip aneurysms (BTA) that were assessed using a diverse array of morphological variables to determine the parameters associated with ruptured aneurysms. Demographic and clinical risk factors of aneurysm rupture were obtained from chart review. CT angiograms (CTA) were evaluated with Slicer, an open source visualization and image analysis software, to generate 3-D models of the aneurysms and surrounding vascular architecture. Morphological parameters examined in each model included aneurysm volume, aspect ratio, size ratio, aneurysm angle, basilar vessel angle, basilar flow angle, and vessel to vessel angles. Univariate and multivariate analyses were performed to determine statistical significance. From 2008–2013, 54 patients with BTA aneurysms were evaluated in a single institution, and CTAs from 33 patients (15 ruptured, 18 unruptured) were available and analyzed. Aneurysms that underwent reoperation, that were associated with arteriovenous malformations, or that lacked preoperative CTA were excluded. Multivariate logistic regression revealed that a larger angle between the posterior cerebral arteries (P1-P1 angle, p = 0.037) was most strongly associated with aneurysm rupture after adjusting for other morphological variables. In this location specific study of BTA aneurysms, the larger the angle formed between posterior cerebral arteries was found to be a new morphological parameter significantly associated with ruptured BTA aneurysms. This is a physically intuitive parameter that can be measured easily and readily applied in the clinical setting.     

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.

     

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.     

Endothelial nitric oxide synthase and calcium production in arterial geometries: an integrated fluid mechanics/cell model

It is well known that atherosclerosis occurs at very specific locations throughout the human vasculature, such as arterial bifurcations and bends, all of which are subjected to low wall shear stress. A key player in the pathology of atherosclerosis is the endothelium, controlling the passage of material to and from the artery wall. Endothelial dysfunction refers to the condition where the normal regulation of processes by the endothelium is diminished. In this paper, the blood flow and transport of the low diffusion coefficient species adenosine triphosphate (ATP) are investigated in a variety of arterial geometries: a bifurcation with varying inner angle, and an artery bend. A mathematical model of endothelial calcium and endothelial nitric oxide synthase cellular dynamics is used to investigate spatial variations in the physiology of the endothelium. This model allows assessment of regions of the artery wall deficient in nitric oxide (NO). The models here aim to determine whether 3D flow fields are important in determining ATP concentration and endothelial function. For ATP transport, the effects of a coronary and carotid wave form on mass transport is investigated for low Womersley number. For the carotid, the Womersley number is then increased to determine whether this is an important factor. The results show that regions of low wall shear stress correspond with regions of impaired endothetial nitric oxide synthase signaling, therefore reduced availability of NO. However, experimental work is required to determine if this level is significant. The results also suggest that bifurcation angle is an important factor and acute angle bifurcations are more susceptible to disease than large angle bifurcations. It has been evidenced that complex 3D flow fields play an important role in determining signaling within endothelial cells. Furthermore, the distribution of ATP in blood is highly dependent on secondary flow features. The models here use ATP concentration simulated under steady conditions. This has been evidenced to reproduce essential features of time-averaged ATP concentration over a cardiac cycle for small Womersley numbers. However, when the Womersley number is increased, some differences are observed. Transient variations are overall insignificant, suggesting that spatial variation is more important than temporal. It has been determined that acute angle bifurcations are potentially more susceptible to atherogenesis and steady-state ATP transport reproduces essential features of time-averaged pulsatile transport for small Womersley number. Larger Womersley numbers appear to be an important factor in time-dependent mass transfer.     

Analysis of Morphological Parameters to Differentiate Rupture Status in Anterior Communicating Artery Aneurysms

In contrast to size, the association of morphological characteristics of intracranial aneurysms with rupture has not been established in a systematic manner. We present an analysis of the morphological variables that are associated with rupture in anterior communicating artery aneurysms to determine site-specific risk variables. One hundred and twenty-four anterior communicating artery aneurysms were treated in a single institution from 2005 to 2010, and CT angiograms (CTAs) or rotational angiography from 79 patients (42 ruptured, 37 unruptured) were analyzed. Vascular imaging was evaluated with 3D Slicer© to generate models of the aneurysms and surrounding vasculature. Morphological parameters were examined using univariate and multivariate analysis and included aneurysm volume, aspect ratio, size ratio, distance to bifurcation, aneurysm angle, vessel angle, flow angle, and parent-daughter angle. Multivariate logistic regression revealed that size ratio, flow angle, and parent-daughter angle were associated with aneurysm rupture after adjustment for age, sex, smoking history, and other clinical risk factors. Simple morphological parameters such as size ratio, flow angle, and parent-daughter angle may thus aid in the evaluation of rupture risk of anterior communicating artery aneurysms.     

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.     

Mathematical models of the flow in the basilar artery

The flow in the basilar artery arises from the merging of the flows from the two vertebral arteries. This study deals with the question whether a parabolic (Poiseuille) profile will have been established before the basilar artery divides into both posterior cerebral arteries. The inlet length (that is, the downstream distance needed for the flow to become approximately equal to the limiting Poiseuille flow) and velocity profiles have been computed from two- and three-dimensional mathematical models in which flow pulsatility and vessel wall distensibility have been neglected and the complex geometry of the junction has been taken into account in a simplified form. The results show that the flow at the end of the basilar artery is far from being parabolic and that an asymmetry in the entrance flow will be carried along towards the end of the basilar artery, thus affecting flows in the circle of Willis.     

Errors in power-law estimations of inflow rates for intracranial aneurysm CFD

Patient-specific inflow rates are rarely available for computational fluid dynamics (CFD) studies of intracranial aneurysms. Instead, inflow rates are often estimated from parent artery diameters via power laws, i.e. Q ∝ Dⁿ, reflecting adaptation of conduit arteries to demanded flow. The present study aimed to validate the accuracy of these power laws. Internal carotid artery (ICA) flow rates were measured from 25 ICA aneurysm patients via 2D phase contrast MRI. ICA diameters, derived from 3D segmentation of rotational angiograms, were used to estimate inflow rates via power laws from the aneurysm CFD literature assuming the same inlet wall shear stress (WSS) (n = 3), velocity (n = 2) or flow rate (n = 0) for all cases. To illustrate the potential impact of errors in flow rate estimates, pulsatile CFD was carried out for four cases having large errors for at least one power law. Flow rates estimated by n = 3 and n = 0 power laws had significant (p < 0.01) mean biases of −22% to +32%, respectively, but with individual errors ranging from −78% to +120%. The n = 2 power law had no significant bias, but had non-negligible individual errors of −58% to +71%. CFD showed similarly large errors for time-averaged sac WSS; however, these were reduced after normalizing by parent artery WSS. High frequency WSS fluctuations, evident in 2/4 aneurysms, were also sensitive to inflow rate errors. Care should therefore be exercised in the interpretation of aneurysm CFD studies that rely on power law estimates of inflow rates, especially if absolute (vs. normalized) WSS, or WSS instabilities, are of interest.     

Numerical simulation of pre- and postsurgical flow in a giant basilar aneurysm

Computational modeling of the flow in cerebral aneurysms is an evolving technique that may play an important role in surgical planning. In this study, we simulated the flow in a giant basilar aneurysm before and after surgical takedown of one vertebral artery. Patient-specific geometry and flowrates obtained from magnetic resonance (MR) angiography and velocimetry were used to simulate the flow prior to and after the surgery. Numerical solutions for steady and pulsatile flows were obtained. Highly three-dimensional flows, with strong secondary flows, were computed in the aneurysm in the presurgical and postsurgical conditions. The computational results predicted that occlusion of a vertebral artery would result in a significant increase of the slow flow region formed in the bulge of the aneurysm, where increased particle residence time and velocities lower than 2.5 cms were computed. The region of slow flow was found to have filled with thrombus following surgery. Predictions of numerical simulation methods are consistent with the observed outcome following surgical treatment of an aneurysm. The study demonstrates that computational models may provide hypotheses to test in future studies, and might offer guidance for the interventional treatment of cerebral aneurysms.     

A study of wall shear stress in 12 aneurysms with respect to different viscosity models and flow conditions

Recent computational fluid dynamics (CFD) studies relate abnormal blood flow to rupture of cerebral aneurysms. However, it is still debated how to model blood flow with sufficient accuracy. Common assumptions made include Newtonian behaviour of blood, traction free outlet boundary conditions and inlet boundary conditions based on available literature. These assumptions are often required since the available patient specific data is usually restricted to the geometry of the aneurysm and the surrounding vasculature. However, the consequences of these assumptions have so far been inadequately addressed.     

Aneurysm Characteristics Associated with the Rupture Risk of Intracranial Aneurysms: A Self-Controlled Study

This study analyzed the rupture risk of intracranial aneurysms (IAs) according to aneurysm characteristics by comparing the differences between two aneurysms in different locations within the same patient. We utilized this self-controlled model to exclude potential interference from all demographic factors to study the risk factors related to IA rupture. A total of 103 patients were diagnosed with IAs between January 2011 and April 2015 and were enrolled in this study. All enrolled patients had two IAs. One IA (the case) was ruptured, and the other (the control) was unruptured. Aneurysm characteristics, including the presence of a daughter sac, the aneurysm neck, the parent artery diameter, the maximum aneurysm height, the maximum aneurysm width, the location, the aspect ratio (AR, maximum perpendicular height/average neck diameter), the size ratio (SR, maximum aneurysm height/average parent diameter) and the width/height ratio (WH ratio, maximum aneurysm width/maximum aneurysm height), were collected and analyzed to evaluate the rupture risks of the two IAs within each patient and to identify the independent risk factors associated with IA rupture. Multivariate, conditional, backward, stepwise logistic regression analysis was performed to identify the independent risk factors associated with IA rupture. The multivariate analysis identified the presence of a daughter sac (odds ratio [OR], 13.80; 95% confidence interval [CI], 1.65–115.87), a maximum aneurysm height ≥7 mm (OR, 4.80; 95% CI, 1.21–18.98), location on the posterior communicating artery (PCOM) or anterior communicating artery (ACOM; OR, 3.09; 95% CI, 1.34–7.11) and SR (OR, 2.13; 95% CI, 1.16–3.91) as factors that were significantly associated with IA rupture. The presence of a daughter sac, the maximum aneurysm height, PCOM or ACOM locations and SR (>1.5±0.7) of unruptured IAs were significantly associated with IA rupture.     

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.