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.

Error analysis of marker-based object localization using a single-plane XRII.

The role of imaging and image guidance is increasing in surgery and therapy, including treatment planning and follow-up. Fluoroscopy is used for two-dimensional (2D) guidance or localization; however, many procedures would benefit from three-dimensional (3D) guidance or localization. Three-dimensional computed tomography (CT) using a C-arm mounted x-ray image intensifier (XRII) can provide high-quality 3D images; however, patient dose and the required acquisition time restrict the number of 3D images that can be obtained. C-arm based 3D CT is therefore limited in applications for x-ray based image guidance or dynamic evaluations. 2D-3D model-based registration, using a single-plane 2D digital radiographic system, does allow for rapid 3D localization. It is our goal to investigate-over a clinically practical range-the impact of x-ray exposure on the resulting range of 3D localization precision. In this paper it is assumed that the tracked instrument incorporates a rigidly attached 3D object with a known configuration of markers. A 2D image is obtained by a digital fluoroscopic x-ray system and corrected for XRII distortions (+/- 0.035 mm) and mechanical C-arm shift (+/- 0.080 mm). A least-square projection-Procrustes analysis is then used to calculate the 3D position using the measured 2D marker locations. The effect of x-ray exposure on the precision of 2D marker localization and on 3D object localization was investigated using numerical simulations and x-ray experiments. The results show a nearly linear relationship between 2D marker localization precision and the 3D localization precision. However, a significant amplification of error, nonuniformly distributed among the three major axes, occurs, and that is demonstrated. To obtain a 3D localization error of less than +/- 1.0 mm for an object with 20 mm marker spacing, the 2D localization precision must be better than +/- 0.07 mm. This requirement was met for all investigated nominal x-ray exposures at 28 cm FOV, and for all but the lowest two at 40 cm FOV. However, even for those two nominal exposures, the expected 3D localization error is less than +/- 1.2 mm. The tracking precision was +/- 0.65 mm for the out-of-plane translations, +/- 0.05 mm for in-plane translations, and +/- 0.05 degrees for the rotations. The root mean square (RMS) difference between the true and projection-Procrustes calculated location was 1.07 mm. It is believed these results show the potential of this technique for dynamic evaluations or real-time image guidance using a single x-ray source and XRII detector.

Comparison of maximum intensity projection and digitally reconstructed radiographic projection for carotid artery stenosis measurement.

Digital subtraction angiography is being supplanted by three-dimensional imaging techniques in many clinical applications, leading to extensive use of maximum intensity projection (MIP) images to depict volumetric vascular data. The MIP algorithm produces intensity profiles that are different than conventional angiograms, and can also increase the vessel-to-tissue contrast-to-noise ratio. We evaluated the effect of the MIP algorithm in a clinical application where quantitative vessel measurement is important: internal carotid artery stenosis grading. Three-dimensional computed rotational angiography (CRA) was performed on 26 consecutive symptomatic patients to verify an internal carotid artery stenosis originally found using duplex ultrasound. These volumes of data were visualized using two different postprocessing projection techniques: MIP and digitally reconstructed radiographic (DRR) projection. A DRR is a radiographic image simulating a conventional digitally subtracted angiogram, but it is derived computationally from the same CRA dataset as the MIP. By visualizing a single volume with two different projection techniques, the postprocessing effect of the MIP algorithm is isolated. Vessel measurements were made, according to the NASCET guidelines, and percentage stenosis grades were calculated. The paired t-test was used to determine if the measurement difference between the two techniques was statistically significant. The CRA technique provided an isotropic voxel spacing of 0.38 mm. The MIPs and DRRs had a mean signal-difference-to-noise-ratio of 30:1 and 26:1, respectively. Vessel measurements from MIPs were, on average, 0.17 mm larger than those from DRRs (P < 0.0001). The NASCET-type stenosis grades tended to be underestimated on average by 2.4% with the MIP algorithm, although this was not statistically significant (P=0.09). The mean interobserver variability (standard deviation) of both the MIP and DRR images was 0.35 mm. It was concluded that the MIP algorithm slightly increased the apparent dimensions of the arteries, when applied to these intra-arterial CRA images. This subpixel increase was smaller than both the voxel size and interobserver variability, and was therefore not clinically relevant.

Virtual angiography for visualization and validation of computational models of aneurysm hemodynamics.

It has recently become possible to simulate aneurysmal blood flow dynamics in a patient-specific manner via the coupling of three-dimensional (3-D) X-ray angiography and cmputational fluid dynamics (CFD). Before such image-based CFD models can be used in a predictive capacity, however, it must be shown that they indeed reproduce the in vivo hemodynamic environment. Motivated by the fact that there are currently no techniques for adequately measuring complex blood velocity fields in vivo, in this paper we describe how cine X-ray angiograms may be simulated for the purpose of indirectly validating patient-sperific CFD models. Mimicking the radiological procedure, a virtual angiogram is constructed by first simulating the time-varying injection of contrast agent into a precomputed, patient-specific CFD model. A time-series of images is then constructed by simulating the attenuation of X-rays through the computed 3-D contrast-agent flow dynamics. Virtual angiographic images and residence time maps, here derived from an image-based CFD model of a giant aneurysm, are shown to be in excellent agreement wiith the corresponding clinical images and residence time maps, but only when the interaction between the quasisteady contrast agent injection and the pulsatile flow are properly accounted for. These virtual angiographic techniques pave the way for validating image-based CFD models against routinely available clinical data, and provide a means of visualizing complex, 3-D blood flow dynamics in a clinically relevant manner. They also clearly show how the contrast agent injection perturbs the noraml blood flow patterns, further highlighting the potential utility of image-based CFD as a window into the true aneurysmal hemodynamics.