Phantom-based experimental validation of fast virtual deployment of self-expandable stents for cerebral aneurysms

Evaluation of Two Fast Virtual Stenting Algorithms for Intracranial Aneurysm Flow Diversion

BACKGROUND

Endovascular treatment of intracranial aneurysms (IAs) by flow diverter (FD) stents depends on flow modification. Patient-specific modeling of FD deployment and computational fluid dynamics (CFD) could enable a priori endovascular strategy optimization. We developed a fast, simplistic, expansion-free balls-weeping algorithm to model FDs in patientspecific aneurysm geometry. However, since such strong simplification could result in less accurate simulations, we also developed a fast virtual stenting workflow (VSW) that explicitly models stent expansion using pseudo-physical forces.

METHODS

To test which of these two fast algorithms more accurately simulates real FDs, we applied them to virtually treat three representative patient-specific IAs. We deployed Pipeline Embolization Device into 3 patient-specific silicone aneurysm phantoms and simulated the treatments using both balls-weeping and VSW algorithms in computational aneurysm models. We then compared the virtually deployed FD stents against experimental results in terms of geometry and post-treatment flow fields. For stent geometry, we evaluated gross configurations and porosity. For post-treatment aneurysmal flow, we compared CFD results against experimental measurements by particle image velocimetry.

RESULTS

We found that VSW created more realistic FD deployments than balls-weeping in terms of stent geometry, porosity and pore density. In particular, balls-weeping produced unrealistic FD bulging at the aneurysm neck, and this artifact drastically increased with neck size. Both FD deployment methods resulted in similar flow patterns, but the VSW had less error in flow velocity and inflow rate.

CONCLUSION

In conclusion, modeling stent expansion is critical for preventing unrealistic bulging effects and thus should be considered in virtual FD deployment algorithms. Also endowed with its high computational efficiency and superior accuracy, the VSW algorithm is a better candidate for implementation into a bedside clinical tool for FD deployment simulation.

Imbalanced flow changes of distal arteries: An important factor in process of delayed ipsilateral parenchymal hemorrhage after flow diversion in patients with cerebral aneurysms

BACKGROUND AND OBJECTIVE

Hemodynamic forces may play a role in symptomatic delayed ipsilateral parenchymal hemorrhage (DIPH) of intracranial aneurysm (IA) after flow diverter placement. We aimed to investigate the hemodynamic risk factors in the postsurgical DIPH process.

METHODS

Six patients with internal carotid artery (ICA) aneurysm developed to DIPH and 12 patients without DIPH (1:2 matched controls) after flow diverter were included between January 2015 to January 2019. Postsurgical hemodynamics of distal arteries (terminal ICA, middle cerebral artery (MCA), anterior cerebral artery (ACA)) were investigated using computational fluid dynamics, as well as the hemodynamic alteration between pre- and post-treatment. The DIPH related and unrelated distal arteries (either MCA or ACA) were discriminated and compared. Definition of imbalance index is the difference in increased velocity post-flow diverter between MCA and ACA and was used to evaluate the blood flow distribution of distal arteries.

RESULTS

The mean and maximum flow velocities in the terminal ICA increased significantly after treatment in both groups. In DIPH group, the increase rate of mean velocity in the DIPH-related artery was significantly higher than that in DIPH-unrelated artery after the treatment (20.98 ± 15.38% vs -6.40 ± 7.74%; p = 0.028). Between the DIPH and control group, the baseline characteristics were well matched. However, a higher imbalance index of mean velocity was found in DIPH group (27.38 ± 13.03% vs 10.85 ± 14.12%; p = 0.031).

CONCLUSION

The mean velocity of DIPH related artery increased more, and the imbalance in increased blood flow distribution of distal arteries might play an important role in DIPH after flow diverter of IAs.

Significant flow velocity reduction at the intracranial aneurysm neck after endovascular treatment leads to favourable angiographic outcome: a prospective study

BACKGROUND

With widely usage of flow diverter in intracranial aneurysm treatment, some previously used predictors may not be effective in evaluating the recurrence risk. We aimed to comprehensively re-evaluate the predictors of intracranial aneurysm outcome with various endovascular treatment methods and devices.

METHODS

This is a prospective registered study. We analysed 6-month and 18-month follow-up angiographic data from the prospective study. Data on patient demographics, aneurysm morphology and type of treatment were recorded. Patient-specific haemodynamic simulations were performed. An unfavourable angiographic outcome was defined as recurrence of aneurysm in cases with coiling or stent-assisted coiling, patency of aneurysm in cases with flow diverters or retreatment during follow-up.

RESULTS

In total, 165 patients (177 intracranial aneurysms) with at least one angiographic follow-up data were analysed. For the short-term (6-month) results, after univariate analysis, the demographic, morphological and treatment-related factors did not achieve significantly statistical differences. The reduction ratio (RR) of velocity at aneurysm neck after embolisation was significantly lower in the unfavourable angiographic group than the favourable angiographic outcome group (p=0.002). After the Cox regression analysis, the RR of velocity at aneurysm neck was the only independent factor associated with favourable angiographic outcome (OR 0.028; p=0.001) and had an acceptable area under the curve (0.714) with a clear cut-off value (46.14%). Similarly, for the analysis of midterm (18-month) results, the RR of velocity at the aneurysm neck was the only independent significant factor for the unfavourable angiographic outcome (OR 0.050; p=0.017). The area under the curve was 0.754 and the cut-off value was 48.20%.

CONCLUSIONS

The haemodynamics showed an independent effect on angiographic follow-up results and may provide helpful suggestions for clinical practice in the future.

What does computational fluid dynamics tell us about intracranial aneurysms? A meta-analysis and critical review

Despite the plethora of published studies on intracranial aneurysms (IAs) hemodynamic using computational fluid dynamics (CFD), limited progress has been made towards understanding the complex physics and biology underlying IA pathophysiology. Guided by 1733 published papers, we review and discuss the contemporary IA hemodynamics paradigm established through two decades of IA CFD simulations. We have traced the historical origins of simplified CFD models which impede the progress of comprehending IA pathology. We also delve into the debate concerning the Newtonian fluid assumption used to represent blood flow computationally. We evidently demonstrate that the Newtonian assumption, used in almost 90% of studies, might be insufficient to describe IA hemodynamics. In addition, some fundamental properties of the Navier-Stokes equation are revisited in supplementary material to highlight some widely spread misconceptions regarding wall shear stress (WSS) and its derivatives. Conclusively, our study draws a roadmap for next-generation IA CFD models to help researchers investigate the pathophysiology of IAs.

Hemodynamics of Cerebral Aneurysms: Connecting Medical Imaging and Biomechanical Analysis

In the last two decades, numerous studies have conducted patient-specific computations of blood flow dynamics in cerebral aneurysms and reported correlations between various hemodynamic metrics and aneurysmal disease progression or treatment outcomes. Nevertheless, intra-aneurysmal flow analysis has not been adopted in current clinical practice, and hemodynamic factors usually are not considered in clinical decision making. This review presents the state of the art in cerebral aneurysm imaging and image-based modeling, discussing the advantages and limitations of each approach and focusing on the translational value of hemodynamic analysis. Combining imaging and modeling data obtained from different flow modalities can improve the accuracy and fidelity of resulting velocity fields and flow-derived factors that are thought to affect aneurysmal disease progression. It is expected that predictive models utilizing hemodynamic factors in combination with patient medical history and morphological data will outperform current risk scores and treatment guidelines. Possible future directions include novel approaches enabling data assimilation and multimodality analysis of cerebral aneurysm hemodynamics.

Relationship between haemodynamic changes and outcomes of intracranial aneurysms after implantation of the pipeline embolisation device: a single centre study

BACKGROUND

Intracranial aneurysms are increasingly being treated by the placement of flow diverters; however, the factors affecting the outcome of aneurysms treated using flow diverters remain unclarified.

METHODS

The present study investigated 94 aneurysms treated with pipeline embolisation device placement, and used a computational fluid dynamics method to explore the factors influencing the outcome of aneurysms.

RESULTS

Seventy-six completely occluded aneurysms and 18 incompletely occluded aneurysms were analysed. Before treatment, inflow jets were found in 13 (72.2%) aneurysms in the incompletely occluded group and 34 (44.7%) in the completely occluded group (P = 0.292). After deployment of the pipeline embolisation device, inflow jets remained in nine (50%) aneurysms in the incompletely occluded group and nine (11.8%) in the completely occluded group (P = 0.001). In the incompletely occluded group, regions with inflow jets after treatment corresponded with the patent areas shown on follow-up digital subtraction angiography. The mean reduction ratios of velocity in the whole aneurysm and on the neck plane were lower in the incompletely occluded than in the completely occluded group (P = 0.003; P = 0.017). Multivariate analysis revealed that the only independent risk factors for incomplete aneurysm occlusion were the reduction ratios of velocity (in the whole aneurysm, threshold 0.362, P = 0.005; on the neck plane, threshold 0.273, P = 0.015).

CONCLUSIONS

After pipeline embolisation device placement, reduction ratios of velocity in the whole aneurysm of less than 0.362 and on the neck plane of less than 0.273 are significantly associated with a greater risk of aneurysm incomplete occlusion. In addition, the persistence of inflow jets in aneurysms is associated with incomplete occlusion in the inflow jet area.

COMPUTATIONAL INVESTIGATION OF THROMBIN CONCENTRATION IN CEREBRAL ANEURYSMS TREATED WITH FLOW-DIVERTING STENTS

Flow-diverting stent is an ongoing embolization device to treat cerebral aneurysms, and it diverts the flow direction to reduce the flow velocity inside the aneurysmal sacs and promote the thrombus formation. However, its effect for aneurysm embolization is controversial. A hemodynamic-biomedical coupling model was constructed to describe the generation and transport of thrombin in arteries, and the model was applied to investigate the variation of thrombin concentration, which plays a key role in thrombus formation, in two patient-specific cerebral aneurysm models when they are treated with Pipeline flow diverting stents. It is observed from computational fluid dynamics simulations that thrombin concentration in the aneurysmal sac without collateral artery increases significantly after Pipeline implantation, however, it has hardly any variation in the aneurysmal sac without collateral artery or in the giant aneurysmal sac after Pipeline implantation. Therefore, we believe that single Pipeline is very effective to embolize a small aneurysm without collateral artery, but cannot embolize a giant aneurysm or a small aneurysm with a collateral artery on its sac effectively.

Virtual stenting of intracranial aneurysms: A pilot study for the prediction of treatment success based on hemodynamic simulations

Endovascular treatment of intracranial aneurysms using flow-diverting devices has revolutionized the treatment of large and complex lesions due to its minimally invasive nature and potential clinical outcomes. However, incomplete or delayed occlusion and persistent intracranial aneurysm growth are still an issue for up to one-third of the patients. We evaluated two patients with intracranial aneurysm located at the internal carotid artery who were treated with flow-diverting devices and had opposite outcomes. Both patients presented with similar aneurysms and were treated with the same device, but after a 1-year follow-up, one case presented with complete occlusion (Case 1) and the other required further treatment (Case 2). To reproduce the interventions, virtual stents were deployed and blood flow simulations were carried out using the respective patient-specific geometries. Afterward, hemodynamic metrics such as aneurysmal inflow reduction, wall shear stresses, oscillatory shear, and inflow concentration indices were quantified. The hemodynamic simulations reveal that for both cases, the neck inflow was clearly reduced due to the therapy (Case 1: 19%, Case 2: 35%). In addition, relevant hemodynamic parameters such as time-averaged wall shear stress (Case 1: 35.6%, Case 2: 57%) and oscillatory shear (Case 1: 33.1%, Case 2: 26.7%) were decreased considerably. However, although stronger relative reductions occurred in the unsuccessful case, the absolute flow values in the successful case were approximately halved. The findings demonstrate that a high relative effect of endovascular devices is not necessarily associated with the desired treatment outcome. Instead, it appears that a successful intracranial aneurysm therapy requires a certain patient-specific inflow threshold.

Efficient simulation of a low-profile visualized intraluminal support device: a novel fast virtual stenting technique

Background

The low-profile visualized intraluminal support (LVIS) stent has become a promising endovascular option for treating intracranial aneurysms. To achieve better treatment of aneurysms using LVIS, we developed a fast virtual stenting technique for use with LVIS (F-LVIS) to evaluate hemodynamic changes in the aneurysm and validate its reliability.

Methods

A patient-specific aneurysm was selected for making comparisons between the real LVIS (R-LVIS) and the F-LVIS. To perform R-LVIS stenting, a hollow phantom based on a patient-specific aneurysm was fabricated using a three-dimensional printer. An R-LVIS was released in the phantom according to standard procedure. F-LVIS was then applied successfully in this aneurysm model. The computational fluid dynamics (CFD) values were calculated for both the F-LVIS and R-LVIS models. Qualitative and quantitative comparisons of the two models focused on hemodynamic parameters.

Results

The hemodynamic characteristics for R-LVIS and F-LVIS were well matched. Representative contours of velocities and wall shear stress (WSS) were consistently similar in both distribution and magnitude. The velocity vectors also showed high similarity, although the R-LVIS model showed faster and more fluid streams entering the aneurysm. Variation tendencies of the velocity in the aneurysm and the WSS on the aneurysm wall were also similar in the two models, with no statistically significant differences in either velocity or WSS.

Conclusions

The results of the computational hemodynamics indicate that F-LVIS is suitable for evaluating hemodynamic factors. This novel F-LVIS is considered efficient, practical, and effective.

Predisposing factors for recanalization of cerebral aneurysms after endovascular embolization: a multivariate study

BACKGROUND

The recanalization of cerebral aneurysms after endovascular embolization (coiling or stent-assisted coiling) has been a matter of concern.

OBJECTIVE

To systematically evaluate the predisposing factors for cerebral aneurysm recanalization using multidimensional analysis in a large patient cohort.

METHODS

In 238 patients with 283 aneurysms, patient baseline characteristics, aneurysm morphological characteristics, treatment-related factors, and changes in flow hemodynamics after endovascular treatment (coiling or stent-assisted coiling) were compared between the recanalization and non-recanalization groups. Multivariate logistic regression analysis was performed to determine independent risk factors correlated with recanalization.

RESULTS

16 aneurysms treated by coiling recanalized, with a recurrence rate of 18.6%, and 24 recanalized in the lesions treated by stent-assisted coiling, with a recanalization rate of 12.2%. Large aneurysms (>10 mm, p=0.002) and a follow-up interval >1 year (p=0.027) were shown to be statistically significant between the recanalization and non-recanalization groups. For flow hemodynamic changes, three parameters (velocity on the neck plane, wall shear stress on the neck wall, and wall shear stress on the whole aneurysm) showed a relatively lower amplitude of decrease after endovascular treatment in the recanalization group. Interestingly, the velocity on the neck plane and wall shear stress on the neck wall may be elevated after treatment. Specifically, the reduction ratio (RR) of velocity on the neck plane showed significant difference between the groups in the multivariate analysis (p=0.013), and was considered an independent risk factor for recanalization.

CONCLUSIONS

The aneurysm size, follow-up interval, and flow hemodynamic changes, especially the RR of velocity on the neck plane, have important roles in aneurysm recanalization.

Association between hemodynamic modifications and clinical outcome of intracranial aneurysms treated using flow diverters

Treatment of intracranial aneurysms (IAs) has been revolutionized by the advent of endovascular Flow Diverters (FDs), which disrupt blood flow within the aneurysm to induce pro-thrombotic conditions, and serves as a scaffold for endothelial ingrowth and arterial remodeling. Despite good clinical success of FDs, complications like incomplete occlusion and post-treatment rupture leading to subarachnoid hemorrhage have been reported. In silico computational fluid dynamic analysis of the pre- and post-treated geometries of IA patients can shed light on the contrasting blood hemodynamics associated with different clinical outcomes. In this study, we analyzed hemodynamic modifications in 15 IA patients treated using a single FD; 10 IAs were completely occluded (successful) and 5 were partially occluded (unsuccessful) at 12-month follow-up. An in-house virtual stenting workflow was used to recapitulate the clinical intervention on these cases, followed by CFD to obtain pre- and post-treatment hemodynamics. Bulk hemodynamic parameters showed comparable reductions in both groups with average inflow rate and aneurysmal velocity reduction of 40.3% and 52.4% in successful cases, and 34.4% and 49.2% in unsuccessful cases. There was a substantial reduction in localized parameter like vortex coreline length and Energy Loss for successful cases, 38.2% and 42.9% compared to 10.1% and 10.5% for unsuccessful cases. This suggest that for successfully treated IAs, the localized complex blood flow is disrupted more prominently by the FD as compared to unsuccessful cases. These localized hemodynamic parameters can be potentially used in prediction of treatment outcome, thus aiding the clinicians in a priori assessment of different treatment strategies.