Effect of flow diverter stent malposition on intracranial aneurysm hemodynamics-An experimental framework using stereoscopic particle image velocimetry.

BACKGROUND: Flow-diverting stents are increasingly used for the minimally-invasive treatment of intracranial aneurysms. However, a correct positioning of such devices can be challenging due to varying vessel diameters as well as the complex anatomy of the neurovasculature. As a consequence, unsuccessful treatment outcomes are increasingly reported requiring an improvement of the understanding of stent-induced flow modification. METHODS: To evaluate the effect of different degrees of flow diverter stent malposition on intra-aneurysmal hemodynamic changes, a controlled hemodynamic configuration was created using an idealized intracranial aneurysms model. Afterwards, four different treatment scenarios were reproduced comprising of 1) the ideal treatment, 2) an insufficient wall apposition in the region of the ostium, 3) a distorted device migrating into the aneurysm sac and 4) an inaccurately deployed stent due to wrong release location. For the assessment of the individual flow modifications, high-resolution stereoscopic particle image velocimetry (PIV) measurements were carried out. RESULTS: The analysis of the precise in-vitro PIV measurements reveals that in all cases a considerable reduction of the cycle-averaged and peak-systolic velocity was obtained. Compared to the untreated aneurysm configuration, the flow reduction ranged from 63% (scenario 4) up to 89% (scenario 3). The ideal treatment reached a reduction of 78%, which is known to be sufficient for a successful therapy. However, inaccurate device positioning leads to increased oscillating flow towards the lateral directions reducing the chances of sufficient thrombus formation. CONCLUSIONS: High-resolution in-vitro PIV measurements enable an accurate quantification of the treatment efficacy for flow-diverting devices. Furthermore, insufficient treatment outcomes can be reproduces allowing for an assessment of intra-aneurysmal hemodynamic changes.

Central moments multiple relaxation time LBM for hemodynamic simulations in intracranial aneurysms: An in-vitro validation study using PIV and PC-MRI.

The lattice Boltzmann method (LBM) has recently emerged as an efficient alternative to classical Navier-Stokes solvers. This is particularly true for hemodynamics in complex geometries. However, in its most basic formulation, i.e. with the so-called single relaxation time (SRT) collision operator, it has been observed to have a limited stability domain in the Courant/Fourier space, strongly constraining the minimum time-step and grid size. The development of improved collision models such as the multiple relaxation time (MRT) operator in central moments space has tremendously widened the stability domain, while allowing to overcome a number of other well-documented artifacts, therefore opening the door for simulations over a wider range of grid and time-step sizes. The present work focuses on implementing and validating a specific collision operator, the central Hermite moments multiple relaxation time model with the full expansion of the equilibrium distribution function, to simulate blood flows in intracranial aneurysms. The study further proceeds with a validation of the numerical model through different test-cases and against experimental measurements obtained via stereoscopic particle image velocimetry (PIV) and phase-contrast magnetic resonance imaging (PC-MRI). For a patient-specific aneurysm both PIV and PC-MRI agree fairly well with the simulation. Finally, low-resolution simulations were shown to be able to capture blood flow information with sufficient accuracy, as demonstrated through both qualitative and quantitative analysis of the flow field while leading to strongly reduced computation times. For instance in the case of the patient-specific configuration, increasing the grid-size by a factor of two led to a reduction of computation time by a factor of 14 with very good similarity indices still ranging from 0.83 to 0.88.

Transient flow prediction in an idealized aneurysm geometry using data assimilation.

Hemodynamic simulations are restricted by modeling assumptions and uncertain initial and boundary conditions, whereas Phase-Contrast Magnetic Resonance Imaging (PC-MRI) data is affected by measurement noise and artifacts. To overcome the limitations of both techniques, the current study uses a Localization Ensemble Transform Kalman Filter (LETKF) to fully incorporate noisy, low-resolution Phase-Contrast MRI data into an ensemble of high-resolution numerical simulations. The analysis output provides an improved state estimate of the three-dimensional blood flow field in an intracranial aneurysm model. Benchmark measurements are carried out in a silicone phantom model of an idealized aneurysm under pulsatile inflow conditions. Validation is ensured with high-resolution Particle Imaging Velocimetry (PIV) obtained in the symmetry plane of the same geometry. Two data assimilation approaches are introduced, which differ in their way to propagate the ensemble members in time. In both cases the velocity noise is significantly reduced over the whole cardiac cycle. Quantitative and qualitative results indicate an improvement of the flow field prediction in comparison to the raw measurement data. Although biased measurement data reveal a systematic deviation from the truth, the LETKF is able to account for stochastically distributed errors. Through the implementation of the data assimilation step, physical constraints are introduced into the raw measurement data. The resulting, realistic high-resolution flow field can be readily used to assess further patient-specific parameters in addition to the velocity distribution, such as wall shear stress or pressure.

Comparison of intracranial aneurysm flow quantification techniques: standard PIV vs stereoscopic PIV vs tomographic PIV vs phase-contrast MRI vs CFD.

Image-based hemodynamic simulations to assess the rupture risk or improve the treatment planning of intracranial aneurysms have become popular recently. However, due to strong modeling assumptions and limitations, the acceptance of numerical approaches remains limited. Therefore, validation using experimental methods is mandatory.In this study, a unique compilation of four in-vitro flow measurements (three particle image velocimetry approaches using a standard (PIV), stereoscopic (sPIV), and tomographic (tPIV) setup, as well as a phase-contrast magnetic resonance imaging (PC-MRI) measurement) were compared with a computational fluid dynamics (CFD) simulation. This was carried out in a patient-specific silicone phantom model of an internal carotid artery aneurysm under steady flow conditions. To evaluate differences between each technique, a similarity index (SI) with respect to the velocity vectors and the average velocity magnitude differences among all involved modalities were computed.The qualitative comparison reveals that all techniques are able to provide a reasonable description of the global flow structures. High quantitative agreement in terms of SI and velocity magnitude differences was found between all PIV methods and CFD. However, quantitative differences were observed between PC-MRI and the other techniques. Deeper analysis revealed that the limited resolution of the PC-MRI technique is a major contributor to the experienced differences and leads to a systematic underestimation of overall velocity magnitude levels inside the vessel. This confirms the necessity of using highly resolving flow measurement techniques, such as PIV, in an in-vitro environment to individually verify the validity of the numerically obtained hemodynamic results.

Tomographic particle image velocimetry for the validation of hemodynamic simulations in an intracranial aneurysm.

Image-based blood flow simulations can provide detailed hemodynamic information in diseased vessels such as intracranial aneurysms. However, validation is essential to evaluate the accuracy of these computations and further improve their acceptance among physicians. In this regard, tomographic particle image velocimetry was used to measure the flow characteristics in a patient specific aneurysm phantom model. Additionally, computational fluid dynamics (CFD) simulations were carried out using a well accepted commercial software package and a clinical research prototype, respectively. The comparison between in-vitro measurement and in-silico computations reveals a good qualitative agreement. Further, computations based on classical CFD agreed well with results from a clinical research prototype. Hence, the results of this study demonstrate the usability of numerical methods to obtain realistic blood flow predictions in a clinical context.

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.