Hemodynamics of Cerebral Aneurysms: Computational Analyses of Aneurysm Progress and Treatment

The progression of a cerebral aneurysm involves degenerative arterial wall remodeling. Various hemodynamic parameters are suspected to be major mechanical factors related to the genesis and progression of vascular diseases. Flow alterations caused by the insertion of coils and stents for interventional aneurysm treatment may affect the aneurysm embolization process. Therefore, knowledge of hemodynamic parameters may provide physicians with an advanced understanding of aneurysm progression and rupture, as well as the effectiveness of endovascular treatments. Progress in medical imaging and information technology has enabled the prediction of flow fields in the patient-specific blood vessels using computational analysis. In this paper, recent computational hemodynamic studies on cerebral aneurysm initiation, progress, and rupture are reviewed. State-of-the-art computational aneurysmal flow analyses after coiling and stenting are also summarized. We expect the computational analysis of hemodynamics in cerebral aneurysms to provide valuable information for planning and follow-up decisions for treatment

An automated multiscale ensemble simulation approach for vascular blood flow

Cerebrovascular diseases such as brain aneurysms are a primary cause of adult disability. The flow dynamics in brain arteries, both during periods of rest and increased activity, are known to be a major factor in the risk of aneurysm formation and rupture. The precise relation is however still an open field of investigation. We present an automated ensemble simulation method for modelling cerebrovascular blood flow under a range of flow regimes. By automatically constructing and performing an ensemble of multiscale simulations, where we unidirectionally couple a 1D solver with a 3D lattice-Boltzmann code, we are able to model the blood flow in a patient artery over a range of flow regimes. We apply the method to a model of a middle cerebral artery, and find that this approach helps us to fine-tune our modelling techniques, and opens up new ways to investigate cerebrovascular flow properties.This work has received funding from the CRESTA project within the EC-FP7 (ICT-2011.9.13) under Grant Agreements no. 287703, and from EPSRC Grants EP/I017909/1 (www.2020science.net) and EP/I034602/1

Estimated pretreatment hemodynamic prognostic factors of aneurysm recurrence after endovascular embolization.

BACKGROUND:Hemodynamic factors play important roles in aneurysm recurrence after endovascular treatment. OBJECTIVE:Predicting the risk of recurrence by hemodynamic analysis using an untreated aneurysm model is important because such prediction is required before treatment. METHODS:We retrospectively analyzed hemodynamic factors associated with aneurysm recurrence from pretreatment models of five recurrent and five stable posterior communicating artery (Pcom) aneurysms with no significant differences in aneurysm volume, coil packing density, or sizes of the dome, neck, or Pcom. Hemodynamic factors of velocity ratio, flow rate, pressure ratio, and wall shear stress were investigated. RESULTS:Among the hemodynamic factors investigated, velocity ratio and flow rate of the Pcom showed significant differences between the recurrence group and stable group (0.630 ± 0.062 and 0.926 ± 0.051, P= 0.016; 56.4 ± 8.9 and 121.6 ± 6.7, P= 0.008, respectively). CONCLUSIONS:Our results suggest that hemodynamic factors may be associated with aneurysm recurrence among Pcom aneurysms. Velocity and flow rate in the Pcom may be a pretreatment prognostic factor for aneurysm recurrence after endovascular treatment

The 'Sphere': A Dedicated Bifurcation Aneurysm Flow-Diverter Device.

We present flow-based results from the early stage design cycle, based on computational modeling, of a prototype flow-diverter device, known as the 'Sphere', intended to treat bifurcation aneurysms of the cerebral vasculature. The device is available in a range of diameters and geometries and is constructed from a single loop of NITINOL(®) wire. The 'Sphere' reduces aneurysm inflow by means of a high-density, patterned, elliptical surface that partially occludes the aneurysm neck. The device is secured in the healthy parent vessel by two armatures in the shape of open loops, resulting in negligible disruption of parent or daughter vessel flow. The device is virtually deployed in six anatomically accurate bifurcation aneurysms: three located at the Basilar tip and three located at the terminus bifurcation of the Internal Carotid artery (at the meeting of the middle cerebral and anterior cerebral arteries). Both steady state and transient flow simulations reveal that the device presents with a range of aneurysm inflow reductions, with mean flow reductions falling in the range of 30.6-71.8% across the different geometries. A significant difference is noted between steady state and transient simulations in one geometry, where a zone of flow recirculation is not captured in the steady state simulation. Across all six aneurysms, the device reduces the WSS magnitude within the aneurysm sac, resulting in a hemodynamic environment closer to that of a healthy vessel. We conclude from extensive CFD analysis that the 'Sphere' device offers very significant levels of flow reduction in a number of anatomically accurate aneurysm sizes and locations, with many advantages compared to current clinical cylindrical flow-diverter designs. Analysis of the device's mechanical properties and deployability will follow in future publications

Modelling the evolution of cerebral aneurysms: biomechanics, mechanobiology and multiscale modelling

Intracranial aneurysms (IAs) are abnormal dilatations of the cerebral vasculature. Computational modelling may shed light on the aetiology of the disease and lead to improved criteria to assist diagnostic decisions. We briefly review models of aneurysm evolution to date and present a novel fluid-solid-growth (FSG) framework for patient-specific modelling of IA evolution. We illustrate its application to 4 clinical cases depicting an IA. The section of arterial geometry containing the IA is removed and replaced with a cylindrical section: this represents an idealised section of healthy artery upon which IA evolution is simulated. The utilisation of patient-specific geometries enables G&R to be explicitly linked to physiologically realistic spatial distributions and magnitudes of haemodynamic stimuli. In this study, we investigate the hypothesis that elastin degradation is driven by locally low wall shear stress (WSS). In 3 out of 4 cases, the evolved model IA geometry is qualitatively similar to the corresponding in vivo IA geometry. This suggests some tentative support for the hypothesis that low WSS plays a role in the mechanobiology of IA evolution

A Rapid and Computationally Inexpensive Method to Virtually Implant Current and Next-Generation Stents into Subject-Specific Computational Fluid Dynamics Models

Computational modeling is often used to quantify hemodynamic alterations induced by stenting, but frequently uses simplified device or vascular representations. Based on a series of Boolean operations, we developed an efficient and robust method for assessing the influence of current and next-generation stents on local hemodynamics and vascular biomechanics quantified by computational fluid dynamics. Stent designs were parameterized to allow easy control over design features including the number, width and circumferential or longitudinal spacing of struts, as well as the implantation diameter and overall length. The approach allowed stents to be automatically regenerated for rapid analysis of the contribution of design features to resulting hemodynamic alterations. The applicability of the method was demonstrated with patient-specific models of a stented coronary artery bifurcation and basilar trunk aneurysm constructed from medical imaging data. In the coronary bifurcation, we analyzed the hemodynamic difference between closed-cell and open-cell stent geometries. We investigated the impact of decreased strut size in stents with a constant porosity for increasing flow stasis within the stented basilar aneurysm model. These examples demonstrate the current method can be used to investigate differences in stent performance in complex vascular beds for a variety of stenting procedures and clinical scenarios

Immersed boundary method predictions of shear stresses for different flow topologies occuring in cerebral aneurysms

A volume-penalizing immersed boundary method is presented that facilitates the computation of incompressible fluid flow in complex flow domains. We apply this method to simulate the flow in cerebral aneurysms, and focus on the accuracy with which the flow field and the corresponding shear stress field are computed. The method is applied to laminar, incompressible flow in curved cylindrical vessels and in a model aneurysm. The time-dependent shear stress distributions over the vessel walls are visualized and interpreted in terms of the flow fields that develop. We compute shear stress levels at two different Reynolds numbers, corresponding to a steady and an unsteady flow. In the latter situation strong fluctuations in the shear stress are observed, that may be connected to raised risk-levels of aneurysm rupture