The Computational Fluid Dynamics Rupture Challenge 2013—Phase I: prediction of rupture status in intracranial aneurysms

Influence of circle of Willis modeling on hemodynamic parameters in anterior communicating artery aneurysms and recommendations for model selection

Computational fluid dynamics (CFD) has been utilized to calculate hemodynamic parameters in anterior communicating artery aneurysm (AComA), which is located at a junction between left and right A1 and A2 segments. However, complete or half circle of Willis (CoW) models are used indiscriminately. This study aims to suggest recommendations for determining suitable CoW model. Five patient-specific CoW models with AComA were used, and each model was divided into complete, left-half, and right-half models. After validating the CFD using a flow experiment, the hemodynamic parameters and flow patterns in five AComAs were compared. In four out of five cases, inflow from one A1 side had a dominant influence on the AComA, while both left and right A1 sides affected the AComA in the remaining case. Also, the average difference in time-averaged wall shear stress between the complete and half models for four cases was 4.6%, but it was 62% in the other case. The differences in the vascular resistances of left and right A1 and A2 segments greatly influenced the flow patterns in the AComA. These results may help to enhance clinicians' understanding of blood flow in the brain, leading to improvements in diagnosis and treatment of cerebral aneurysms.

Analysis of Intracranial Aneurysm Haemodynamics Altered by Wall Movement

Computational fluid dynamics is intensively used to deepen our understanding of aneurysm growth and rupture in an attempt to support physicians during therapy planning. Numerous studies assumed fully rigid vessel walls in their simulations, whose sole haemodynamics may fail to provide a satisfactory criterion for rupture risk assessment. Moreover, direct in vivo observations of intracranial aneurysm pulsation were recently reported, encouraging the development of fluid-structure interaction for their modelling and for new assessments. In this work, we describe a new fluid-structure interaction functional setting for the careful evaluation of different aneurysm shapes. The configurations consist of three real aneurysm domes positioned on a toroidal channel. All geometric features, employed meshes, flow quantities, comparisons with the rigid wall model and corresponding plots are provided for the sake of reproducibility. The results emphasise the alteration of flow patterns and haemodynamic descriptors when wall deformations were taken into account compared with a standard rigid wall approach, thereby underlining the impact of fluid-structure interaction modelling.

Image-based hemodynamic simulations for intracranial aneurysms: the impact of complex vasculature

PURPOSE

Hemodynamics play an important role in the assessment of intracranial aneurysm (IA) development and rupture risk. The purpose of this study was to examine the impact of complex vasculatures onto the intra-vessel and intra-aneurysmal blood flow.

METHODS

Complex segmentation of a subject-specific, 60-outlet and 3-inlet circle of Willis model captured with 7T magnetic resonance imaging was performed. This model was trimmed to a 10-outlet model version. Two patient-specific IAs were added onto both models yielding two pathological versions, and image-based blood flow simulations of the four resulting cases were carried out. To capture the differences between complex and trimmed model, time-averaged and centerline velocities were compared. The assessment of intra-saccular blood flow within the IAs involved the evaluation of wall shear stresses (WSS) at the IA wall and neck inflow rates (NIR).

RESULTS

Lower flow values are observed in the majority of the complex model. However, at specific locations (left middle cerebral artery 0.5 m/s, left posterior cerebral artery 0.25 m/s), higher flow rates were visible when compared to the trimmed counterpart. Furthermore, at the centerlines the total velocity values reveal differences up to 0.15 m/s. In the IAs, the reduction in the neck inflow rate and WSS in the complex model was observed for the first IA (IA-A δNIRmean =  - 0.07ml/s, PCA.l δWSSmean =  - 0.05 Pa). The second IA featured an increase in the neck inflow rate and WSS (IA-B δNIRmean = 0.04 ml/s, PCA.l δWSSmean = 0.07 Pa).

CONCLUSION

Both the magnitude and shape of the flow distribution vary depending on the model's complexity. The magnitude is primarily influenced by the global vessel model, while the shape is determined by the local structure. Furthermore, intra-aneurysmal flow strongly depends on the location in the vessel tree, emphasizing the need for complex model geometries for realistic hemodynamic assessment and rupture risk analysis.

StrokeNet: An automated approach for segmentation and rupture risk prediction of intracranial aneurysm

Intracranial Aneurysms (IA) present a complex challenge for neurosurgeons as the risks associated with surgical intervention, such as Subarachnoid Hemorrhage (SAH) mortality and morbidity, may outweigh the benefits of aneurysmal occlusion in some cases. Hence, there is a critical need for developing techniques that assist physicians in assessing the risk of aneurysm rupture to determine which aneurysms require treatment. However, a reliable IA rupture risk prediction technique is currently unavailable. To address this issue, this study proposes a novel approach for aneurysm segmentation and multidisciplinary rupture prediction using 2D Digital Subtraction Angiography (DSA) images. The proposed method involves training a fully connected convolutional neural network (CNN) to segment aneurysm regions in DSA images, followed by extracting and fusing different features using a multidisciplinary approach, including deep features, geometrical features, Fourier descriptor, and shear pressure on the aneurysm wall. The proposed method also adopts a fast correlation-based filter approach to drop highly correlated features from the set of fused features. Finally, the selected fused features are passed through a Decision Tree classifier to predict the rupture severity of the associated aneurysm into four classes: Mild, Moderate, Severe, and Critical. The proposed method is evaluated on a newly developed DSA image dataset and on public datasets to assess its generalizability. The system's performance is also evaluated on DSA images annotated by expert neurosurgeons for the rupture risk assessment of the segmented aneurysm. The proposed system outperforms existing state-of-the-art segmentation methods, achieving an 85 % accuracy against annotated DSA images for the risk assessment of aneurysmal rupture.

Association of local solid mechanical, hemodynamic and morphological characteristics with ruptured intracranial aneurysm

The rupture of intracranial aneurysms (IAs) is a complicated phenomenon of which the mechanism is not fully understood. The purpose of this study is to associate local solid mechanical, hemodynamic, and morphological characteristics with rupture regions through statistical means, in an attempt to identify the parameters that are indicative of rupture propensity for IAs. Twenty patient-specific ruptured IA models were reconstructed from digital subtraction angiography (DSA), and applied in the analysis of wall tension, wall shear stress (WSS) and curvature. The precise rupture locations were marked out through intraoperative videos. Pearson correlation analysis was employed to investigate the correlations of these three parameters with patient characteristics and global geometric features. Univariate and multivariate logistic regression analysis were further performed on wall tension, WSS and curvature with regards to rupture and nonrupture regions. Receiver operating characteristic (ROC) analysis defining area under the curve (AUC) was performed on these three parameters. The univariate model of wall tension (AUC, 0.9750), WSS (AUC, 0.9300), curvature (0.8150) and their combined multivariate model (AUC, 0.9875) all present high AUC values. The wall tension, WSS and curvature are acceptable parameters relating to rupture regions. The rupture odd is more sensitive to the wall tension and WSS than curvature. Each logistic model is capable in discriminating ruptures from nonrupture regions, while the multivariate model is the most efficient.

Influence of blood viscosity models and boundary conditions on the computation of hemodynamic parameters in cerebral aneurysms using computational fluid dynamics

BACKGROUND

Computational fluid dynamics (CFD) is widely used to calculate hemodynamic parameters that are known to influence cerebral aneurysms. However, the boundary conditions for CFD are chosen without any specific criteria. Our objective is to establish the recommendations for setting the analysis conditions for CFD analysis of the cerebral aneurysm.

METHOD

The plug and the Womersley flow were the inlet boundary conditions, and zero and pulsatile pressures were the outlet boundary conditions. In addition, the difference in the assumption of viscosity was analyzed with respect to the flow rate. The CFD process used in our research was validated using particle image velocimetry experiment data from Tupin et al.'s work to ensure the accuracy of the simulations.

RESULTS

It was confirmed that if the entrance length was sufficiently secured, the inlet and outlet boundary conditions did not affect the CFD results. In addition, it was observed that the difference in the hemodynamic parameter between Newtonian and non-Newtonian fluid decreased as the flow rate increased. Furthermore, it was confirmed that similar tendencies were evaluated when these recommendations were utilized in the patient-specific cerebral aneurysm models.

CONCLUSIONS

These results may help conduct standardized CFD analyses regardless of the research group.

Hemodynamics of anterior circulation intracranial aneurysms with daughter blebs: investigating the multidirectionality of blood flow fields

Recent advances in diagnostic neuroradiological imaging, allowed the detection of unruptured intracranial aneurysms (IAs). The shape - irregular or multilobular - of the aneurysmal dome, is considered as a possible rupture risk factor, independently of the size, the location and patient medical background. Disturbed blood flow fields in particular is thought to play a key role in IAs progression. However, there is an absence of widely-used hemodynamic indices to quantify the extent of a multi-directional disturbed flow. We simulated blood flow in twelve patient-specific anterior circulation unruptured intracranial aneurysms with daughter blebs utilizing the spectral/hp element framework Nektar++. We simulated three cardiac cycles using a volumetric flow rate waveform while we considered blood as a Newtonian fluid. To investigate the multidirectionality of the blood flow fields, besides the time-averaged wall shear stress (TAWSS), we calculated the oscillatory shear index (OSI), the relative residence time (RRT) and the time-averaged cross flow index (TACFI). Our CFD simulations suggest that in the majority of our vascular models there is a formation of complex intrasaccular flow patterns, resulting to low and highly oscillating WSS, especially in the area of the daughter blebs. The existence of disturbed multi-directional blood flow fields is also evident by the distributions of the RRT and the TACFI. These findings further support the theory that IAs with daughter blebs are linked to a potentially increased rupture risk.

Detection and analysis of cerebral aneurysms based on X-ray rotational angiography - the CADA 2020 challenge

The Cerebral Aneurysm Detection and Analysis (CADA) challenge was organized to support the development and benchmarking of algorithms for detecting, analyzing, and risk assessment of cerebral aneurysms in X-ray rotational angiography (3DRA) images. 109 anonymized 3DRA datasets were provided for training, and 22 additional datasets were used to test the algorithmic solutions. Cerebral aneurysm detection was assessed using the F2 score based on recall and precision, and the fit of the delivered bounding box was assessed using the distance to the aneurysm. The segmentation quality was measured using the Jaccard index and a combination of different surface distance measures. Systematic errors were analyzed using volume correlation and bias. Rupture risk assessment was evaluated using the F2 score. 158 participants from 22 countries registered for the CADA challenge. The U-Net-based detection solutions presented by the community show similar accuracy compared to experts (F2 score 0.92), with a small number of missed aneurysms with diameters smaller than 3.5 mm. In addition, the delineation of these structures, based on U-Net variations, is excellent, with a Jaccard score of 0.92. The rupture risk estimation methods achieved an F2 score of 0.71. The performance of the detection and segmentation solutions is equivalent to that of human experts. The best results are obtained in rupture risk estimation by combining different image-based, morphological, and computational fluid dynamic parameters using machine learning methods. Furthermore, we evaluated the best methods pipeline, from detecting and delineating the vessel dilations to estimating the risk of rupture. The chain of these methods achieves an F2-score of 0.70, which is comparable to applying the risk prediction to the ground-truth delineation (0.71).

Integral Curve Clustering and Simplification for Flow Visualization: A Comparative Evaluation

Unsupervised clustering techniques have been widely applied to flow simulation data to alleviate clutter and occlusion in the resulting visualization. However, there is an absence of systematic guidelines for users to evaluate (both quantitatively and visually) the appropriate clustering technique and similarity measures for streamline and pathline curves. In this work, we provide an overview of a number of prevailing curve clustering techniques. We then perform a comprehensive experimental study to qualitatively and quantitatively compare these clustering techniques coupled with popular similarity measures used in the flow visualization literature. Based on our experimental results, we derive empirical guidelines for selecting the appropriate clustering technique and similarity measure given the requirements of the visualization task. We believe our work will inform the task of generating meaningful reduced representations for large-scale flow data and inspire the continuous investigation of a more refined guidance on clustering technique selection.

The biophysical role of hemodynamics in the pathogenesis of cerebral aneurysm formation and rupture

The pathogenesis of intracranial aneurysms remains complex and multifactorial. While vascular, genetic, and epidemiological factors play a role, nascent aneurysm formation is believed to be induced by hemodynamic forces. Hemodynamic stresses and vascular insults lead to additional aneurysm and vessel remodeling. Advanced imaging techniques allow us to better define the roles of aneurysm and vessel morphology and hemodynamic parameters, such as wall shear stress, oscillatory shear index, and patterns of flow on aneurysm formation, growth, and rupture. While a complete understanding of the interplay between these hemodynamic variables remains elusive, the authors review the efforts that have been made over the past several decades in an attempt to elucidate the physical and biological interactions that govern aneurysm pathophysiology. Furthermore, the current clinical utility of hemodynamics in predicting aneurysm rupture is discussed.

How patient specific are patient-specific computational models of cerebral aneurysms? An overview of sources of error and variability

Computational modeling of cerebral aneurysms, derived from clinical 3D angiography, has become widespread over the past 15 years. While such "image-based" or "patient-specific" models have shown promise for the assessment of rupture risk, much debate remains about their reliability in light of necessary modeling assumptions and incomplete or uncertain model input parameters derived from the clinic. The aims of this review were to walk through the various steps of this so-called patient-specific modeling pipeline and to highlight evidence supporting those steps that we can or cannot rely on. The relative importance of the different sources of error and variability on hemodynamic predictions is summarized, with recommendations to standardize for those that can be avoided and to pay closer attention those to that cannot.

A review on the reliability of hemodynamic modeling in intracranial aneurysms: why computational fluid dynamics alone cannot solve the equation

Computational blood flow modeling in intracranial aneurysms (IAs) has enormous potential for the assessment of highly resolved hemodynamics and derived wall stresses. This results in an improved knowledge in important research fields, such as rupture risk assessment and treatment optimization. However, due to the requirement of assumptions and simplifications, its applicability in a clinical context remains limited.This review article focuses on the main aspects along the interdisciplinary modeling chain and highlights the circumstance that computational fluid dynamics (CFD) simulations are embedded in a multiprocess workflow. These aspects include imaging-related steps, the setup of realistic hemodynamic simulations, and the analysis of multidimensional computational results. To condense the broad knowledge, specific recommendations are provided at the end of each subsection.Overall, various individual substudies exist in the literature that have evaluated relevant technical aspects. In this regard, the importance of precise vessel segmentations for the simulation outcome is emphasized. Furthermore, the accuracy of the computational model strongly depends on the specific research question. Additionally, standardization in the context of flow analysis is required to enable an objective comparison of research findings and to avoid confusion within the medical community. Finally, uncertainty quantification and validation studies should always accompany numerical investigations.In conclusion, this review aims for an improved awareness among physicians regarding potential sources of error in hemodynamic modeling for IAs. Although CFD is a powerful methodology, it cannot provide reliable information, if pre- and postsimulation steps are inaccurately carried out. From this, future studies can be critically evaluated and real benefits can be differentiated from results that have been acquired based on technically inaccurate procedures.

How patient-specific do internal carotid artery inflow rates need to be for computational fluid dynamics of cerebral aneurysms?

BACKGROUND

Computational fluid dynamics (CFD) has become a popular tool for studying 'patient-specific' blood flow dynamics in cerebral aneurysms; however, rarely are the inflow boundary conditions patient-specific. We aimed to test the impact of widespread reliance on generalized inflow rates.

METHODS

Internal carotid artery (ICA) flow rates were measured via 2D cine phase-contrast MRI for 24 patients scheduled for endovascular therapy of an ICA aneurysm. CFD models were constructed from 3D rotational angiography, and pulsatile inflow rates imposed as measured by MRI or estimated using an average older-adult ICA flow waveform shape scaled by a cycle-average flow rate (Q_avg) derived from the patient's ICA cross-sectional area via an assumed inlet velocity.

RESULTS

There was good overall qualitative agreement in the magnitudes and spatial distributions of time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and spectral power index (SPI) using generalized versus patient-specific inflows. Sac-averaged quantities showed moderate to good correlations: R²=0.54 (TAWSS), 0.80 (OSI), and 0.68 (SPI). Using patient-specific Q_avg to scale the generalized waveform shape resulted in near-perfect agreement for TAWSS, and reduced bias, but not scatter, for SPI. Patient-specific waveform had an impact only on OSI correlations, which improved to R²=0.93.

CONCLUSIONS

Aneurysm CFD demonstrates the ability to stratify cases by nominal hemodynamic 'risk' factors when employing an age- and vascular-territory-specific recipe for generalized inflow rates. Q_avg has a greater influence than waveform shape, suggesting some improvement could be achieved by including measurement of patient-specific Q_avg into aneurysm imaging protocols.

Computational Hemodynamic Modeling of Arterial Aneurysms: A Mini-Review

Arterial aneurysms are pathological dilations of blood vessels, which can be of clinical concern due to thrombosis, dissection, or rupture. Aneurysms can form throughout the arterial system, including intracranial, thoracic, abdominal, visceral, peripheral, or coronary arteries. Currently, aneurysm diameter and expansion rates are the most commonly used metrics to assess rupture risk. Surgical or endovascular interventions are clinical treatment options, but are invasive and associated with risk for the patient. For aneurysms in locations where thrombosis is the primary concern, diameter is also used to determine the level of therapeutic anticoagulation, a treatment that increases the possibility of internal bleeding. Since simple diameter is often insufficient to reliably determine rupture and thrombosis risk, computational hemodynamic simulations are being developed to help assess when an intervention is warranted. Created from subject-specific data, computational models have the potential to be used to predict growth, dissection, rupture, and thrombus-formation risk based on hemodynamic parameters, including wall shear stress, oscillatory shear index, residence time, and anomalous blood flow patterns. Generally, endothelial damage and flow stagnation within aneurysms can lead to coagulation, inflammation, and the release of proteases, which alter extracellular matrix composition, increasing risk of rupture. In this review, we highlight recent work that investigates aneurysm geometry, model parameter assumptions, and other specific considerations that influence computational aneurysm simulations. By highlighting modeling validation and verification approaches, we hope to inspire future computational efforts aimed at improving our understanding of aneurysm pathology and treatment risk stratification.

Population-specific modelling of between/within-subject flow variability in the carotid arteries of the elderly

Computational fluid dynamics models are increasingly proposed for assisting the diagnosis and management of vascular diseases. Ideally, patient-specific flow measurements are used to impose flow boundary conditions. When patient-specific flow measurements are unavailable, mean values of flow measurements across small cohorts are used as normative values. In reality, both the between-subjects and within-subject flow variabilities are large. Consequently, neither one-shot flow measurements nor mean values across a cohort are truly indicative of the flow regime in a given person. We develop models for both the between-subjects and within-subject variability of internal carotid flow. A log-linear mixed effects model is combined with a Gaussian process to model the between-subjects flow variability, while a lumped parameter model of cerebral autoregulation is used to model the within-subject flow variability in response to heart rate and blood pressure changes. The model parameters are identified from carotid ultrasound measurements in a cohort of 103 elderly volunteers. We use the models to study intracranial aneurysm flow in 54 subjects under rest and exercise and conclude that OSI, a common wall shear-stress derived quantity in vascular CFD studies, may be too sensitive to flow fluctuations to be a reliable biomarker.

4D-CT angiography versus 3D-rotational angiography as the imaging modality for computational fluid dynamics of cerebral aneurysms

Background and purpose

Computational fluid dynamics (CFD) can provide valuable information regarding intracranial hemodynamics. Patient-specific models can be segmented from various imaging modalities, which may influence the geometric output and thus hemodynamic results. This study aims to compare CFD results from aneurysm models segmented from three-dimensional rotational angiography (3D-RA) versus novel four-dimensional CT angiography (4D-CTA).

Methods

Fourteen patients with 16 cerebral aneurysms underwent novel 4D-CTA followed by 3D-RA. Endoluminal geometries were segmented from each modality using an identical workflow, blinded to the other modality, to produce 28 'original' models. Each was then minimally edited a second time to match length of branches, producing 28 additional 'matched' models. CFD simulations were performed using estimated flow rates for 'original' models (representing real-world experience) and patient-specific flow rates from 4D-CTA for 'matched' models (to control for influence of modality alone).

Results

Overall, geometric and hemodynamic results were consistent between models segmented from 3D-RA and 4D-CTA, with correlations improving after matching to control for operator-introduced variability. Despite smaller 4D-CTA parent artery diameters (3.49±0.97 mm vs 3.78±0.92 mm for 3D-RA; p=0.005) and sac volumes (157 (37–750 mm3) vs 173 (53–770 mm3) for 3D-RA; p=0.0002), sac averages of time-averaged wall shear stress (TAWSS), oscillatory shear (OSI), and high frequency fluctuations (measured by spectral power index, SPI) were well correlated between 3D-RA and 4D-CTA 'matched' control models (TAWSS, R2=0.91; OSI, R2=0.79; SPI, R2=0.90).

Conclusions

Our study shows that CFD performed using 4D-CTA models produces reliable geometric and hemodynamic information in the intracranial circulation. 4D-CTA may be considered as a follow-up imaging tool for hemodynamic assessment of cerebral aneurysms.

Classification of Blood Flow Patterns in Cerebral Aneurysms

We present a Cerebral Aneurysm Vortex Classification (CAVOCLA) that allows to classify blood flow in cerebral aneurysms. Medical studies assume a strong relation between the progression and rupture of aneurysms and flow patterns. To understand how flow patterns impact the vessel morphology, they are manually classified according to predefined classes. However, manual classifications are time-consuming and exhibit a high inter-observer variability. In contrast, our approach is more objective and faster than manual methods. The classification of integral lines, representing steady or unsteady blood flow, is based on a mapping of the aneurysm surface to a hemisphere by calculating polar-based coordinates. The lines are clustered and for each cluster a representative is calculated. Then, the polar-based coordinates are transformed to the representative as basis for the classification. Classes are based on the flow complexity. The classification results are presented by a detail-on-demand approach using a visual transition from the representative over an enclosing surface to the associated lines. Based on seven representative datasets, we conduct an informal interview with five domain experts to evaluate the system. They confirmed that CAVOCLA allows for a robust classification of intra-aneurysmal flow patterns. The detail-on-demand visualization enables an efficient exploration and interpretation of flow patterns.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-Phase Ib: Effect of morphology on hemodynamics

Background

Image-based blood flow simulations have been increasingly applied to investigate intracranial aneurysm (IA) hemodynamics. However, the acceptance among physicians remains limited due to the high variability in the underlying assumptions and quality of results.

Methods

To evaluate the vessel segmentation as one of the most important sources of error, the international Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was announced. 26 research groups from 13 different countries segmented three datasets, which contained five IAs in total. Based on these segmentations, 73 time-dependent blood flow simulations under consistent conditions were carried out. Afterwards, relevant flow and shear parameters (e.g., neck inflow rate, parent vessel flow rate, spatial mean velocity, and wall shear stress) were analyzed both qualitatively and quantitatively.

Results

Regarding the entire vasculature, the variability of the segmented vessel radius is 0.13 mm, consistent and independent of the local vessel radius. However, the centerline velocity shows increased variability in more distal vessels. Focusing on the aneurysms, clear differences in morphological and hemodynamic parameters were observed. The quantification of the segmentation-induced variability showed approximately a 14% difference among the groups for the parent vessel flow rate. Regarding the mean aneurysmal velocity and the neck inflow rate, a variation of 30% and 46% was observed, respectively. Finally, time-averaged wall shear stresses varied between 28% and 51%, depending on the aneurysm in question.

Conclusions

MATCH reveals the effect of state-of-the-art segmentation algorithms on subsequent hemodynamic simulations for IA research. The observed variations may lead to an inappropriate interpretation of the simulation results and thus, can lead to inappropriate conclusions by physicians. Therefore, accurate segmentation of the region of interest is necessary to obtain reliable and clinically helpful flow information.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-phase II: rupture risk assessment

PURPOSE

Assessing the rupture probability of intracranial aneurysms (IAs) remains challenging. Therefore, hemodynamic simulations are increasingly applied toward supporting physicians during treatment planning. However, due to several assumptions, the clinical acceptance of these methods remains limited.

METHODS

To provide an overview of state-of-the-art blood flow simulation capabilities, the Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was conducted. Seventeen research groups from all over the world performed segmentations and hemodynamic simulations to identify the ruptured aneurysm in a patient harboring five IAs. Although simulation setups revealed good similarity, clear differences exist with respect to the analysis of aneurysm shape and blood flow results. Most groups (12/71%) included morphological and hemodynamic parameters in their analysis, with aspect ratio and wall shear stress as the most popular candidates, respectively.

RESULTS

The majority of groups (7/41%) selected the largest aneurysm as being the ruptured one. Four (24%) of the participating groups were able to correctly select the ruptured aneurysm, while three groups (18%) ranked the ruptured aneurysm as the second most probable. Successful selections were based on the integration of clinically relevant information such as the aneurysm site, as well as advanced rupture probability models considering multiple parameters. Additionally, flow characteristics such as the quantification of inflow jets and the identification of multiple vortices led to correct predictions.

CONCLUSIONS

MATCH compares state-of-the-art image-based blood flow simulation approaches to assess the rupture risk of IAs. Furthermore, this challenge highlights the importance of multivariate analyses by combining clinically relevant metadata with advanced morphological and hemodynamic quantification.

Real-World Variability in the Prediction of Intracranial Aneurysm Wall Shear Stress: The 2015 International Aneurysm CFD Challenge

Purpose

Image-based computational fluid dynamics (CFD) is widely used to predict intracranial aneurysm wall shear stress (WSS), particularly with the goal of improving rupture risk assessment. Nevertheless, concern has been expressed over the variability of predicted WSS and inconsistent associations with rupture. Previous challenges, and studies from individual groups, have focused on individual aspects of the image-based CFD pipeline. The aim of this Challenge was to quantify the total variability of the whole pipeline.

Methods

3D rotational angiography image volumes of five middle cerebral artery aneurysms were provided to participants, who were free to choose their segmentation methods, boundary conditions, and CFD solver and settings. Participants were asked to fill out a questionnaire about their solution strategies and experience with aneurysm CFD, and provide surface distributions of WSS magnitude, from which we objectively derived a variety of hemodynamic parameters.

Results

A total of 28 datasets were submitted, from 26 teams with varying levels of self-assessed experience. Wide variability of segmentations, CFD model extents, and inflow rates resulted in interquartile ranges of sac average WSS up to 56%, which reduced to < 30% after normalizing by parent artery WSS. Sac-maximum WSS and low shear area were more variable, while rank-ordering of cases by low or high shear showed only modest consensus among teams. Experience was not a significant predictor of variability.

Conclusions

Wide variability exists in the prediction of intracranial aneurysm WSS. While segmentation and CFD solver techniques may be difficult to standardize across groups, our findings suggest that some of the variability in image-based CFD could be reduced by establishing guidelines for model extents, inflow rates, and blood properties, and by encouraging the reporting of normalized hemodynamic parameters.

Modeling hemodynamics in intracranial aneurysms: Comparing accuracy of CFD solvers based on finite element and finite volume schemes

Image-based computational fluid dynamics (CFD) has shown potential to aid in the clinical management of intracranial aneurysms, but its adoption in the clinical practice has been missing, partially because of lack of accuracy assessment and sensitivity analysis. To numerically solve the flow-governing equations, CFD solvers generally rely on 2 spatial discretization schemes: finite volume (FV) and finite element (FE). Since increasingly accurate numerical solutions are obtained by different means, accuracies and computational costs of FV and FE formulations cannot be compared directly. To this end, in this study, we benchmark 2 representative CFD solvers in simulating flow in a patient-specific intracranial aneurysm model: (1) ANSYS Fluent, a commercial FV-based solver, and (2) VMTKLab multidGetto, a discontinuous Galerkin (dG) FE-based solver. The FV solver's accuracy is improved by increasing the spatial mesh resolution (134k, 1.1m, 8.6m, and 68.5m tetrahedral element meshes). The dGFE solver accuracy is increased by increasing the degree of polynomials (first, second, third, and fourth degree) on the base 134k tetrahedral element mesh. Solutions from best FV and dGFE approximations are used as baseline for error quantification. On average, velocity errors for second-best approximations are approximately 1 cm/s for a [0,125] cm/s velocity magnitude field. Results show that high-order dGFE provides better accuracy per degree of freedom but worse accuracy per Jacobian nonzero entry as compared with FV. Cross-comparison of velocity errors demonstrates asymptotic convergence of both solvers to the same numerical solution. Nevertheless, the discrepancy between underresolved velocity fields suggests that mesh independence is reached following different paths.

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.

Better Than Nothing: A Rational Approach for Minimizing the Impact of Outflow Strategy on Cerebrovascular Simulations

BACKGROUND AND PURPOSE

Computational fluid dynamics simulations of neurovascular diseases are impacted by various modeling assumptions and uncertainties, including outlet boundary conditions. Many studies of intracranial aneurysms, for example, assume zero pressure at all outlets, often the default ("do-nothing") strategy, with no physiological basis. Others divide outflow according to the outlet diameters cubed, nominally based on the more physiological Murray's law but still susceptible to subjective choices about the segmented model extent. Here we demonstrate the limitations and impact of these outflow strategies, against a novel "splitting" method introduced here.

MATERIALS AND METHODS

With our method, the segmented lumen is split into its constituent bifurcations, where flow divisions are estimated locally using a power law. Together these provide the global outflow rate boundary conditions. The impact of outflow strategy on flow rates was tested for 70 cases of MCA aneurysm with 0D simulations. The impact on hemodynamic indices used for rupture status assessment was tested for 10 cases with 3D simulations.

RESULTS

Differences in flow rates among the various strategies were up to 70%, with a non-negligible impact on average and oscillatory wall shear stresses in some cases. Murray-law and splitting methods gave flow rates closest to physiological values reported in the literature; however, only the splitting method was insensitive to arbitrary truncation of the model extent.

CONCLUSIONS

Cerebrovascular simulations can depend strongly on the outflow strategy. The default zero-pressure method should be avoided in favor of Murray-law or splitting methods, the latter being released as an open-source tool to encourage the standardization of outflow strategies.

Initial Clinical Experience with AView-A Clinical Computational Platform for Intracranial Aneurysm Morphology, Hemodynamics, and Treatment Management

BACKGROUND

The management of intracranial aneurysm (IA) is challenging. Clinicians often rely on varied and intuitively disparate ways of evaluating rupture risk that may only partially take into account complex hemodynamic and morphologic factors. We developed a prototype of a clinically oriented, streamlined, computational platform, AView, for rapid assessment of hemodynamics and morphometrics in clinical settings. To show the potential clinical utility of AView, we report our initial multicenter experience highlighting the possible advantages of morphologic and hemodynamic analysis of IAs.

METHODS

AView software was deployed across 8 medical centers (6 in the United States, 2 in Japan). Eight clinicians were trained and used the AView software between September 2012 and January 2013.

RESULTS

We present 12 illustrative cases that show the potential clinical utility of AView. For all, morphology and hemodynamics, flow visualization, and rupture resemblance score (a surrogate for rupture risk) were provided. In 3 cases, AView could confirm the clinicians' decision to treat; in 3 cases, it could suggest which aneurysms may be at greater risk among multiple aneurysms; in 5 cases, AView could provide additional information for use during treatment decisions for ambiguous situations. In one stent-assisted coiling case, flow visualization predicted that the intuitive choice for stent placement could have resulted in sacrifice of an anterior cerebral artery due to blockage by coils and led clinicians to reconsider treatment plans.

CONCLUSIONS

AView has the potential to confirm decisions to treat IAs, suggest which among multiple aneurysms to treat, and guide treatment decisions. Furthermore, the flow visualization it affords can inform aneurysm treatment planning and potentially avoid poor outcomes.

Multiple intracranial aneurysms: a direct hemodynamic comparison between ruptured and unruptured vessel malformations

PURPOSE

Despite numerous studies addressing the rupture risk of intracranial aneurysms that have been published, the assessment thereof still remains challenging. Image-based simulations enable a precise prediction of patient-specific blood flow information. However, those approaches normally consider only small segments of the complete cerebral vasculature.

METHODS

To test the validity of the consideration of single aneurysms in one computational setup, domains of the complete anterior and posterior circulations with multiple intracranial aneurysms (MIA) were simulated. Six patients with MIA were investigated, while 3D surfaces of eleven unruptured and six ruptured aneurysms were segmented. The segmentations were used for the determination of morphological parameters and also for image-based blood flow simulations used to characterize the hemodynamic properties of each aneurysm.

RESULTS

In the geometric comparison, neck aspect ratios of unruptured and ruptured aneurysms did not differ significantly. In contrast, size ratios, aspect ratios, surface areas, volumes, and non-sphericity indices were significantly higher in the ruptured cases. The analysis of hemodynamic parameters demonstrated that in each patient, the ruptured aneurysm exhibited the lowest averaged wall shear stresses and highest oscillatory shears. Unstable flow was also detected in ruptured aneurysms based on increased oscillatory velocity.

CONCLUSION

In this small study involving patients with MIA, different morphologies and flow patterns were observed between ruptured and unruptured aneurysms. The analysis of the hemodynamics in such patients revealed a good agreement with studies that only considered single malformations. Additionally, complex flow patterns are detected in ruptured cases, which require deeper investigation.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre-operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient-specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error-free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post-treatment intra-aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random-effects meta-analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta-analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.1385 For further resources related to this article, please visit the WIREs website.

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.

Impact of Using Conventional Inlet/Outlet Boundary Conditions on Haemodynamic Metrics in a Subject-Specific Rabbit Aorta

Computational fluid dynamics is a tool capable of accurately measuring metrics currently used to predict the behaviour of cardiovascular diseases. This study quantifies the impact various commonly used inlet and outlet boundary conditions have on various shear rate–based haemodynamic metrics currently used for predicting the localisation of cardiovascular diseases. Simulations are conducted on an accurately represented rabbit aorta configuration and comparison has been made against available in vivo data. The boundary conditions studied include two different inlet profiles, three outlet boundary conditions, and steady-state versus pulsatile flow cases. Large variations were found in the results, particularly when using different outlet boundary conditions, and the discrepancies were evaluated both quantitatively and qualitatively. The results clearly highlight the importance of the type of boundary condition used when simulating complex cardiovascular models. By restricting the attention to the flow within the aorta and the intercostal branches, the results suggest that prescribing transient simulation and fully developed flow at the inlet are not required. Furthermore, assuming the widely accepted low wall shear stress theory of Caro, it was found that Murray’s law–based outlet boundary condition returns the most physiologically accurate results when compared to in vivo data.

Combined Visualization of Vessel Deformation and Hemodynamics in Cerebral Aneurysms

We present the first visualization tool that combines patient-specific hemodynamics with information about the vessel wall deformation and wall thickness in cerebral aneurysms. Such aneurysms bear the risk of rupture, whereas their treatment also carries considerable risks for the patient. For the patient-specific rupture risk evaluation and treatment analysis, both morphological and hemodynamic data have to be investigated. Medical researchers emphasize the importance of analyzing correlations between wall properties such as the wall deformation and thickness, and hemodynamic attributes like the Wall Shear Stress and near-wall flow. Our method uses a linked 2.5D and 3D depiction of the aneurysm together with blood flow information that enables the simultaneous exploration of wall characteristics and hemodynamic attributes during the cardiac cycle. We thus offer medical researchers an effective visual exploration tool for aneurysm treatment risk assessment. The 2.5D view serves as an overview that comprises a projection of the vessel surface to a 2D map, providing an occlusion-free surface visualization combined with a glyph-based depiction of the local wall thickness. The 3D view represents the focus upon which the data exploration takes place. To support the time-dependent parameter exploration and expert collaboration, a camera path is calculated automatically, where the user can place landmarks for further exploration of the properties. We developed a GPU-based implementation of our visualizations with a flexible interactive data exploration mechanism. We designed our techniques in collaboration with domain experts, and provide details about the evaluation.

Association of Hemodynamic Factors With Intracranial Aneurysm Formation and Rupture: Systematic Review and Meta-analysis

BACKGROUND

Recent evidence suggests a link between the magnitude and distribution of hemodynamic factors and the formation and rupture of intracranial aneurysms. However, there are many conflicting results.

OBJECTIVE

To quantify the effect of hemodynamic factors on aneurysm formation and their association with ruptured aneurysms.

METHODS

We performed a systematic review and meta-analysis through October 2014. Analysis of the effects of hemodynamic factors on aneurysm formation was performed by pooling the results of studies that compared geometrical models of intracranial aneurysms and "preaneurysm" models where the aneurysm was artificially removed. Furthermore, we calculated pooled standardized mean differences between ruptured and unruptured aneurysms to quantify the association of hemodynamic factors with ruptured aneurysms. Standard PRISMA guidelines were followed.

RESULTS

The hemodynamic factors that showed high positive correlations with location of aneurysm formation were high wall shear stress (WSS) and high gradient oscillatory number, with pooled proportions of 78.8% and 85.7%, respectively. Positive correlations were largely seen in bifurcation aneurysms, whereas negative correlations were seen in sidewall aneurysms. Mean and normalized WSS were significantly lower and low shear area significantly higher in ruptured aneurysms.

CONCLUSION

Pooled analyses of computational fluid dynamics models suggest that an increase in WSS and gradient oscillatory number may contribute to aneurysm formation, whereas low WSS is associated with ruptured aneurysms. The location of the aneurysm at the bifurcation or sidewall may influence the correlation of these hemodynamic factors.

Blood Flow Into Basilar Tip Aneurysms

Background and Purpose—

Hemodynamic forces may play a role in the recanalization of coiled aneurysms. The purpose of this study was to investigate the influence of presurgical hemodynamics on the efficacy of coil embolization for basilar tip aneurysms.

Methods—

We identified 82 patients who underwent endovascular coil embolization for basilar tip aneurysms with a follow-up of >1 year. Presurgical hemodynamics were investigated using computational fluid dynamics with 3-dimensional data derived from rotational angiography. During postprocessing, we quantified the rate of net flow entering the aneurysm through its neck and calculated the proportion of the aneurysmal inflow rate to the basilar artery flow rate. In addition, we investigated the correlation between the basilar bifurcation configuration and the hemodynamics.

Results—

Twenty-five of the 82 patients were excluded because of difficult vascular geometry reconstruction. Among the 57 examined patients, angiographic recanalization was observed in 19 patients (33.3%). The proportion of the aneurysmal inflow rate to the basilar artery flow rate and a coil packing density <30% were independent and significant predictors for the recanalization of coiled aneurysms. Additional investigation revealed that a small branch angle formed by the basilar artery and the posterior cerebral artery increased blood flow into the aneurysm.

Conclusions—

The proportion of the aneurysmal inflow rate to the basilar artery flow rate, influenced by the basilar bifurcation configuration, was an independent and significant predictor for recanalization after coil embolization in basilar tip aneurysms.

Bringing hemodynamic simulations closer to the clinics: a CFD prototype study for intracranial aneurysms

Computational Fluid Dynamics enables the investigation of patient-specific hemodynamics for rupture predictions and treatment support of intracranial aneurysms. However, due to numerous simplifications to decrease the computations effort, clinical applicability is limited until now. To overcome this situation a clinical research software prototype was tested that can be easily operated by attending physicians. In order to evaluate the accuracy of this prototype, four patient-specific intracranial aneurysms were investigated using four different spatial resolutions. The results demonstrate that physicians were able to generate hemodynamic predictions within several minutes at low spatial resolution. However, depending on the parameter of interest and the desired accuracy, higher resolutions are required, which will lead to an increase of computational times that still look very attractive towards clinical usability. The study shows that the next step towards an applicable individualized therapy for patients harboring intracranial aneurysms can be done. However, further in vivo validations are required to guarantee realistic predictions in future studies.

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.