Philipp Berg

Variability of Computational Fluid Dynamics Solutions for Pressure and Flow in a Giant Aneurysm: The ASME 2012 Summer Bioengineering Conference CFD Challenge

Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

Cerebral Blood Flow in a Healthy Circle of Willis and Two Intracranial Aneurysms: Computational Fluid Dynamics Versus Four-Dimensional Phase-Contrast Magnetic Resonance Imaging

Computational fluid dynamics (CFD) opens up multiple opportunities to investigate the hemodynamics of the human vascular system. However, due to numerous assumptions the acceptance of CFD among physicians is still limited in practice and validation through comparison is mandatory. Time-dependent quantitative phase-contrast magnetic resonance imaging PC-MRI measurements in a healthy volunteer and two intracranial aneurysms were carried out at 3 and 7 Tesla. Based on the acquired images, three-dimensional (3D) models of the aneurysms were reconstructed and used for the numerical simulations. Flow information from the MR measurements were applied as boundary conditions. The four-dimensional (4D) velocity fields obtained by CFD and MRI were qualitatively as well as quantitatively compared including cut planes and vector analyses. For all cases a high similarity of the velocity patterns was observed. Additionally, the quantitative analysis revealed a good agreement between CFD and MRI. Deviations were caused by minor differences between the reconstructed vessel models and the actual lumen. The comparisons between diastole and systole indicate that relative differences between MRI and CFD are intensified with increasing velocity. The findings of this study lead to the conclusion that CFD and MRI agree well in predicting intracranial velocities when realistic geometries and boundary conditions are provided. Due to the considerably higher temporal and spatial resolution of CFD compared to MRI, complex flow patterns can be further investigated in order to evaluate their role with respect to aneurysm formation or rupture. Nevertheless, special care is required regarding the vessel reconstruction since the geometry has a major impact on the subsequent numerical results.

The Computational Fluid Dynamics Rupture Challenge 2013—Phase I: Prediction of Rupture Status in Intracranial Aneurysms

BACKGROUND AND PURPOSE: Rupture risk assessment for intracranial aneurysms remains challenging, and risk factors, including wall shear stress, are discussed controversially. The primary purpose of the presented challenge was to determine how consistently aneurysm rupture status and rupture site could be identified on the basis of computational fluid dynamics. MATERIALS AND METHODS: Two geometrically similar MCA aneurysms were selected, 1 ruptured, 1 unruptured. Participating computational fluid dynamics groups were blinded as to which case was ruptured. Participants were provided with digitally segmented lumen geometries and, for this phase of the challenge, were free to choose their own flow rates, blood rheologies, and so forth. Participants were asked to report which case had ruptured and the likely site of rupture. In parallel, lumen geometries were provided to a group of neurosurgeons for their predictions of rupture status and site. RESULTS: Of 26 participating computational fluid dynamics groups, 21 (81%) correctly identified the ruptured case. Although the known rupture site was associated with low and oscillatory wall shear stress, most groups identified other sites, some of which also experienced low and oscillatory shear. Of the 43 participating neurosurgeons, 39 (91%) identified the ruptured case. None correctly identified the rupture site. CONCLUSIONS: Geometric or hemodynamic considerations favor identification of rupture status; however, retrospective identification of the rupture site remains a challenge for both engineers and clinicians. A more precise understanding of the hemodynamic factors involved in aneurysm wall pathology is likely required for computational fluid dynamics to add value to current clinical decision-making regarding rupture risk.

Endothelialization of over- and undersized flow-diverter stents at covered vessel side branches: An in vivo and in silico study.

Although flow-diverting devices are promising treatment options for intracranial aneurysms, jailed side branches might occlude leading to insufficient blood supply. Especially differences in the local stent strut compression may have a drastic influence on subsequent endothelialization. To investigate the outcome of different treatment scenarios, over- and undersized stent deployments were realized experimentally and computationally. Two Pipeline Embolization Devices were placed in the right common carotid artery of large white swine, crossing the right ascending pharyngeal artery. DSA and PC-MRI measurements were acquired pre- and post-stenting and after three months. To evaluate the stent strut endothelialization and the corresponding ostium patency, the swine were sacrificed and scanning electron microscopy measurements were carried out. A more detailed analysis of the near-stent hemodynamics was enabled by a realistic virtual stenting in combination with highly resolved Computational Fluid Dynamics simulations using case-specific boundary conditions. The oversizing resulted in an elongated stent deployment with more open stent pores, while for the undersized case a shorter deployment with more condensed pores was present. In consequence, the side branch of the first case remained patent after three months and the latter almost fully occluded. The virtual investigation confirmed the experimental findings by identifying differences between the individual velocities as well as stent shear stresses at the distal part of the ostia. The choice of flow-diverting device and the subsequent deployment strategy strongly influences the patency of jailed side branches. Therefore, careful treatment planning is required, to guarantee sufficient blood supply in the brain territories supplied those branches.

An automatic CFD-based flow diverter optimization principle for patient-specific intracranial aneurysms

The optimal treatment of intracranial aneurysms using flow diverting devices is a fundamental issue for neuroradiologists as well as neurosurgeons. Due to highly irregular manifold aneurysm shapes and locations, the choice of the stent and the patient-specific deployment strategy can be a very difficult decision. To support the therapy planning, a new method is introduced that combines a three-dimensional CFD-based optimization with a realistic deployment of a virtual flow diverting stent for a given aneurysm. To demonstrate the feasibility of this method, it was applied to a patient-specific intracranial giant aneurysm that was successfully treated using a commercial flow diverter. Eight treatment scenarios with different local compressions were considered in a fully automated simulation loop. The impact on the corresponding blood flow behavior was evaluated qualitatively as well as quantitatively, and the optimal configuration for this specific case was identified. The virtual deployment of an uncompressed flow diverter reduced the inflow into the aneurysm by 24.4% compared to the untreated case. Depending on the positioning of the local stent compression below the ostium, blood flow reduction could vary between 27.3% and 33.4%. Therefore, a broad range of potential treatment outcomes was identified, illustrating the variability of a given flow diverter deployment in general. This method represents a proof of concept to automatically identify the optimal treatment for a patient in a virtual study under certain assumptions. Hence, it contributes to the improvement of virtual stenting for intracranial aneurysms and can support physicians during therapy planning in the future.

Recommendations for Accurate Numerical Blood Flow Simulations of Stented Intracranial Aneurysms

Abstract The number of scientific publications dealing with stented intracranial aneurysms is rapidly increasing. Powerful computational facilities are now available; an accurate computational modeling of hemodynamics in patient-specific configurations is, however, still being sought. Furthermore, there is still no general agreement on the quantities that should be computed and on the most adequate analysis for intervention support. In this article, the accurate representation of patient geometry is first discussed, involving successive improvements. Concerning the second step, the mesh required for the numerical simulation is especially challenging when deploying a stent with very fine wire structures. Third, the description of the fluid properties is a major challenge. Finally, a founded quantitative analysis of the simulation results is obviously needed to support interventional decisions. In the present work, an attempt has been made to review the most important steps for a high-quality computational fluid dynamics computation of virtually stented intracranial aneurysms. In consequence, this leads to concrete recommendations, whereby the obtained results are not discussed for their medical relevance but for the evaluation of their quality. This investigation might hopefully be helpful for further studies considering stent deployment in patient-specific geometries, in particular regarding the generation of the most appropriate computational model.

AmniVis - a system for qualitative exploration of near-wall hemodynamics in cerebral aneurysms

The qualitative exploration of near‐wall hemodynamics in cerebral aneurysms provides important insights for risk assessment. For instance, a direct relation between complex flow patterns and aneurysm formation could be observed. Due to the high complexity of the underlying time‐dependent flow data, the exploration is challenging, in particular for medical researchers not familiar with such data. We present the AmniVis‐Explorer, a system that is designed for the preparation of a qualitative medical study. The provided features were developed in close collaboration with medical researchers involved in the study. This comprises methods for a purposeful selection of surface regions of interest and a novel approach to provide a 2D overview of flow patterns that are represented by streamlines at these regions. Furthermore, we present a specialized interface that supports binary classification of patterns and temporal exploration as well as methods for selection, highlighting and automatic 3D navigation to particular patterns. Based on eight representative datasets, we conducted informal interviews with two bord‐certified radiologists and a flow expert to evaluate the system. It was confirmed that the AmniVis‐Explorer allows for an easy selection, qualitative exploration and classification of near‐wall flow patterns that are represented by streamlines.

Role of terminal and anastomotic circulation in the patency of arteries jailed by flow-diverting stents: animal flow model evaluation and preliminary results

OBJECTIVE The authors describe herein the creation of an animal model capable of producing quantifiable data regarding blood flow rate and velocity modifications in terminal and anastomotic types of cerebrofacial circulation. They also present the preliminary results of a translational study aimed at investigating the role of terminal and anastomotic types of circulation in arterial branches jailed by flow-diverting stents as factors contributing to arterial patency or occlusion. METHODS Two Large White swine were used to validate a terminal-type arterial model at the level of the right ascending pharyngeal artery (APhA), created exclusively by endovascular means. Subsequently 4 Large White swine, allocated to 2 groups corresponding to the presence (Group B) or absence (Group A) of terminal-type flow modification, underwent placement of flow-diverting stents. Blood flow rates and velocities were quantified using a dedicated time-resolved 3D phase-contrast MRA sequence before and after stenting. Three months after stent placement, the stented arteries were evaluated with digital subtraction angiography (DSA) and scanning electron microscopy (SEM). Patent (circulating) ostia quantification was performed on the SEM images. RESULTS Terminal-type flow modification was feasible; an increase of 75.8% in mean blood velocities was observed in the right APhAs. The mean blood flow rate for Group A was 0.31 ± 0.19 ml/sec (95% CI -1.39 to 2.01) before stenting and 0.21 ± 0.07 ml/sec (95% CI -0.45 to 0.87) after stenting. The mean blood flow rate for Group B was 0.87 ± 0.32 ml/sec (95% CI -1.98 to 3.73) before stenting and 0.76 ± 0.13 ml/sec (95% CI -0.41 to 1.93) after stenting. Mean flow rates after stenting showed a statistically significant difference between Groups A and B (Welch test). Mean and maximal blood velocities were reduced in Group A cases and did not decrease in Group B cases. Control DSA and SEM findings showed near occlusion of the jailed APhAs in both cases of anastomotic circulation (mean patent ostium surface 32,776 μm2) and patency in both cases of terminal-type circulation (mean patent ostium surface 422,334 μm2). CONCLUSIONS Terminal-type arterial modification in swine APhAs is feasible. Sufficient data were acquired to perform an a priori analysis for further research. Flow diversion at the level of the APhA ostium resulted in significant stenosis in cases of anastomotic circulation, while sufficient patency was observed in terminal-type circulation.

Does the DSA reconstruction kernel affect hemodynamic predictions in intracranial aneurysms? An analysis of geometry and blood flow variations

Background Computational fluid dynamics (CFD) blood flow predictions in intracranial aneurysms promise great potential to reveal patient-specific flow structures. Since the workflow from image acquisition to the final result includes various processing steps, quantifications of the individual introduced potential error sources are required. Methods Three-dimensional (3D) reconstruction of the acquired imaging data as input to 3D model generation was evaluated. Six different reconstruction modes for 3D digital subtraction angiography (DSA) acquisitions were applied to eight patient-specific aneurysms. Segmentations were extracted to compare the 3D luminal surfaces. Time-dependent CFD simulations were carried out in all 48 configurations to assess the velocity and wall shear stress (WSS) variability due to the choice of reconstruction kernel. Results All kernels yielded good segmentation agreement in the parent artery; deviations of the luminal surface were present at the aneurysm neck (up to 34.18%) and in distal or perforating arteries. Observations included pseudostenoses as well as noisy surfaces, depending on the selected reconstruction kernel. Consequently, the hemodynamic predictions show a mean SD of 11.09% for the aneurysm neck inflow rate, 5.07% for the centerline-based velocity magnitude, and 17.83%/9.53% for the mean/max aneurysmal WSS, respectively. In particular, vessel sections distal to the aneurysms yielded stronger variations of the CFD values. Conclusions The choice of reconstruction kernel for DSA data influences the segmentation result, especially for small arteries. Therefore, if precise morphology measurements or blood flow descriptions are desired, a specific reconstruction setting is required. Furthermore, research groups should be encouraged to denominate the kernel types used in future hemodynamic studies.

Role of terminal and anastomotic circulation in the patency of arteries jailed by flow-diverting stents: from hemodynamic changes to ostia surface modifications

OBJECTIVE The outcome for jailing arterial branches that emerge near intracranial aneurysms during flow-diverting stent (FDS) deployment remains controversial. In this animal study, the authors aimed to elucidate the role of collateral supply with regard to the hemodynamic changes and neointimal modifications that occur from jailing arteries with FDSs. To serve this purpose, the authors sought to quantify 1) the hemodynamic changes that occur at the jailed arterial branches immediately after stent placement and 2) the ostia surface values at 3 months after stenting; both parameters were investigated in the presence or absence of collateral arterial flow. METHODS After an a priori power analysis, 2 groups (Group A and Group B) were created according to an animal flow model for terminal and anastomotic arterial circulation; each group contained 7 Large White swine. Group A animals possessed an anastomotic-type arterial configuration to supply the territory of the right ascending pharyngeal artery (APhA), while Group B animals possessed a terminal-type arterial configuration to supply the right APhA territory. Subsequently, all animals underwent FDS placement, thereby jailing the right APhAs. Mean flow rates and velocities inside the jailed branches were quantified using time-resolved 3D phase-contrast MR angiography before and after stenting. Three months after stent placement, the jailed ostia surface values were quantified on scanning electron micrographs. The data were analyzed using descriptive statistics and group comparisons with parametric and nonparametric tests. RESULTS The endovascular procedures were feasible, and there were no findings of in situ thrombus formation on postprocedural optical coherence tomography or ischemia on postprocedural diffusion-weighted imaging. In Group A, the mean flow rate values at the jailed right APhAs were reduced immediately following stent placement as compared with values obtained before stent placement (p = 0.02, power: 0.8). In contrast, the mean poststenting flow rates for Group B remained similar to those obtained before stent placement. Three months after stent placement, the mean ostia surface values were significantly higher for Group B (527,911 ± 306,229 μm2) than for Group A (89,329 ± 59,762 μm2; p < 0.01, power: 1.00), even though the initial dimensions of the jailed ostia were similar between groups. A statistically significant correlation was found between groups (A or B), mean flow rates after stent placement, and ostia surface values at 3 months. CONCLUSIONS When an important collateral supply was present, the jailing of side arteries with flow diverters resulted in an immediate and significant reduction in the flow rate inside these arteries as compared with the prestenting values. In contrast, when competitive flow was absent, jailing did not result in significant flow rate reductions inside the jailed arteries. Ostium surface values at 3 months after stent placement were significantly higher in the terminal group of jailed arteries (Group B) than in the anastomotic group (Group A) and strongly correlated with poststenting reductions in the velocity value.