The Computational Fluid Dynamics Rupture Challenge 2013--Phase II: Variability of Hemodynamic Simulations in Two Intracranial Aneurysms.

With the increased availability of computational resources, the past decade has seen a rise in the use of computational fluid dynamics (CFD) for medical applications. There has been an increase in the application of CFD to attempt to predict the rupture of intracranial aneurysms, however, while many hemodynamic parameters can be obtained from these computations, to date, no consistent methodology for the prediction of the rupture has been identified. One particular challenge to CFD is that many factors contribute to its accuracy; the mesh resolution and spatial/temporal discretization can alone contribute to a variation in accuracy. This failure to identify the importance of these factors and identify a methodology for the prediction of ruptures has limited the acceptance of CFD among physicians for rupture prediction. The International CFD Rupture Challenge 2013 seeks to comment on the sensitivity of these various CFD assumptions to predict the rupture by undertaking a comparison of the rupture and blood-flow predictions from a wide range of independent participants utilizing a range of CFD approaches. Twenty-six groups from 15 countries took part in the challenge. Participants were provided with surface models of two intracranial aneurysms and asked to carry out the corresponding hemodynamics simulations, free to choose their own mesh, solver, and temporal discretization. They were requested to submit velocity and pressure predictions along the centerline and on specified planes. The first phase of the challenge, described in a separate paper, was aimed at predicting which of the two aneurysms had previously ruptured and where the rupture site was located. The second phase, described in this paper, aims to assess the variability of the solutions and the sensitivity to the modeling assumptions. Participants were free to choose boundary conditions in the first phase, whereas they were prescribed in the second phase but all other CFD modeling parameters were not prescribed. In order to compare the computational results of one representative group with experimental results, steady-flow measurements using particle image velocimetry (PIV) were carried out in a silicone model of one of the provided aneurysms. Approximately 80% of the participating groups generated similar results. Both velocity and pressure computations were in good agreement with each other for cycle-averaged and peak-systolic predictions. Most apparent "outliers" (results that stand out of the collective) were observed to have underestimated velocity levels compared to the majority of solutions, but nevertheless identified comparable flow structures. In only two cases, the results deviate by over 35% from the mean solution of all the participants. Results of steady CFD simulations of the representative group and PIV experiments were in good agreement. The study demonstrated that while a range of numerical schemes, mesh resolution, and solvers was used, similar flow predictions were observed in the majority of cases. To further validate the computational results, it is suggested that time-dependent measurements should be conducted in the future. However, it is recognized that this study does not include the biological aspects of the aneurysm, which needs to be considered to be able to more precisely identify the specific rupture risk of an intracranial aneurysm.

Rupture point is associated with divergent hemodynamics in intracranial aneurysms.

Background: Understanding the risk factors leading to intracranial aneurysm (IA) rupture have still not been fully clarified. They are vital for proper medical guidance of patients harboring unruptured IAs. Clarifying the hemodynamics associated with the point of rupture could help could provide useful information about some of the risk factors. Thus far, few studies have studied this issue with often diverging conclusions. Methods: We identified a point of rupture in patients operated for an IAs during surgery, using a combination of preoperative computed tomography (CT) and computed tomography angiography (CTA). Hemodynamic parameters were calculated both for the aneurysm sac as a whole and the point of rupture. In two cases, the results of CFD were compared with those of the experiment using particle image velocimetry (PIV). Results: We were able to identify 6 aneurysms with a well-demarcated point of rupture. In four aneurysms, the rupture point was near the vortex with low wall shear stress (WSS) and high oscillatory shear index (OSI). In one case, the rupture point was in the flow jet with high WSS. In the last case, the rupture point was in the significant bleb and no specific hemodynamic parameters were found. The CFD results were verified in the PIV part of the study. Conclusion: Our study shows that different hemodynamic scenarios are associated with the site of IA rupture. The numerical simulations were confirmed by laboratory models. This study further supports the hypothesis that various pathological pathways may lead to aneurysm wall damage resulting in its rupture.

Morphological and Hemodynamic Changes during Cerebral Aneurysm Growth.

Computational fluid dynamics (CFD) has grown as a tool to help understand the hemodynamic properties related to the rupture of cerebral aneurysms. Few of these studies deal specifically with aneurysm growth and most only use a single time instance within the aneurysm growth history. The present retrospective study investigated four patient-specific aneurysms, once at initial diagnosis and then at follow-up, to analyze hemodynamic and morphological changes. Aneurysm geometries were segmented via the medical image processing software Mimics. The geometries were meshed and a computational fluid dynamics (CFD) analysis was performed using ANSYS. Results showed that major geometry bulk growth occurred in areas of low wall shear stress (WSS). Wall shape remodeling near neck impingement regions occurred in areas with large gradients of WSS and oscillatory shear index. This study found that growth occurred in areas where low WSS was accompanied by high velocity gradients between the aneurysm wall and large swirling flow structures. A new finding was that all cases showed an increase in kinetic energy from the first time point to the second, and this change in kinetic energy seems correlated to the change in aneurysm volume.