Virtual stenting workflow with vessel-specific initialization and adaptive expansion for neurovascular stents and flow diverters

Endovascular intervention using traditional neurovascular stents and densely braided flow diverters (FDs) have become the preferred treatment strategies for traditionally challenging intracranial aneurysms. Modeling stent and FD deployment in patient-specific aneurysms and its flow modification results prior to the actual intervention can potentially predict the patient outcome and treatment optimization. We present a clinically focused, streamlined virtual stenting workflow that efficiently simulates stent and FD treatment in patient-specific aneurysms based on expanding a simplex mesh structure. The simplex mesh is generated using an innovative vessel-specific initialization technique, which uses the patient’s parent artery diameter to identify the initial position of the simplex mesh inside the artery. A novel adaptive expansion algorithm enables the acceleration of deployment process by adjusting the expansion forces based on the distance of the simplex mesh from the parent vessel. The virtual stenting workflow was tested by modeling the treatment of two patient-specific aneurysms using the Enterprise stent and the Pipeline Embolization Device (commercial FD). Both devices were deployed in the aneurysm models in a few seconds. Computational fluid dynamics analyses of pre- and post-treatment aneurysmal hemodynamics show flow reduction in the aneurysmal sac in treated aneurysms, with the FD diverting more flow than the Enterprise stent. The test results show that this workflow can rapidly simulate clinical deployment of stents and FDs, hence paving the way for its future clinical implementation.     

Impact of stents and flow diverters on hemodynamics in idealized aneurysm models

Cerebral aneurysms constitute a major medical challenge as treatment options are limited and often associated with high risks. Statistically, up to 3% of patients with a brain aneurysm may suffer from bleeding for each year of life. Eight percent of all strokes are caused by ruptured aneurysms. In order to prevent this rupture, endovascular stenting using so called flow diverters is increasingly being regarded as an alternative to the established coil occlusion method in minimally invasive treatment. Covering the neck of an aneurysm with a flow diverter has the potential to alter the hemodynamics in such a way as to induce thrombosis within the aneurysm sac, stopping its further growth, preventing its rupture and possibly leading to complete resorption. In the present study the influence of different flow diverters is quantified considering idealized patient configurations, with a spherical sidewall aneurysm placed on either a straight or a curved parent vessel. All important hemodynamic parameters (exchange flow rate, velocity, and wall shear stress) are determined in a quantitative and accurate manner using computational fluid dynamics when varying the key geometrical properties of the aneurysm. All simulations are carried out using an incompressible, Newtonian fluid with steady conditions. As a whole, 72 different cases have been considered in this systematic study. In this manner, it becomes possible to compare the efficiency of different stents and flow diverters as a function of wire density and thickness. The results show that the intra-aneurysmal flow velocity, wall shear stress, mean velocity, and vortex topology can be considerably modified thanks to insertion of a suitable implant. Intra-aneurysmal residence time is found to increase rapidly with decreasing stent porosity. Of the three different implants considered in this study, the one with the highest wire density shows the highest increase of intra-aneurysmal residence time for both the straight and the curved parent vessels. The best hemodynamic modifications are always obtained for a small aneurysm diameter.     

Statistical hemodynamics: a tool for evaluating the effect of fluid dynamic forces on vascular biology in vivo

BACKGROUND: In vivo experimentation is the most realistic approach for exploring the vascular biological response to the hemodynamic stresses that are present in life. Post-mortem vascular casting has been used to define the in vivo geometry for hemodynamic simulation; however, this procedure damages or destroys the tissue and cells on which biological assays are to be performed. METHOD OF APPROACH: Two statistical approaches, regional (RSH) and linear (LSH) statistical hemodynamics, are proposed and illustrated, in which flow simulations from one series of experiments are used to define a best estimate of the hemodynamic environment in a second series. As an illustration of the technique, RSH is used to compare the gene expression profiles of regions of the proximal external iliac arteries of swine exposed to different levels of time-average shear stress. RESULTS: The results indicate that higher shears promote a more atheroprotective expression phenotype in porcine arterial endothelium. CONCLUSION: Statistical hemodynamics provides a realistic estimate of the hemodynamic stress on vascular tissue that can be correlated against biological response.     

Finite element modeling of embolic coil deployment: Multifactor characterization of treatment effects on cerebral aneurysm hemodynamics

Endovascular coiling is the most common treatment for cerebral aneurysms. During the treatment, a sequence of embolic coils with different stiffness, shapes, sizes, and lengths is deployed to fill the aneurysmal sac. Although coil packing density has been clinically correlated with treatment success, many studies have also reported success at low packing densities, as well as recurrence at high packing densities. Such reports indicate that other factors may influence treatment success. In this study, we used a novel finite element approach and computational fluid dynamics (CFD) to investigate the effects of packing density, coil shape, aneurysmal neck size, and parent vessel flow rate on aneurysmal hemodynamics. The study examines a testbed of 80 unique CFD simulations of post-treatment flows in idealized basilar tip aneurysm models. Simulated coil deployments were validated against in vitro and in vivo deployments. Among the investigated factors, packing density had the largest effect on intra-aneurysmal velocities. However, multifactor analysis of variance showed that coil shape can also have considerable effects, depending on packing density and neck size. Further, linear regression analysis showed an inverse relationship between mean void diameter in the aneurysm and mean intra-aneurysmal velocities, which underscores the importance of coil distribution and thus coil shape. Our study suggests that while packing density plays a key role in determining post-treatment hemodynamics, other factors such as coil shape, aneurysmal geometry, and parent vessel flow may also be very important.     

Combined Effects of Flow Diverting Strategies and Parent Artery Curvature on Aneurysmal Hemodynamics: A CFD Study

Purpose

Flow diverters (FD) are increasingly being considered for treating large or giant wide-neck aneurysms. Clinical outcome is highly variable and depends on the type of aneurysm, the flow diverting device and treatment strategies. The objective of this study was to analyze the effect of different flow diverting strategies together with parent artery curvature variations on altering intra-aneurysmal hemodynamics.

Methods

Four ideal intracranial aneurysm models with different parent artery curvature were constructed. Computational fluid dynamics (CFD) simulations of the hemodynamics before and after applying five types of flow diverting strategies (single FD, single FD with 5% and 10% packing density of coils, two FDs with 25% and 50% overlapping rate) were performed. Changes in pressure, wall shear stress (WSS), relative residence time (RRT), inflow velocity and inflow volume rate were calculated and compared.

Results

Each flow diverting strategy resulted in enhancement of RRT and reduction of normalized mean WSS, inflow volume rate and inflow velocity in various levels. Among them, 50% overlapped FD induced most effective hemodynamic changes in RRT and inflow volume rate. The mean pressure only slightly decreased after treatment. Regardless of the kind of implantation of FD, the mean pressure, inflow volume rate and inflow velocity increased and the RRT decreased as the curvature of the parent artery increased.

Conclusions

Of all flow diverting strategies, overlapping FDs induced most favorable hemodynamic changes. Hemodynamics alterations post treatment were substantially influenced by parent artery curvature. Our results indicate the need of an individualized flow diverting strategy that is tailored for a specific aneurysm.     

Mapping vascular response to in vivo Hemodynamics: application to increased flow at the basilar terminus

Hemodynamic forces play critical roles in vascular pathologies such as atherosclerosis, aneurysms, and stenosis. However, detailed relationships between the specific in vivo hemodynamic microenvironment and vascular responses leading to the triggering or exacerbation of pathological remodeling of the vessel remain elusive. We have developed a hemodynamics-biology co-mapping technique that enables in situ correlation between the in vivo blood flow field and vascular changes secondary to hemodynamic insult. The hemodynamics profile is obtained from computational fluid dynamics simulation within the vascular geometry reconstructed from three-dimensional in vivo images, whereas the vascular response is obtained from histology or immunohistochemistry on harvested vascular tissue. The hemodynamics field is virtually sectioned in the histological slicing planes and digitally co-mapped with the histological images, thereby enabling correlation of the specific local vascular responses with the inciting hemodynamic stresses. We demonstrate application of this technique to rabbit basilar terminus subjected to elevated flow. Morphological changes at the basilar terminus 5 days after the flow increase were co-mapped with the initial wall shear stress and wall shear stress gradient distributions, from which localization of destructive remodeling in a specific hemodynamic zone was noticed. This method paves the way for further investigations to determine the connection between in vivo mechanical stimuli and biological responses, such as initiation of aneurysmal remodeling.     

Stratification of a population of intracranial aneurysms using blood flow metrics

Indices of the intra-aneurysm hemodynamic environment have been proposed as potentially indicative of their longitudinal outcome. To be useful, the indices need to be used to stratify large study populations and tested against known outcomes. The first objective was to compile the diverse hemodynamic indices reported in the literature. Furthermore, as morphology is often the only patient-specific information available in large population studies, the second objective was to assess how the ranking of aneurysms in a population is affected by the use of steady flow simulation as an approximation to pulsatile flow simulation, even though the former is clearly non-physiological. Sixteen indices of aneurysmal hemodynamics reported in the literature were compiled and refined where needed. It was noted that, in the literature, these global indices of flow were always time-averaged over the cardiac cycle. Steady and pulsatile flow simulations were performed on a population of 198 patient-specific and 30 idealised aneurysm models. All proposed hemodynamic indices were estimated and compared between the two simulations. It was found that steady and pulsatile flow simulations had a strong linear dependence (r ≥ 0.99 for 14 indices; r ≥ 0.97 for 2 others) and rank the aneurysms in an almost identical fashion (ρ ≥ 0.99 for 14 indices; ρ ≥ 0.96 for other 2). When geometry is the only measured piece of information available, stratification of aneurysms based on hemodynamic indices reduces to being a physically grounded substitute for stratification of aneurysms based on morphology. Under such circumstances, steady flow simulations may be just as effective as pulsatile flow simulation for estimating most key indices currently reported in the literature.     

False discovery rate control in two-stage designs

ABSTRACT: BACKGROUND: For gene expression or gene association studies with a large number of hypotheses the number of measurements per marker in a conventional single-stage design is often low due to limited resources. Two-stage designs have been proposed where in a first stage promising hypotheses are identified and further investigated in the second stage with larger sample sizes. For two types of two-stage designs proposed in the literature we derive multiple testing procedures controlling the False Discovery Rate (FDR) demonstrating FDR control by simulations: designs where a fixed number of top-ranked hypotheses are selected and designs where the selection in the interim analysis is based on an FDR threshold. In contrast to earlier approaches which use only the second-stage data in the hypothesis tests (pilot approach), the proposed testing procedures are based on the pooled data from both stages (integrated approach). Results: For both selection rules the multiple testing procedures control the FDR in the considered simulation scenarios. This holds for the case of independent observations across hypotheses as well as for certain correlation structures. Additionally, we show that in scenarios with small effect sizes the testing procedures based on the pooled data from both stages can give a considerable improvement in power compared to tests based on the second-stage data only. Conclusion: The proposed hypothesis tests provide a tool for FDR control for the considered two-stage designs. Comparing the integrated approaches for both selection rules with the corresponding pilot approaches showed an advantage of the integrated approach in many simulation scenarios.     

In silico simulation of biological network dynamics

Realistic simulation of biological networks requires stochastic simulation approaches because of the small numbers of molecules per cell. The high computational cost of stochastic simulation on conventional microprocessor-based computers arises from the intrinsic disparity between the sequential steps executed by a microprocessor program and the highly parallel nature of information flow within biochemical networks. This disparity is reduced with the Field Programmable Gate Array (FPGA)-based approach presented here. The parallel architecture of FPGAs, which can simulate the basic reaction steps of biological networks, attains simulation rates at least an order of magnitude greater than currently available microprocessors.     

Influence of stent configuration on cerebral aneurysm fluid dynamics

Embolic coiling is the most popular endovascular treatment available for cerebral aneurysms. Nevertheless, the embolic coiling of wide-neck aneurysms is challenging and, in many cases, ineffective. Use of highly porous stents to support coiling of wide-neck aneurysms has become a common procedure in recent years. Several studies have also demonstrated that high porosity stents alone can significantly alter aneurysmal hemodynamics, but differences among different stent configurations have not been fully characterized. As a result, it is usually unclear which stent configuration is optimal for treatment. In this paper, we present a flow study that elucidates the influence of stent configuration on cerebral aneurysm fluid dynamics in an idealized wide-neck basilar tip aneurysm model. Aneurysmal fluid dynamics for three different stent configurations (half-Y, Y and, cross-bar) were first quantified using particle image velocimetry and then compared. Computational fluid dynamics (CFD) simulations were also conducted for selected stent configurations to facilitate validation and provide more detailed characterizations of the fluid dynamics promoted by different stent configurations. In vitro results showed that the Y stent configuration reduced cross-neck flow most significantly, while the cross-bar configuration reduced velocity magnitudes within the aneurysmal sac most significantly. The half-Y configuration led to increased velocity magnitudes within the aneurysmal sac at high parent-vessel flow rates. Experimental results were in strong agreement with CFD simulations. Simulated results indicated that differences in fluid dynamic performance among the different stent configurations can be attributed primarily to protruding struts within the bifurcation region.     

Hemodynamics of the mouse abdominal aortic aneurysm

The abdominal aortic aneurysm (AAA) is a significant cause of death and disability in the Western world and is the subject of many clinical and pathological studies. One of the most commonly used surrogates of the human AAA is the angiotensin II (Ang II) induced model used in mice. Despite the widespread use of this model, there is a lack of knowledge concerning its hemodynamics; this study was motivated by the desire to understand the fluid dynamic environment of the mouse AAA. Numerical simulations were performed using three subject-specific mouse models in flow conditions typical of the mouse. The numerical results from one model showed a shed vortex that correlated with measurements observed in vivo by Doppler ultrasound. The other models had smaller aneurysmal volumes and did not show vortex shedding, although a recirculation zone was formed in the aneurysm, in which a vortex could be observed, that elongated and remained attached to the wall throughout the systolic portion of the cardiac cycle. To link the hemodynamics with aneurysm progression, the remodeling that occurred between week one and week two of the Ang II infusion was quantified and compared with the hemodynamic wall parameters. The strongest correlation was found between the remodeled distance and the oscillatory shear index, which had a correlation coefficient greater than 0.7 for all three models. These results demonstrate that the hemodynamics of the mouse AAA are driven by a strong shear layer, which causes the formation of a recirculation zone in the aneurysm cavity during the systolic portion of the cardiac waveform. The recirculation zone results in areas of quiescent flow, which are correlated with the locations of the aneurysm remodeling.     

Marker-based parentage analysis facilitates low input ‘breeding without breeding’ strategies for forest trees

Controlled crosses and progeny testing are important components of tree breeding programs. Recently, a new approach, ‘breeding without breeding’ (BWB), has been proposed to obtain structured pedigrees for the breeding and testing from wind-pollinated progenies; the technique uses high-resolution markers instead of conventional approaches. The BWB approach is attractive for outcrossing conifers where the breeding cycle is long and reproductive maturity is often delayed. For the application of the BWB approach, the accurate assignment of parents is essential. The first step in the application of the BWB approach is to conduct tests into the reliability of paternity and parent-pair exclusion procedures in complex situations. In this study we conducted multi-locus-based empirical and simulation analyses for a Scots pine seed orchard crop, in order to develop the appropriate parentage assignment techniques for BWB and to investigate whether pedigrees from a seed orchard crop can be determined with sufficient accuracy. The results were promising, as 93–97% of the seeds generated by mating among the seed orchard clones were correctly assigned to a parent-pair. When mating occurred with foreign pollen, the success rate of identifying mother clones was 60–78%. Based on these results, we propose three novel low input breeding strategies for forest trees. The principles of the different options and their potential benefits and costs are described and discussed using Scots pine as a case study.     

Numerical Simulation of Physiological Blood Flow in 2-way Coronary Artery Bypass Grafts

The Coronary Artery Bypass Graft (CABG) yields excellent results and remains the modern standard of care for treatment of occlusive disease in the cardiovascular system. However, the development of anastomotic Intimal Hyperplasia (IH) and restenosis can compromise the medium-and-long term effects of the CABG. This problem can be correlated with the geometric configuration and hemodynamics of the bypass graft. A novel geometric configuration was proposed for the CABG with two symmetrically implanted grafts for the purpose of improving the hemodynamics. Physiological blood flows in two models of bypass grafts were simulated using numerical methods. One model was for the conventional bypass configuration with a single graft (1-way model); the other model was for the proposed bypass configuration with two grafts (2-way model). The temporal and spatial distributions of hemodynamics, such as flow patterns and Wall Shear Stress (WSS) in the vicinity of the distal anastomoses, were analyzed and compared. Calculation results showed that the 2-way model possessed favorable hemodynamics with uniform longitudinal flow patterns and WSS distributions, which could decrease the probability of restenosis and improve the effect of the surgical treatment. Concerning the limitations of the 2-way bypass grafts, it is necessary to perform animal experiments to verify the viability of this novel idea for the CABG.     

Reliable prediction of adsorption isotherms via genetic algorithm molecular simulation

Conventional molecular simulation techniques such as grand canonical Monte Carlo (GCMC) strictly rely on purely random search inside the simulation box for predicting the adsorption isotherms. This blind search is usually extremely time demanding for providing a faithful approximation of the real isotherm and in some cases may lead to non-optimal solutions. A novel approach is presented in this article which does not use any of the classical steps of the standard GCMC method, such as displacement, insertation, and removal. The new approach is based on the well-known genetic algorithm to find the optimal configuration for adsorption of any adsorbate on a structured adsorbent under prevailing pressure and temperature. The proposed approach considers the molecular simulation problem as a global optimization challenge. A detailed flow chart of our so-called genetic algorithm molecular simulation (GAMS) method is presented, which is entirely different from traditions molecular simulation approaches. Three real case studies (for adsorption of CO₂ and H₂ over various zeolites) are borrowed from literature to clearly illustrate the superior performances of the proposed method over the standard GCMC technique. For the present method, the average absolute values of percentage errors are around 11% (RHO-H₂), 5% (CHA-CO₂), and 16% (BEA-CO₂), while they were about 70%, 15%, and 40% for the standard GCMC technique, respectively.     

PIV-validated numerical modeling of pulsatile flows in distal coronary end-to-side anastomoses

This study employed particle image velocimetry (PIV) to validate a numerical model in a complementary approach to quantify hemodynamic factors in distal coronary anastomoses and to gain more insights on their relationship with anastomotic geometry. Instantaneous flow fields and wall shear stresses (WSS) were obtained from PIV measurement in a modified life-size silastic anastomosis model adapted from a conventional geometry by incorporating a smooth graft-artery transition. The results were compared with those predicted by a concurrent numerical model. The numerical method was then used to calculate cycle-averaged WSS (WSS(cyc)) and spatial wall shear stress gradient (SWSSG), two critical hemodynamic factors in the pathogenesis of intimal thickening (IT), to compare the conventional and modified geometries. Excellent qualitative agreement and satisfactory quantitative agreement with averaged normalized error in WSS between 0.8% and 8.9% were achieved between the PIV experiment and numerical model. Compared to the conventional geometry, the modified geometry produces a more uniform WSS(cyc) distribution eliminating both high and low WSS(cyc) around the toe, critical in avoiding IT. Peak SWSSG on the artery floor of the modified model is less than one-half that in the conventional case, and high SWSSG at the toe is eliminated. The validated numerical model is useful for modeling unsteady coronary anastomotic flows and elucidating the significance of geometry regulated hemodynamics. The results suggest the clinical relevance of constructing smooth graft-artery transition in distal coronary anastomoses to improve their hemodynamic performance.     

Realistic non-Newtonian viscosity modelling highlights hemodynamic differences between intracranial aneurysms with and without surface blebs

Most computational fluid dynamic (CFD) simulations of aneurysm hemodynamics assume constant (Newtonian) viscosity, even though blood demonstrates shear-thinning (non-Newtonian) behavior. We sought to evaluate the effect of this simplifying assumption on hemodynamic forces within cerebral aneurysms, especially in regions of low wall shear stress, which are associated with rupture. CFD analysis was performed for both viscosity models using 3D rotational angiography volumes obtained for 26 sidewall aneurysms (12 with blebs, 12 ruptured), and parametric models incorporating blebs at different locations (inflow/outflow zone). Mean and lowest 5% values of time averaged wall shear stress (TAWSS) computed over the dome were compared using Wilcoxon rank-sum test. Newtonian modeling not only resulted in higher aneurysmal TAWSS, specifically in areas of low flow and blebs, but also showed no difference between aneurysms with or without blebs. In contrast, for non-Newtonian analysis, bleb-bearing aneurysms showed significantly lower 5% TAWSS compared to those without (p=0.005), despite no significant difference in mean dome TAWSS (p=0.32). Non-Newtonian modeling also accentuated the differences in dome TAWSS between ruptured and unruptured aneurysms (p<0.001). Parametric models further confirmed that realistic non-Newtonian viscosity resulted in lower bleb TAWSS and higher focal viscosity, especially when located in the outflow zone. The results show that adopting shear-thinning non-Newtonian blood viscosity in CFD simulations of intracranial aneurysms uncovered hemodynamic differences induced by bleb presence on aneurysmal surfaces, and significantly improved discriminant statistics used in risk stratification. These findings underline the possible implications of using a realistic model of blood viscosity in predictive computational hemodynamics.     

Blood flow patterns and pressure loss in the ascending aorta: A comparative study on physiological and aneurysmal conditions

An aortic aneurysm is defined as a balloon-shaped bulging of all three histologic components of the aortic vessel walls (intima, media and adventitia). This dilation results from vessel weakening owing to remodeling, i.e. due to cystic degeneration of the Tunica media (Marfan), progression of atherosclerosis or presence of a bicuspid aortic valve. The growth rate of the aortic diameter varies from patient to patient and may progress until the aneurysm ultimately ruptures. The role of hemodynamics, i.e. blood flow patterns, and shear stresses that are supposed to intensify during aneurysm growth are not yet fully understood, but thought to play a key role in the enlargement process. The aim of this study is to characterize the aortic blood flow in a silicone model of a pathological aorta with ascending aneurysm, to analyze the differences in the blood flow pattern compared to a healthy aortic model, and to single out possible blood flow characteristics measurable using phase contrast magnetic resonance imaging (MRI) that could serve as indicators for aneurysm severity. MRI simulations were performed under physiological, pulsatile flow conditions using data obtained from optical three dimensional particle tracking measurements. In comparison to the healthy geometry, elevated turbulence intensity and pressure loss are measured in the diseased aorta, which we propose as a complimentary indicator for assessing the aneurysmal severity. Our results shed a light on the interplay between the blood flow dynamics and their contribution to the pathophysiology and possible role for future risk assessment of ascending aortic aneurysms.     

Influence of head position on carotid hemodynamics in young adults

Studies in adults have shown marked changes in geometry and relative positions of the carotid arteries when rotating the head. The aim of this study was to quantify the change in geometry and analyze its effect on carotid hemodynamics as a result of head rotation. The right carotid arteries of nine young adult subjects were investigated in supine position with straight and left turned head positions, respectively. The three-dimensional (3D) carotid geometry was reconstructed by using 3D ultrasound (3D US), and the carotid hemodynamics were calculated by combining 3D US with computational fluid dynamics. It was observed that cross-sectional areas and shapes did not change markedly with head rotation, but carotid vessel center lines altered with planarification of the common carotid artery as a main feature (P < 0.05). Measured common carotid flow rates changed significantly at the individual level when the head was turned, but on the average, the change in mean common carotid flow rate was relatively small (0.37 +/- 1.11 ml/s). The effect of the altered center lines and flow rates on the atherogenic nature of the carotid bifurcation was evaluated by using calculated hemodynamic wall parameters, such as wall shear stress (WSS) and oscillatory shear index (OSI). It was found that WSS and OSI patterns changed significantly with head rotation, but the variations were very subject dependent and could not have been predicted without assessing the altered geometry and flow of the carotid bifurcation for individual cases. This study suggests that there is a need for standardization of the choice of head position in the 3D US scan protocol, and that carotid stents and emboli diverters should be studied in different head positions.     

Anomalous hemodynamic effects of a self-expanding intracranial stent: Comparing in-vitro and ex-vivo models using ultra-high resolution MicroCT based CFD

Previous research on the effects of intracranial stents on arterial hemodynamics has involved computational hemodynamics (CHD) simulations applied to artificially generated stent models. In this study, accurate geometric reconstructions of in-vitro (PTFE tube) and ex-vivo (canine artery) deployed stents based on ultra-high resolution MicroCT imaging were used. The primary goal was to compare the hemodynamic effects of deployment in these two different models and to identify flow perturbations due to deployment anomalies such as stent malapposition and strut prolapse, important adverse mechanics occurring in clinical practice, but not considered in studies using idealized stent models.Ultra-high resolution MicroCT data provided detailed visualization of deployment characteristics allowing for accurate in-stent flow simulation. For stent cells that are regularly and symmetrically deployed, the near wall flow velocities and wall shear stresses were similar to previously published results derived from idealized models. In-stent hemodynamics were significantly altered by misaligned or malapposed stent cells, important effects not realistically captured in previous models. This research shows the feasibility and value of an ex-vivo stent model for MicroCT based CHD studies. It validates previous in-vitro studies and further contributes to the understanding of in-stent hemodynamics associated with adverse mechanics of self-expanding intracranial stents.