Analysis of Hemodynamics and Aneurysm Occlusion after Flow-Diverting Treatment in Rabbit Models =============================================================================================== * J.R. Cebral * F. Mut * M. Raschi * S. Hodis * Y.-H. Ding * B.J. Erickson * R. Kadirvel * D.F. Kallmes ## Abstract **BACKGROUND AND PURPOSE:** Predicting the outcome of flow diversion treatment of cerebral aneurysms remains challenging. Our aim was to investigate the relationship between hemodynamic conditions created immediately after flow diversion and subsequent occlusion of experimental aneurysms in rabbits. **MATERIALS AND METHODS:** The hemodynamic environment before and after flow-diversion treatment of elastase-induced aneurysms in 20 rabbits was modeled by using image-based computational fluid dynamics. Local aneurysm occlusion was quantified by using a voxelization technique on 3D images acquired 8 weeks after treatment. Global and local voxel-by-voxel hemodynamic variables were used to statistically compare aneurysm regions that later thrombosed to regions that remained patent. **RESULTS:** Six aneurysms remained patent at 8 weeks, while 14 were completely or nearly completely occluded. Patent aneurysms had statistically larger neck sizes (*P* = .0015) and smaller mean transit times (*P* = .02). The velocity, vorticity, and shear rate were approximately 2.8 times (*P* < .0001) larger in patent regions—that is, they had larger “flow activity” than regions that progressed to occlusion. Statistical models based on local hemodynamic variables were capable of predicting local occlusion with good precision (84% accuracy), especially away from the neck (92%–94%). Predictions near the neck were poorer (73% accuracy). **CONCLUSIONS:** These results suggests that the dominant healing mechanism of occlusion within the aneurysm dome is related to slow-flow-induced thrombosis, while near the neck, other processes could be at play simultaneously. ## ABBREVIATIONS: CFD : computational fluid dynamics FD : flow-diverting 3DRA : 3D rotational angiography Despite the increased use of flow-diverting (FD) devices for treating intracranial aneurysms, the exact effects and processes responsible for the evolution of the aneurysm, and the ultimate outcome of these procedures are poorly understood.1⇓⇓–4 Intrasaccular thrombosis and endothelial cell growth across the aneurysm neck have been observed after flow diversion and have been proposed as the main mechanisms driving the healing process.5⇓–7 However, the interaction, dominance, relative importance, and time precedence of these processes are not well-understood. Knowledge of these mechanisms is important to facilitate the development of future devices and therapies. The purpose of this study was to investigate the relationship between the hemodynamic conditions created immediately after deployment of FD devices and the subsequent aneurysm occlusion in a rabbit model of saccular aneurysms. ## Materials and Methods ### Animal Models and Imaging Twenty-three elastase-induced aneurysms were created in New Zealand white rabbits, following the approach described in Altes et al.8 Four weeks after their creation, the aneurysms were treated with a FD device (Pipeline Embolization Device; Covidien, Irvine, California). Two days before treatment, the subjects were premedicated with aspirin (10 mg/kg by mouth) and clopidogrel (10 mg/kg by mouth) and were continued on these for 1 month after treatment. Immediately before treatment, 3D rotational angiography (3DRA) images were acquired, and velocities in the surrounding vessels were measured with Doppler sonography. Six animals were sacrificed before 1 week after treatment; all others, after 8 weeks. Immediately before sacrifice, 3DRA imaging was repeated. Some of the rabbits used in this study were part of another investigation in which we analyzed the mechanism of the endothelialization after flow diverter implantation. This article is entirely unrelated to the previous study. ### Hemodynamics Modeling Subject-specific computational fluid dynamics (CFD) models were constructed from pretreatment 3DRA images.9 Unstructured grids were generated with a resolution of 0.2 mm. Models of the FD devices used to treat the aneurysms were created and virtually deployed within the reconstructed vascular models.10 Blood flows were modeled by solving the unsteady 3D incompressible Navier-Stokes equations.11 Physiologic flow conditions were derived from the Doppler sonography velocity measurements and imposed as boundary conditions in the computational models. Wall compliance was not included in the model. Blood attenuation was set to ρ = 1.0 g/cm3, and blood viscosity, to μ = 0.04 poise. The governing equations were numerically solved by using an efficient finite-element solver with a time-step of 0.01 seconds.12 After deployment of the FD device, the mesh was adaptively refined to resolve the stent wires and a new CFD simulation was performed by using immersed grids.13,14 ### Occlusion Modeling Regions of the aneurysm that remained open or patent and regions that were occluded at the time of sacrifice were identified by constructing a second vascular model from the 3DRA image acquired before sacrifice. The pretreatment model and this premortem or follow-up model were manually aligned by using rigid registration. A new grid was generated filling the volume of the follow-up model. The aneurysm neck was interactively delineated on the pretreatment model by connecting selected points along the path of minimum geodesic distance.15 The aneurysm orifice defined by the delineated neck was triangulated, and grid points on each side of this surface were labeled as belonging to the aneurysm or the vessel.15 The bounding box of the aneurysm was found and voxelized with isotropic voxels of the same resolution as the 3DRA images. Mesh points in the aneurysm region were then labeled as “open” or “occluded,” depending on whether they were inside or outside the aligned follow-up model. The methodology is illustrated in Fig 1. ![Fig 1.](http://www.ajnr.org/https://www.ajnr.org/content/ajnr/35/8/1567/F1.medium.gif) [Fig 1.](http://www.ajnr.org/content/35/8/1567/F1) Fig 1. Imaging, modeling, and data analysis. ### Global (per Aneurysm) Data Analysis The degree of aneurysm occlusion at the time of sacrifice was quantified by the percentage of occluded voxels identified as described above. Aneurysms were then classified into 2 groups: 1) patent group: <80% of the voxels occluded at follow-up; 2) occluded group: >80% of the voxels occluded at follow-up. This 80% threshold was chosen because aneurysms with >80% occlusion of their volume looked almost completely closed on DSA images with perhaps a small “bump” in the parent artery—that is, they did not exhibit remnant necks of persistent filling of the sac. Using the 3D meshes, we calculated geometric variables, including volume, maximum size, neck area, and maximum neck size for each aneurysm. Similarly, the following hemodynamic quantities were calculated over the aneurysm region and averaged over the cardiac cycle: mean aneurysm inflow rate , mean aneurysm velocity , mean aneurysm wall shear stress <τ>, and mean aneurysm shear rate <γ̇>. These quantities were computed before and after treatment. Additionally, a new variable denoted “mean aneurysm transit time” (MATT) was defined as ![Formula][1] where *V* is the aneurysm volume, *A* is the neck area, and <*v*> is the mean aneurysm velocity. The rationale is that this variable, which has units of time, should do the following: 1. Increase with the aneurysm volume (ie, it takes longer to traverse a large aneurysm) 2. Decrease with the neck area (ie, a larger neck implies more chances of getting out of the aneurysm, while a smaller neck implies more chances of staying within the aneurysm and recirculating) 3. Decrease with the blood velocity (ie, slower velocities imply longer transit times). The values of geometric and hemodynamic variables were then averaged over the patent and occluded groups. Cases belonging to the incompletely occluded group sacrificed before 1 week after treatment (between brackets in the Tables) were excluded from the statistical analysis. The Wilcoxon nonparametric test was used to test whether the means were statistically different between the occluded and patent groups. ### Local (Voxel-by-Voxel) Data Analysis Velocity (v), vorticity (ω), and shear rate (γ̇) were calculated on the CFD mesh and interpolated and averaged over each aneurysm voxel. The mean values of these quantities over open or occluded voxels of each aneurysm were calculated and compared. Predictive statistical models of regional occlusion were then constructed. We defined 2 classes or outcomes: 1) occluded voxels (+1), and 2) open voxels (−1). Seven attributes or features were considered: 1) pre, 2) <ω>pre, 3) <γ̇>pre, 4) post, 5) <ω>post, 6) <γ̇>post, and 7) d = path distance to the neck. These values were calculated for each voxel from the pretreatment and posttreatment CFD simulations. A logistic regression model16 was trained by using a “leave-one-out” methodology. Specifically, it was trained with all open and occluded voxels from all aneurysms except 1. The aneurysm left out was then used to test the accuracy of the model. The probabilities of each voxel of this aneurysm belonging to each class were calculated, and the voxel was “predicted” to belong to the class of the highest probability. The accuracy of the prediction was calculated by counting the number of correctly and incorrectly classified voxels. The process was repeated by leaving each aneurysm out one at a time, and the total predictive accuracy of the model was computed. Next, the voxels of each aneurysm were divided into 3 groups according to their distance to the aneurysm neck. Denoting by d the distance from a voxel to the neck and dmax the maximum distance from any voxel to the neck, we assigned voxels to the following 3 groups: a) dome if d/dmax > 2/3, b) body if 1/3 ≤ d/dmax ≤ 2/3, and c) neck if d/dmax < 1/3. The predictive statistical analysis described above by using the logistic regression model was repeated for voxels in each of these 3 regions. The corresponding predictive accuracies were calculated and compared. Finally, we computed and compared the predictive power of the following statistical models typically used in machine-learning studies 17: 1) logistic regression, 2) neural network, and 3) support vector machine. Each of these models was trained with data from all aneurysms except 1 and was tested on the one left out. The process was repeated by leaving each aneurysm out in turn. The total predictive accuracy of each model was computed and compared. ## Results ### Hemodynamics Flow visualizations revealed that velocities within the aneurysm are substantially reduced after treatment, that flow structures may change and, in particular, become smoother and simpler (ie, less swirling), and that the location of the inflow stream may shift from the distal part of the neck to the proximal part. To illustrate these observations, we present 2 examples in Fig 2, corresponding to incomplete and complete occlusions. The placement of FD devices thus causes both quantitative and qualitative changes in the hemodynamic environment within the aneurysm. ![Fig 2.](http://www.ajnr.org/https://www.ajnr.org/content/ajnr/35/8/1567/F2.medium.gif) [Fig 2.](http://www.ajnr.org/content/35/8/1567/F2) Fig 2. Visualization of peak systole flow structures pre- and post-treatment in an incompletely occluded aneurysm (case 1, top) and in a completely occluded aneurysm (case 10, bottom). Each panel shows 20 cm/s velocity isosurfaces (top) and flow streamlines (bottom) before (left) and after (right) treatment. ### Global (per Aneurysm) Characteristics Geometric characteristics of the aneurysms, the parent arteries, and the FD devices are presented in Tables 1 and 2. Aneurysms have been ordered by percentage occlusion. The aneurysm volume was not statistically different between the open and occluded aneurysm groups. The neck area (*P* = .0006), maximum neck size (*P* = .0015), and maximum aneurysm size (*P* = .02) were significantly larger in the patent than in the occluded groups. View this table: [Table 1:](http://www.ajnr.org/content/35/8/1567/T1) Table 1: Aneurysm occlusion, geometry, and device dataa View this table: [Table 2:](http://www.ajnr.org/content/35/8/1567/T2) Table 2: Statistics of aneurysm geometrical data Results from the global hemodynamics analysis are presented in the On-line Table. Only the mean inflow rate and mean aneurysm transit time were significantly different between the patent and occluded groups. The inflow rates before (*P* = .0015) and after treatment (*P* = .0305) were approximately 2.8 times larger in the patent group. The mean aneurysm transit time after treatment (*P* = .02) was approximately 3 times larger in occluded aneurysms than in patent aneurysms. The mean posttreatment shear rate achieved marginal significance (*P* = .0757) and was approximately 1.6 times larger in the patent group. All other quantities were not statistically different between the open and occluded groups. To assess the possibility of predicting whether the aneurysm will be completely occluded at follow-up, we performed a receiving operating curve analysis. The area under the curve was calculated for each geometric and hemodynamic variable (listed in Tables 1 and 2 and the On-line Table). The best predictor was neck area and size, with an accuracy of approximately 70%–73%. However, our sample included a small number of aneurysms. and it is likely that this accuracy could be improved with larger samples. ### Regional (Voxel-by-Voxel) Characteristics Results from the local hemodynamics analysis are presented in Fig 3. This figure shows the mean value of posttreatment hemodynamics variables (velocity, vorticity, and shear rate) averaged over the open (red) and occluded (green) voxels of each aneurysm. Aneurysms in the x-axis are ordered by percentage occlusion. Open voxels have larger posttreatment velocity, vorticity, and shear rates than occluded voxels. These differences are less pronounced in aneurysms that had large patent regions (toward the left of Fig 3). The differences seem to increase with the percentage of occluded voxels (to the right of Fig 3). Ratios of posttreatment hemodynamic variables in open voxels over occluded voxels are presented in Table 3 for aneurysms that were not completely occluded. These results indicate that velocity, vorticity, and shear rate are on average larger by a factor of approximately 2.6 in open regions than in occluded regions. The ratios, however, are larger in aneurysms that were nearly completely occluded than in incompletely occluded ones. The *P* values indicate that these differences are statistically significant. ![Fig 3.](http://www.ajnr.org/https://www.ajnr.org/content/ajnr/35/8/1567/F3.medium.gif) [Fig 3.](http://www.ajnr.org/content/35/8/1567/F3) Fig 3. Posttreatment hemodynamic variables averaged over open and occluded regions. View this table: [Table 3:](http://www.ajnr.org/content/35/8/1567/T3) Table 3: Ratio of mean posttreatment hemodynamic variables in open-over-occluded regionsa ### Predictive Models Table 4 presents the accuracy and number of correctly predicted voxels over the total number of voxels for logistic regression models created for the entire aneurysm (sac), the dome, the body, and the neck regions. In addition, a range of accuracy values is provided for each region. These ranges were computed by calculating the accuracy of all possible combinations of the 7 features used to build the models (255 models for each region). Overall, the logistic regression model is capable of predicting which voxels will be occluded with very good accuracy and that voxels near the dome are more easily predicted than voxels near the neck. View this table: [Table 4:](http://www.ajnr.org/content/35/8/1567/T4) Table 4: Accuracy of logistic regression models of local occlusion A comparison of different statistical classifiers typically used in machine learning is presented in Table 5. Overall, all statistical models had similar accuracies. The best performance was obtained with the neural network. The cases that were most difficult to predict were the ones with patent regions that were approximately 60%–70% of the aneurysm volume. The predictive accuracy for aneurysms that were approximately 30%–46% occluded ranged from approximately 40% to 50% (data not shown). View this table: [Table 5:](http://www.ajnr.org/content/35/8/1567/T5) Table 5: Accuracy of different statistical models of local occlusion ## Discussion Flow diverters cause a significant disruption of the aneurysm inflow stream, which induces both qualitative and quantitative changes in the hemodynamic environment within the aneurysm. Aneurysms completely or almost completely occluded after treatment had statistically smaller necks than aneurysms that remained patent. They also had larger mean aneurysm transit times (related to the aneurysm residence time). Immediately after flow diverter implantation, the local velocity, vorticity (speed of rotation or swirling of fluid elements), and shear rate (rate of deformation of fluid elements) were all statistically smaller in regions later observed to occlude than in regions that remained patent. Statistical models based on local hemodynamic conditions were able to predict local occlusion with a very good accuracy. For a given aneurysm, regions close to the neck were more difficult to predict than regions that were close to the aneurysm dome. The mechanism of intrasaccular thrombosis and subsequent aneurysm occlusion is intimately related to the hemodynamic environment created after treatment. This environment depends mainly on the geometry of the aneurysm and parent artery and on the hemodynamic effects of the FD device. Aneurysm morphology, in particular the location of the orifice on the parent artery and the neck size, determines the inflow stream into the aneurysm. The characteristics of the FD device, in turn, determine the disruption of this inflow stream and subsequent recirculation within the aneurysm sac. Thus, associations between both geometric and hemodynamic variables with outcome are expected. However, in this series, all aneurysms were created at the same location and were treated with similar FD devices. Thus, it is not surprising that the best predictor of complete occlusion was the neck size. Further studies with larger samples are needed to better understand the relationships among geometry, hemodynamics, and outcomes. Previous studies have largely focused on the changes of hemodynamic quantities from pre- to posttreatment. Pereira et al18 used dynamic DSA images to estimate the mean aneurysm velocity (mean aneurysm flow amplitude) before and after FD treatment in a pilot clinical study (*n* = 21), and found that at a later follow-up examination, aneurysms that were completely occluded (*n* = 18) had larger reductions of the mean aneurysm flow amplitude than aneurysms still patent (*n* = 3). Zhang et al19 used CFD porous media to model the flow in 2 stented aneurysms, one remaining patent and the other completely occluded. They reported a larger reduction of aneurysm inflow velocity in the occluded than in the patent aneurysm. Huang et al20 analyzed the flow dynamics in 14 rabbit aneurysm models treated with flow diverters. They used micro-CT images of harvested arteries with the flow diverter to create a model of the stent and place it in CFD models and found larger reductions in aneurysm inflow and larger increases of relative residence time in occluded aneurysms (*n* = 10) than in patent aneurysms (*n* = 4), though these differences did not reach statistical significance. They also observed changes in the location of the inflow stream, from the distal part of the neck to a more central/proximal position. In our study, we used CFD models with subject-specific flow conditions and FD models that reproduced the geometric characteristics of the implanted devices and compared global and local hemodynamic quantities in both occluded and patent aneurysms as well as between occluded and patent regions of each aneurysm. Our study focused on hemodynamic conditions created immediately after treatment rather than on changes from pre- to posttreatment. The results of our study suggest that global hemodynamic conditions may not be sufficient to understand or predict aneurysm occlusion following flow-diversion treatment. They also suggest that thrombosis induced by the local hemodynamic conditions seems to be the main driver of aneurysm occlusion at the dome and body of the aneurysm, while other mechanisms such as endothelial tissue growth may be important near the neck. It is also possible that the thrombosis starts at the dome and progresses toward the neck, altering the intra-aneurysmal hemodynamic environment and inducing thrombus formation near the neck later followed by endothelialization. The regional analysis described in this work provides local information that could potentially be used to establish thresholds for thrombosis initiation and develop quantitative predictive models of local aneurysm occlusion. The current study has a number of limitations that should be considered when interpreting the results. Unlike human cerebral arteries, flow reversals were observed in the parent arteries of our animal models. This could affect the values of the time-averaged hemodynamic quantities investigated. Alignment of follow-up and pretreatment vascular models was performed manually. This could result in inaccuracies of the quantification of the patent and occluded regions of the aneurysm. Our analysis did not distinguish between voxels that were close to the wall or near the middle of the aneurysm, which may have different hemodynamic environments. Likewise, the changes in the hemodynamic environment as the aneurysm progressively thromboses were not taken into account, only the hemodynamic conditions immediately posttreatment were studied. More refined models incorporating thrombus formation and proximity to the wall should be considered in future studies. The virtual deployment, the estimation of flow conditions from Doppler sonography images, as well as the many assumptions and approximations made during the CFD modeling process may affect the hemodynamic results. However, previous studies have indicated that these models could reproduce in vivo velocity measurements and observed flow characteristics.21 Finally, the number of animal models included in the study is limited. Studies with a larger number of subjects are needed to confirm the trends observed in this work and to achieve stronger statistical significance. Despite these limitations, this study confirms that local hemodynamic conditions created immediately after flow diverter implantation play an important role in determining whether different regions of the aneurysm will thrombose and occlude or remain patent. The next set of studies should focus on the connection between hemodynamic conditions and biologic processes responsible for not only thrombosis and aneurysm occlusion but also endothelial tissue growth and parent artery reconstruction. ## Conclusions Aneurysms completely occluded at follow-up had statistically smaller necks and longer mean transit times than aneurysms that remained patent. Posttreatment velocity, vorticity, and shear rate were significantly larger in regions that remained patent (ie, had larger “flow activity”) than in regions that occluded. Machine-learning models based on local hemodynamic variables were capable of predicting local occlusion with good precision, especially away from the neck. This suggests that the dominant healing mechanism at the dome and body of the aneurysm is related to flow-induced thrombosis, while near the neck, other processes such as flow-modulated endothelialization may also play a fundamental role. ## Acknowledgments We thank Covidien Inc for generously providing flow diverters. ## Footnotes * Disclosures: Juan R. Cebral—*RELATED*: *Grant*: NIH,* *Comments*: research grant, *UNRELATED*: *Grants/Grants Pending*: National Institutes of Health,* Philips Healthcare,* Fernando Mut—*RELATED*: *Grant*: National Institutes of Health.* David F. Kallmes—*RELATED*: *Grant*: ev3,* *Comments*: supplied devices for the study, *UNRELATED*: *Consultancy*: ev3,* Medtronic,* *Comments*: Consultation for design and implementation of clinical trials, *Grants/Grants Pending*: ev3,* Sequent,* SurModics,* MicroVention,* Cordis,* NeuroSigmas,* *Comments*: preclinical research and clinical trials implementation, *Royalties*: University of Virginia Patent Foundation, *Comments*: Spine fusion, *Travel/Accommodations/Meeting Expenses Unrelated to Activities Listed*: ev3,* *Comments*: participation in proctoring activities. *Money paid to the institution. * This work was supported by National Institutes of Health grant NS076491. Indicates open access to non-subscribers at [www.ajnr.org](http://www.ajnr.org) ## REFERENCES 1. 1. Lylyk P, Miranda C, Ceratto R, et al. Curative endovascular reconstruction of cerebral aneurysms with the Pipeline embolization device: the Buenos Aires experience. Neurosurgery 2009;64:632–642 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1227/01.NEU.0000339109.98070.65&link_type=DOI) [PubMed](http://www.ajnr.org/lookup/external-ref?access_num=19349825&link_type=MED&atom=%2Fajnr%2F35%2F8%2F1567.atom) [Web of Science](http://www.ajnr.org/lookup/external-ref?access_num=000264783500018&link_type=ISI) 2. 2. Nelson PK, Lylyk P, Szikora I, et al. The Pipeline embolization device for the intracranial treatment of aneurysms trial. AJNR. AJNR Am J Neuroradiol 2011;32:34–40 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo3OiIzMi8xLzM0IjtzOjQ6ImF0b20iO3M6MjA6Ii9ham5yLzM1LzgvMTU2Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. 3. Seong J, Wakhloo AK, Lieber BB. In vitro evaluation of flow divertors in an elastase-induced saccular aneurysm model in rabbit. J Biomech Eng 2007;129:863–72 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1115/1.2800787&link_type=DOI) [PubMed](http://www.ajnr.org/lookup/external-ref?access_num=18067390&link_type=MED&atom=%2Fajnr%2F35%2F8%2F1567.atom) 4. 4. Kulcsár Z, Houdart E, Bonafe A, et al. Intra-aneurysmal thrombosis as a possible cause of delayed aneurysm rupture after flow-diversion treatment. AJNR Am J Neuroradiol 2011;32:20–25 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo3OiIzMi8xLzIwIjtzOjQ6ImF0b20iO3M6MjA6Ii9ham5yLzM1LzgvMTU2Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. 5. Kallmes D, Ding YH, Dai D, et al. A new endoluminal, flow disrupting device for treatment of saccular aneurysms. Stroke 2007;38:2346–52 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToic3Ryb2tlYWhhIjtzOjU6InJlc2lkIjtzOjk6IjM4LzgvMjM0NiI7czo0OiJhdG9tIjtzOjIwOiIvYWpuci8zNS84LzE1NjcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 6. 6. Kallmes D, Ding YH, Dai D, et al. Low-disrupting device for treatment of saccular aneurysms. AJNR Am J Neuroradiol 2009;30:1153–58 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo5OiIzMC82LzExNTMiO3M6NDoiYXRvbSI7czoyMDoiL2FqbnIvMzUvOC8xNTY3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 7. 7. Szikora I, Berentei Z, Kulcsar Z, et al. Treatment of intracranial aneurysms by functional reconstruction of the parent artery: the Budapest experience with the Pipeline embolization device. AJNR Am J Neuroradiol 2010;31:1139–47 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo5OiIzMS82LzExMzkiO3M6NDoiYXRvbSI7czoyMDoiL2FqbnIvMzUvOC8xNTY3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 8. 8. Altes TA, Cloft HJ, Short JG, et al. 1999 ARRS Executive Council Award: creation of saccular aneurysms in the rabbit—a model suitable for testing endovascular devices. AJR Am J Roentgenol 2000;174:349–54 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.2214/ajr.174.2.1740349&link_type=DOI) [PubMed](http://www.ajnr.org/lookup/external-ref?access_num=10658703&link_type=MED&atom=%2Fajnr%2F35%2F8%2F1567.atom) [Web of Science](http://www.ajnr.org/lookup/external-ref?access_num=000084885000013&link_type=ISI) 9. 9. Cebral JR, Castro MA, Appanaboyina S, et al. Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity. IEEE Trans Med Imaging 2005;24:457–67 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1109/TMI.2005.844159&link_type=DOI) [PubMed](http://www.ajnr.org/lookup/external-ref?access_num=15822804&link_type=MED&atom=%2Fajnr%2F35%2F8%2F1567.atom) [Web of Science](http://www.ajnr.org/lookup/external-ref?access_num=000228051200005&link_type=ISI) 10. 10. Mut F, Cebral JR. Effects of flow-diverting device oversizing on hemodynamics alteration in cerebral aneurysms. AJNR Am J Neuroradiol 2012;33:2010–16 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czoxMDoiMzMvMTAvMjAxMCI7czo0OiJhdG9tIjtzOjIwOiIvYWpuci8zNS84LzE1NjcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 11. 11. Kundu PK, Cohen LM, Dowling DL. Fluid Mechanics. 5th ed. New York: Academic Press (Elsevier); 2011 12. 12. Mut F., Aubry R, Löhner R, et al. Fast numerical solutions in patient-specific blood flows in 3D arterial systems. Int J Num Method Biomed Eng 2010;26:73–85 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1002/cnm.1235&link_type=DOI) 13. 13. Appanaboyina S, Mut F, Löhne R, et al. Computational fluid dynamics of stented intracranial aneurysms using adaptive embedded unstructured grids. Int J Numer Methods Fluids 2008;57:457–93 14. 14. Appanaboyina S, Mut F, Löhne R, et al. Simulation of intracranial aneurysm stenting: techniques and challenges. Comput Methods Appl Mech Eng 2009;198:3567–82 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1016/j.cma.2009.01.017&link_type=DOI) [Web of Science](http://www.ajnr.org/lookup/external-ref?access_num=000270702500006&link_type=ISI) 15. 15. Mut F, Löhner R, Chien A, et al. Computational hemodynamics framework for the analysis of cerebral aneurysms. Int J Numer Meth Biomed Eng 2011;27:822–39 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1002/cnm.1424&link_type=DOI) 16. 16. Fab RE, Chang KW, Hsieh CJ, et al. LIBLINEAR: a library for large linear classification. J Machine Learning Res 2008;9:1871–74 17. 17. Hall M, Frank E, Holmes G, et al. The WEKA data mining software: an update. SIGKDD Explorations 2009;11:10–18 18. 18. Pereira VM, Bonnefous O, Ouared R, et al. A DSA-based method using contrast-motion estimation for the assessment of the intra-aneurysmal flow changes induced by flow-diverter stents. AJNR Am J Neuroradiol 2013;34:805–15 19. 19. Zhang Y, Chong W, Qian Y. Investigation of intracranial aneurysm hemodynamics following flow diverter stent treatment. Med Eng Phys 2013;35:608–15 [CrossRef](http://www.ajnr.org/lookup/external-ref?access_num=10.1016/j.medengphy.2012.07.005&link_type=DOI) [PubMed](http://www.ajnr.org/lookup/external-ref?access_num=22884174&link_type=MED&atom=%2Fajnr%2F35%2F8%2F1567.atom) 20. 20. Huang Q, Xu J, Cheng J, et al. Hemodynamic changes by flow diverters in rabbit aneurysm models: a computational fluid dynamic study based on micro-computed tomography reconstruction. Stroke 2013;44:1936–41 [Abstract/FREE Full Text](http://www.ajnr.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToic3Ryb2tlYWhhIjtzOjU6InJlc2lkIjtzOjk6IjQ0LzcvMTkzNiI7czo0OiJhdG9tIjtzOjIwOiIvYWpuci8zNS84LzE1NjcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 21. 21. Cebral JR, Mut F, Raschi M, et al. Flow diversion in rabbit aneurysm models. In: Proceedings of the American Society of Mechanical Engineering Summer Bioengineering Conference, Sunriver, Oregon; June 26–29, 2013 * Received September 16, 2013. * Accepted after revision January 9, 2014. * © 2014 by American Journal of Neuroradiology [1]: /embed/graphic-2.gif