Computational fluid dynamic (CFD) simulations of blood flow provide new insights into the hemodynamics of vascular pathologies such as cerebral aneurysms. Understanding the relations between hemodynamics and aneurysm initiation, progression, and risk of rupture is crucial in diagnosis and treatment. Recent studies link the existence of vortices in the blood flow pattern to aneurysm rupture and report observations of embedded vortices -a larger vortex encloses a smaller one flowing in the opposite direction -whose implications are unclear. We present a clustering-based approach for the visual analysis of vortical flow in simulated cerebral aneurysm hemodynamics. We show how embedded vortices develop at saddle-node bifurcations on vortex core lines and convey the participating flow at full manifestation of the vortex by a fast and smart grouping of streamlines and the visualization of group representatives. The grouping result may be refined based on spectral clustering generating a more detailed visualization of the flow pattern, especially further off the core lines. We aim at supporting CFD engineers researching the biological implications of embedded vortices.

Rupture of intracranial saccular aneurysms is the most common cause of spontaneous subarachnoid hemorrhage, which has significant morbidity and mortality. Although there is still controversy regarding the decision on which unruptured aneurysms should be treated, this is based primarily on their size. Nonetheless, many large lesions do not rupture whereas some small ones do. It is commonly accepted that hemodynamical factors are important to better understand the natural history of cerebral aneurysms. However, it might not always be practical to carry out a detailed computational analysis of such factors if a prompt assessment is required. Since shape is likely to be dependent on the balance between hemodynamic forces and the aneurysmal surrounding environment, an appropriate morphological 3-D characterization is likely to provide a practical surrogate to quickly evaluate the risk of rupture. In this paper, an efficient and novel methodology for 3-D shape characterization of cerebral aneurysms is described. The aneurysms are isolated by taking into account a portion of their adjacent vessels. Two methods to characterize the morphology of the aneurysms models using moment invariants have been considered: geometrical moment invariants (GMI) and Zernike moment invariants (ZMI). The results have been validated in a database containing 53 patients with a total of 31 ruptured aneurysms and 24 unruptured aneurysms. It has been found that ZMI indices are more robust than GMI, and seem to provide a reliable way to discriminate between ruptured and unruptured aneurysms. Correct rupture prediction rates of approximately equal to 80% were achieved in contrast to 66% that is found when the aspect ratio index is considered.

Hemodynamic factors are thought to be implicated in the progression and rupture of intracranial aneurysms. Current efforts aim to study the possible associations of hemodynamic characteristics such as complexity and stability of intra-aneurysmal flow patterns, size and location of the region of flow impingement with the clinical history of aneurysmal rupture. However, there are no reliable methods for measuring blood flow patterns in vivo. In this paper, an efficient methodology for patient-specific modeling and characterization of the hemodynamics in cerebral aneurysms from medical images is described. A sensitivity analysis of the hemodynamic characteristics with respect to variations of several variables over the expected physiologic range of conditions is also presented. This sensitivity analysis shows that although changes in the velocity fields can be observed, the characterization of the intra-aneurysmal flow patterns is not altered when the mean input flow, the flow division, the viscosity model, or mesh resolution are changed. It was also found that the variable that has the greater impact on the computed flow fields is the geometry of the vascular structures. We conclude that with the proposed modeling pipeline clinical studies involving large numbers cerebral aneurysms are feasible.

The simulation of blood flow past endovascular devices such as coils and stents is a challenging problem due to the complex geometry of the devices. Traditional unstructured grid computational fluid dynamics relies on the generation of finite element grids that conform to the boundary of the computational domain. However, the generation of such grids for patient-specific modeling of cerebral aneurysm treatment with coils or stents is extremely difficult and time consuming. This paper describes the application of an adaptive grid embedding technique previously developed for complex fluid structure interaction problems to the simulation of endovascular devices. A hybrid approach is used: the vessel walls are treated with body conforming grids and the endovascular devices with an adaptive mesh embedding technique. This methodology fits naturally in the framework of image-based computational fluid dynamics and opens the door for exploration of different therapeutic options and personalization of endovascular procedures.

Virtual bronchoscopy reconstructions of the airway noninvasively provide useful morphologic information of structural abnormalities such as stenoses and masses. In this paper, we show how virtual bronchoscopy can be used to perform aerodynamic calculations in anatomically realistic models. Pressure and flow patterns in a human airway were computed noninvasively. These showed decreased pressure and increased shear stress in the region of a stenosis.

Three-dimensional (3-D) angiographic methods are gaining acceptance for evaluation of atherosclerotic disease. However, measurement of vessel stenosis from 3-D angiographic methods can be problematic due to limited image resolution and contrast. We present a method for reconstructing vessel surfaces from 3-D angiographic methods that allows for objective measurement of vessel stenosis. The method is a deformable model that employs a tubular coordinate system. Vertex merging is incorporated into the coordinate system to maintain even vertex spacing and to avoid problems of self-intersection of the surface. The deformable model was evaluated on clinical magnetic resonance (MR) images of the carotid (n = 6) and renal (n = 2) arteries, on an MR image of a physical vascular phantom and on a digital vascular phantom. Only one gross error occurred for all clinical images. All reconstructed surfaces had a realistic, smooth appearance. For all segments of the physical vascular phantom, vessel radii from the surface reconstruction had an error of less than 0.2 of the average voxel dimension. Variability of manual initialization of the deformable model had negligible effect on the measurement of the degree of stenosis of the digital vascular phantom.