Ignacio Larrabide

Fast virtual deployment of self-expandable stents: Method and in vitro evaluation for intracranial aneurysmal stenting

INTRODUCTION Minimally invasive treatment approaches, like the implantation of percutaneous stents, are becoming more popular every day for the treatment of intracranial aneurysms. The outcome of such treatments is related to factors like vessel and aneurysm geometry, hemodynamic conditions and device design. For this reason, having a tool for assessing stenting alternatives beforehand is crucial. METHODOLOGY The Fast Virtual Stenting (FVS) method, which provides an estimation of the configuration of intracranial stents when released in realistic geometries, is proposed in this paper. This method is based on constrained simplex deformable models. The constraints are used to account for the stent design. An algorithm for its computational implementation is also proposed. The performance of the proposed methodology was contrasted with real stents released in a silicone phantom. RESULTS In vitro experiments were performed on the phantom where a contrast injection was performed. Subsequently, corresponding Computational Fluid Dynamics (CFD) analyzes were carried out on a digital replica of the phantom with the virtually released stent. Virtual angiographies are used to compare in vitro experiments and CFD analysis. Contrast time-density curves for in vitro and CFD data were generated and used to compare them. CONCLUSIONS Results of both experiments resemble very well, especially when comparing the contrast density curves. The use of FVS methodology in the clinical environment could provide additional information to clinicians before the treatment to choose the therapy that best fits the patient.

Computational hemodynamics in cerebral aneurysms: the effects of modeled versus measured boundary conditions.

Modeling of flow in intracranial aneurysms (IAs) requires flow information at the model boundaries. In absence of patient-specific measurements, typical or modeled boundary conditions (BCs) are often used. This study investigates the effects of modeled versus patient-specific BCs on modeled hemodynamics within IAs. Computational fluid dynamics (CFD) models of five IAs were reconstructed from three-dimensional rotational angiography (3DRA). BCs were applied using in turn patient-specific phase-contrast-MR (pc-MR) measurements, a 1D-circulation model, and a physiologically coherent method based on local WSS at inlets. The Navier–Stokes equations were solved using the Ansys®-CFX™ software. Wall shear stress (WSS), oscillatory shear index (OSI), and other hemodynamic indices were computed. Differences in the values obtained with the three methods were analyzed using boxplot diagrams. Qualitative similarities were observed in the flow fields obtained with the three approaches. The quantitative comparison showed smaller discrepancies between pc-MR and 1D-model data, than those observed between pc-MR and WSS-scaled data. Discrepancies were reduced when indices were normalized to mean hemodynamic aneurysmal data. The strong similarities observed for the three BCs models suggest that vessel and aneurysm geometry have the strongest influence on aneurysmal hemodynamics. In absence of patient-specific BCs, a distributed circulation model may represent the best option when CFD is used for large cohort studies.

How do coil configuration and packing density influence intra-aneurysmal hemodynamics?

BACKGROUND AND PURPOSE: Endovascular coiling is a well-established therapy for treating intracranial aneurysms. Nonetheless, postoperative hemodynamic changes induced by this therapy remain not fully understood. The purpose of this work is to assess the influence of coil configuration and packing density on intra-aneurysmal hemodynamics. MATERIALS AND METHODS: Three 3D rotational angiography images of 3 intracranial aneurysms before and after endovascular coiling were used. For each aneurysm, a 3D representation of the vasculature was obtained after the segmentation of the images. Afterward, a virtual coiling technique was used to treat the aneurysm geometries with coil models. The aneurysms were coiled with 5 packing densities, and each was generated by using 3 coil configurations. Computational fluid dynamics analyses were carried out in both untreated and treated aneurysm geometries. Statistical tests were performed to evaluate the relative effect of coil configuration on local hemodynamics. RESULTS: The intra-aneurysmal blood flow velocity and wall shear stress were diminished as packing density increased. Aneurysmal flow velocity was reduced >50% due to the first inserted coils (packing density <12%) but with a high dependency on coil configuration. Nonsignificant differences (P > .01) were found in the hemodynamics due to coil configuration for high packing densities (near 30%). A damping effect was observed on the intra-aneurysmal blood flow waveform after coiling. CONCLUSIONS: Intra-aneurysmal hemodynamics are altered by coils. Coil configuration might reduce its influence on intra-aneurysmal hemodynamics as the packing density increases until an insignificant influence could be achieved for high packing densities.

Patient-specific simulations of stenting procedures in coronary bifurcations: Two clinical cases

Computational simulations of stenting procedures in idealized geometries can only provide general guidelines and their use in the patient-specific planning of percutaneous treatments is inadequate. Conversely, image-based patient-specific tools that are able to realistically simulate different interventional options might facilitate clinical decision-making and provide useful insights on the treatment for each individual patient. The aim of this work is the implementation of a patient-specific model that uses image-based reconstructions of coronary bifurcations and is able to replicate real stenting procedures following clinical indications. Two clinical cases are investigated focusing the attention on the open problems of coronary bifurcations and their main treatment, the provisional side branch approach. Image-based reconstructions are created combining the information from conventional coronary angiography and computed tomography angiography while structural finite element models are implemented to replicate the real procedure performed in the patients. First, numerical results show the biomechanical influence of stents deployment in the coronary bifurcations during and after the procedures. In particular, the straightening of the arterial wall and the influence of two overlapping stents on stress fields are investigated here. Results show that a sensible decrease of the vessel tortuosity occurs after stent implantation and that overlapping devices result in an increased stress state of both the artery and the stents. Lastly, the comparison between numerical and image-based post-stenting configurations proved the reliability of such models while replicating stent deployment in coronary arteries.

Computational fluid dynamic simulations of image-based stented coronary bifurcation models.

One of the relevant phenomenon associated with in-stent restenosis in coronary arteries is an altered haemodynamics in the stented region. Computational fluid dynamics (CFD) offers the possibility to investigate the haemodynamics at a level of detail not always accessible within experimental techniques. CFD can quantify and correlate the local haemodynamics structures which might lead to in-stent restenosis. The aim of this work is to study the fluid dynamics of realistic stented coronary artery models which replicate the complete clinical procedure of stent implantation. Two cases of pathologic left anterior descending coronary arteries with their bifurcations are reconstructed from computed tomography angiography and conventional coronary angiography images. Results of wall shear stress and relative residence time show that the wall regions more prone to the risk of restenosis are located next to stent struts, to the bifurcations and to the stent overlapping zone for both investigated cases. Considering a bulk flow analysis, helical flow structures are generated by the curvature of the zone upstream from the stent and by the bifurcation regions. Helical recirculating microstructures are also visible downstream from the stent struts. This study demonstrates the feasibility to virtually investigate the haemodynamics of patient-specific coronary bifurcation geometries.

Intra-Aneurysmal Pressure and Flow Changes Induced by Flow Diverters: Relation to Aneurysm Size and Shape

BACKGROUND AND PURPOSE: Effects of blood flow modification by flow diverters are observed to lead often to aneurysm thrombosis and reverse remodeling. For this process, to further understand the potential roles of intra-aneurysmal blood pressure changes and aneurysm morphologies, 23 patients were studied by numeric simulation. MATERIALS AND METHODS: 3D imaging of aneurysms of different sizes and shapes, all located at the supraclinoid segment of the ICA (n = 23), was prepared for CFD simulations. Hemodynamic variables were calculated for conditions before and after virtual FD implantation, reconstituting a vessel wall scaffold across the aneurysm neck. WSS, velocity, residence time, turnover time, and intra-aneurysmal pressure were assessed statistically. RESULTS: After placement of FDs, significant reductions inside the aneurysm were observed for most hemodynamic variables (P < .01) except mean intra-aneurysmal pressures. For minimum/maximum intra-aneurysmal pressure values, small but significant changes were found; however, they were considered too small to be of relevance. CONCLUSIONS: Calculations in 23 cases did not reveal significant intra-aneurysmal mean or peak pressure changes, indicating a minor role of pressure changes in the rare event of secondary ruptures after FD use. Other hemodynamic variables (WSS and velocity) exhibited more significant changes, indicating their role in intra-aneurysmal thrombus formation. Size-dependent, significantly higher reduction in WSS (P = .069) and velocity (P = .013) was observed in small aneurysms compared with larger ones. When it came to shape, there were significantly higher reductions in WSS (P = .055) and velocity (P = .065) and a significantly higher increase in turnover time in fusiform aneurysms compared with saccular aneurysms.

Patient-specific computational hemodynamics of intracranial aneurysms from 3D rotational angiography and CT angiography: an in vivo reproducibility study

BACKGROUND AND PURPOSE: Patient-specific simulations of the hemodynamics in intracranial aneurysms can be constructed by using image-based vascular models and CFD techniques. This work evaluates the impact of the choice of imaging technique on these simulations. MATERIALS AND METHODS: Ten aneurysms, imaged with 3DRA and CTA, were analyzed to assess the reproducibility of geometric and hemodynamic variables across the 2 modalities. RESULTS: Compared with 3DRA models, we found that CTA models often had larger aneurysm necks (P = .05) and that most of the smallest vessels (between 0.7 and 1.0 mm in diameter) could not be reconstructed successfully with CTA. With respect to the values measured in the 3DRA models, the flow rate differed by 14.1 ± 2.8% (mean ± SE) just proximal to the aneurysm and 33.9 ± 7.6% at the aneurysm neck. The mean WSS on the aneurysm differed by 44.2 ± 6.0%. Even when normalized to the parent vessel WSS, a difference of 31.4 ± 9.9% remained, with the normalized WSS in most cases being larger in the CTA model (P = .04). Despite these substantial differences, excellent agreement (κ ≥ 0.9) was found for qualitative variables that describe the flow field, such as the structure of the flow pattern and the flow complexity. CONCLUSIONS: Although relatively large differences were found for all evaluated quantitative hemodynamic variables, the main flow characteristics were reproduced across imaging modalities.

neurIST complex information processing toolchain for the integrated management of cerebral aneurysms

Cerebral aneurysms are a multi-factorial disease with severe consequences. A core part of the European project @neurIST was the physical characterization of aneurysms to find candidate risk factors associated with aneurysm rupture. The project investigated measures based on morphological, haemodynamic and aneurysm wall structure analyses for more than 300 cases of ruptured and unruptured aneurysms, extracting descriptors suitable for statistical studies. This paper deals with the unique challenges associated with this task, and the implemented solutions. The consistency of results required by the subsequent statistical analyses, given the heterogeneous image data sources and multiple human operators, was met by a highly automated toolchain combined with training. A testimonial of the successful automation is the positive evaluation of the toolchain by over 260 clinicians during various hands-on workshops. The specification of the analyses required thorough investigations of modelling and processing choices, discussed in a detailed analysis protocol. Finally, an abstract data model governing the management of the simulation-related data provides a framework for data provenance and supports future use of data and toolchain. This is achieved by enabling the easy modification of the modelling approaches and solution details through abstract problem descriptions, removing the need of repetition of manual processing work.

Newtonian and non-Newtonian blood flow in coiled cerebral aneurysms

Endovascular coiling aims to isolate the aneurysm from blood circulation by altering hemodynamics inside the aneurysm and triggering blood coagulation. Computational fluid dynamics (CFD) techniques have the potential to predict the post-operative hemodynamics and to investigate the complex interaction between blood flow and coils. The purpose of this work is to study the influence of blood viscosity on hemodynamics in coiled aneurysms. Three image-based aneurysm models were used. Each case was virtually coiled with a packing density of around 30%. CFD simulations were performed in coiled and untreated aneurysm geometries using a Newtonian and a Non-Newtonian fluid models. Newtonian fluid slightly overestimates the intra-aneurysmal velocity inside the aneurysm before and after coiling. There were numerical differences between fluid models on velocity magnitudes in coiled simulations. Moreover, the non-Newtonian fluid model produces high viscosity (>0.007 [Pas]) at aneurysm fundus after coiling. Nonetheless, these local differences and high-viscous regions were not sufficient to alter the main flow patterns and velocity magnitudes before and after coiling. To evaluate the influence of coiling on intra-aneurysmal hemodynamics, the assumption of a Newtonian fluid can be used.

Influence of different computational approaches for stent deployment on cerebral aneurysm haemodynamics

Cerebral aneurysms are abnormal focal dilatations of artery walls. The interest in virtual tools to help clinicians to value the effectiveness of different procedures for cerebral aneurysm treatment is constantly growing. This study is focused on the analysis of the influence of different stent deployment approaches on intra-aneurysmal haemodynamics using computational fluid dynamics (CFD). A self-expanding stent was deployed in an idealized aneurysmatic cerebral vessel in two initial positions. Different cases characterized by a progression of simplifications on stent modelling (geometry and material) and vessel material properties were set up, using finite element and fast virtual stenting methods. Then, CFD analysis was performed for untreated and stented vessels. Haemodynamic parameters were analysed qualitatively and quantitatively, comparing the cases and the two initial positions. All the cases predicted a reduction of average wall shear stress and average velocity of almost 50 per cent after stent deployment for both initial positions. Results highlighted that, although some differences in calculated parameters existed across the cases based on the modelling simplifications, all the approaches described the most important effects on intra-aneurysmal haemodynamics. Hence, simpler and faster modelling approaches could be included in clinical workflow and, despite the adopted simplifications, support clinicians in the treatment planning.