Figures
Abstract
Cerebral aneurysms are a major risk factor for intracranial bleeding with devastating consequences for the patient. One recently established treatment is the implantation of flow-diverters (FD). Methods to predict their treatment success before or directly after implantation are not well investigated yet. The aim of this work was to quantitatively study hemodynamic parameters in patient-specific models of treated cerebral aneurysms and its correlation with the clinical outcome. Hemodynamics were evaluated using both computational fluid dynamics (CFD) and phase contrast (PC) MRI. CFD simulations and in vitro MRI measurements were done under similar flow conditions and results of both methods were comparatively analyzed. For preoperative and postoperative distribution of hemodynamic parameters, CFD simulations and PC-MRI velocity measurements showed similar results. In both cases where no occlusion of the aneurysm was observed after six months, a flow reduction of about 30-50% was found, while in the clinically successful case with complete occlusion of the aneurysm after 6 months, the flow reduction was about 80%. No vortex was observed in any of the three models after treatment. The results are in agreement with recent studies suggesting that CFD simulations can predict post-treatment aneurysm flow alteration already before implantation of a FD and PC-MRI could validate the predicted hemodynamic changes right after implantation of a FD.
Citation: Sindeev S, Arnold PG, Frolov S, Prothmann S, Liepsch D, Balasso A, et al. (2018) Phase-contrast MRI versus numerical simulation to quantify hemodynamical changes in cerebral aneurysms after flow diverter treatment. PLoS ONE 13(1): e0190696. https://doi.org/10.1371/journal.pone.0190696
Editor: Josué Sznitman, Technion Israel Institute of Technology, ISRAEL
Received: May 18, 2017; Accepted: December 19, 2017; Published: January 5, 2018
Copyright: © 2018 Sindeev et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data can be found at the Technical University of Munich University Library repository with the following link: https://mediatum.ub.tum.de/1359655.
Funding: This work was supported by Russian Science Foundation (Project 16-15-10327). This work was supported by the German Research Foundation (DFG) and the Technical University of Munich (TUM) in the framework of the Open Access Publishing Program. We thank Acandis GmbH & Co. KG for providing the DERIVO Embolisation Device and silicone models of aneurysms for experimental studies. The reported study was supported by the Supercomputing Center of Lomonosov Moscow State University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: We have the following interests: Acandis GmbH & Co. KG provided the DERIVO Embolisation Device and silicone models of aneurysms for experimental studies. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials, as detailed online in the guide for authors.
Introduction
A flow-diverter (FD) is a promising tool for the treatment of wide-necked and fusiform aneurysms which has a low complication rate [1–4]. Nevertheless the prediction of the treatment success is complicated since the flow parameters which influence the clinical outcome are still not well understood. Commonly in clinical practice a post-treatment contrast agent stasis is considered to assess the flow diverting effect immediately after the FD placement. However the observed flow stasis could not be a predictor of a complete aneurysm occlusion in some cases [5]. Therefore the clinical decision-making process is largely dependent on global risk factors and on experience and qualification of the neurosurgeon or neuroradiologist [6].
A number of studies, both experimental and numerical, have been done during recent years to find reliable hemodynamic parameters that could evaluate the FD performance and predict the clinical outcome [7–9]. Tsang et al [10] found a correlation between a significant change in turnover time after stenting and aneurysm occlusion, suggesting that it could be one of the criteria to evaluate the flow diversion effect and clinical success. Similarly Janiga et al [11] used the turnover time (residence time) to analyze different treatment scenarios and select an appropriate stent model and position leading to maximal increase in residence time.
Additionally some other parameters taking into account the flow pattern inside an aneurysm sac could be considered. In an experimental study by Balasso et al [12] the FD performance was assessed by complex analysis of a change of maximum velocity at the inflow zone, in the dome and at the outflow zone of the aneurysm. Moreover in the study by Suzuki et al [13] the relative reduction of the spatial-averaged and maximum velocity within the aneurysm sac as well as the maximum wall shear stress (WSS) were considered to analyze the flow reduction effects after the treatment. Another alternative, proposed by Lieber et al, is an assessment of the stenting effect according to kinetic energy reduction [14] and change of the vorticity inside an aneurysm [15].
Despite recent achievements of numerical and experimental methods for evaluation of the FD performance, their application in clinical practice is still constrained due to their complexity and time restrictions. For example, evaluation of post-treatment blood flow immediately after the stent placement using CFD methods is problematic, since typical numerical simulation with FD could take several weeks depending on the size of the computational mesh. In the same time several feasibility studies showed that MRI could serve as a tool for the clinical assessment of cerebral hemodynamics which could rapidly evaluate the intra-aneurysmal flow field with sufficient precision [16–18]. However Pereira et al [19] reported technical issues in the measurement of low velocities and the need for more validation work in a feasibility study regarding post-treatment evaluation of blood flow in aneurysms after FD placement using MRI. Further MacDonald et al [20] assessed the pre- and postoperative hemodynamic parameters for a giant cerebral aneurysm and found that the intra-aneurysmal pressure derived from MRI measured velocity field differed only by 6.1% from the pressure measured by invasive techniques. However there is a lack of studies which evaluate hemodynamic changes in cerebral aneurysms after a FD placement using both experimental and numerical methods.
The purpose of this work was to quantitatively study the distribution of hemodynamic parameters in patient-specific models of cerebral aneurysms before and after treatment with flow diverters using PC-MRI and CFD and to correlate these results with the clinical outcome.
Materials and methods
Ethics committee approval
As our experiments were conducted retrospectively with permanently anonymized patient data, our local ethics committee deemed the study exempt from the requirement for approval.
Clinical cases
Twenty three patients treated with a FD at the department of Neuroradiology at the Klinikum rechts der Isar between 2008 and 2014 were screened for availability of pre-treatment high spatial resolution 3D angiographic data. These datasets were available for three patients, who were included in this study and subgrouped regarding their clinical treatment outcome. This study was HIPAA compliant and in line with the local ethical and legislative requirements. Patient consent was not required due to the fully anonymous retrospective analysis.
Group A consisted of one patient (A01) with a clinically successful treatment. This aneurysm showed no perfusion 6 months after FD placement. Group B included two patients (B01, B02) who showed persistent perfusion of the aneurysm for more than 6 months after treatment. A detailed information about the studied cases is presented in Table 1.
Experimental phantoms
Realistic silicone phantoms were produced for each patient. The patient’s pre-intervention CT angiographic data was used to segment the aneurysm and adjacent vessels to generate STL models. These STL models served as template for a high precision 3D wax printer. The silicone Elastosil 601 was molded around the wax which was later removed using a dissolvent. For each case 2 silicone phantoms were manufactured by Acandis (Pforzheim, Germany).
One phantom of each aneurysm was used to study preoperative hemodynamics. The other one was used to study hemodynamic changes after the placement of the corresponding FD. The aneurysm phantoms are presented in Fig 1. These phantoms were used for MRI flow measurements. The dimensions of the studied phantoms are presented in Table 2.
Experimental setup
The measurement setup (Fig 2) included a computer controlled piston pump (8). The valve (7) mimicked the aortic valve. The compliance chamber (6) simulated the windkessel function of the aorta. The experimental fluid was a 58% aqueous glycerol mixture, which exhibits Newtonian fluid behavior. The flow rate (5), the inlet pressure (4) and the outlet pressure (10) were measured with the corresponding sensors and recorded by the measurement computer (13). The aneurysm phantom (3) was placed inside a MRI hand wrist coil (2) which allowed for better signal to contrast ratio and thus spatial resolution. The MRI (1) was synchronized with the piston pump (8) using a photo diode (9) which provided a trigger signal. The pressure chamber (12) was used to provide a constant diastolic pressure. The realistic inlet flow rate of 7 l/h was used, which corresponded to the average velocity of 0.25 m/s for B01, 0.14 m/s for B02 and 0.19 m/s for A01 at the expanded inlet segment during the systolic peak.
1—MRI Philips Ingenia 3T; 2—MRI handwrist-coil; 3—aneurysm phantom; 4—inlet pressure sensor; 5—flow sensor; 6—compliance chamber; 7—valve; 8—piston pump; 9—photo diode; 10—outlet pressure sensor; 11—LED; 12—pressure chamber; 13—measurement computer; 14—MRI computer.
The measurement protocol is available as supporting information file (S1 Protocol) or at protocols.io (http://dx.doi.org/10.17504/protocols.io.kekctcw).
MRI settings
A 3T Philips Ingenia (Philips Healthcare, Netherlands) was used as MRI velocity measuring system. The aneurysm phantoms were placed inside a hand wrist coil to maximize spatial resolution and signal-to-noise ratio (SNR). A 4D phase contrast sequence with a multi-shot 3D Gradient-Echo sequence with a turbo-factor of 5 without SENSE was used for image acquisition. The acquisition parameters were set to: TE = 2.9 ms, TR = 6.1 ms, α = 10°. Flow in all 3 spatial directions (ap, rl, cc) were acquired consecutively in one single scan. Total scan duration was 30:44 min.
The quantitative velocity field was reconstructed with corrections for magnetic field inhomogeneities by static measurements. The spatial resolution was RLxAPxFH = 0.53x0.53x0.60 mm3 with a matrix size of 112x112x94 voxels. The temporal resolution was 71 ms. After inspection of the CFD velocity data velocity encoding level (VENC) was set to 100 cm/s for each spatial direction to prevent aliasing in the aneurysm region. Partial volume effects close to the vessel wall were neglected as the drop in magnitude signal for all cases with and without FD was steep and abrupt.
CFD settings
The Navier-Stokes equations for incompressible fluid were utilized for CFD simulations. The inlet velocity curves were obtained from the data of the MRI velocity measurements. The fluid density was set to ρ = 1141 kg/m3 which corresponds to the density of the experimental fluid. A Newtonian model was used to represent the fluid behavior. The viscosity was set to η = 4.1 mPa ⋅ s, which corresponded to the used experimental fluid.
For the numerical studies the geometrical models of the aneurysms were created using corresponding STL-models of the silicone phantoms. The geometrical models of Acandis DERIVO FDs were created according to the manufacturer specifications. Fast virtual stenting technique was used for the virtual deployment of the FDs in the aneurysm models. Here, nominal diameters, initial lengths, strut diameters as well as strut angles were respected. Hence, arbitrary configurations can be reproduced. This explicit method considers the individual stent pores and is therefore clearly superior compared to simplified approaches such as the assumption of a porous medium to account for flow-diverting effects. The virtually deployed FDs were visually compared with the available in-vitro deployments and showed sufficient correlation. A detailed description of the virtual stenting technique was presented in the studies by Berg et al [21, 22]. Furthermore, a validation of the fast virtual stenting approach can be found in the study by Janiga et al. [11].
The hexahedral computational meshes were generated using snappyHexMesh tool from OpenFOAM CFD Toolbox (CFD Direct, Caversham, England). Additional refinement of the cells around the FD was used to properly represent the minimal elements of FD braid. A mesh independence test showed that generated meshes were sufficient to capture all significant features of the flow pattern in the studied aneurysms. Five cardiac cycles were simulated to omit the initial perturbations of the flow field. The results for the last cardiac cycle were used for analysis.
Analysis of hemodynamic parameters
The moment of systolic peak was selected to analyze the hemodynamic parameters in the aneurysms. This moment was used since hemodynamic parameters reach their maximum values at this time. The flow pattern in the central cross-section of each aneurysm was analyzed in detail as well as the WSS distribution, streamlines, kinetic energy and recirculation. Additionally, the hemodynamic changes after FD placement were analyzed and correlated with the available clinical outcome.
The WSS was computed to analyze the shear stress distribution over the aneurysm sac. To derive the WSS from the MRI-measured velocity field, the wall shear rate was computed, which then was multiplied by dynamic viscosity η. The normalized WSS (nWSS) was used to compare the WSS distribution between CFD and MRI: where WSS(x, y, z) is the WSS magnitude at the point with coordinates (x, y, z); WSSmax is the maximum WSS magnitude over the aneurysm sac.
To evaluate a kinetic energy in the aneurysm sac we used a method proposed by Seong et al. [14, 23], according to which a measure of the intra-aneurysmal kinetic energy E could be estimated as a sum of squares of velocity magnitudes ui at every i-th measured point (x, y, z) in the aneurysm sac, i.e.: (1) where N is a number of measured points.
Mean kinetic energy over the cardiac cycle was employed to quantify hemodynamic changes after the FD placement, which is an average of sum of instantaneous kinetic energies over the cardiac cycle. Since CFD results have a much higher spatial and temporal resolution a comparison of absolute values of kinetic energy (1) between MRI and CFD is not valid. Therefore only relative kinetic energy reductions ER for CFD and MRI results were considered in this study: where EnoFD and EFD are the intra-aneurysmal kinetic energies before and after FD placement respectively.
Another hemodynamic parameters used to assess FD performance were peak and mean recirculation [15]. The instantaneous recirculation R could be obtained by computing vorticity ω of the intra-aneurysmal velocity field: (2)
The mean recirculation Rm could be found by integrating Eq (2) over the cardiac cycle: where T is the length of the cardiac cycle.
Results
Preoperative analysis
The central cross-sections of the aneurysm models were used to analyze the flow pattern in the aneurysm sac. The cross-sections are presented in Figs 3, 4 and 5 for cases B01, B02 and A01 respectively. The CFD simulations and MRI measurements predicted a similar main flow pattern for the all three cases. The region which belongs to the aneurysm sac itself, excluding the parent artery, was used to compute an average intra-aneurysmal blood flow. It was done in order to correctly compare the average intra-aneurysmal velocity before and after the treatment. The obtained average intra-aneurysmal velocity is presented in Table 3. The average difference between CFD and MRI results was 9.47% for B01; 6.1% for B02 and 6.21% for A01.
a) CFD simulation; b) MRI measurement.
a) CFD simulation; b) MRI measurement.
a) CFD simulation; b) MRI measurement.
noFD—before the treatment; FD—after the treatment; U—an average intra-aneurysmal velocity during the systolic peak (m/s); E—mean kinematic energy (m2/s2); R—mean recirculation (mm2/s).
The nWSS is shown in Fig 6 both for CFD and MRI. Similar general nWSS distributions were observed for the aneurysm models, but detailed variations in WSS were not adequately revealed by MRI. Regions of high nWSS (≥ 0.75) were found at the proximal and distal segments of the parent artery. Contrary, the aneurysm sac and aneurysm bleb were characterized by low values of nWSS (≤ 0.25). However the distribution of nWSS for MRI was influenced by noise and low spatial resolution.
Streamlines were computed from CFD data to detect vortices in the aneurysm sac. The computed streamlines for the studied aneurysms are presented in Fig 7. The preoperative blood flow in the aneurysms was characterized by the presence of a vortex in all three models. The swirling flow in the aneurysm sac produced the region of low velocities in the center of the aneurysms.
Postoperative analysis
The flow-patterns in the aneurysms were decidedly changed after FD placement. The central cross-sections for the treated aneurysms are presented in Figs 8, 9 and 10 for cases B01, B02 and A01 respectively. Both CFD and MRI showed a general similarity of flow-patterns for the treated cases. However the presence of metallic braiding of the FD led to non-visualizable regions in the central cross-sections for MRI data. The maximum velocity region was jailed by the braid of FD, which led to reduction of blood flow in the aneurysm sac. The obtained reduced intra-aneurysmal velocities for studied models are presented in Table 4. The observed flow reduction was about 52%, 28% and 80% for cases B01, B02 and A01 respectively.
a) CFD simulation; b) MRI measurement.
a) CFD simulation; b) MRI measurement.
a) CFD simulation; b) MRI measurement.
Also nWSS was computed to determine the changes in WSS distribution after treatment. No significant changes in distribution pattern were found (Fig 11). However the average WSS magnitude was decidedly reduced in all three cases. Streamlines were computed to determine the changes of flow patterns after FD placement. The computed streamlines are presented in Fig 12. The major part of the flow was redirected along the FD, while the intra-aneurysmal flow was decidedly reduced. The vortex in the aneurysm sac disappeared in all three models.
Discussion
In the present study we analyzed the hemodynamic changes in cerebral aneurysms of 3 patients before and after FD treatment and correlated them with their known clinical outcome. Similar velocities and flow patterns were obtained for both CFD and MRI methods while substantial differences were found for calculated WSS. We observed a larger flow reduction for the successfully treated aneurysm, while a smaller reduction was found for the two cases where additional treatment was needed for occlusion.
A relative flow reduction is considered a major parameter correlated with the clinical outcome of the treatment. Xiang et al [24] found a reduction of average aneurysmal velocity by 76% for the aneurysm which was occluded within the first three months and a reduction by 40% for the aneurysms that occluded after six months. It was proposed that the post-treatment reduction of average intra-aneurysmal velocity could be correlated with aneurysm occlusion. Also the flow reduction is used to evaluate the flow-diversion efficiency and is considered as the objective function for optimization of FDs [25]. However Cebral et al [26] demonstrated that the reduction of intra-aneurysmal velocity was not sufficient for a successful treatment of aneurysms with a FD. Despite the decreased velocity in the aneurysm sac, an increase of intra-aneurysmal pressure after the FD placement was observed for the ruptured cases.
The results of the present study are in agreement with recent studies, which proposed that for a successful treatment with a FD the flow reduction should be at least one-third of the preoperative state [27–29]. On the other side, a flow reduction of 52% was observed for the the giant fusiform aneurysm, where an immediate clinical success was not achieved. It demonstrates that the clinical outcome is determined not only by the relative flow reduction, but it depends on a complex set of parameters, which characterize the geometry of the aneurysm [30, 31], the hemodynamics in the aneurysm sac [32], the aneurysm wall structure and blood coagulation properties [33–35].
The preoperative state of all three aneurysms was characterized by the presence of a vortex in the aneurysm sac. A zone of low velocities was observed in the center of the vortex, which could contribute to thrombus formation [36–38]. The placement of FD led to significant changes in flow-patterns for the all three cases. No vortex was observed in the aneurysms sac after the treatment. The observed relative reduction of recirculation in the studied cases is in agreement with results of Lieber et al [14] suggesting that clinical success cannot be determined by using only a recirculation as a criterion.
The virtual deployment of a FD enables a precise description of stent-induced hemodynamics at no risk for the individual patient. In the frame of this study, a fast virtual stenting approach was chosen based on Janiga et al. [39]. Its advantage clearly lies in the clinical applicability, since realistic deformations are performed within seconds. Furthermore, the complete geometric appearance of the stent (e.g., each single stent strut and individual pores) is considered. In contrast, other approaches reported in literature are either highly computationally expensive or over-simplify those devices due to the consideration of porous media. In this regard, the recent virtual FD deployments are well suited for the comparisons with the real stentings of this study.
WSS distribution is another important parameter which is commonly analyzed for evaluation of pre- and postoperative hemodynamics in cerebral aneurysms. There is still no clear theory whether high or low WSS values are prone to aneurysm growth and rupture [40–42]. However in some cases a correlation of high WSS with the rupture of saccular aneurysms is observed [9, 43]. In our study only the distribution of nWSS was considered, since a precise calculation of WSS using the PC-MRI measurements was problematic due to the lack of spatial resolution, especially close to the aneurysm wall, which correlates with results reported by Petersson et al [44]. Therefore the absolute values of WSS were different for CFD and MRI results and small alterations in WSS could not be revealed by MRI, however the distribution of normalized WSS showed a similar pattern for both methods, which is in agreement with the results of the recent study by Cibis et al [45].
This study has some limitations. A Newtonian fluid was used for investigation of the intra-aneurysmal hemodynamics. However recent studies showed that non-Newtonian fluid behavior could play a significant role in distribution of hemodynamic parameters in the aneurysm, especially at zones of recirculating flow [46]. The use of rigid phantoms also could lead to overestimation of hemodynamic parameters such as peak velocity magnitude and WSS [47]. Additionally no patient-specific boundary conditions were available for the studied cases, which could have a substantial effect on velocity calculations [48, 49]. However, the comparison of CFD and MRI should not be affected by these shortcomings as similar parameters were used for both methods. For better comparison within our simulations, we implanted similar FDs (DERIVO) in all the different models. As sizing of the in-vitro FD was determined based on an angiography of the respective phantom, slight differences in FD dimensions occurred compared to the patients. This may limit the comparability with the in-vivo cases, however little differences may be expected from different manufacturers and clinical corresponded well to our experimental studies. Furthermore, only a single successful patient and two unsuccessful ones were considered, limiting the generalizability of obtained individual results to other aneurysms with anatomic/geometric variability due to the limited sample size of selected patients. Future studies will address the listed limitations.
Conclusion
In the present work we quantitatively studied the hemodynamics in three patient-specific aneurysms before and after FD treatment using both experimental and numerical methods. For preoperative and postoperative distribution of hemodynamic parameters similarity between results obtained by CFD simulations and MRI velocity measurements was shown. In cases without aneurysm occlusion after FD placement a flow reduction of about 30-50% was found, while for the clinically successful case the flow reduction was 80%. These results are in agreement with recent studies of post-treatment complications after FD placement. Both CFD and MRI methods could be used for evaluation of cerebral hemodynamics before and after treatment and could be used complementary in supporting a physician during an intervention planning. While CFD may be advantageous in a pre-treatment planning to predict the hemodynamic changes after a FD implementation, in-vivo MRI measurements are not dependent on virtual stenting and thus could assess the blood flow right after FD deployment.
While velocity reduction is clearly one aspect of the mechanism of action, it is not solely responsible for treatment success. The treatment result is more complicated and likely involves an interaction between velocity reduction, hemodynamic changes and coagulation. While no simulation can realistically incorporate all of these variables due to many unknowns, MRI measurements could provide velocity information that may factor into treatment outcome.
Supporting information
S1 Protocol. MRI velocity measurement protocol.
https://doi.org/10.1371/journal.pone.0190696.s001
(PDF)
Acknowledgments
We thank Acandis GmbH & Co. KG for providing the DERIVO Embolisation Device and silicone models of aneurysms for experimental studies.
The reported study was supported by the Supercomputing Center of Lomonosov Moscow State University.
References
- 1. Byrne JV, Beltechi R, Yarnold JA, Birks J, Kamran M. Early experience in the treatment of intra-cranial aneurysms by endovascular flow diversion: a multicentre prospective study. PLoS ONE. 2010;5(9):e12492. pmid:20824070
- 2. Arrese I, Sarabia R, Pintado R, Delgado-Rodriguez M. Flow-diverter devices for intracranial aneurysms: systematic review and meta-analysis. Neurosurgery. 2013;73(2):193–200. pmid:23624409
- 3. Lubicz B, der Elst OV, Collignon L, Mine B, Alghamdi F. Silk flow-diverter stent for the treatment of intracranial aneurysms: a series of 58 patients with emphasis on long-term results. AJNR Am J Neuroradiol. 2014;36(3):542–546. pmid:25376806
- 4. Lv X, Yang H, Liu P, Li Y. Flow-diverter devices in the treatment of intracranial aneurysms: a meta-analysis and systematic review. Neuroradiol J. 2016;29(1):66–71. pmid:26838174
- 5. Tsang ACO, Fung AMY, Tsang FCP, Leung GKK, Lee R, Lui WM. Failure of flow diverter treatment of intracranial aneurysms related to the fetal-type posterior communicating artery. Neurointervention. 2015;10(2):60–66. pmid:26389008
- 6. Etminan N, Brown RD, Beseoglu K, Juvela S, Raymond J, Morita A, et al. The unruptured intracranial aneurysm treatment score: a multidisciplinary consensus. Neurology. 2015;85(10):881–889. pmid:26276380
- 7. Raschi M, Mut F, Byrne G, Putman CM, Tateshima S, Viñuela F, et al. CFD and PIV analysis of hemodynamics in a growing intracranial aneurysm. Int J Numer Method Biomed Eng. 2011;28(2):214–228.
- 8. Fukazawa K, Ishida F, Umeda Y, Miura Y, Shimosaka S, Matsushima S, et al. Using computational fluid dynamics analysis to characterize local hemodynamic features of middle cerebral artery aneurysm rupture points. World Neurosurg. 2015;83(1):80–86. pmid:23403347
- 9. Can A, Du R. Association of hemodynamic factors with intracranial aneurysm formation and rupture. Neurosurgery. 2016;78(4):510–520. pmid:26516819
- 10. Tsang AO, Tang AS, Chung W, Leung GK, Chow K. Correlating hemodynamic changes and occlusion time after flow diverter treatment of bilateral large internal carotid artery aneurysms. Clin Neuroradiol. 2016;26(4):477–480. pmid:26869443
- 11. Janiga G, Rossl C, Skalej M, Thevenin D. Realistic virtual intracranial stenting and computational fluid dynamics for treatment analysis. J Biomech. 2013;46(1):7–12. pmid:23063770
- 12. Balasso A, Bauer JS, Liebig T, Dorn F, Zimmer C, Liepsch D, et al. Evaluation of intra-aneurysmal hemodynamics after flow diverter placement in a patient-specific aneurysm model. Biorheology. 2014;51(6):341–354. pmid:25777299
- 13. Suzuki T, Takao H, Fujimura S, Dahmani C, Ishibashi T, Mamori H, et al. Selection of helical braided flow diverter stents based on hemodynamic performance and mechanical properties. J Neurointerv Surg. 2016; p. neurintsurg–2016. pmid:27646987
- 14. Seong J, Wakhloo AK, Lieber BB. In vitro evaluation of flow divertors in an elastase-induced saccular aneurysm model in rabbit. J Biomech Eng. 2007;129(6):863–872. pmid:18067390
- 15. Lieber BB, Livescu V, Hopkins L, Wakhloo AK. Particle image velocimetry assessment of stent design influence on intra-aneurysmal flow. Ann Biomed Eng. 2002;30(6):768–777. pmid:12220077
- 16. Meckel S, Stalder AF, Santini F, Radü EW, Rüfenacht DA, Markl M, et al. In vivo visualization and analysis of 3-D hemodynamics in cerebral aneurysms with flow-sensitized 4-D MR imaging at 3 T. Neuroradiology. 2008;50(6):473–484. pmid:18350286
- 17. Berg P, Stucht D, Janiga G, Beuing O, Speck O, Thévenin D. Cerebral blood flow in a healthy Circle of Willis and two intracranial aneurysms: computational fluid dynamics versus four-dimensional phase-contrast magnetic resonance imaging. J Biomech Eng. 2014;136(4):041003.
- 18. Schnell S, Wu C, Ansari SA. 4D MRI flow examinations in cerebral and extracerebral vessels. Ready for clinical routine? Curr Opin Neurol. 2016;29(4):419–428. pmid:27262148
- 19. Pereira VM, Brina O, Delattre BMA, Ouared R, Bouillot P, Erceg G, et al. Assessment of intra-aneurysmal flow modification after flow diverter stent placement with four-dimensional flow MRI: a feasibility study. J Neurointerv Surg. 2014;7(12):913–919. pmid:25280571
- 20. MacDonald M, Dolati P, Mitha A, Eesa M, Wong J, Frayne R. Hemodynamic alterations measured with phase-contrast MRI in a giant cerebral aneurysm treated with a flow-diverting stent. Radiol Case Rep. 2016;10(2):1109. pmid:27398123
- 21. Berg P, Iosif C, Ponsonnard S, Yardin C, Janiga G, Mounayer C. Endothelialization of over- and undersized flow-diverter stents at covered vessel side branches: an in vivo and in silico study. J Biomech. 2016;49(1):4–12. pmid:26607220
- 22.
Berg P, Daróczy L, Janiga G. Virtual stenting for intracranial aneurysms. In: Computing and visualization for intravascular imaging and computer-assisted stenting. Elsevier BV; 2017. p. 371–411. Available from: https://doi.org/10.1016%2Fb978-0-12-811018-8.00014-x.
- 23.
Seong J, Lieber BB, Wakhloo AK. Reduction of intra-aneurysmal kinetic energy by intralumenal flow diverting devices. In: ASME 2007 summer bioengineering conference. American Society of Mechanical Engineers; 2007. p. 195–196.
- 24. Xiang J, Damiano RJ, Lin N, Snyder KV, Siddiqui AH, Levy EI, et al. High-fidelity virtual stenting: modeling of flow diverter deployment for hemodynamic characterization of complex intracranial aneurysms. J Neurosurg. 2015;123(4):832–840. pmid:26090829
- 25. Zhang M, Anzai H, Chopard B, Ohta M. Towards the patient-specific design of flow diverters made from helix-like wires: an optimization study. Biomed Eng Online. 2016;15(S2).
- 26. Cebral JR, Mut F, Raschi M, Scrivano E, Ceratto R, Lylyk P, et al. Aneurysm rupture following treatment with flow-diverting stents: computational hemodynamics analysis of treatment. AJNR Am J Neuroradiol. 2010;
- 27. Ouared R, Larrabide I, Brina O, Bouillot P, Erceg G, Yilmaz H, et al. Computational fluid dynamics analysis of flow reduction induced by flow-diverting stents in intracranial aneurysms: a patient-unspecific hemodynamics change perspective. J Neurointerv Surg. 2016;8(12):1288–1293.
- 28. Chong W, Zhang Y, Qian Y, Lai L, Parker G, Mitchell K. Computational hemodynamics analysis of intracranial aneurysms treated with flow diverters: correlation with clinical outcomes. AJNR Am J Neuroradiol. 2013;35(1):136–142. pmid:24287091
- 29. Kulcsár Z, Augsburger L, Reymond P, Pereira VM, Hirsch S, Mallik AS, et al. Flow diversion treatment: intra-aneurismal blood flow velocity and WSS reduction are parameters to predict aneurysm thrombosis. Acta Neurochir (Wien). 2012;154(10):1827–1834.
- 30. Ujiie H, Tamano Y, Sasaki K, Hori T. Is the aspect ratio a reliable index for predicting the rupture of a saccular aneurysm? Neurosurgery. 2001;48(3):495–503. pmid:11270538
- 31. Zanaty M, Chalouhi N, Tjoumakaris SI, Gonzalez LF, Rosenwasser RH, Jabbour PM. Aneurysm geometry in predicting the risk of rupture. A review of the literature. Neurol Res. 2014;36(4):308–313. pmid:24617931
- 32. Jeong W, Rhee K. Hemodynamics of cerebral aneurysms: computational analyses of aneurysm progress and treatment. Comput Math Methods Med. 2012;2012:1–11.
- 33. Holzapfel GA. Modeling of saccular aneurysm growth in a human middle cerebral artery. J Biomech Eng. 2008;130(5):051012. pmid:19045519
- 34. Watton PN, Ventikos Y, Holzapfel GA. Modelling the growth and stabilization of cerebral aneurysms. Math Med Biol. 2009;26(2):133–164. pmid:19234094
- 35. Watton PN, Selimovic A, Raberger NB, Huang P, Holzapfel GA, Ventikos Y. Modelling evolution and the evolving mechanical environment of saccular cerebral aneurysms. Biomech Model Mechanobiol. 2010;10(1):109–132. pmid:20496095
- 36. Wang X, Li X. Biomechanical behaviour of cerebral aneurysm and its relation with the formation of intraluminal thrombus: a patient-specific modelling study. Comput Methods Biomech Biomed Engin. 2013;16(11):1127–1134. pmid:22292428
- 37. Rayz VL, Boussel L, Lawton MT, Acevedo-Bolton G, Ge L, Young WL, et al. Numerical modeling of the flow in intracranial aneurysms: prediction of regions prone to thrombus formation. Ann Biomed Eng. 2008;36(11):1793–1804. pmid:18787954
- 38. Rayz VL, Boussel L, Ge L, Leach JR, Martin AJ, Lawton MT, et al. Flow residence time and regions of intraluminal thrombus deposition in intracranial aneurysms. Ann Biomed Eng. 2010;38(10):3058–3069. pmid:20499185
- 39. Janiga G, Daróczy L, Berg P, Thévenin D, Skalej M, Beuing O. An automatic CFD-based flow diverter optimization principle for patient-specific intracranial aneurysms. J Biomech. 2015;48(14):3846–3852. pmid:26472308
- 40. Meng H, Tutino VM, Xiang J, Siddiqui A. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol. 2013;35(7):1254–1262. pmid:23598838
- 41. Zhang Y, Takao H, Murayama Y, Qian Y. Propose a wall shear stress divergence to estimate the risks of intracranial aneurysm rupture. ScientificWorldJournal. 2013;2013:1–8.
- 42. Tanoue T, Tateshima S, Villablanca JP, Vinuela F, Tanishita K. Wall shear stress distribution inside growing cerebral aneurysm. AJNR Am J Neuroradiol. 2011;32(9):1732–1737. pmid:21984256
- 43. Cebral JR, Vazquez M, Sforza DM, Houzeaux G, Tateshima S, Scrivano E, et al. Analysis of hemodynamics and wall mechanics at sites of cerebral aneurysm rupture. J Neurointerv Surg. 2014;7(7):530–536. pmid:24827066
- 44. Petersson S, Dyverfeldt P, Ebbers T. Assessment of the accuracy of MRI wall shear stress estimation using numerical simulations. J Magn Reson Imaging. 2012;36(1):128–138. pmid:22336966
- 45. Cibis M, Potters WV, Gijsen FJH, Marquering H, vanBavel E, van der Steen AFW, et al. Wall shear stress calculations based on 3D cine phase contrast MRI and computational fluid dynamics: a comparison study in healthy carotid arteries. NMR Biomed. 2014;27(7):826–834. pmid:24817676
- 46. Frolov SV, Sindeev SV, Liepsch D, Balasso A. Experimental and CFD flow studies in an intracranial aneurysm model with Newtonian and non-Newtonian fluids. Technol Health Care. 2016;24(3):317–333. pmid:26835725
- 47. Xu L, Sugawara M, Tanaka G, Ohta M, Liu H, Yamaguchi R. Effect of elasticity on wall shear stress inside cerebral aneurysm at anterior cerebral artery. Technol Health Care. 2016;24(3):349–357. pmid:26835728
- 48. Venugopal P, Valentino D, Schmitt H, Villablanca JP, Viñuela F, Duckwiler G. Sensitivity of patient-specific numerical simulation of cerebal aneurysm hemodynamics to inflow boundary conditions. J Neurosurg. 2007;106(6):1051–1060. pmid:17564178
- 49. McGah PM, Levitt MR, Barbour MC, Morton RP, Nerva JD, Mourad PD, et al. Accuracy of computational cerebral aneurysm hemodynamics using patient-specific endovascular measurements. Ann Biomed Eng. 2014;42(3):503–514. pmid:24162859