Towards the patient-specific design of flow diverters made from helix-like wires: an optimization study
- Mingzi Zhang^{1},
- Hitomi Anzai^{2},
- Bastien Chopard^{3} and
- Makoto Ohta^{4}Email author
https://doi.org/10.1186/s12938-016-0257-z
© The Author(s) 2016
Published: 28 December 2016
Abstract
Background
Flow diverter (FD) intervention is an emerging endovascular technique for treating intracranial aneurysms. High flow-diversion efficiency is desired to accelerate thrombotic occlusion inside the aneurysm; however, the risk of post-stenting stenosis in the parent artery is posed when flow-diversion efficiency is pursued by simply decreasing device porosity. For improving the prognosis of FD intervention, we develop an optimization method for the design of patient-specific FD devices that maintain high levels of porosity.
Methods
An automated structure optimization method for FDs with helix-like wires was developed by applying a combination of lattice Boltzmann fluid simulation and simulated annealing procedure. Employing intra-aneurysmal average velocity as the objective function, the proposed method tailored the wire structure of an FD to a given vascular geometry by rearranging the starting phase of the helix wires.
Results
FD optimization was applied to two idealized (S and C) vascular models and one realistic (R) model. Without altering the original device porosity of 80%, the flow-reduction rates of optimized FDs were improved by 5, 2, and 28% for the S, C, and R models, respectively. Furthermore, the aneurysmal flow patterns after optimization exhibited marked alterations. We confirmed that the disruption of bundle of inflow is of great help in blocking aneurysmal inflow. Axial displacement tests suggested that the optimal FD implanted in the R model possesses good robustness to tolerate uncertain axial positioning errors.
Conclusions
The optimization method developed in this study can be used to identify the FD wire structure with the optimal flow-diversion efficiency. For a given vascular geometry, custom-designed FD structure can maximally reduce the aneurysmal inflow with its porosity maintained at a high level, thereby lowering the risk of post-stenting stenosis. This method facilitates the study of patient-specific designs for FD devices.
Keywords
Intracranial aneurysm Flow diverter Design optimization Computational fluid dynamicsBackground
Flow diverter (FD) intervention has become increasingly attractive for the treatment of wide-neck and fusiform intracranial aneurysms (IAs), which has been studied intensively by many groups [1–4] in recent years. Meanwhile, clinical follow-ups revealed that FD recipients might incur post-stenting complications such as delayed aneurysm ruptures and post-stenting stenosis [5–7].
The conventional, commercially available FD devices such as Pipeline embolization device (PED; Irvine, CA, USA) and SILK (Balt, Montmorency, France) are constructed by homogeneous helix-like wires. The porosity of FD device is associated with post-stenting stenosis as suggested by prior studies [8], since a high metal-to-arterial tissue ratio resulted from low device porosity may pose the risk of vascular injury. Animal experiments, on the other hand, have confirmed that a lower device porosity promotes a complete thrombotic occlusion of an aneurysm [9, 10]. Therefore, simply modifying FD wires by decreasing or increasing the porosity may result in an increased risk of post-stenting stenosis or long-term thrombosis formation. To accelerate thrombotic occlusion and avoid post-stenting stenosis, a possible solution may be the application of patient-specifically tailored FDs with porosity maintained at a high level. Attempts of introducing optimization to FD structures have been made to improve the flow-diversion efficiency [11–13]. However, a practical optimization strategy that can be feasibly applied to the conventional FDs has not yet been developed.
In this study, we demonstrate an automated optimization method on a conventional, homogeneous, helix-like FD to adapt its wire structure to an assigned aneurysm. The proposed optimization was designed to rearrange the starting phases of FD wires, so that the original device porosity was kept to maintain the metal-to-arterial tissue ratio. After the FD structure with the highest flow-diversion efficiency was identified, its robust performance was then investigated by axial displacement test.
Methods
Vascular model
Aneurysmal local hemodynamics is sensitive to the morphological characteristics of the parent artery. Thus, three vascular geometries were used to investigate the proposed optimization method under various hemodynamic conditions.
Idealized aneurysm geometry
Two idealized aneurysm models—the Straight (S, Fig. 1a) model and the Curved (C, Fig. 1a) model—were constructed; for both models, the aneurysmal diameter (D) was 4.8 mm, the neck diameter (N) was 2.8 mm, and the arterial diameter (d _{1}) was 3.5 mm. The curvature radius (r) of the C model was 6.0 mm [14].
Realistic aneurysm geometry
The 3D patient-specific geometry of a human internal carotid artery (ICA) with an aneurysm was reconstructed (R model, Fig. 1a). The parent artery had an inlet diameter (d _{2}) of 3.8 mm [15].
FD model
Random modification
In this manner, the modification resulted in the axial displacement of the arbitrarily selected helix along the centerline of the parent artery, whereas the device porosity remained unchanged (Fig. 1).
Hemodynamic simulation
A previous study suggested that the peak value of hemodynamic parameters computed for pulsatile flow matches those of the corresponding steady flow [23]. Therefore, we performed steady flow simulation with standard D3Q19 lattice topology [22] in consideration of hundreds of CFD simulation steps in the following optimization procedure.
The bounce-back rule was used to impose the no-slip boundary conditions as well as to define FD wire structures. The spatial discretization \( (\Delta r) \) was set at 0.05 mm after sensitivity tests, in which we found no obvious differences in the intra-aneurysmal average velocity and flow patterns after doubling the lattice grid resolution. The numbers of fluid cells for the S, C, and R models were 3.57 × 10^{6}, 3.01 × 10^{6}, and 4.06 × 10^{6}, respectively.
The blood flow was assumed to be an incompressible Newtonian fluid. To reach the same Reynolds number (Re) of 200, velocity was defined as parabolic profiles at inlets of 0.23, 0.23, and 0.21 m/s for the S, C, and R models, respectively. A constant pressure boundary was imposed at the outlets. The density and kinematic viscosity were assumed to be constant at 1040 kg/m^{3} and 4.0 × 10^{−6} m^{2}/s, respectively. The kinetic viscosity of the lattice (ν _{ LB }) was chosen as 0.012, giving a relaxation time τ of \( (6\nu_{LB} + 1)/2\, = \,0.536 \). We assumed that the simulation had reached a convergent state when the change in the average energy of fluid cells was less than 10^{−6} kg m^{2}/s^{2}.
Simulated annealing
To control the random modifications progressing towards the optimal solution, a SA procedure (Fig. 1) was implemented to identify the FD structure with the lowest intra-aneurysmal average velocity within a certain range of temperature drop [24, 25]. We selected intra-aneurysmal average velocity as the objective function of SA because of its possible correlation to thrombotic occlusion [26].
Initial temperatures and lower temperature limits for the S, C, and R models
S model | C model | R model | |
---|---|---|---|
Initial temperature | 6.02 × 10^{−5} | 1.27 × 10^{−4} | 1.01 × 10^{−3} |
Lower temperature limit | 1.08 × 10^{−7} | 2.28 × 10^{−7} | 1.81 × 10^{−6} |
We used the scalar parallel computing system (SGI UV2000) at the Institute of Fluid Science, Tohoku University. The computational times for one stage of CFD simulation (using 256 cores) were approximately 30, 45, and 70 min for the S, C, and R models, respectively.
Axial displacement test (ADT)
The manual operation of FD delivery and deployment leads to an uncertain axial displacement during the interventional procedure. An axial displacement test (ADT) is designed to investigate the robustness of a given FD structure. For a given FD structure, we sequentially added the variable \( \theta_{\Delta } \), ranging from −π to π, to the starting phase θ of each helix to mimic the axial displacement along the centerline of the parent artery. A CFD simulation was subsequently performed to calculate the difference in velocity resulting from the displacement.
Flow reduction (FR) rate
Results
Intra-aneurysmal average velocity and the R _{ f } under non-stent, initial stent, and optimal stent placements of S, C, and R models, respectively
Non-stent | Initial stent placement | Optimal stent placement | |||
---|---|---|---|---|---|
V _{ w/o } (mm/s) | V _{ initial } (mm/s) | R _{ f } (%) | V _{ optimal } (mm/s) | R _{ f } (%) | |
S model | 2.827 | 0.269 | 90.48 | 0.123 | 95.65 |
C model | 10.629 | 0.528 | 95.04 | 0.249 | 97.65 |
R model | 39.309 | 13.071 | 66.75 | 2.025 | 94.85 |
Aneurysmal inflow is observed as a flow bundle entering from the aneurysm orifice. The concept of bundle of inflow (BOI) area describes the inflow feature of an aneurysm [13]. Likewise, the bundle of outflow (BOO) area indicates the region(s) where the bloodstream exits an aneurysm. To demonstrate the unique features of different flow patterns, the red and blue arrows/circles are used in Fig. 3b, c to indicate the BOI and BOO areas, respectively; the yellow circles with dotted lines in Fig. 3a, b, d depict the areas where BOI concentrates. The black arrows in all figures are used to identify the flow directions.
S model
Before FD stent implantation, the flow enters the aneurysm from the distal end and exits through the proximal end. After FD implantation with the initial structure, the symmetric flow distribution was disrupted; the flow circulation inside of the aneurysm sac became irregular and sparse. The BOI area shifted from the central distal end to the proximal end with a negative Y offset; meanwhile, the outflow zone switched to the distal end (Fig. 3c). A marked reduction in both the breadth and magnitude of the BOI was observed for the optimized FD structure (Fig. 3b). Rotational flow circulations were found inside the aneurysm. The optimized wire structure showed a concentration of FD wires inside the BOI area, resulting in a denser strut distribution in the proximal orifice end (Fig. 3a).
C model
Under the non-stent condition, a strong inflow jet and two outflow jets were observed (Fig. 3d). The BOI was located in the proximal area of the neck, whereas the BOOs were symmetrically dispersed along both sides of the aneurysmal neck. After the initial FD implantation, the magnitude of the velocity of the BOI decreased considerably (Fig. 3b), and the outflow bundles were disrupted. For the optimized FD, the velocity magnitude of the BOI was further reduced. A denser strut distribution inside the BOI region was observed in the optimized case (Fig. 3a).
R model
Compared with the S and C models, the non-stent R model has a strong and sharp inflow jet, as seen in Fig. 3d. The bloodstream flows into the deep aneurysm sac through the orifice and circulates inside the cavity, finally flowing out of the aneurysm as a wide and strong BOO. After the FD implantation of the initial structure, the width and velocity magnitude of BOI were reduced; both the size and volume of the isovelocity surface decreased (Fig. 3d). After the implantation of the optimized FD, the flow circulation was drastically modified. The depth traveled by the inflow jet into the aneurysm sac decreased, and rotational flow circulation was clearly observed inside the sac (Fig. 3b). The sac circulation and parent flow were further separated, and both BOI and BOO were split into two weak and thin streams (Fig. 3c). Similarly, the wire concentration inside the BOI area of the optimized FD structure can be seen in Fig. 3a.
Discussion
In this study, we have demonstrated an optimization approach to improve the flow-diversion efficiency of conventional homogeneous FD stents. Our optimization involves rearranging the starting phases of homogeneous helix wires, thereby modifying FD structure without altering its porosity, which finally enables the optimal wire configuration to maximally block the aneurysmal inflow.
Selecting the modification parameter for optimization
In consideration of the possible correlation between porosity and post-stenting stenosis, we chose the starting phase θ of each FD helix as the modification parameter.
Previous studies on FD structure optimization have introduced various modification parameters to improve flow reduction effectiveness. Anzai et al. [13] modified the configurations of the strut segments within the aneurysmal neck domain and found that the flow reduction effectiveness was improved when the stent porosity was strictly maintained at 80%, while the isolated and disconnected strut segments denied the manufactural possibility of the optimal FD structure. Lee et al. [12] used strut size and gap spacing as modification parameters and obtained relatively optimal designs. However, the porosity of FD devices cannot be precisely controlled during modification.
Theta is one of the parameters required in FD design and manufacturing process. For design optimization, using θ as modification parameter could maintain the device porosity at a pre-defined value. Our results demonstrate (Fig. 2) that modifying θ effectively improves the flow reduction rate of a conventional homogeneous FD structure.
Objective function of optimization and SA procedure
We chose intra-aneurysmal average velocity as the objective function. Anzai et al. [13, 27] used average velocity, while Srinivas et al. [11, 12] used average velocity and vorticity, and Janiga et al. [28] used wall shear stress for optimization.
It should be noted that the selection of objective function remains a controversial issue. Corbett et al. [29] reported an in vitro study that thrombosis could occur in a specific threshold of shear stress or shear rate using bovine blood. A review by Moiseyev et al. [30] revealed that the shear-induced activation of platelets is a basic element for blood coagulation. On the other hand, Janiga et al. [28] described that the flow reduction within an aneurysm is relevant to wall shear stress, and recently, Chung et al. [26] showed that the average velocity in the aneurysm may be related to a shorter period of aneurysm occlusion. Revealed by these results, the average velocity seems to be correlated with thrombotic occlusion. However, further in vivo study is still indispensable for confirming its relevance in clinical practice.
We employed SA procedure to identify a global optimum for the average velocity as shown in Fig. 2. To prevent optimization from resulting in a local optimum, modifications to FD structures with inferior flow diversion performance might also be accepted according to our pre-defined cooling schedule. Our optimization approach could accept other objective functions, as long as the initial temperature and cooling schedule are well established [24, 25].
FD design with inhomogeneous wire structure
Our optimization generated modified FD designs with inhomogeneous wire structure. Similarly, previous studies also reported the benefits of inhomogeneous and asymmetric FD designs. Rudin et al. [31] showed that a local low-porosity design can decrease flow velocity inside an in vitro model. By animal experiments, Ionita et al. [32] showed the good performance of asymmetric stents in occluding rabbit elastase aneurysms. These studies suggest the possibility of using inhomogeneous device to achieve favorite treatment outcomes. In this study, we demonstrated how a conventional homogeneous FD can be tailored to an inhomogeneous one by merely changing the values of ‘θ’, which is applicable to be modified without affecting the manufacturing possibility.
In addition, our optimization results have revealed a practical approach for the conventional homogeneous FD devices to improve its flow diversion efficiency, that is, compacting FD wires into BOI areas during deployment may achieve a marked difference in blocking the aneurysmal inflow.
Robust performance of the optimal devices
We performed ADTs to investigate the robustness of the wire structures. The flow reductions achieved by the optimized FDs were greater than those obtained using the homogeneous FDs within displacement ranges of −0.25 to 0.25, −0.5 to 0.25, and −0.75 to 0.75 mm for the S, C, and R models, respectively (Fig. 4). It is indicated that a homogeneous wire configuration can nonspecifically prevent a strong inflow jet from passing through an aneurysm orifice; however, its flow-diversion efficiency is inferior to that of the optimal wire configuration when the device was desirably deployed.
The robustness of an optimized FD is associated with the wire coverage of an aneurysm orifice. After homogeneous FD implantation, if the BOI areas are axially distributed along the orifice (e.g. R model, Fig. 3c), the robustness of an optimized FD might be superior. In contrast, if the BOI areas are radially distributed (e.g. S and C models, Fig. 3c), the optimized FD robustness could be inferior. The optimized FD wires concentrate in BOI regions, whereas only a small number of wires are assigned in the remaining areas where large holes can be found. When axial displacement occurs, the inflow jet may cross the orifice through the holes, causing considerable fluctuations of R _{ f }.
The robustness of an optimized FD may also relate to the shape of an BOI area. When BOI is strong and concentrates in a small region (e.g. C model, Fig. 3c), axial displacements may result in jet flow entering aneurysm cavity through areas with less wires. It is implied that the BOI characteristics of the FD recipient needs to be investigated before an optimized device can be applied.
Limitations
This study has several limitations. We applied steady flow and Newtonian fluid assumptions to reduce the computational cost, since the objective of this study is to develop a feasible and manufacture-oriented optimization approach for FD device. When computational cost is no longer a problem, the setting is readily changed by adopting time-dependent boundary conditions and a non-Newtonian rheology in the LB solver.
Another limitation from the viewpoint of clinical practice is that sending the optimized FD to the aneurysm location could be a challenge for interventionists, since a larger device displacement after deployment may considerably increase the aneurysmal inflow. This might be solved in the future by embedding a reliability test into the optimization loop to achieve improved stability.
It should be noted that we addressed merely the hemodynamic factors that may affect the FD performance, while the mechanical and material properties of the modified FDs have not been investigated. In future work, we plan to include these parameters as a part of objective functions for optimization to improve FD’s hemodynamic compatibility.
Conclusions
A practical optimization method for commercially available helix-wire FDs was developed in this study. By rearranging the starting phase of each helix subset, the structure of FD can be tailored to efficiently block the inflow for a patient-specific aneurysm. Using this optimization method, three optimized FD structures with unchanged device porosity were obtained corresponding to three different vascular geometries. The developed method potentially enhances the study of the patient-specific design of FD devices.
Declarations
Declarations
Authors’ contributions
MO, MZ, and HA were responsible for the design, data analysis and overall investigation. MZ, HA, and BC were responsible for the computational modeling and numerical simulation. MZ and HA were responsible for constructing simulation programs. All authors (1) have made substantial contributions to conception and design, or acquisition of data, or analysis and interpretation of data; (2) have been involved in drafting the manuscript or revising it critically for important intellectual content; and (3) have given final approval of the version to be published. Each author has participated sufficiently in the work to take public responsibility for appropriate portions of the content. All authors read and approved the final manuscript.
Acknowledgements
This work was supported by the Grant-in-Aid for Scientific Research (B), Japan Ministry of Education, Science, Sports and Culture, 2013–2015 (25282140, Makoto Ohta) and in part by the Excellent Graduate Schools program, Tohoku University, Japan.
Competing interests
The authors declare that they have no competing interests.
About this supplement
This article has been published as part of BioMedical Engineering OnLine Volume 15 Supplement 2, 2016. Computational and experimental methods for biological research: cardiovascular diseases and beyond. The full contents of the supplement are available online http://biomedical-engineering-online.biomedcentral.com/articles/supplements/volume-15-supplement-2.
Availability of data and materials
The datasets of this study are available from the corresponding author on reasonable request.
Funding
Publication charges for this article have been funded by the Grand-in-Aid for scientific Research (B), Japan Ministry of Education, Science, Sports and Culture.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Sforza DM, Putman CM, Cebral JR. Hemodynamics of cerebral aneurysms. Annu Rev Fluid Mech. 2009;41:91–107.View ArticleMATHGoogle Scholar
- Appanaboyina S, Mut F, Löhner R, Putman C, Cebral J. Simulation of intracranial aneurysm stenting: techniques and challenges. Comput Methods Appl Mech Eng. 2009;198:3567–82.View ArticleMATHGoogle Scholar
- D’Urso PI, Lanzino G, Cloft HJ, Kallmes DF. Flow diversion for intracranial aneurysms: a review. Stroke. 2011;42:2363–8.View ArticleGoogle Scholar
- Chung B, Cebral JR. CFD for evaluation and treatment planning of aneurysms: review of proposed clinical uses and their challenges. Ann Biomed Eng. 2015;43(1):122–38.View ArticleGoogle Scholar
- Lubicz B, Collignon L, Raphaeli G, Pruvo J-P, Bruneau M, DeWitte O, et al. Flow-diverter stent for the endovascular treatment of intracranial aneurysms: a prospective study in 29 patients with 34 aneurysms. Stroke. 2010;41:2247–53.View ArticleGoogle Scholar
- Berge J, Biondi A, Machi P, Brunel H, Pierot L, Gabrillargues J, et al. Flow-diverter silk stent for the treatment of intracranial aneurysms: 1-year follow-up in a multicenter study. Am J Neuroradiol. 2012;33:1150–5.View ArticleGoogle Scholar
- Cohen JE, Gomori JM, Moscovici S, Leker RR, Itshayek E. Delayed complications after flow-diverter stenting: reactive in-stent stenosis and creeping stents. J Clin Neurosci. 2014;21:1116–22.View ArticleGoogle Scholar
- Tominaga R, Harasaki H, Sutton C, Emoto H, Kambic H, Hollman J. Effects of stent design and serum cholesterol level on the restenosis rate in atherosclerotic rabbits. Am Heart J. 1993;126:1049–58.View ArticleGoogle Scholar
- Lieber BB, Stancampiano AP, Wakhloo AK. Alteration of hemodynamics in aneurysm models by stenting: influence of stent porosity. Ann Biomed Eng. 1997;25:460–9.View ArticleGoogle Scholar
- Rhee K, Han MH, Cha SH. Changes of flow characteristics by stenting in aneurysm models: influence of aneurysm geometry and stent porosity. Ann Biomed Eng. 2002;30:894–904.View ArticleGoogle Scholar
- Srinivas K, Townsend S, Lee CJ, Nakayama T, Ohta M, Obayashi S, et al. Two-dimensional optimization of a stent for an aneurysm. J Med Devices. 2010;4:21003.View ArticleGoogle Scholar
- Lee C-J, Srinivas K, Qian Y. Three-dimensional hemodynamic design optimization of stents for cerebral aneurysms. Proc Inst Mech Eng. 2014;228:213–24.View ArticleGoogle Scholar
- Anzai H, Falcone J-L, Chopard B, Hayase T, Ohta M. Optimization of strut placement in flow diverter stents for four different aneurysm configurations. J Biomech Eng. 2014;136:61006.View ArticleGoogle Scholar
- Imai Y, Sato K, Ishikawa T, Yamaguchi T. Inflow into saccular cerebral aneurysms at arterial bends. Ann Biomed Eng. 2008;36:1489–95.View ArticleGoogle Scholar
- Radaelli AG, Augsburger L, Cebral JR, Ohta M, Rüfenacht DA, Balossino R, et al. Reproducibility of hemodynamical simulations in a subject-specific stented aneurysm model—a report on the virtual intracranial stenting challenge 2007. J Biomech. 2008;41:2069–81.View ArticleGoogle Scholar
- www.vmtk.org. Accessed 8 July 2015.
- Lieber BB, Sadasivan C. Endoluminal scaffolds for vascular reconstruction and exclusion of aneurysms from the cerebral circulation. Stroke. 2010;41:S21–5.View ArticleGoogle Scholar
- www.palabos.org. Accessed 26 Oct 2015.
- Chopard B, Lagrava D, Malaspinas O, Ouared R, Latt J, Lovblad K-O, et al. A lattice Boltzmann modeling of bloodflow in cerebral aneurysms. V Eur Conf Comput Fluid Dyn ECCOMAS CFD. 2010;453:12.Google Scholar
- Harrison SE, Smith SM, Bernsdorf J, Hose DR, Lawford PV. Application and validation of the lattice Boltzmann method for modelling flow-related clotting. J Biomech. 2007;40:3023–8.View ArticleGoogle Scholar
- Kim YH, Xu X, Lee JS. The effect of stent porosity and strut shape on saccular aneurysm and its numerical analysis with lattice Boltzmann method. Ann Biomed Eng. 2010;38:2274–92.View ArticleGoogle Scholar
- Succi S. The lattice-Boltzmann equation: for fluid dynamics and beyond. Oxford university press; 2001.Google Scholar
- Seo T, Schachter LG, Barakat AI. Computational study of fluid mechanical disturbance induced by endovascular stents. Ann Biomed Eng. 2005;33:444–56.View ArticleGoogle Scholar
- Kirkpatrick S. Optimization by simulated annealing: quantitative studies. J Stat Phys. 1984;34:975–86.MathSciNetView ArticleGoogle Scholar
- Johnson DS, Aragon CR, McGeoch LA, Schevon C. Optimization by simulated annealing: an experimental evaluation; part II, graph coloring and number partitioning. Oper Res. 1991;39:378–406.View ArticleMATHGoogle Scholar
- Chung B, Mut F, Kadirvel R, Lingineni R, Kallmes DF, Cebral JR. Hemodynamic analysis of fast and slow aneurysm occlusions by flow diversion in rabbits. J Neurointerv Surg. 2018;7(12):931–5.View ArticleGoogle Scholar
- Anzai H, Ohta M, Falcone J-L, Chopard B. Optimization of flow diverters for cerebral aneurysms. J Comput Sci. 2012;3:1–7.View ArticleGoogle Scholar
- 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–52.View ArticleGoogle Scholar
- Corbett SC, Ajdari A, Coskun AU, N-Hashemi H. In vitro and computational thrombosis on artificial surfaces with shear stress. Artif Organs. 2010;34:561–9.View ArticleGoogle Scholar
- Moiseyev G, Bar-Yoseph PZ. Computational modeling of thrombosis as a tool in the design and optimization of vascular implants. J Biomech. 2013;46:248–52.View ArticleGoogle Scholar
- Rudin S, Wang Z, Kyprianou I, Hoffmann KR, Wu Y, Meng H, et al. Measurement of flow modification in phantom aneurysm model: comparison of coils and a longitudinally and axially asymmetric stent—initial findings 1. Radiology. 2004;231:272–6.View ArticleGoogle Scholar
- Ionita CN, Paciorek AM, Dohatcu A, Hoffmann KR, Bednarek DR, Kolega J, et al. The asymmetric vascular stent: efficacy in a rabbit aneurysm model. Stroke. 2009;40:959–65.View ArticleGoogle Scholar