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.

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.

Closing the loop: modelling of heart failure progression from health to end-stage using a meta-analysis of left ventricular pressure-volume loops.

INTRODUCTION: The American Heart Association (AHA)/American College of Cardiology (ACC) guidelines for the classification of heart failure (HF) are descriptive but lack precise and objective measures which would assist in categorising such patients. Our aim was two fold, firstly to demonstrate quantitatively the progression of HF through each stage using a meta-analysis of existing left ventricular (LV) pressure-volume (PV) loop data and secondly use the LV PV loop data to create stage specific HF models. METHODS AND RESULTS: A literature search yielded 31 papers with PV data, representing over 200 patients in different stages of HF. The raw pressure and volume data were extracted from the papers using a digitising software package and the means were calculated. The data demonstrated that, as HF progressed, stroke volume (SV), ejection fraction (EF%) decreased while LV volumes increased. A 2-element lumped parameter model was employed to model the mean loops and the error was calculated between the loops, demonstrating close fit between the loops. The only parameter that was consistently and statistically different across all the stages was the elastance (Emax). CONCLUSIONS: For the first time, the authors have created a visual and quantitative representation of the AHA/ACC stages of LVSD-HF, from normal to end-stage. The study demonstrates that robust, load-independent and reproducible parameters, such as elastance, can be used to categorise and model HF, complementing the existing classification. The modelled PV loops establish previously unknown physiological parameters for each AHA/ACC stage of LVSD-HF, such as LV elastance and highlight that it this parameter alone, in lumped parameter models, that determines the severity of HF. Such information will enable cardiovascular modellers with an interest in HF, to create more accurate models of the heart as it fails.

Accuracy vs. computational time: translating aortic simulations to the clinic.

State of the art simulations of aortic haemodynamics feature full fluid-structure interaction (FSI) and coupled 0D boundary conditions. Such analyses require not only significant computational resource but also weeks to months of run time, which compromises the effectiveness of their translation to a clinical workflow. This article employs three computational fluid methodologies, of varying levels of complexity with coupled 0D boundary conditions, to simulate the haemodynamics within a patient-specific aorta. The most comprehensive model is a full FSI simulation. The simplest is a rigid walled incompressible fluid simulation while an alternative middle-ground approach employs a compressible fluid, tuned to elicit a response analogous to the compliance of the aortic wall. The results demonstrate that, in the context of certain clinical questions, the simpler analysis methods may capture the important characteristics of the flow field.

Computational modelling and evaluation of cardiovascular response under pulsatile impeller pump support.

This study presents a numerical simulation of cardiovascular response in the heart failure condition under the support of a Berlin Heart INCOR impeller pump-type ventricular assist device (VAD). The model is implemented using the CellML modelling language. To investigate the potential of using the Berlin Heart INCOR impeller pump to produce physiologically meaningful arterial pulse pressure within the various physiological constraints, a series of VAD-assisted cardiovascular cases are studied, in which the pulsation ratio and the phase shift of the VAD motion profile are systematically changed to observe the cardiovascular responses in each of the studied cases. An optimization process is proposed, including the introduction of a cost function to balance the importance of the characteristic cardiovascular variables. Based on this cost function it is found that a pulsation ratio of 0.35 combined with a phase shift of 200° produces the optimal cardiovascular response, giving rise to a maximal arterial pulse pressure of 12.6 mm Hg without inducing regurgitant pump flow while keeping other characteristic cardiovascular variables within appropriate physiological ranges.