Finite element modeling of endovascular coiling and flow diversion enables hemodynamic prediction of complex treatment strategies for intracranial aneurysm

Optimizing distal and proximal splenic artery embolization with patient-specific computational fluid dynamics

Splenic artery embolization (SAE) has become a favored alternative to splenectomy, offering a less invasive intervention for injured spleens while preserving spleen function. However, our understanding of the role that hemodynamics plays during embolization remains limited. In this study, we utilized patient-specific computational fluid dynamics (CFD) simulations to study distal and proximal embolization strategies commonly used in SAE. Detailed 3D computer models were constructed considering the descending aorta, various major visceral arteries, and the iliac arteries. Subsequently, the blood flow and pressure associated with different coil placement locations in proximal embolization were studied considering the collateral vessels. Coil induced variations in pressure fields were quantified and compared to baseline. The coil induced flow stagnation was also quantified with particle residence time. Distal embolization was modeled with Lagrangian particle tracking and the effect of particle size, release location, and timing on embolization outcome was studied. Our findings highlight the crucial role of collateral vessels in maintaining blood supply to the spleen following proximal embolization. It was demonstrated that coil location can affect distal pressure and that strategic coil placement guided by patient-specific CFD simulations can further reduce this pressure as desired. Additionally, the results point to the critical roles that particle size, release timing, and location play in distal embolization. Our study provides an early attempt to use patient-specific computer modeling for optimizing embolization strategies and ultimately improving patient outcomes during SAE procedures.

In vitro and in silico assessment of flow modulation after deploying the Contour Neurovascular System in intracranial aneurysm models

BACKGROUND

The novel Contour Neurovascular System (Contour) has been reported to be efficient and safe for the treatment of intracranial, wide-necked bifurcation aneurysms. Flow in the aneurysm and posterior cerebral arteries (PCAs) after Contour deployment has not been analyzed in detail yet. However, this information is crucial for predicting aneurysm treatment outcomes.

METHODS

Time-resolved three-dimensional velocity maps in 14 combinations of patient-based basilar tip aneurysm models with and without Contour devices (sizes between 5 and 14 mm) were analyzed using four-dimensionsal (4D) flow MRI and numerical/image-based flow simulations. A complex virtual processing pipeline was developed to mimic the experimental shape and position of the Contour together with the simulations.

RESULTS

On average, the Contour significantly reduced intra-aneurysmal flow velocity by 67% (mean w/ = 0.03m/s; mean w/o = 0.12m/s; p-value=0.002), and the time-averaged wall shear stress by more than 87% (mean w/ = 0.17Pa; mean w/o = 1.35Pa; p-value=0.002), as observed by numerical simulations. Furthermore, a significant reduction in flow (P<0.01) was confirmed by the neck inflow rate, kinetic energy, and inflow concentration index after Contour deployment. Notably, device size has a stronger effect on reducing flow than device positioning. However, positioning affected flow in the PCAs, while being robust in effectively reducing flow.

CONCLUSIONS

This study showed the high efficacy of the Contour device in reducing flow within aneurysms regardless of the exact position. However, we observed an effect on the flow in PCAs, which needs to be investigated further.

Influence of framing coil orientation and its shape on the hemodynamics of a basilar aneurysm model

Aneurysms are bulges of an artery, which require clinical management solutions. Due to the inherent advantages, endovascular coil filling is emerging as the treatment of choice for intracranial aneurysms (IAs). However, after successful treatment of coil embolization, there is a serious risk of recurrence. It is well known that optimal packing density will enhance treatment outcomes. The main objective of endovascular coil embolization is to achieve flow stasis by enabling significant reduction in intra-aneurysmal flow and facilitate thrombus formation. The present study numerically investigates the effect of framing coil orientation on intra-aneurysmal hemodynamics. For the purpose of analysis, actual shape of the embolic coil is used, instead of simplified ideal coil shape. Typically used details of the framing coil are resolved for the analysis. However, region above the framing coil is assumed to be filled with a porous medium. Present simulations have shown that orientation of the framing coil loop (FCL) greatly influences the intra-aneurysmal hemodynamics. The FCLs which were placed parallel to the outlets of basilar tip aneurysm (Coil A) were found to reduce intra-aneurysmal flow velocity that facilitates thrombus formation. Involving the coil for the region is modeled using a porous medium model with a packing density of 20 % . The simulations indicate that the framing coil loop (FCL) has a significant influence on the overall outcome.

Hemodynamic study on the therapeutic effects of varying diameter embolic coils in the treatment of intracranial aneurysms

Endovascular coiling is the predominant method for treating cerebral aneurysms. Extensive reports on selecting coil length, hardness, and material are available. However, the impact of coil diameter on postoperative outcomes remains unclear. This study enrolled six personalized geometric models of intracranial aneurysms: three bifurcation aneurysms and three sidewall aneurysms. Four coil models were constructed by changing the coil diameter. Coil embolization was simulated using the finite element method. Computational fluid dynamics was used to characterize hemodynamics in the aneurysms after embolization. Evaluation parameters included velocity reduction, wall shear stress (WSS), low WSS (LWSS), oscillatory shear index (OSI), relative residence time (RRT), and residual flow volume in the aneurysms. At the peak time (t = 0.17 s), the proportion of LWSS area in bifurcation aneurysms increase with the rise in coil diameter: 0.8D, 71.28 ± 12.62% versus 1D, 74.97 ± 19.17% versus 1.2D, 78.88 ± 18.56% versus 1.4D, 84.00 ± 11.53% (mean ± SD). The proportion of high OSI area decreases as the coil diameter increases: 0.8D, 4.41% ± 2.82% versus 1.0D, 3.78 ± 3.33% versus 1.2D, 2.28% ± 1.77% versus 1.4D, 1.58% ± 1.11% (mean ± SD). The proportion of high RRT area increases as the coil diameter rises: 0.8D, 3.40% ± 1.68% versus 1.0D, 7.67 ± 4.12% versus 1.2D, 9.84% ± 9.50% versus 1.4D, 22.29% ± 14.28% (mean ± SD). Side wall aneurysms do not exhibit the aforementioned trend. Bifurcation aneurysms plugged with a coil of 1.4 times the diameter have the largest RFVs (<10 mm/s) within the group. Aforementioned patterns are not found in sidewall aneurysms. In the treatment of aneurysms with coiling, varying coil diameters can result in different hemodynamic environments within the aneurysm. Larger coil diameters have improved hemodynamic performance for bifurcation aneurysms. However, coil diameter and embolization effectiveness have no significant relationship for sidewall aneurysms.

Hemodynamic alterations of flow diverters on aneurysms at the fetal posterior communicating artery: A simulation study using CFD to compare the surpass streamline, pipeline flex, and tubridge devices

PURPOSE

Traditional flow diverters (FDs) for treating aneurysms at the fetal posterior communicating artery are unsatisfactory. Surpass Streamline is a novel FD with different mesh characteristics; however, the outcomes for such aneurysms remain unclear. This study aimed to compare hemodynamic alterations induced by Surpass Streamline, Pipeline Flex, and Tubridge devices and explore possible strategies for aneurysms at the fetal posterior communicating artery.

METHODS

Two simulated aneurysms (Case 1, Case 2) were constructed from digital subtraction angiography (DSA). The three FDs were virtually deployed, and hemodynamic analysis based on computational fluid dynamics was performed. Hemodynamic parameters, including the sac-averaged velocity magnitude (Velocity), high-flow volume (HFV), and wall shear stress (WSS), were compared between each FD and the untreated model (control). Surpass Streamline was performed in real life for two aneurysms and the clinical outcomes were collected for analysis.

RESULTS

Compared to the control, the Surpass resulted in the most significant reduction in flow. In Case 1, the Velocity, HFV, and WSS were reduced by 51.6%, 78.1%, and 64.3%, respectively. In Case 2, the Velocity, HFV, and WSS were reduced by 48.0%, 81.1%, and 65.3%, respectively. Tubridge showed slightly larger changes in hemodynamic parameters than Pipeline. In addition, our analysis suggested that metal coverage was correlated with the WSS, Velocity, and HFV. The postoperative DSA showed that the aneurysm was nearly occluded in Case 1 and decreased in Case 2.

CONCLUSION

Compared to that with the Pipeline and Tubridge, the Surpass resulted in the greatest reduction in hemodynamic parameters and might be effective for aneurysms at the fetal posterior communicating artery. Virtual FD deployment and computational fluid dynamics analysis may be used to predict the treatment outcomes.

Numerical Study of a Thrombus Migration Risk in Aneurysm After Coil Embolization in Patient Cases: FSI Modelling

PURPOSE

There are still many challenges for modelling a thrombus migration process in aneurysms. The main novelty of the present research lies in the modelling of aneurysm clot migration process in a realistic cerebral aneurysm, and the analysis of forces suffered by clots inside an aneurysm, through transient FSI simulations.

METHODS

The blood flow has been modelled using a Womersley velocity profile, and following the Carreau viscosity model. Hyperelastic Ogden model has been used for clot and isotropic linear elastic model for the artery walls. The FSI coupled model was implemented in ANSYS® software. The hemodynamic forces suffered by the clot have been quantified using eight different clot sizes and positions inside a real aneurysm.

RESULTS

The obtained results have shown that it is almost impossible for clots adjacent to aneurysm walls, to leave the aneurysm. Nevertheless, in clots positioned in the centre of the aneurysm, there is a real risk of clot migration. The risk of migration of a typical post-coiling intervention clot in an aneurysm, in contact with the wall and occupying a significant percentage of its volume is very low in the case studied, even in the presence of abnormally intense events, associated with sneezes or impacts.

CONCLUSIONS

The proposed methodology allows evaluating the clot migration risk, vital for evaluating the progress after endovascular interventions, it is a step forward in the personalized medicine, patient follow-up, and helping the medical team deciding the optimal treatment.

Computational fluid dynamics (CFD) analysis in a ruptured vertebral artery dissecting aneurysm implanted by Pipeline when recurrent after LVIS-assisted coiling treatment: Case report and review of the literatures

BACKGROUNDS

Hemodynamics plays an important role in the natural history of the process of rupture and recurrence of intracranial aneurysms. This study aimed to investigate the role of hemodynamics for recurrence in a vertebral artery dissecting aneurysm (VADA).

METHODS

A patient with a ruptured VADA firstly treated by low-profile visualized intraluminal support (LVIS)-assisted coiling, and was implanted with a Pipeline Embolization Device (PED) after aneurysm recurrence. Finite element analysis and computational fluid dynamics simulations were conducted in 6 serial imaging procedures, and the calculated hemodynamics was correlated with aneurysm recurrence.

RESULTS

Wall shear stress (WSS) was not effectively suppressed, resulting in aneurysm recurrence with initial entry tear to occur above the protuberance after 7 months of LVIS stent-assisted coiling. With the implantation of PED, WSS, inflow stream and velocity at the aneurysm neck significantly decreased. During the 3-month follow-up after PED deployment, there was significant shrinkage of the sac and the blood flow in the sac was reduced considerably. The 27-month follow-up after PED deployment indicated the aneurysm was stable.

CONCLUSIONS

The present case study suggests that insufficient suppression of high WSS and high inflow velocity at the neck of the parent artery, especially near the posterior inferior cerebellar artery, might be associated with aneurysm recurrence.

Improving the accuracy of computational fluid dynamics simulations of coiled cerebral aneurysms using finite element modeling

Cerebral aneurysms are a serious clinical challenge, with ∼half resulting in death or disability. Treatment via endovascular coiling significantly reduces the chances of rupture, but the techniquehas failure rates of ∼20 %. This presents a pressing need to develop a method fordetermining optimal coildeploymentstrategies. Quantification of the hemodynamics of coiled aneurysms using computational fluid dynamics (CFD) has the potential to predict post-treatment outcomes, but representing the coil mass in CFD simulations remains a challenge. We use the Finite Element Method (FEM) for simulating patient-specific coil deployment for n = 4 ICA aneurysms for which 3D printed in vitro models were also generated, coiled, and scanned using ultra-high resolution synchrotron micro-CT. The physical and virtual coil geometries were voxelized onto a binary structured grid and porosity maps were generated for geometric comparison. The average binary accuracy score is 0.8623 and the average error in porosity map is 4.94 %. We then conduct patient-specific CFD simulations of the aneurysm hemodynamics using virtual coils geometries, micro-CT generated oil geometries, and using the porous medium method to represent the coil mass. Hemodynamic parameters including Neck Inflow Rate (Qneck) and Wall Shear Stress (WSS) were calculated for each of the CFD simulations. The average relative error in Qneck and WSS from CFD using FEM geometry were 6.6 % and 21.8 % respectively, while the error from CFD using a porous media approximation resulted in errors of 55.1 % and 36.3 % respectively; demonstrating a marked improvement in the accuracy of CFD simulations using FEM generated coil geometries.

An Inverse Method to Determine Mechanical Parameters of Porcine Vitreous Bodies Based on the Indentation Test

The vitreous body keeps the lens and retina in place and protects these tissues from physical insults. Existing studies have reported that the mechanical properties of vitreous body varied after liquefaction, suggesting mechanical properties could be effective parameters to identify vitreous liquefaction process. Thus, in this work, we aimed to propose a method to determine the mechanical properties of vitreous bodies. Fresh porcine eyes were divided into three groups, including the untreated group, the 24 h liquefaction group and the 48 h liquefaction group, which was injected collagenase and then kept for 24 h or 48 h. The indentation tests were carried out on the vitreous body in its natural location while the posterior segment of the eye was fixed in the container. A finite element model of a specimen undertaking indentation was constructed to simulate the indentation test with surface tension of vitreous body considered. Using the inverse method, the mechanical parameters of the vitreous body and the surface tension coefficient were determined. For the same parameter, values were highest in the untreated group, followed by the 24 h liquefaction group and the lowest in the 48 h liquefaction group. For C₁₀ in the neo-Hookean model, the significant differences were found between the untreated group and liquefaction groups. This work quantified vitreous body mechanical properties successfully using inverse method, which provides a new method for identifying vitreous liquefactions related studies.

How to perform intra-aneurysmal coil embolization after Pipeline deployment: a study from a hemodynamic viewpoint

BACKGROUND

Pipeline embolization device (PED) deployment combined with coil therapy for large complex intracranial aneurysms is effective and considered superior to PED deployment alone. However, the optimal strategy for use of coils remains unclear. We used patient-specific aneurysm models and finite element analysis to determine the ideal packing density of coils after PED placement.

METHODS

Finite element analysis was used to provide a higher-fidelity model for accurate post-treatment computational fluid dynamics analysis to simulate the real therapeutic process of PED and all coils. We then calculated and analyzed the reduction ratio of velocity to identify the hemodynamic change during PED deployment and each coil embolization.

RESULTS

Sixteen consecutive patients underwent PED plus coil procedures to treat internal carotid artery intracranial aneurysms. After PED deployment, the intra-aneurysmal flow velocity significantly decreased (15.3 vs 10.0 cm/s; p<0.001). When the first coil was inserted, the flow velocity in the aneurysm further decreased and the reduction was significant (10.0 vs 5.3 cm/s; p<0.001). Analysis of covariance showed that the effect of the reduction ratio of velocity of the second coil was significantly lower than that of the first coil (p<0.001)-that is, when the packing density increased to 7.06%, the addition of coils produced no further hemodynamic effect.

CONCLUSION

Adjunct coiling could improve the post-PED hemodynamic environment in treated intracranial aneurysms. However, dense packing is not necessary because the intra-aneurysmal hemodynamics tend to stabilize as the packing density reaches an average of 7.06% or after insertion of the second coil.

Modeling Flow in Cerebral Aneurysm After Coils Embolization Treatment: A Realistic Patient-Specific Porous Model Approach

PURPOSE

Computational fluid dynamics (CFD) has been used to evaluate the efficiency of endovascular treatment in coiled cerebral aneurysms. The explicit geometry of the coil mass cannot typically be incorporated into CFD simulations since the coil mass cannot be reconstructed from clinical images due to its small size and beam hardening artifacts. The existing methods use imprecise porous medium representations. We propose a new porous model taking into account the porosity heterogeneity of the coils deployed in the aneurysm.

METHODS

The porosity heterogeneity of the coil mass deployed inside two patients' cerebral aneurysm phantoms is first quantified based on 3D X-ray synchrotron images. These images are also used to compute the permeability and the inertial factor arising in porous models. A new homogeneous porous model (porous crowns model), considering the coil's heterogeneity, is proposed to recreate the flow within the coiled aneurysm. Finally, the validity of the model is assessed through comparisons with coil-resolved simulations.

RESULTS

The strong porosity gradient of the coil measured close to the aneurysmal wall is well captured by the porous crowns model. The permeability and the inertial factor values involved in this model are closed to the ideal homogeneous porous model leading to a mean velocity in the aneurysmal sac similar as in the coil-resolved model.

CONCLUSION

The porous crowns model allows for an accurate description of the mean flow within the coiled cerebral aneurysm.

Endovascular treatment of challenging aneurysms with FRED Jr flow diverter stents: a single-center experience

PURPOSE

To analyze clinical safety and efficacy of flow re-direction endoluminal device (FRED) Jr flow diverter for treatment of unruptured, ruptured, or recanalyzed aneurysms.

MATERIALS AND METHODS

Between October 2019 and February 2022, 25 patients with 31 aneurysms treated with FRED Jr were included in the study. Clinical and radiological records, procedural details, and follow-up outcomes were retrospectively evaluated. Eighteen (72%) patients were female. Median age was 48.8 (age range 9-85). Mean follow-up was 21 months (6-28 months). Location of the aneurysms were as follows; 13 in middle cerebral artery (MCA), 7 in anterior cerebral artery (ACA), 4 in posterior cerebral artery (PCA), 3 in true posterior communicating artery (PCom), 2 in anterior communicating artery (ACom), 1 in superior cerebellar artery (SCA), 1 in true ophthalmic artery. Five patients (20%) presented with acute subarachnoid hemorrhage (aSAH).

RESULTS

In all procedures, FRED Jr was successfully deployed. Three true Pcom aneurysms and a true ophthalmic aneurysm were treated with FRED Jr. Three patients with two adjacent aneurysms were treated with a single FRED Jr. In two (8%) patients in-stent thrombosis occurred intraoperatively, they were treated with iv tirofiban and thrombectomy without any sequelae. Post-discharge 2 weeks later, intraparenchymal hemorrhage occurred in a patient. He was treated with surgical drainage, the clinical course was modified Rankin score (mRS) 2. Digital subtraction angiography (DSA) was performed on 16 (64%) patients with 21 (67%) aneurysms. Near complete-complete occlusion (O'Kelly-Morata grading scale (OKM C-D) was documented in 15/16 (93.7%) patients, 20/21 (95.2%) aneurysms. In nine (36%) patients, no residual filling was observed in the magnetic resonance angiography (MRA). Good clinical outcome (mRS 0-1) was achieved in 24/25 (96%) of patients.

CONCLUSION

Endovascular treatment of small cerebral aneurysms with FRED Jr is safe and effective even in complex and challenging morphologies allowing high rates of aneurysm occlusion with low periprocedural complications. Our cohort, consisting of a rate 20% acute ruptured aneurysms, is the major additive data to the published literature.

Predicting flow diverter sizing using the AneuGuideTM software: a validation study

BACKGROUND

Stent sizing remains a challenging task for flow diverter implantation because of stent foreshortening. In this study, we aimed to quantify the change in length after implantation and assess the error in length prediction using AneuGuide^TM software.

METHODS

In a retrospective cohort of 101 patients with 102 aneurysms undergoing treatment with a pipeline embolization device (PED; Covidien, Irvine, California, USA), we used AneuGuide^TM software to obtain measured lengths (ML) and calculated lengths (CL) after stent implantation. Stent elongation was defined as the ratio of ML-LL to the labeled length (LL). Simulation error was defined as the ratio of the absolute value of CL-ML to ML. The correlation and consistency between ML and LL and between ML and CL were analyzed using Pearson's correlation test and the Bland-Altman plot. Statistical significance was set at p<0.05.

RESULTS

The mean elongation of ML was 32.6% (range 26.3-109.2%). Moderate consistency was observed between LL and ML (ρ=0.74, p<0.001). With the AneuGuide^TM software, the mean simulation error was 6.6% (range 0.32-21.2%). Pearson's correlation test and the Bland-Altman plot showed a high correlation and consistency between ML and CL (ρ=0.96, p<0.001).

CONCLUSION

Labeled length provides only a low reference value for predicting the actual length of the flow diverter after implantation. The high consistency between ML and CL obtained from AneuGuide^TM software shows its great potential for the optimization of the flow diverter sizing process.

Hemodynamic mechanism of pulsatile tinnitus caused by venous diverticulum treated with coil embolization

BACKGROUND AND OBJECTIVE

Coil embolization has become a new treatment method for pulsatile tinnitus (PT) caused by sigmoid sinus diverticulum (SSD). Although this therapy has achieved good results in clinical reports, the hemodynamic mechanism of coils in the treatment of PT in SSD remained unclear.

METHODS

Finite element method (FEM) and computational fluid dynamics (CFD) were combined to explore the hemodynamic mechanism of coil embolization in SSD treatment. Three personalized geometric models of sigmoid sinus were established according to the CTA data of patients. Coil model were established by FEM, and the hemodynamic differences of SSD before and after coiling were compared by transient CFD method.

RESULTS

Velocity streamlines disappeared in the SSD after coiling. At the peak time (t1 = 0.22 s), the SSD-average velocity decreased in every patient. The average value of the decreased in three patients was 0.154 ± 0.028 m/s (mean ± SD). Wall average pressure (P_avg) also showed a decline in every patient. Average of decrements of three patients was 17.69 ± 4.91 Pa (mean ± SD). Average WSS (WSS_avg) was also reduced in every patient. The average value of WSS drop was 9.74 ± 3.02 Pa (mean ± SD). After coiling, the proportion of low-velocity region in the sigmoid sinus cortical plate dehiscence (SSCPD) area increased. Average of increments was 22.1 ± 5.36% (mean ± SD).

CONCLUSIONS

A reduction in SSD-average velocity, wall pressure, and WSS were the short-term hemodynamic mechanism of coil embolization for PT. Coil embolization increased the proportion of low-velocity region in the SSCPD area, thereby creating a hemodynamic environment that easily produced thrombus and protects blood vessels from the impact of blood flow. This phenomenon was the long-term effect of coil embolization.

Hemodynamic Effects of Stent-Induced Straightening of Parent Artery vs. Stent Struts for Intracranial Bifurcation Aneurysms

Objective

This study aims to compare the hemodynamic impact of stent-mesh and stent-induced straightening of the parent artery in intracranial bifurcation aneurysms using finite element method simulation.

Material and Methods

Three intracranial bifurcation aneurysms treated with different stent-assisted coil embolization were evaluated. Simulation using the finite element method was conducted for Solitaire, LVIS and Neuroform stents. Four models of each stent were established, including a pre-treatment baseline, stenting without parent artery straightening (presented as stent-mesh effect), no-stent with parent artery reconstruction (to reveal the straightening impact), and stenting with straightening (categorized as Models I–IV respectively). Hemodynamic characteristics of the four models for each stent were compared.

Results

In the Neuroform stent, compared with the pre-treatment model (100%), the mean WSS decreased to 82.3, 71.4, and 57.0% in Models II-IV, velocity to 88.3, 74.4, and 62.8%, and high flow volume (HFV, >0.3 m/s) to 77.7, 44.0, and 19.1%. For the LVIS stent, the mean WSS changed to 105.0, 40.2, and 39.8% in Models II to IV; velocity to 91.2, 58.1, and 52.5%, and HFV to 92.0, 56.1, and 43.9%. For the Solitaire stent, compared with the pre-treatment model (100%), the mean WSS of Models II-IV changed altered by 105.7, 42.6, and 39.4%, sac-averaged velocity changed to 111.3, 46.6, and 42.8%, and HFV 115.6, 15.1, and 13.6%.

Conclusion

The hemodynamic effect of straightening the parent artery of intracranial bifurcation aneurysms by stenting was noticeably improved over stent mesh diversion in all three stents tested. Therefore stent-induced remodeling of the parent artery appears to be the best method of decreasing recurrence in intracranial bifurcation aneurysms.

Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle

Interventional neuroradiology is characterized by engineering- and experience-driven device development with design improvements every few months. However, clinical validation of these new devices requires lengthy and expensive randomized controlled trials. This contribution proposes a machine learning-based in silico study design to evaluate new devices more quickly with a small sample size. Acute diffusion- and perfusion-weighted MRI, segmented one-week follow-up imaging, and clinical variables were available for 90 acute ischemic stroke patients. Three treatment option-specific random forest models were trained to predict the one-week follow-up lesion segmentation for (1) patients successfully recanalized using intra-arterial mechanical thrombectomy, (2) patients successfully recanalized using intravenous thrombolysis, and (3) non-recanalizing patients as an analogue for conservative treatment for each patient in the sample, independent of the true group membership. A repeated-measures analysis of the three predicted follow-up lesions for each patient revealed significantly larger lesions for the non-recanalizing group compared to the successful intravenous thrombolysis treatment group, which in turn showed significantly larger lesions compared to the successful mechanical thrombectomy treatment group (p < 0.001). A groupwise comparison of the true follow-up lesions for the three treatment options showed the same trend but did not reach statistical significance (p = 0.19). We conclude that the proposed machine learning-based in silico trial design leads to clinically feasible results and can support new efficacy studies by providing additional power and potential early intermediate results.

FSI modeling on the effect of artery-aneurysm thickness and coil embolization in patient cases

BACKGROUND AND OBJECTIVE

The attainment of a methodology to simulate the hemodynamic in patient-specific cerebral vessels with aneurysms is still a challenge. The novelty of this work is focused on the effect of coil embolization in a realistic cerebral aneurysm, according to the vessel wall thickness and aneurysm thickness, through transient FSI simulations.

METHODS

The quality of the mesh for simulations was checked with a specific mesh convergence study; and the numerical methodology was validated using numerical research data of the literature. The model was implemented in ANSYS® software. The total deformation and equivalent stress evolution in the studied cases, before and after coil embolization, were compared. More than 20 different models were employed due to different arterial wall thickness and aneurysm wall thickness combinations.

RESULTS

The obtained results have showed that deformation and stress values are highly influenced with the sac thickness. The thinner sac aneurysm thickness is, the greater deformation and stress are. The results after coil embolization process have highlighted that considering typical values of arterial wall thickness and aneurysm thickness 0.3 mm and 0.15 mm respectively, a deformation reduction around 50% and a stress reduction around 70% can be achieved.

CONCLUSIONS

The proposed methodology is a step forward in the personalized medicine, quantifying the aneurysm rupture risk reduction, and helping the medical team in the preoperative planning, or to deciding the optimal treatment.

Structural Design of Vascular Stents: A Review

Percutaneous Coronary Intervention (PCI) is currently the most conventional and effective method for clinically treating cardiovascular diseases such as atherosclerosis. Stent implantation, as one of the ways of PCI in the treatment of coronary artery diseases, has become a hot spot in scientific research with more and more patients suffering from cardiovascular diseases. However, vascular stent implanted into vessels of patients often causes complications such as In-Stent Restenosis (ISR). The vascular stent is one of the sophisticated medical devices, a reasonable structure of stent can effectively reduce the complications. In this paper, we introduce the evolution, performance evaluation standards, delivery and deployment, and manufacturing methods of vascular stents. Based on a large number of literature pieces, this paper focuses on designing structures of vascular stents in terms of “bridge (or link)” type, representative volume unit (RVE)/representative unit cell (RUC), and patient-specific stent. Finally, this paper gives an outlook on the future development of designing vascular stents.

Association Between Aneurysmal Haemodynamics and Device Microstructural Characteristics After Flow-Diversion Treatments With Dual Stents of Different Sizes: A Numerical Study

Objectives

Treating intracranial aneurysms with flow-diverting stents sometimes requires deployment of a second device. Herein we quantify the sizing effects of devices in dual-stent treatments upon the final stent microstructure and the post-treatment aneurysmal haemodynamics.

Methods

Fifteen sidewall ICA aneurysm geometries were included. Using a virtual stenting technique, we implanted either one or two stents for each aneurysm treatment considered, with each stent specified as one of two different sizes, yielding a total of two single-stent and fouir dual-stent treatment scenarios for each aneurysm. Three stent microstructural parameters and nine aneurysmal haemodynamic parameters were quantified and systematically compared across the 90 treatment scenarios.

Results

Deployment of a second stent further reduced the aneurysmal inflow rate (IR) and energy loss (EL) by, respectively, 14 ± 11% (p = 0.001) and 9 ± 12% (p = 0.056), relative to the untreated condition. Sizing effects of the earlier-deployed stent led to largest differences of 6.9% for the final IR reduction and 11.1% for the EL, whereas sizing effects from the later-deployed stent were minor (≤2.1%). The change in stent pore size was the only microstructural parameter demonstrating a strong correlation with the reduction in the post-treatment aneurysmal haemodynamics, in terms of the IR (r = 0.50, p < 0.001) and pressure drop (r = 0.63, p < 0.001).

Conclusion

Size of the earlier-deployed stent has substantial effects on the final haemodynamic outcomes after dual-stent treatment. The average pore size of stent wires at the aneurysm orifice shows promise as a potential index for predicting the efficacy of flow-diversion treatments.

The effect of Dean, Reynolds and Womersley numbers on the flow in a spherical cavity on a curved round pipe. Part 2. The haemodynamics of intracranial aneurysms treated with flow-diverting stents

The flow in a spherical cavity on a curved round pipe is a canonical flow that describes well the flow inside a sidewall aneurysm on an intracranial artery. Intracranial aneurysms are often treated with a flow-diverting stent (FDS), a low-porosity metal mesh that covers the entrance to the cavity, to reduce blood flow into the aneurysm sac and exclude it from mechanical stresses imposed by the blood flow. Successful treatment is highly dependent on the degree of reduction of flow inside the cavity, and the resulting altered fluid mechanics inside the aneurysm following treatment. Using stereoscopic particle image velocimetry, we characterize the fluid mechanics in a canonical configuration representative of an intracranial aneurysm treated with a FDS: a spherical cavity on the side of a curved round pipe covered with a metal mesh formed by an actual medical FDS. This porous mesh coverage is the focus of Part 2 of the paper, characterizing the effects of parent vessel Re, De and pulsatility, Wo, on the fluid dynamics, compared with the canonical configuration with no impediments to flow into the cavity that is described in Part 1 (Chassagne et al., J. Fluid Mech., vol. 915, 2021, A123). Coverage with a FDS markedly reduces the flow Re in the aneurysmal cavity, creating a viscous-dominated flow environment despite the parent vessel Re > 100. Under steady flow conditions, the topology that forms inside the cavity is shown to be a function of the parent vessel De. At low values of De, flow enters the cavity at the leading edge and remains attached to the wall before exiting at the trailing edge, a novel behaviour that was not found under any conditions of the high-Re, unimpeded cavity flow described in Part 1. Under these conditions, flow in the cavity co-rotates with the direction of the free-stream flow, similar to Stokes flow in a cavity. As De increases, the flow along the leading edge begins to separate, and the recirculation zone grows with increasing De, until, above De ≈ 180, the flow inside the cavity is fully recirculating, counter-rotating with respect to the free-stream flow. Under pulsatile flow conditions, the vortex inside the cavity progresses through the same cycle - switching from attached and co-rotating with the free-stream flow at the beginning of the cycle (low velocity and positive acceleration) to separated and counter-rotating as De reaches a critical value. The location of separation within the harmonic cycle is shown to be a function of both De and Wo. The values of aneurysmal cavity Re based on both the average velocity and the circulation inside the cavity are shown to increase with increasing values of De, while Wo is shown to have little influence on the time-averaged metrics. As De increases, the strength of the secondary flow in the parent vessel grows, due to the inertial instability in the curved pipe, and the flow rate entering the cavity increases. Thus, the effectiveness of FDS treatment to exclude the aneurysmal cavity from the haemodynamic stresses is compromised for aneurysms located on high-curvature arteries, i.e. vessels with high De, and this can be a fluid mechanics criterion to guide treatment selection.

Fast virtual coiling algorithm for intracranial aneurysms using pre-shape path planning

To aid in predicting and improving treatment outcome of endovascular coiling of intracranial aneurysms, simulation of patient-specific coil deployment should be both accurate and fast. We developed a fast virtual coiling algorithm called Pre-shape Path Planning (P3). It captures the mechanical propensity of a released coil to restore its pre-shape for bending energy minimization, producing coils without unrealistic kinks and bends. A coil is discretized into finite-length segments and extruded from the delivery catheter segment-by-segment following a generic coil pre-shape. With the release of each segment, coil-wall and coil-coil collisions are detected and resolved. Modeling of each case took seconds to minutes. To test the algorithm, we evaluated its output against the literature, experiments, and patient angiograms. The periphery-to-core ratio of coils deployed by P3 decreased with increasing coil packing density, consistent with observations in the literature. Coils deployed by P3 compared well with in vitro experiments, free from unphysical kinks and loops that arose from previous virtual coiling algorithms. Simulations of coiling in four patient-specific aneurysms agreed well with the patient angiograms. To test the influence of coil pre-shape on P3, we performed hemodynamic simulations in aneurysms with coils deployed by P3 using the generic pre-shape, P3 using a coil-specific pre-shape, and full finite-element-method simulation. We found that the generic pre-shape was sufficient to produce results comparable to virtual coiling by finite element modeling. Based on these findings, P3 can rapidly simulate coiling in patient-specific aneurysms with good accuracy and is thus a potential candidate for clinical treatment planning.

Hemodynamic effects of intracranial aneurysms from stent-induced straightening of parent vessels by stent-assisted coiling embolization

BACKGROUND

Straightening of parent vessels happens for stent-assisted coiling embolization (SACE) treatment of intracranial aneurysms. This study aims to investigate aneurysmal hemodynamic modifications caused by stent-induced vessel straightening.

METHODS

Stent and coil deployments of a SACE-treated distal bifurcation aneurysm by finite element method were performed first with the preoperative (not straightened, NS) and postoperative (straightened, S) vessel models respectively. Computational fluid dynamics were then performed for eight models, including (I) NS only model, (II) NS+stent model, (III) NS+coils model, (IV) NS+stent+coils model, (V) S only model, (VI) S+stent model, (VII) S+coils model, and (VIII) S+stent+coils model. Finally, changes in aneurysmal flow velocity, isovelocity surface and wall shear stress (WSS) were analyzed qualitatively and quantitatively.

RESULTS

The flow was less in the S models than that in the corresponding NS models. Coils blocked most of the flow into the aneurysm sac in both NS models and S models and vessel straightening had more profound effect on the high aneurysmal flow volume reduction than coiling, while stenting generated adverse effect on flow reduction. Taking the NS only model as baseline (100%), the sac-averaged velocities of models II to VIII were 112%, 36%, 42%, 45%, 39%, 12%, 13%, and high flow volumes were 119%, 21%, 30%, 10%, 8%, 3%, 3%, while the sac-averaged WSSs were 106%, 37%, 44%, 41%, 35%, 17% and 24%, respectively.

CONCLUSIONS

Stent-induced vessel straightening combined coil embolization has the best performance in hemodynamic modifications and may reduce the recurrence rate, whereas stenting may generate adverse effect on hemodynamic alterations.

[Finite element simulation of stent implantation and its applications in the interventional planning for hemorrhagic cardio-cerebrovascular diseases]

Numerical simulation of stent deployment is very important to the surgical planning and risk assess of the interventional treatment for the cardio-cerebrovascular diseases. Our group developed a framework to deploy the braided stent and the stent graft virtually by finite element simulation. By using the framework, the whole process of the deployment of the flow diverter to treat a cerebral aneurysm was simulated, and the deformation of the parent artery and the distributions of the stress in the parent artery wall were investigated. The results provided some information to improve the intervention of cerebral aneurysm and optimize the design of the flow diverter. Furthermore, the whole process of the deployment of the stent graft to treat an aortic dissection was simulated, and the distributions of the stress in the aortic wall were investigated when the different oversize ratio of the stent graft was selected. The simulation results proved that the maximum stress located at the position where the bare metal ring touched the artery wall. The results also can be applied to improve the intervention of the aortic dissection and the design of the stent graft.

Implementation of computer simulation to assess flow diversion treatment outcomes: systematic review and meta-analysis

Introduction

Despite a decade of research into virtual stent deployment and the post-stenting aneurysmal hemodynamics, the hemodynamic factors which correlate with successful treatment remain inconclusive. We aimed to examine the differences in various post-treatment hemodynamic parameters between successfully and unsuccessfully treated cases, and to quantify the additional flow diversion achievable through stent compaction or insertion of a second stent.

Methods

A systematic review and meta-analysis were performed on eligible studies published from 2000 to 2019. We first classified cases according to treatment success (aneurysm occlusion) and then calculated the pooled standardized mean differences (SMD) of each available parameter to examine their association with clinical outcomes. Any additional flow diversion arising from the two common strategies for improving the stent wire density was quantified by pooling the results of such studies.

Results

We found that differences in the aneurysmal inflow rate (SMD −6.05, 95% CI −10.87 to −1.23, p=0.01) and energy loss (SMD −5.28, 95% CI −7.09 to −3.46, p<0.001) between the successfully and unsuccessfully treated groups were indicative of statistical significance, in contrast to wall shear stress (p=0.37), intra-aneurysmal average velocity (p=0.09), vortex core-line length (p=0.46), and shear rate (p=0.09). Compacting a single stent could achieve additional flow diversion comparable to that by dual-stent implantation.

Conclusions

Inflow rate and energy loss have shown promise as identifiers to discriminate between successful and unsuccessful treatment, pending future research into their diagnostic performance to establish optimal cut-off values.

Computational fluid dynamics as a risk assessment tool for aneurysm rupture

OBJECTIVE

The authors reviewed the clinical role of computational fluid dynamics (CFD) in assessing the risk of intracranial aneurysm rupture.

METHODS

A literature review was performed to identify reports on CFD assessment of aneurysms using PubMed. The usefulness of various hemodynamic parameters, such as wall shear stress (WSS) and the Oscillatory Shear Index (OSI), and their role in aneurysm rupture risk analysis, were analyzed.

RESULTS

The authors identified a total of 258 published articles evaluating rupture risk, growth, and endovascular device assessment. Of these 258 articles, 113 matching for CFD and hemodynamic parameters that contribute to the risk of rupture (such as WSS and OSI) were identified. However, due to a lack of standardized methodology, controversy remains on each parameter's role.

CONCLUSIONS

Although controversy continues to exist on which risk factors contribute to predict aneurysm rupture, CFD can provide additional parameters to assess this rupture risk. This technology can contribute to clinical decision-making or evaluation of efficacy for endovascular methods and devices.

Flow diversion of fusiform intracranial aneurysms

Fusiform aneurysms are less common than saccular aneurysms, but have higher associated mortality and rebleeding rates. Recently, flow diversion has emerged as a possible treatment option. The purpose of this study was to determine the safety and efficacy of the Pipeline Embolization Device (PED) for the treatment of ruptured and unruptured fusiform aneurysms. This was a retrospective analysis of patients with fusiform intracranial aneurysms treated with a PED at a quaternary care center between January 2012 and September 2019. Occlusion rates, neurologic morbidity/mortality, and other clinical variables were analyzed. Twenty-nine patients with 30 fusiform aneurysms were treated with a PED. Sixteen aneurysms (53%) were located in the anterior circulation and 14 aneurysms (47%) were in the posterior circulation. The mean maximal diameter of the aneurysms was 10.1 ± 5.6 mm (range 2.3-25 mm). Angiographic and clinical follow-up were available for 28 aneurysms (93%). The median follow-up was 17.4 months (IQR 4.8 to 28 months) and occlusion rates were graded according to the O'Kelly-Marotta (OKM) scale. Of patients with DSA follow-up, 15 aneurysms (60%) were completely occluded (OKM D) and 19 aneurysms (76%) had a favorable occlusion result (OKM C1-3 and D). The overall complication rate was 26.7% with a neurological morbidity rate of 6.7% and neurological mortality rate of 3.4%. Flow diversion can be an effective treatment for both ruptured and unruptured fusiform aneurysms. Nevertheless, complete occlusion rates are lower than for saccular aneurysms. Therefore, flow diversion should be considered only if other more direct treatment options, such as clipping or stent/coiling are not applicable. Flow diversion should be used cautiously in patients presenting with rupture.

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.

Improving accuracy for finite element modeling of endovascular coiling of intracranial aneurysm

BACKGROUND

Computer modeling of endovascular coiling intervention for intracranial aneurysm could enable a priori patient-specific treatment evaluation. To that end, we previously developed a finite element method (FEM) coiling technique, which incorporated simplified assumptions. To improve accuracy in capturing real-life coiling, we aimed to enhance the modeling strategies and experimentally test whether improvements lead to more accurate coiling simulations.

METHODS

We previously modeled coils using a pre-shape based on mathematical curves and mechanical properties based on those of platinum wires. In the improved version, to better represent the physical properties of coils, we model coil pre-shapes based on how they are manufactured, and their mechanical properties based on their spring-like geometric structures. To enhance the deployment mechanics, we include coil advancement to the aneurysm in FEM simulations. To test if these new strategies produce more accurate coil deployments, we fabricated silicone phantoms of 2 patient-specific aneurysms in duplicate, deployed coils in each, and quantified coil distributions from intra-aneurysmal cross-sections using coil density (CD) and lacunarity (L). These deployments were simulated 9 times each using the original and improved techniques, and CD and L were calculated for cross-sections matching those in the experiments. To compare the 2 simulation techniques, Euclidean distances (dMin, dMax, and dAvg) between experimental and simulation points in standardized CD-L space were evaluated. Univariate tests were performed to determine if these distances were significantly different between the 2 simulations.

RESULTS

Coil deployments using the improved technique agreed better with experiments than the original technique. All dMin, dMax, and dAvg values were smaller for the improved technique, and the average values across all simulations for the improved technique were significantly smaller than those from the original technique (dMin: p = 0.014, dMax: p = 0.013, dAvg: p = 0.045).

CONCLUSION

Incorporating coil-specific physical properties and mechanics improves accuracy of FEM simulations of endovascular intracranial aneurysm coiling.

Aneurysm characteristics, coil packing, and post-coiling hemodynamics affect long-term treatment outcome

BACKGROUND

Recurrence of intracranial aneurysms after endovascular coiling is a serious clinical concern.

OBJECTIVE

We hypothesized that recurrence is associated with aneurysm morphology and flow, as well as the coil intervention and the induced flow modifications.

METHODS

We collected 52 primary-coiling aneurysm cases that were either occluded (n=34) or recurrent (n=18) at >1 year follow-up. We created aneurysm models from pre-coiling digital subtraction angiographic images, calculated aneurysm morphology, simulated pre-coiling hemodynamics, modeled coil deployment, and obtained post-coiling hemodynamics for each case. We performed univariable analysis on 26 morphologic, treatment-specific, and hemodynamic parameters to distinguish between recurrent and occluded groups, and multivariable analysis to identify independently significant parameters associated with recurrence. Univariable analysis was also performed on ruptured and unruptured aneurysm subcohorts separately to investigate if they shared specific significant parameters.

RESULTS

Recurrence was associated with pre-coiling aneurysm morphologic and flow parameters including larger size (maximum dimension and volume), larger neck (diameter, area, and neck-to-parent-artery ratio), and higher flow momentum and kinetic energy. Recurrence was also associated with lower coil packing (packing density and uncoiled volume), higher post-treatment flow (velocity, momentum, and kinetic energy), lower post-treatment washout time, and higher post-treatment impingement force at the neck. Multivariable analysis identified two aneurysmal characteristics (neck diameter and pre-coiling flow kinetic energy), one coil packing parameter (uncoiled volume), and one post-treatment hemodynamic parameter (flow momentum) that were independently associated with recurrence. In ruptured aneurysms, recurrence was associated with larger neck (diameter and area), whereas in unruptured aneurysms, recurrence was associated with larger size (maximum dimension and volume). In both subcohorts, recurrence was associated with higher post-coiling flow momentum and kinetic energy.

CONCLUSION

Recurrence at >1 year after coil treatment is associated with intrinsic aneurysm characteristics, coiling itself, and flow changes induced by coiling. Larger aneurysm size and neck, less coil packing, and higher intra-aneurysmal flow before and after coiling predict recurrence.

Variation of Mass Effect After Using a Flow Diverter With Adjunctive Coil Embolization for Symptomatic Unruptured Large and Giant Intracranial Aneurysms

Background: Mass effect associated with large or giant aneurysms is an intractable problem for traditional endovascular treatments. Preventing recurrence of aneurysms requires dense coiling, which may aggravate the mass effect. However, the flow diverter (FD) is a new device that avoids the need for dense coiling. This study was performed to investigate whether use of FDs with adjunctive coil embolization can relieve the aneurysmal mass effect and to explore the factors that affect the variation of compressional symptoms.

Methods: We retrospectively evaluated patients with compressional symptoms caused by unruptured aneurysms who underwent endovascular treatment with an FD with adjunctive coil embolization at our center from January 2015 to December 2017. Imaging follow-up included digital subtraction angiography (DSA) ranging from 11 to 14 months and magnetic resonance imaging (MRI) ranging from 24 to 30 months; the former was used to evaluate the intracavitary volume, and the latter was used to measure the variation of the mass effect. Follow-up physical examinations were performed to observe variations of symptoms.

Results: In total, 22 patients with 22 aneurysms were treated by an FD combined with coil embolization. All 22 patients underwent the last clinical follow-up. Regarding compressional symptoms, 12 (54.54%) patients showed improvement, 6 (27.27%) were fully recovered, and 6 (27.27%) showed improvement but with incomplete cranial palsy. However, five (22.72%) patients showed no change, four (18.18%) showed worsening symptoms compared with their preoperative state, and one (4.55%) died of delayed rupture. Seventeen of the 22 patients underwent MRI. Of these 17 patients, the aneurysm shrank in 13 (76.47%) and no significant change occurred in 4 (23.53%). In the multivariate analysis, a short duration from symptom occurrence to treatment (p = 0.03) and younger patient age (p = 0.038) were statistically significant factors benefiting symptom improvement, and shrinkage of the aneurysm was associated with favorable clinical outcomes (p = 0.006).

Conclusions: Use of the FD with adjunctive loose coil embolization might help to alleviate the compressional symptoms caused by intracranial aneurysms. Shrinkage of the aneurysm, a short duration of symptoms, and younger patient age might contribute to favorable outcomes of mass effect-related symptoms.

Numerical studies of hemodynamic alterations in pre- and post-stenting cerebral aneurysms using a multiscale modeling

The aim of this work was to use a multiscale modeling to study the influence of stent deployment, with generic stents, on flow distributions within the vascular network and the hemodynamic alterations within the cerebral aneurysms pre- and post-stenting. To achieve this goal, two image-based anatomical cerebral aneurysm models were reconstructed along with the respective aneurysms post-stenting models after deploying a 16- or 24-wire stent. The investigation results revealed that the stent may increase the local pressure resistance resulting in flow alterations. The hemodynamic parameters demonstrated stent placement can reduce the intra-aneurysmal pressure, decrease wall shear stress (WSS) at the neck region, and increase blood turnover time for aneurysm case I (sidewall aneurysm). These findings are consistent with the trends of hemodynamic changes reported previously. However, aneurysm case II (bifurcation aneurysm) showed gradually increased intra-aneurysmal pressure and the pressure at the neck region, decreased WSS over the sac surface, and enhanced flow vortices within the aneurysm. When simulating the hemodynamics of pre- and post-stenting aneurysms for a patient using measured flow waveforms, the flow alteration induced by the stent deployment may affect the hemodynamic predictions for the post-stenting aneurysm. Thus, the remeasurement of boundary conditions once the morphology of the aneurysm is deformed is needed in follow-up studies with a focus on aneurysm growth and stent deployment.

Numerical simulation of stent deployment within patient-specific artery and its validation against clinical data

BACKGROUND AND OBJECTIVE

One of the most widely adopted endovascular treatment procedures is the stent implantation. The effectiveness of the treatment depends on the appropriate stent expansion. However, it is difficult to accurately predict the outcome of such an endovascular intervention. Numerical simulations represent a useful tool to study the complex behavior of the stent during deployment. This study presents a numerical model capable of simulating this process.

METHODS

The numerical model consists of three parts: modeling of stent expansion, modeling the interaction of the stent with the arterial wall and the deformation of the arterial wall. The model is able to predict the shapes of both stent and arterial wall during the entire deployment process. Simulations are performed using patient-specific clinical data that ensures more realistic results.

RESULTS

The numerical simulations of stent deployment are performed using the extracted geometry of the coronary arteries of two patients. The obtained results are validated against clinical data from the follow up examination and both quantitative and qualitative analysis of the results is presented. The areas of several slices of the arterial wall are calculated for all the three states (before, after and follow up) and the standard error of the area when comparing simulation and follow up examination is 5.27% for patient #1 and 4.5% for patient #2.

CONCLUSIONS

The final goal of numerical simulations in stent deployment should be to provide a clinical tool that is capable of reliably predicting the treatment outcome. This study showed through the good agreement of results of the numerical simulations and clinical data that the presented numerical model represents a step towards this final goal. These simulations can also provide valuable information about distribution of forces and stress in the arterial wall that can improve pre-operative planning and treatment optimization.

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.

Differences in hemodynamic characteristics under high packing density between the porous media model and finite element analysis in computational fluid dynamics of intracranial aneurysm virtual treatment

Objective

This study aimed to compare the hemodynamic differences among no sac (NOS), porous media (POM) and finite element analysis (FEA) models to investigate the recurrence-related risks for coiled intracranial aneurysms (IAs).

Methods

The study enrolled 10 patients with 11 IAs who received simple coiling treatment and hemodynamic simulations were performed for all IAs using the above three models. Velocity, wall shear stress (WSS) and residual flow volume (RFV) were calculated and compared in order to assess the model differences for both aneurysm sac and parent vessel regions.

Results

For parent artery regions, all three models produced similar flow patterns and quantification analysis did not indicate differences in velocity and WSS (p>0.05). For aneurysm sac regions, the FEA model resulted in higher sac-maximized (0.18 m/s vs 0.06 m/s) and sac-averaged velocity (0.013 m/s vs 0.007 m/s), and higher sac-averaged (0.55 Pa vs 0.36 Pa, p=0.006) and sac-maximized WSS (12.1 Pa vs 6.6 Pa) than the POM model. The differences in RFV between the POM and FEA models under 11 different isovelocity thresholds (0.0001 m/s, 0.001 m/s, 0.002 m/s, 0.005 m/s, 0.01 m/s, 0.02 m/s, 0.05 m/s, 0.1 m/s, 0.2 m/s, 0.5 m/s, and 1 m/s) showed that the POM RFV was generally larger than those of the FEA model.

Conclusions

Compared with the FEA model, the POM model provides a lower velocity and WSS and higher RFV for the aneurysm sac, which could lead to incorrect estimates of the recurrent risk of coiled IAs under high packing density.

Ostium Ratio and Neck Ratio Could Predict the Outcome of Sidewall Intracranial Aneurysms Treated with Flow Diverters

BACKGROUND AND PURPOSE

Incompletely occluded flow diverter treated aneurysms remain at risk of rupture and thromboembolic complications. Our aim was to identify the potential for incomplete occlusion of intracranial aneurysms treated by flow diverters. We investigated whether aneurysm ostium size in relation to parent artery size affects angiographic outcomes of flow diverter-treated sidewall aneurysms.

MATERIALS AND METHODS

Flow diverter-treated sidewall aneurysms were divided into "occluded" and "residual" (incomplete occlusion) groups based on 6-month angiographic follow-up. We calculated the ostium ratio, a new parameter defined as the aneurysm ostium surface area versus the circumferential surface area of the parent artery. We also calculated the neck ratio, defined as clinical aneurysm neck diameter versus parent artery diameter from pretreatment 2D DSA, as a 2D surrogate. We compared the performance of these ratios with existing aneurysm morphometrics (size, neck diameter, volume, aspect ratio, size ratio, undulation index, nonsphericity index, ellipticity index, bottleneck factor, aneurysm angle, and parent vessel angle) and flow diverter-related parameters (metal coverage rate and pore density). Statistical tests and receiver operating characteristic analyses were performed to identify significantly different parameters between the 2 groups and test their predictive performances.

RESULTS

We included 63 flow diverter-treated aneurysms, 46 occluded and 17 residual. The ostium ratio and neck ratio were significantly higher in the residual group than in the occluded group (P < .001 and P = .02, respectively), whereas all other parameters showed no statistical difference. As discriminating parameters for occlusion, ostium ratio and neck ratio achieved areas under the curve of 0.912 (95% CI, 0.838-0.985) and 0.707 (95% CI, 0.558-0.856), respectively.

CONCLUSIONS

High ostium ratios and neck ratios could predict incomplete occlusion of flow diverter-treated sidewall aneurysms. Neck ratio can be easily calculated by interventionists to predict flow-diverter treatment outcomes.

Numerical simulation of patient-specific endovascular stenting and coiling for intracranial aneurysm surgical planning

Background

In this study, we develop reliable and practical virtual coiling and stenting methods for intracranial aneurysm surgical planning. Since the purpose of deploying coils and stents is to provide device geometries for subsequent accurate post-treatment computational fluid dynamics analysis, we do not need to accurately capture all the details such as the stress and force distribution for the devices and vessel walls. Our philosophy for developing these methods is to balance accuracy and practicality.

Methods

We consider the mechanical properties of the devices and recapitulate the clinical practice using a finite element method (FEM) approach. At the same time, we apply some simplifications for FEM modeling to make our methods efficient. For the virtual coiling, the coils are modeled as 3D Euler–Bernoulli beam elements, which is computationally efficient and provides good geometry representation. During the stent deployment process, the stent–catheter system is transformed according to the centerline of the parent vessel since the final configuration of the stent is not dependent of the deployment history. The aneurysm and vessel walls are assumed to be rigid and are fully constrained during the simulation. All stent–catheter system and coil–catheter system are prepared and packaged as a library which contains all types of stents, coils and catheters, which improves the efficiency of surgical planning process.

Results

The stent was delivered to the suitable position during the clinical treatment, achieving good expansion and apposition of the stent to the arterial wall. The coil was deployed into the aneurysm sac and deformed to different shapes because of the stored strain energy during coil package process and the direction of the microcatheter.

Conclusions

The method which we develop here could become surgical planning for intracranial aneurysm treatment in the clinical workflow.

The Combined Use of Intraluminal and Intrasaccular Flow Diversion for the Treatment of Intracranial Aneurysms: Report of 25 Cases

PURPOSE

The Medina Embolic Device (MED) is a new intrasaccular device with promising early results. Previously we documented our initial experience of this device both alone and in combination with other devices including flow diverter stents (FDS). We sought to determine the effect of the MED + FDS strategy for the treatment of selected aneurysms.

MATERIALS AND METHODS

We performed a retrospective analysis of prospectively collected data to identify all patients with aneurysms treated using both the MED and intraluminal FDS. We present our technical success rate, early and mid-term angiographic follow-up, and clinical outcome data.

RESULTS

We identified 25 non-consecutive patients. The treatment was staged in 9 patients and in a single session 16 patients. The average age was 61±12.8 years (range 40-82). The average fundus height was 11±3.6 mm and average fundus width was 10.1±3.4 mm. In the staged cohort (n=9) at delayed angiography (mean 10 mths) 8 aneurysms (89%) showed complete exclusion (mRRC 1) and in one patient there was a parent vessel occlusion. In the simultaneous cohort delayed angiography (n=10, mean 8.1 months) demonstrated complete occlusion (mRRC 1) in 6 aneurysms (60%), 3 neck remnants (mRRC 2) (30%) and 1 patient (10%) showed persistent aneurysmal filling (mRRC 3a). There were 5 complications with permanent morbidity (mRS >2) in two patients. There were no mortalities.

CONCLUSION

The MED can be successfully used in combination with intraluminal FDS and in selected aneurysms this may represent an alternative to FDS and adjunctive coiling.

Effect of Local Coil Density on Blood Flow Stagnation in Densely Coiled Cerebral Aneurysms: A Computational Study Using a Cartesian Grid Method

Aneurysm recurrence is the most critical concern following coil embolization of a cerebral aneurysm. Adequate packing density (PD) and coil uniformity are believed necessary to achieve sufficient flow stagnation, which decreases the risk of aneurysm recurrence. The effect of coil distribution on the extent of flow stagnation, however, especially in cases of dense packing (high PD), has received less attention. Thus, the cause of aneurysm recurrence despite dense packing is still an open question. The primary aim of this study is to evaluate the effect of local coil density on the extent of blood flow stagnation in densely coiled aneurysms. For this purpose, we developed a robust computational framework to determine blood flow using a Cartesian grid method, by which the complex fluid pathways in coiled aneurysms could be flexibly treated using an implicit function. This tool allowed us to conduct blood flow analyses in two patient-specific geometries with 50 coil distribution patterns in each aneurysm at clinically adequate PD. The results demonstrated that dense packing in the aneurysm may not necessarily block completely the inflow into the aneurysm and local flow that formed in the neck region, whose strength was inversely related to this local PD. This finding suggests that local coil density in the neck region still plays an important role in disturbing the remaining local flow, which possibly prevents thrombus formation in a whole aneurysm sac, increasing the risk of aneurysm regrowth and subsequent recurrence.

Hemodynamics and coil distribution with changing coil stiffness and length in intracranial aneurysms

Purpose

The purpose of this study was to investigate hemodynamics and coil distribution with changing coil stiffness and length using the finite element method (FEM) and computational fluid dynamics (CFD) analysis.

Methods

Basic side-wall and bifurcation type aneurysm models were used. Six types of coil models were generated by changing the coil stiffness and length, based on commercially available embolic coils. Coil embolization was simulated using FEM. CFD was performed to characterize the hemodynamics in the aneurysms after embolization. Coil distribution and velocity reduction in the aneurysms were evaluated.

Results

The median value of radial coil distribution was shifted from the center to the outer side of the aneurysmal dome by changing coil stiffness: harder coils entered the outer side of the aneurysmal dome more easily. Short coils were more distributed at the neck region, since their small size made it easy for them to enter the tighter area. CFD results also indicated that velocity in the aneurysm was effectively reduced when the coils were more distributed at the neck region and the outer side of the aneurysmal dome because of the disturbance in blood inflow.

Conclusions

It is easier for coils to enter the outer side of the aneurysmal sphere when they are harder. If coils are short, they can enter tighter areas more easily. In addition, high coil density at the outer side of the aneurysmal dome and at the neck region is important to achieve effective velocity reduction.

A new combined parameter predicts re-treatment for coil-embolized aneurysms: a computational fluid dynamics multivariable analysis study

Purpose

Coil embolization is a minimally invasive method used to treat cerebral aneurysms. Although this endovascular treatment has a high success rate, aneurysmal re-treatment due to recanalization remains a major problem of this method. The purpose of this study was to determine a combined parameter that can be useful for predicting aneurysmal re-treatment due to recanalization.

Methods

Patient-specific geometries were used to retrospectively analyze the blood flow for 26 re-treated and 74 non-retreated aneurysms. Post-operatively aneurysms were evaluated at 12-month follow-up. The hemodynamic differences between the re-treatment and non-retreatment aneurysms were analyzed before and after coil embolization using computation fluid dynamics. Basic fluid characteristics, rates of change, morphological factors of aneurysms and patient-specific clinical information were examined. Multivariable analysis and logistic regression analysis were performed to determine a combined parameter—re-treatment predictor (RP).

Results

Among examined hemodynamic, morphological, and clinical parameters, slight reduction of blood flow velocity rate in the aneurysm, slight increase of pressure rate at the aneurysmal neck and neck area, and hypertension were the main factors contributing to re-treatment. Notably, hemodynamic parameters between re-treatment and non-retreatment groups before embolization were similar: however, we observed significant differences between the groups in the post-embolization average velocity and the rate of reduction in this velocity in the aneurysmal dome.

Conclusions

The combined parameter, RP, which takes into consideration hemodynamic, morphological, and clinical parameters, accurately predicts aneurysm re-treatment. Calculation of RP before embolization may be able to predict the aneurysms that will require re-treatment.

Compacting a Single Flow Diverter versus Overlapping Flow Diverters for Intracranial Aneurysms: A Computational Study

BACKGROUND AND PURPOSE

Locally compacting the mesh of a flow diverter by a dynamic push-pull technique can accelerate intracranial aneurysm healing. We asked how this deployment strategy compares with overlapping 2 flow diverters for aneurysmal flow reduction.

MATERIALS AND METHODS

Using a high-fidelity virtual stent placement method, we simulated 3 flow-diverter strategies (single noncompacted, 2 overlapped, and single compacted) in 3 aneurysms (fusiform, large saccular, and medium saccular). Computational fluid dynamics analysis provided posttreatment hemodynamic parameters, including time-averaged inflow rate, aneurysm-averaged velocity, wall shear stress, total absolute circulation, and turnover time. We examined the relationship between the achieved degree of compaction and aneurysm orifice area.

RESULTS

Flow-diverter compaction resulted in a compaction coverage of 57%, 47%, and 22% over the orifice of the fusiform, large, and medium saccular aneurysm, respectively. Compaction coverage increased linearly with orifice area. In the fusiform aneurysm, the single compacted flow diverter accomplished more aneurysmal flow reduction than the other 2 strategies, as indicated by all 5 hemodynamic parameters. In the 2 saccular aneurysms, the overlapped flow diverters achieved the most flow reduction, followed by the single compacted and the noncompacted flow diverter.

CONCLUSIONS

Compacting a single flow diverter can outperform overlapping 2 flow diverters in aneurysmal flow reduction, provided that the compaction produces a mesh denser than 2 overlapped flow diverters and this denser mesh covers a sufficient portion of the aneurysm orifice area, for which we suggest a minimum of 50%. This strategy is most effective for aneurysms with large orifices, especially fusiform aneurysms.

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

BACKGROUND

Endovascular intervention using a stent is a mainstream treatment for cerebral aneurysms. To assess the effect of intervention strategies on aneurysm hemodynamics, we have developed a fast virtual stenting (FVS) technique to simulate stent deployment in patient-specific aneurysms. However, quantitative validation of the FVS against experimental data has not been fully addressed. In this study, we performed in vitro analysis of a patient-specific model to illustrate the realism and usability of this novel FVS technique.

METHODS

We selected a patient-specific aneurysm and reproduced it in a manufactured realistic aneurismal phantom. Three numerical simulation models of the aneurysm with an Enterprise stent were constructed. Three models were constructed to obtain the stented aneurysms: a physical phantom scanned by micro-CT, fast virtual stenting technique and finite element method. The flow in the three models was simulated using a computational fluid dynamics software package, and the hemodynamics parameters for the three models were calculated and analyzed.

RESULTS

The computational hemodynamics in the patient-specific aneurysm of the three models resembled the very well. A qualitative comparison revealed high similarity in the wall shear stress, streamline, and velocity plane among the three different methods. Quantitative comparisons revealed that the difference ratios of the hemodynamic parameters were less than 10%, with the difference ratios for area average of wall shear stress in the aneurysm being very low.

CONCLUSIONS

In conclusion, the results of the computational hemodynamics indicate that FVS is suitable for evaluation of the hemodynamic factors that affect treatment outcomes.

Safety and efficacy of flow re-direction endoluminal device (FRED) in the treatment of cerebral aneurysms: a single center experience

BACKGROUND

Experience with the endovascular treatment of cerebral aneurysms by the Flow Re-Direction Endoluminal Device (FRED) is still limited. The aim of this study is to discuss the results and complications of this new flow diverter device (FDD).

METHODS

Between November 2013 and April 2015, 20 patients (15 female and five male) harboring 24 cerebral aneurysms were treated with FRED FDD in a single center.

RESULTS

Complete occlusion was obtained in 20/24 aneurysms (83 %) and partial occlusion in four (17 %). Intraprocedural technical complication occurred in one case (4 %) and post-procedural complications in three (12 %). None reported neurological deficits (mRS = 0). All FRED were patent at follow-up. No early or delayed aneurysm rupture, no subarachnoid (SAH) or intraparenchymal hemorrhage (IPH) no ischemic complications and no deaths occurred.

CONCLUSIONS

Endovascular treatment with FRED FDD is a safe treatment for unruptured cerebral aneurysms, resulting in a high rate of occlusion. The FRED is substantially equivalent to the other known FDDs, which show similar functions and technical profiles.

Computational fluid dynamics of cerebral aneurysm coiling using high-resolution and high-energy synchrotron X-ray microtomography: comparison with the homogeneous porous medium approach

BACKGROUND

Computational modeling of intracranial aneurysms provides insights into the influence of hemodynamics on aneurysm growth, rupture, and treatment outcome. Standard modeling of coiled aneurysms simplifies the complex geometry of the coil mass into a homogeneous porous medium that fills the aneurysmal sac. We compare hemodynamics of coiled aneurysms modeled from high-resolution imaging with those from the same aneurysms modeled following the standard technique, in an effort to characterize sources of error from the simplified model.

MATERIALS

Physical models of two unruptured aneurysms were created using three-dimensional printing. The models were treated with coil embolization using the same coils as those used in actual patient treatment and then scanned by synchrotron X-ray microtomography to obtain high-resolution imaging of the coil mass. Computational modeling of each aneurysm was performed using patient-specific boundary conditions. The coils were modeled using the simplified porous medium or by incorporating the X-ray imaged coil surface, and the differences in hemodynamic variables were assessed.

RESULTS

X-ray microtomographic imaging of coils and incorporation into computational models were successful for both aneurysms. Porous medium calculations of coiled aneurysm hemodynamics overestimated intra-aneurysmal flow, underestimated oscillatory shear index and viscous dissipation, and over- or underpredicted wall shear stress (WSS) and WSS gradient compared with X-ray-based coiled computational fluid dynamics models.

CONCLUSIONS

Computational modeling of coiled intracranial aneurysms using the porous medium approach may inaccurately estimate key hemodynamic variables compared with models incorporating high-resolution synchrotron X-ray microtomographic imaging of complex aneurysm coil geometry.