Abstract
Endovascular treatment has become the standard therapy for cerebral aneurysms, while the effective treatment for middle cerebral artery (MCA) bifurcation aneurysms remains a challenge. Current flow-diverting techniques with endovascular coils cover the aneurysm orifice as well as adjacent vessel branches, which may lead to branch occlusion. Novel endovascular flow disruptors, such as the Contour device (Cerus Endovascular), are of great potential to eliminate the risk of branch occlusion. However, there is a lack of valid comparison between novel flow disruptors and conventional (intraluminal) flow-diverters. In this study, two in silico MCA bifurcation aneurysm models were treated by specific Contour devices and flow-diverters using fast-deployment algorithms. Computational fluid dynamic simulations were used to examine the performance and efficiency of deployed devices. Hemodynamic parameters, including aneurysm inflow and wall shear stress, were compared among each Contour device, conventional flow-diverter, and untreated condition. Our results show that the placement of devices can effectively reduce the risk of aneurysm rupture, while the deployment of a Contour device causes more flow reduction than using flow-diverters (e.g. Silk Vista Baby). Besides, the Contour device presents the flow diversion capability of targeting the aneurysm neck without occluding the daughter vessel. In summary, the in silico aneurysm models presented in this study can serve as a powerful pre-planning tool for testing new treatment techniques, optimising device deployment, and predicting the performance in patient-specific aneurysm cases. Contour device is proved to be an effective treatment of MCA bifurcation aneurysms with less daughter vessel occlusion.
Similar content being viewed by others
1 Introduction
The effective treatment of middle cerebral artery (MCA) bifurcation aneurysms remains challenging, as several lesions present in such sites (Alfano et al. 2013; Zhang et al. 2018). Current treatment options focus either on aneurysm coil retention supported by a stent-like device positioned in the parent vascular lumen or on conventional (intraluminal) flow-diverters to disrupt blood flow within the aneurysm sac (Zhang et al. 2017; Mokin et al. 2020). One limitation of using conventional flow-diverters for such bifurcation aneurysms is that one daughter needs to be jailed, with evidence of causing branch occlusion and mandating antiplatelet medication (Raychev et al. 2016; Narata et al. 2018). Significant issues regarding jailed vessel occlusion in clinical cases have been reported in the literature. A study by (Saatchi et al. 2012) covered 46 ‘uncoilable’ aneurysms that originated at vessel bifurcations, with successful treatment by the Pipeline™ Embolisation Device (PED) reported in 41 (89.1%) cases. Among the remaining 5 cases (10.9% of the total), 1 patient died due to complete occlusion of the jailed daughter vessel. Similar results were reported by (Saleme et al. 2014). A retrospective analysis of 37 cases of bifurcation aneurysms treated with flow-diverters. At 6-month follow-up, complete aneurysm occlusion was observed in 97.3% of cases, while the new permanent neurological deficit was observed in 9.4% of patients. During follow-up, jailed daughter vessel occlusion was found in 32.4% of aneurysms and 40.5% of cases showed significant vessel stenosis. However, the effects of vessel occlusion or narrowing were only symptomatic in 13.5% of all aneurysms.
As an alternative, intravascular devices create a barrier at the neck inside the aneurysm, which potentially eliminates the risk of branch occlusion and the need for long-term anticoagulation therapy (Pierot et al. 2012; Liebig et al. 2015; Arthur et al. 2019; Youssef et al. 2021). Endovascular flow disruptors include the Woven-EndoBridge (WEB, MicroVention) (Goyal et al. 2020), the LUNA-AED (Medtronic) (Piotin et al. 2018), and, most recently reported, the Contour (Cerus Endovascular) (Akhunbay-Fudge et al. 2020; Bhogal et al. 2021). Both the WEB and LUNA-AED devices feature a spherical braid that fills the aneurysm dome, while the Contour is built on a planar braid (Bhogal et al. 2019), entirely positioned at the neck (see Fig. 1).
The Contour Neurovascular System™ (Cerus Endovascular Inc., Fremont, CA) is designed to target the neck of aneurysms by providing flow disruption within the aneurysm sac combined flow diversion at the neck of the aneurysm (Thormann et al. 2022). This Contour is a circular, dual-layered structure of 2 × 72 nitinol wires (Fig. 1a) with 1 platinum marker (Fig. 1b). It is radiopaque and retrievable until electrolytically detached. On deployment, the device adapts to the lower half of the aneurysm, covering the neck (Fig. 1c). The Contour is deployed through a 0.027 microcatheter (″). Sizing is performed relative to the aneurysm neck width and diameter, disregarding height, with 5-, 7-, 9-, and 11-mm-diameter devices available (Liebig et al. 2022). Moreover, the device is designed to support neo-intimal growth across the mesh at the neck after thrombosis of the aneurysm.
Recent data from the CERUS study (the first multicentre clinical trial evaluating safety and effectiveness outcomes with the Contour neurovascular system) have indicated promising results in the treatment of unruptured intracranial bifurcation aneurysms (Lawson et al. 2017; Biondi et al. 2023). To further evaluate the performance of Contour and its differences compared to conventional flow-diverters when treating bifurcation aneurysms, patient-specific computational fluid dynamics (CFD) models can be utilised to evaluate the effects of novel Contour device.
In this work, the Contour devices were deployed in two MCA bifurcation aneurysm in silico models and then compared with conventional endovascular flow-diverter devices deployed in the same model. The structure of the Contour device was acquired via micro-computed tomography (CT) scanning. This study evaluates two important hemodynamic parameters in the context of aneurysms, such as wall shear stress (WSS) and aneurysm inflow rate (Q).
2 Methods
2.1 In silico aneurysm models and clinical approach
Two in silico MCA bifurcation aneurysms models were manually segmented in OsiriX (OsiriX v.4.1.1, Freeware) from CT angiography imaging data. Both pre-intervention models were further reconstructed and trimmed (Fig. 2a). The structure of the Contour device was obtained from a Micro-CT scan (Fig. 2b). A conventional flow-diverter device was reconstructed in a range of sizes as per the methodology previously detailed in our previous study (Lyu et al. 2021). In the clinical approach, the selection of the sizing process of the Contour device is based on analysing aneurysm neck width (dn) and maximum aneurysm width (da).
As can be seen in Fig. 3, the aneurysm neck and aneurysm widths are 8.58 mm and 9.37 mm, respectively, in Case I; therefore, a CNS014-15 (diameter: 14 mm) Contour device and Silk Vista Baby (SVB) (2.25 × 15 mm) were chosen for the deployment. Both sizes are in line with manufacturer instructions for use. As for Case II, the aneurysm neck and aneurysm width are 3.11 mm and 5.56 mm, respectively; therefore, a CNS21007-15 (diameter: 7 mm) Contour device and 3 × 15 mm SVB were chosen for the deployment inside the sac of the aneurysm.
2.2 Virtual deployment
To virtually model the deployment of both the Contour and SVB devices, four variants of the devices were created (Fig. 4). Stents I and II mimic the SVB devices (2.25 × 15 mm) and (3 × 15 mm) deployed as per manufacturer instructions with a free expansion of up to 2.25 mm and 3 mm in diameter. Contours I and II mimic the Contour device deployed in the neck of the aneurysm, and the size of the Contour device is based on the relative size of the aneurysm and specifications of the Contour device (Table 1).
A fast virtual deployment algorithm was used for device placement. Deployment algorithms were reported in our previous study (Peach et al. 2014a; Spranger and Ventikos 2014; Lyu et al. 2021). For SVB, a representative device wireframe is compressed to minimise the diameters and then aligned with the vessel centreline to simulate the device sheathing process. The device unsheathing process is accomplished by the SVB releasing progressively along its length. Within the limit of the vessel wall, the device extends to a stress-free shape. As for the Contour device, the wireframe is obtained from a Micro-CT scan and then compressed to reduce the diameter and align with a deployment path into an aneurysm to mimic the sheathing process. The device unsheathing is again achieved by relaxing the device, which then expands to its original shape inside the aneurysm. By relaxing the device, the unsheathing process is achieved by the device expanding to its original shape inside the aneurysm. Device sizing, porosity, and pore density are shown in Table 2 for the SVB and Contour device as per manufacturers’ guidance (Eker et al. 2022).
2.3 CFD simulation
Blood was modelled as an incompressible fluid, and the 3D Navier–Stokes equations were solved in CFD-ACE + (ESI Group, Paris, France) using the finite volume method with a central differencing scheme for spatial interpolations. We have used a second-order geometry-weighted interpolation for the convection terms of the momentum conservation equation, and we used pure 2nd order central finite differences both for the diffusion terms and for all relevant derivatives in the pressure correction equation. The SIMPLE Consistent (SIMPLEC) pressure correction method (Van Doormaal and Raithby 1984; Lonsdale 1993) and an algebraic multigrid method for convergence acceleration were used (Ni and Abdou 2007). Although previous studies (Valencia and Solis 2006; Liu et al. 2021) have confirmed that the non-Newtonian effect of blood in the brain is small in the cerebral circulation, blood is modelled as a Newtonian fluid with a density of 1000 kg/m3 and a dynamic viscosity of 0.004 Pa·s as many well-known CFD studies of cerebrovascular blood flow studies (Xu et al. 2018; Hossain et al. 2021; Boite et al. 2023). Rigid wall assumption was applied to the arterial walls (Xu et al. 2018; Hossain et al. 2021; Boite et al. 2023), with the effect of such an assumption on flow patterns having been shown to be small for the cerebral vasculature. Both the device struts and vessel walls were applied as no-slip boundary conditions. Convergence criteria of absolute or relative residual reduction to 1 × 10−8 and 1 × 10−5 were employed for all variables. The solution was considered to be converged if the change of residual is smaller than 1 × 10−5 (relative residual reduction) or the residual is below 1 × 10−8 (absolute residual). The convergence is typically reached in fewer than 100 iterations.
A radially symmetric inlet velocity boundary condition was applied to each geometry scaling the corresponding velocity to a mean internal carotid artery (ICA) flow rate of 230 ml/min (Valencia et al. 2006). A fixed pressure outlet boundary condition was applied to all geometry outlets; more complex outflow conditions incorporating windkessel models were considered but rejected given very little variation in daughter vessel flowrates under a constant pressure condition when compared to physiological values (to within 5% of mean flow rates reported in the literature). One of the objectives of this study is to find a balance between time and accuracy, therefore steady-state simulations were run on 32 × 2.80 GHz CPUs 128 GB RAM for 1–1.5 h depending on mesh size. Our previous study (Peach et al. 2017, 2019) reported that around 75 h is required for a transient simulation.
2.4 Mesh Independence
The middle cerebral artery (MCA) bifurcation aneurysm models, both without and with deployed devices, were meshed in CFD-VisCART (ESI Group, Paris, France). The meshing of each geometry is completed with a Projected Single Domain con-conforming mesh, an Omnitree Cartesian tree type, and three near-wall Cartesian layers to give a smooth and well-resolved boundary definition (Peach et al. 2016). The detail of mesh around the Contour device and the fluid domain of Case II is shown in Fig. 5.
As reported previously by our group (Peach et al. 2016), simulations run for mesh independence achieved with a typical minimum cell density finer than 4000 elements/mm3 for Cases I and II; such mesh densities meet the requirements discussed in the literature for smooth velocity and WSS resolution (Stuhne and Steinman 2004; Valen-Sendstad and Steinman 2014) and are more than the mesh-independent values calculated for cerebral aneurysms both with and without devices deployed previously obtained by the authors (Peach et al. 2014a, b). Additionally, at this level, the number of mesh elements in each geometry measurement plane exceeds the recommendation of Jou et al. that is required to fully resolve flow features (Jou and Mawad 2011). The minimum cell densities of meshes are generated from 1 to 50,000 elements/mm3 to increase the levels of fineness. Aneurysm inflow is measured through a plane defined at the aneurysm neck. Mesh independence to within 1% for both aneurysm inflow (Q) and wall shear stress (WSS) magnitude was achieved, by meshing the geometries with a mesh density > 4000 elements per mm3.
CFD simulation results were post-processed in CFD-VIEW (ESI Group, Paris, France), by visualising the velocity streamlines and WSS distributions inside the aneurysm dome. Two planes are used to evaluate to monitor inflow: One is a neck plane placed proximal to the devices deployed that define the aneurysm boundary, and another plane is placed on the proximal daughter vessel perpendicular to the daughter vessel wall (Peach et al. 2017, 2019). WSS max is defined as the largest value of WSS in the region of interest, referring to the aneurysm sac or the aneurysm neck.
3 Results and discussion
3.1 Virtual deployment of devices
The virtual fast stenting algorithm is applied to all devices (Stents I–II and Contours I–II) in two patient-specific aneurysm models, as shown in Table 2, Figs. 6 and 7. Similar to the clinical procedure, placement of the flow-diverter begins from distal side to proximal of the device with respect to the aneurysm sac. Figures 6 and 7 show the overall configurations of the SVB and the Contour device for Case I and Case II at the three stages of the deployment process.
The fast virtual stenting algorithms are developed in Visual Studio 2019 (Microsoft, Albuquerque, New Mexico, USA) running on a single 2.60 GHz core without any parallelisation (i.e. multithreading). After approximately 48 iterations, each device in each model reached the deployed position, with a computation time of less than 1 min per case. The device deployed made good contact with the vessel wall, and the vessel wall was considered fixed during the deployment process.
3.2 Post- and pre-intervention hemodynamics
Calculations of inflow entering (Qin) through the aneurysm neck with no device deployed in each case are shown in Table 3 with values of 126.4 and 51.8 ml/min, respectively, representing approximately 100% and 69.63%, respectively, of the parent vessel (MCA) average flowrate. Velocity streamlines for No Device cases are shown in the top rows of Figs. 8 and 9.
In Case I, a relatively fast flow (\(\sim \) 0.5 ms−1) entered deeply into the aneurysm dome, whereas the flow at a slower speed (\(\sim \) 0.24 ms−1) was found at the entrance of the aneurysm dome in Case II. Most prominent in Case I, the jets of flow cause a region of increased WSS magnitude inside the dome of the aneurysm by the impact of the jet at the vessel wall of the aneurysm tip and the impingement flow leaving the neck of the aneurysm dome. The risk of aneurysm rupture is correlated with concentrated jet inflow, high aneurysm inflow rates, and regions of flow impingement and confirms the fragile nature of these particular aneurysm cases prior to intervention (Hassan et al. 2005; Castro et al. 2009; Doddasomayajula et al. 2017).
In Case I, a significant reduction in inflow entering the aneurysm sac by around 89.8% after the deployment of both devices (Table 3). The inflow reduction has been accomplished by eliminating the inflow jet almost completely and by decreasing the flow velocity that does enter the aneurysm substantially (second row, Fig. 8).
Compared to the performance of the SVB, the Contour device shows a slightly better inflow reduction ability with a reduction of 91.1% inflow entering the dome of the aneurysm (Table 3). Also, this device almost eliminates the aneurysm inflow jet with only isolated flow streamlines entering the left-hand side of the aneurysm (third row, Fig. 8). Besides, the velocity streamlines of Case I (A and B) further demonstrate that the right-hand daughter vessel serves as the primary outlet for the majority of the flow entering the aneurysm dome for both devices.
The change of flow pattern has effectively resolved a significant portion of the flow impingement observed in the No Device case, where a portion of the flow exits aneurysm through the left-hand daughter vessel resulting in elevated WSS max (\(\sim \) 8 Pa) at the neck of the aneurysm. The WSS magnitude plot for Case I-A in Fig. 8 confirms that the introduction of the device has eliminated the flow jet and impingement region entirely. Consequently, the WSS throughout the entire aneurysm dome remains at approximately 2 Pa, which is a typical value for healthy blood vessels.
There are similarities in the performance of both devices for Case II as summarised in Table 3. Deployment of the SVB and Contour devices again results in an inflow reduction of around 90%. When comparing the first and second rows of Fig. 9, it can be observed that the results obtained from the SVB are similar to the previous findings. The presence of jets (> 0.25 ms−1) entering the aneurysm is reduced. The complex and impinging flow within the aneurysm dome is now transformed into a relatively simple circulating flow pattern, where blood enters the aneurysm from the right-hand side and exits through the left-hand side. Notably, the high WSS (max magnitude > 8 Pa) concentrated region seen close to the amorphous neck of the aneurysm in the No Device case is reduced dramatically after the device placements of SVB and the Contour device. However, the WSS max magnitude remains high in the left-hand daughter vessel after placement of the Contour device, whereas the WSS max magnitude is reduced to an average of 2 Pa in the region of the SVB placement location in the left-hand daughter vessel. Similar to Case I, the elimination of increased WSS region and impinging flow would suggest that both devices substantially reduce the rupture risk of aneurysm and facilitate aneurysm stabilisation and thrombosis.
3.3 The effect on the daughter vessel
Table 4 shows a 20.60% increase in the daughter vessel flow rate after Contour device deployment, whereas there is a 15.38% decrease in the daughter vessel flow rate when SVB is deployed. As for Case II, there is a 3.37% increase in the daughter vessel flow rate after Contour deployment compared to an appropriately 1.81% decrease in the daughter vessel flow rate. The average daughter vessel flow rate increased by 11.99% for the Contour device and the average flow reduction in daughter vessel flow rate is 8.60% for the SVB device.
An average change of 2.59% is observed in the non-stented daughter vessel flow rate after the deployment of both the Contour device and SVB in Case II, whereas there is a significant change (average change of 17.99%) in the non-stented daughter vessel flow rate after the deployment of both devices in Case I. This indicates that the aneurysm morphology might play a vital role in affecting fluid alternation to the device implantation. The large non-stented daughter vessel diameter ratios (SDRs; ratios between non-stented daughter vessel and main vessel diameters), small bifurcating angles (angle between main and non-stented daughter vessel) in Case I might contribute to the significant change in the non-stented daughter vessel flow rate.
Additionally, the average Relative Change (%) after the deployment of the Contour device is around 11.99%, whereas the average Relative Change (%) is reduced by 8.60% after the implantation of SVB. One possible cause is that the Contour device tends to “cover” the neck of an aneurysm which leads to more flow into the non-stented daughter vessel. On the other hand, SVB tends to lead more flow into the stented daughter vessel. This result illustrates that the Contour device offers the possibility of a flow diversion effect at the aneurysm neck, without the drawbacks of daughter vessel occlusion. However, without incorporating biochemical models to capture device endothelialisation the actual occlusion of the daughter vessel cannot be predicted.
There are several limitations to this study. First, the number of cases used in this proof-of-concept study is modest, with the future improvement in the automated 3D reconstruction algorithm and simulation pipeline more cases are warranted. Second, the movement of vessels during the cardiac cycle was not considered. (4D simulation) with time makes the simulation realistic at the expense of computational time. Third, the use of average flow rates and fixed pressure outlet boundary condition may be considered unrealistic, and certainly it reduces the patient-specific nature of the simulations. The Windkessel model can be introduced to facilitate the "correction" of left-to-right vessel flow imbalance which can improve the accuracy of the flow model. A time-dependent solution can provide more hemodynamic parameters such as Relative Residence Time (RRT) for the analysis of successful aneurysm occlusion, but as a trade-off, the computation time will increase significantly. Extending this study to include proper pulsatile flow simulations will certainly improve the accuracy and relevance of the implantation effect predictions, and it is a straightforward step, as it has been demonstrated in the past. We would like to emphasise again the focus of this paper, which is the comparison between scenarios and the speed of the interactive part of the planning cycle, i.e. the actual device implantation phase. Fourth, mechanical contact analysis should be carried out in the future for more accurate modelling of the vessel wall in response to the deployment of the stent. Fifth, the virtual stenting process can be modelled within a few seconds, but the steady-state CFD computation typically requires 1–1.5 h to complete. A deep learning method such as a physics-informed neural network (PINN) can be integrated into this work to further reduce the CFD computation time. As for more complicated cases such as fusiform aneurysms treated with a peripheral stent using spiral flow technology, transient simulations can be carried out in the future to provide more realistic flow predictions and to facilitate better understanding of potential treatment strategies.
4 Conclusions
In this work, a fast virtual stenting algorithm developed in our previous study (Peach et al. 2014a; Spranger and Ventikos 2014; Lyu et al. 2021) was used to deploy different types of endovascular devices into the patient-specific middle cerebral artery (MCA) bifurcation aneurysms. The algorithms proposed previously achieved fast and accurate device deployment which allows for the device positioning and size optimisation. Different deployment strategies and positions may also be tested in future work.
Regardless of the device type, aneurysm inflow is decreased in both of the patient-specific MCA bifurcation aneurysms by at least ~ 85%, and areas of high WSS from flow jetting and impingement are eliminated, all features associated with successful aneurysm isolation. These results suggest that the performance of SVB and Contour device is sufficient for the treatment of these patient-specific bifurcation aneurysms.
The comparison of SVB and Contour device in two cases suggests that the Contour device shows a slightly better inflow reduction ability with an average of 91.25% compared to 87.95% of SVB. However, the placement of SVB shows better side branch WSS magnitude reduction around the placement region. The study on the daughter vessel flow change suggests Contour device deployment increases the daughter vessel flow rate by an average of 11.99%, whereas SVB deployment reduces the daughter vessel flow rate by an average of 8.60%.
Data accessibility
All data generated or analysed during this study are included in the manuscript and supplementary tables and figures. The raw data that support the findings of this study are available from the corresponding author, [YV], upon request.
References
Akhunbay-Fudge CY, Deniz K, Tyagi AK, Patankar T (2020) Endovascular treatment of wide-necked intracranial aneurysms using the novel contour neurovascular system: a single-center safety and feasibility study. J NeuroInterv Surg 12:987–992. https://doi.org/10.1136/neurintsurg-2019-015628
Alfano JM, Kolega J, Natarajan SK et al (2013) Intracranial aneurysms occur more frequently at bifurcation sites that typically experience higher hemodynamic stresses. Neurosurgery 73:497–505. https://doi.org/10.1227/NEU.0000000000000016
Arthur AS, Molyneux A, Coon AL et al (2019) The safety and effectiveness of the Woven EndoBridge (WEB) system for the treatment of wide-necked bifurcation aneurysms: final 12-month results of the pivotal WEB Intrasaccular Therapy (WEB-IT) Study. J NeuroInterv Surg 11:924–930. https://doi.org/10.1136/neurintsurg-2019-014815
Bhogal P, Lylyk I, Chudyk J et al (2021) The Contour—early human experience of a novel aneurysm occlusion device. Clin Neuroradiol 31:147–154. https://doi.org/10.1007/s00062-020-00876-4
Bhogal P, Udani S, Cognard C et al (2019) Endosaccular flow disruption: where are we now? J NeuroInterv Surg 11:1024–1025. https://doi.org/10.1136/neurintsurg-2018-014623
Biondi A, Primikiris P, Vitale G, Charbonnier G (2023) Endosaccular flow disruption with the contour neurovascular system: angiographic and clinical results in a single-center study of 60 unruptured intracranial aneurysms. J NeuroInterv Surg 15:838–843. https://doi.org/10.1136/jnis-2022-019271
Boite Y, Suaiden Klein T, de Andrade MR, Wajnberg E (2023) Numerical simulation of flow-diverting stent: comparison between branches in bifurcation brain aneurysm. Biomech Model Mechanobiol 22:1801–1814. https://doi.org/10.1007/s10237-023-01733-2
Castro MA, Putman CM, Sheridan MJ, Cebral JR (2009) Hemodynamic patterns of anterior communicating artery aneurysms: a possible association with rupture. Am J Neuroradiol 30:297–302. https://doi.org/10.3174/ajnr.A1323
Doddasomayajula R, Chung B, Hamzei-Sichani F et al (2017) Differences in hemodynamics and rupture rate of aneurysms at the bifurcation of the basilar and internal carotid arteries. Am J Neuroradiol 38:570–576. https://doi.org/10.3174/ajnr.A5088
Eker OF, Levrier O, Houdart E et al (2022) Endovascular treatment of complex intracranial aneurysms with LEO+ stents: the LEO II cohort study. Front Neurol 13:848683. https://doi.org/10.3389/fneur.2022.848683
Goyal N, Hoit D, DiNitto J et al (2020) How to WEB: a practical review of methodology for the use of the Woven EndoBridge. J NeuroInterv Surg 12:512–520. https://doi.org/10.1136/neurintsurg-2019-015506
Hassan T, Timofeev EV, Saito T et al (2005) A proposed parent vessel geometry—based categorization of saccular intracranial aneurysms: computational flow dynamics analysis of the risk factors for lesion rupture. J Neurosurg 103:662–680. https://doi.org/10.3171/jns.2005.103.4.0662
Hossain SS, Starosolski Z, Sanders T et al (2021) Image-based patient-specific flow simulations are consistent with stroke in pediatric cerebrovascular disease. Biomech Model Mechanobiol 20:2071–2084. https://doi.org/10.1007/s10237-021-01495-9
Jou L-D, Mawad ME (2011) Hemodynamic effect of Neuroform stent on intimal hyperplasia and thrombus formation in a carotid aneurysm. Med Eng Phys 33:573–580. https://doi.org/10.1016/j.medengphy.2010.12.013
Lawson A, Goddard T, Ross S et al (2017) Endovascular treatment of cerebral aneurysms using the Woven EndoBridge technique in a single center: preliminary results. J Neurosurg 126:17–28. https://doi.org/10.3171/2015.4.JNS142456
Liebig T, Kabbasch C, Strasilla C et al (2015) Intrasaccular flow disruption in acutely ruptured aneurysms: a multicenter retrospective review of the use of the WEB. AJNR Am J Neuroradiol 36:1721–1727. https://doi.org/10.3174/ajnr.A4347
Liebig T, Killer-Oberpfalzer M, Gal G et al (2022) The safety and effectiveness of the contour neurovascular system (Contour) for the treatment of bifurcation aneurysms: the CERUS study. Neurosurgery 90:270. https://doi.org/10.1227/NEU.0000000000001783
Liu H, Lan L, Abrigo J et al (2021) Comparison of newtonian and non-newtonian fluid models in blood flow simulation in patients with intracranial arterial stenosis. Front Physiol 12:718540. https://doi.org/10.3389/fphys.2021.718540
Lonsdale RD (1993) An algebraic multigrid solver for the Navier–Stokes equations on unstructured meshes. Int J Numer Methods Heat Fluid Flow 3:3–14. https://doi.org/10.1108/eb017512
Lyu M, Ventikos Y, Peach TW et al (2021) Virtual Flow-T Stenting for two patient-specific bifurcation aneurysms. Front Neurol 12:726980. https://doi.org/10.3389/fneur.2021.726980
Mokin M, Primiani CT, Ren Z et al (2020) Stent-assisted coiling of cerebral aneurysms: multi-center analysis of radiographic and clinical outcomes in 659 patients. J NeuroInterv Surg 12:289–297. https://doi.org/10.1136/neurintsurg-2019-015182
Narata AP, de Moura FS, Larrabide I et al (2018) The role of hemodynamics in intracranial bifurcation arteries after aneurysm treatment with flow-diverter Stents. Am J Neuroradiol 39:323–330. https://doi.org/10.3174/ajnr.A5471
Ni M-J, Abdou MA (2007) A bridge between projection methods and SIMPLE type methods for incompressible Navier–Stokes equations. Int J Numer Meth Eng 72:1490–1512. https://doi.org/10.1002/nme.2054
Peach TW, Ngoepe M, Spranger K et al (2014a) Personalizing flow-diverter intervention for cerebral aneurysms: from computational hemodynamics to biochemical modeling. Int J Numer Methods Biomed Eng 30:1387–1407. https://doi.org/10.1002/cnm.2663
Peach TW, Cornhill JF, Nguyen A et al (2014b) The ‘Sphere’: a dedicated bifurcation aneurysm flow-diverter device. Cardiovasc Eng Tech 5:334–347. https://doi.org/10.1007/s13239-014-0188-4
Peach TW, Ricci D, Ventikos Y (2019) A virtual comparison of the eCLIPs device and conventional flow-diverters as treatment for cerebral bifurcation aneurysms. Cardiovasc Eng Tech 10:508–519. https://doi.org/10.1007/s13239-019-00424-3
Peach TW, Spranger K, Ventikos Y (2016) Towards predicting patient-specific flow-diverter treatment outcomes for bifurcation aneurysms: from implantation rehearsal to virtual angiograms. Ann Biomed Eng 44:99–111. https://doi.org/10.1007/s10439-015-1395-3
Peach TW, Spranger K, Ventikos Y (2017) Virtual flow-diverter treatment planning: the effect of device placement on bifurcation aneurysm haemodynamics. Proc Inst Mech Eng H 231:432–443. https://doi.org/10.1177/0954411916673674
Pierot L, Liebig T, Sychra V et al (2012) Intrasaccular flow-disruption treatment of intracranial aneurysms: preliminary results of a multicenter clinical study. Am J Neuroradiol 33:1232–1238. https://doi.org/10.3174/ajnr.A3191
Piotin M, Biondi A, Sourour N et al (2018) The LUNA aneurysm embolization system for intracranial aneurysm treatment: short-term, mid-term and long-term clinical and angiographic results. J NeuroInterv Surg 10:e34–e34. https://doi.org/10.1136/neurintsurg-2018-013767
Raychev R, Tateshima S, Vinuela F et al (2016) Predictors of thrombotic complications and mass effect exacerbation after pipeline embolization: the significance of adenosine diphosphate inhibition, fluoroscopy time, and aneurysm size. Interv Neuroradiol 22:34–41. https://doi.org/10.1177/1591019915609125
Saatci I, Yavuz K, Ozer C et al (2012) Treatment of intracranial aneurysms using the pipeline flow-diverter embolization device: a single-center experience with long-term follow-up results. AJNR Am J Neuroradiol 33:1436–1446. https://doi.org/10.3174/ajnr.A3246
Saleme S, Iosif C, Ponomarjova S et al (2014) Flow-diverting stents for intracranial bifurcation aneurysm treatment. Neurosurgery 75:623–631. https://doi.org/10.1227/NEU.0000000000000522
Spranger K, Ventikos Y (2014) Which spring is the best? Comparison of methods for virtual stenting. IEEE Trans Biomed Eng 61:1998–2010. https://doi.org/10.1109/TBME.2014.2311856
Stuhne GR, Steinman DA (2004) Finite-element modeling of the hemodynamics of stented aneurysms. J Biomech Eng 126:382–387. https://doi.org/10.1115/1.1762900
Thormann M, Mpotsaris A, Behme D (2022) Treatment of a middle cerebral artery bifurcation aneurysm with the novel contour neurovascular system compatible with 0.021″ catheters. Neuroradiol J 35:396–398. https://doi.org/10.1177/19714009211041523
Valencia A, Solis F (2006) Blood flow dynamics and arterial wall interaction in a saccular aneurysm model of the basilar artery. Comput Struct 84:1326–1337. https://doi.org/10.1016/j.compstruc.2006.03.008
Valencia AA, Guzmán AM, Finol EA, Amon CH (2006) Blood flow dynamics in saccular aneurysm models of the basilar artery. J Biomech Eng 128:516–526. https://doi.org/10.1115/1.2205377
Valen-Sendstad K, Steinman DA (2014) Mind the gap: Impact of computational fluid dynamics solution strategy on prediction of intracranial aneurysm hemodynamics and rupture status indicators. Am J Neuroradiol 35:536–543. https://doi.org/10.3174/ajnr.A3793
Van Doormaal JP, Raithby GD (1984) Enhancements of the simple method for predicting incompressible fluid flows. Numer Heat Transf 7:147–163. https://doi.org/10.1080/01495728408961817
Xu P, Liu X, Zhang H et al (2018) Assessment of boundary conditions for CFD simulation in human carotid artery. Biomech Model Mechanobiol 17:1581–1597. https://doi.org/10.1007/s10237-018-1045-4
Youssef PP, Iii DD, Peterson J et al (2021) Woven EndoBridge (WEB) device in the treatment of ruptured aneurysms. J NeuroInterv Surg 13:443–446. https://doi.org/10.1136/neurintsurg-2020-016405
Zhang X-J, Gao B-L, Hao W-L et al (2018) Presence of anterior communicating artery aneurysm is associated with age, bifurcation angle, and vessel diameter. Stroke 49:341–347. https://doi.org/10.1161/STROKEAHA.117.019701
Zhang Y, Yang M, Zhang H et al (2017) Stent-assisted coiling may prevent the recurrence of very small ruptured intracranial aneurysms: a multicenter study. World Neurosur 100:22–29. https://doi.org/10.1016/j.wneu.2016.12.107
Acknowledgements
The authors would like to acknowledge the financial and technical support provided by the ESI Group in association with the use of its CFD-ACE+ software package throughout this research. The authors would like to thank Shukun (Beijing) Network Technology Co., Ltd. for providing additional computing resources and technical support.
Funding
Open Access funding enabled and organized by CAUL and its Member Institutions. This study is supported by the National Key R&D Program of China (2018AAA0102600).
Author information
Authors and Affiliations
Contributions
ML involved in conceptualisation, data curation, formal analysis, funding acquisition, investigation, methodology, writing original manuscript, reviewing, and editing; RT took part in investigation, methodology, reviewing, and editing; CL participated in formal analysis, investigation, methodology, reviewing, and editing; TP involved in formal analysis, data analysis, reviewing, and editing; PB involved in data curation, study design, reviewing, and editing; LM took part in data curation, reviewing, and editing; QL involved in formal analysis, writing original manuscript, reviewing, and editing; DC took part in conceptualisation, funding acquisition, formal analysis, investigation, supervision, reviewing, and editing; YV involved in conceptualisation, funding acquisition, formal analysis, investigation, supervision, reviewing, and editing.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethical approval
Ethical review and approval were not required for the study on human participants in accordance with the local legislation and institutional requirements. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lyu, M., Torii, R., Liang, C. et al. Treatment for middle cerebral artery bifurcation aneurysms: in silico comparison of the novel Contour device and conventional flow-diverters. Biomech Model Mechanobiol (2024). https://doi.org/10.1007/s10237-024-01829-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10237-024-01829-3