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.     

A multi-scale mechanobiological model of in-stent restenosis: deciphering the role of matrix metalloproteinase and extracellular matrix changes

Since their first introduction, stents have revolutionised the treatment of atherosclerosis; however, the development of in-stent restenosis still remains the Achilles' heel of stent deployment procedures. Computational modelling can be used as a means to model the biological response of arteries to different stent designs using mechanobiological models, whereby the mechanical environment may be used to dictate the growth and remodelling of vascular cells. Changes occurring within the arterial wall due to stent-induced mechanical injury, specifically changes within the extracellular matrix, have been postulated to be a major cause of activation of vascular smooth muscle cells and the subsequent development of in-stent restenosis. In this study, a mechanistic multi-scale mechanobiological model of in-stent restenosis using finite element models and agent-based modelling is presented, which allows quantitative evaluation of the collagen matrix turnover following stent-induced arterial injury and the subsequent development of in-stent restenosis. The model is specifically used to study the influence of stent deployment diameter and stent strut thickness on the level of in-stent restenosis. The model demonstrates that there exists a direct correlation between the stent deployment diameter and the level of in-stent restenosis. In addition, investigating the influence of stent strut thickness using the mechanobiological model reveals that thicker strut stents induce a higher level of in-stent restenosis due to a higher extent of arterial injury. The presented mechanobiological modelling framework provides a robust platform for testing hypotheses on the mechanisms underlying the development of in-stent restenosis and lends itself for use as a tool for optimisation of the mechanical parameters involved in stent design.     

Hemodynamic study of overlapping bare-metal stents intervention to aortic aneurysm

To investigate the hemodynamic performance of overlapping bare-metal stents intervention treatment to thoracic aortic aneurysms (TAA), three simplified TAA models, representing, no stent, with a single stent and 2 overlapped stents deployed in the aneurismal sac, were studied and compared in terms of flow velocity, wall shear stress (WSS) and pressure distributions by means of computational fluid dynamics. The results showed that overlapping stents intervention induced a flow field of slow velocity near the aneurismal wall. Single stent deployment in the sac reduced the jet-like flow formed prior to the proximal neck of the aneurysm, which impinged on the internal wall of the aneurysm. This jet-like flow vanished completely in the overlapping double stents case. Overlapping stents intervention led to an evident decrease in WSS; meanwhile, the pressure acting on the wall of the aneurysm was reduced slightly and presented more uniform distribution. The results therefore indicated that overlapping stents intervention may effectively isolate the thoracic aortic aneurysm, protecting it from rupture. In conclusion, overlapping bare-metal stents may serve a purpose similar to that of the multilayer aneurysm repair system (MARS) manufactured by Cardiatis SA (Isnes, Belgium).     

Influence of overlapping pattern of multiple overlapping uncovered stents on the local mechanical environment: A patient-specific parameter study

BackgroundMultiple overlapping uncovered stents (MOUS) system has shown potentials in managing complex aortic aneurysms with side branches involvement. It promotes the development of thrombus by modulating local flow pattern that reduces the wall tension, while maintaining patency of side branches. However the modulation of local hemodynamic parameters depends on various factors that have not been assessed comprehensively.MethodsAneurysm 3D geometry was reconstructed based on CT images. One-way fluid-structure interaction analysis was performed to quantify structural stress concentration in the wall, and changes of blood velocity, wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT) and pressure in the sac due to the stent deployment.ResultsHigh structural stress concentration due to stent deployment was found in the landing zone and it increased linearly with the number of stents deployed. The wall tension in the sac was unaffected by the stent deployment. Stress within the wall was insensitive to the different overlapping pattern. After one stent was deployed, the mean flow velocity in the sac reduced by 36.4%. The deployment of the 2nd stent further reduced the mean sac velocity by 10%. WSS decreased while both OSI and RRT increased after stent deployment, however pressure in the sac remained nearly unchanged. Except for the cases with complete stents struts alignment, different overlapping pattern had little effect on flow parameters.ConclusionsMechanical parameters modulated by the MOUS are insensitive to different overlapping pattern suggesting that endovascular procedure can be performed with less attention to the overlapping pattern.     

Sequential structural and fluid dynamic numerical simulations of a stented bifurcated coronary artery

Despite their success, stenting procedures are still associated to some clinical problems like sub-acute thrombosis and in-stent restenosis. Several clinical studies associate these phenomena to a combination of both structural and hemodynamic alterations caused by stent implantation. Recently, numerical models have been widely used in the literature to investigate stenting procedures but always from either a purely structural or fluid dynamic point of view. The aim of this work is the implementation of sequential structural and fluid dynamic numerical models to provide a better understanding of stenting procedures in coronary bifurcations. In particular, the realistic geometrical configurations obtained with structural simulations were used to create the fluid domains employed within transient fluid dynamic analyses. This sequential approach was applied to investigate the final kissing balloon (FKB) inflation during the provisional side branch technique. Mechanical stresses in the arterial wall and the stent as well as wall shear stresses along the arterial wall were examined before and after the FKB deployment. FKB provoked average mechanical stresses in the arterial wall almost 2.5 times higher with respect to those induced by inflation of the stent in the main branch only. Results also enlightened FKB benefits in terms of improved local blood flow pattern for the side branch access. As a drawback, the FKB generates a larger region of low wall shear stress. In particular, after FKB the percentage of area characterized by wall shear stresses lower than 0.5?Pa was 79.0%, while before the FKB it was 62.3%. For these reasons, a new tapered balloon dedicated to bifurcations was proposed. The inclusion of the modified balloon has reduced the mechanical stresses in the proximal arterial vessel to 40% and the low wall shear stress coverage area to 71.3%. In conclusion, these results show the relevance of the adopted sequential approach to study the wall mechanics and the hemodynamics created by stent deployment.     

Quantification of thrombus formation in malapposed coronary stents deployed in vitro through imaging analysis

Stent thrombosis is a major complication of coronary stent and scaffold intervention. While often unanticipated and lethal, its incidence is low making mechanistic examination difficult through clinical investigation alone. Thus, throughout the technological advancement of these devices, experimental models have been indispensable in furthering our understanding of device safety and efficacy. As we refine model systems to gain deeper insight into adverse events, it is equally important that we continue to refine our measurement methods. We used digital signal processing in an established flow loop model to investigate local flow effects due to geometric stent features and ultimately its relationship to thrombus formation. A new metric of clot distribution on each microCT slice termed normalized clot ratio was defined to quantify this distribution. Three under expanded coronary bare-metal stents were run in a flow loop model to induce clotting. Samples were then scanned in a MicroCT machine and digital signal processing methods applied to analyze geometric stent conformation and spatial clot formation. Results indicated that geometric stent features play a significant role in clotting patterns, specifically at a frequency of 0.6225 Hz corresponding to a geometric distance of 1.606 mm. The magnitude-squared coherence between geometric features and clot distribution was greater than 0.4 in all samples. In stents with poor wall apposition, ranging from 0.27 mm to 0.64 mm maximum malapposition (model of real-world heterogeneity), clots were found to have formed in between stent struts rather than directly adjacent to struts. This early work shows how the combination of tools in the areas of image processing and signal analysis can advance the resolution at which we are able to define thrombotic mechanisms in in vitro models, and ultimately, gain further insight into clinical performance.     

Mechanistic evaluation of long-term in-stent restenosis based on models of tissue damage and growth

Development and application of advanced mechanical models of soft tissues and their growth represent one of the main directions in modern mechanics of solids. Such models are increasingly used to deal with complex biomedical problems. Prediction of in-stent restenosis for patients treated with coronary stents remains a highly challenging task. Using a finite element method, this paper presents a mechanistic approach to evaluate the development of in-stent restenosis in an artery following stent implantation. Hyperelastic models with damage, verified with experimental results, are used to describe the level of tissue damage in arterial layers and plaque caused by such intervention. A tissue-growth model, associated with vessel damage, is adopted to describe the growth behaviour of a media layer after stent implantation. Narrowing of lumen diameter with time is used to quantify the development of in-stent restenosis in the vessel after stenting. It is demonstrated that stent designs and materials strongly affect the stenting-induced damage in the media layer and the subsequent development of in-stent restenosis. The larger the artery expansion achieved during balloon inflation, the higher the damage introduced to the media layer, leading to an increased level of in-stent restenosis. In addition, the development of in-stent restenosis is directly correlated with the artery expansion during the stent deployment. The correlation is further used to predict the effect of a complex clinical procedure, such as stent overlapping, on the level of in-stent restenosis developed after percutaneous coronary intervention.

     

Effects of different stent designs on local hemodynamics in stented arteries

Following the deployment of a coronary stent and disruption of an atheromatous plaque, the deformation of the arterial wall and the presence of the stent struts create a new fluid dynamic field, which can cause an abnormal biological response. In this study 3D computational models were used to analyze the fluid dynamic disturbances induced by the placement of a stent inside a coronary artery. Stents models were first expanded against a simplified arterial plaque, with a solid mechanics analysis, and then subjected to a fluid flow simulation under pulsatile physiological conditions. Spatial and temporal distribution of arterial wall shear stress (WSS) was investigated after the expansion of stents of different designs and different strut thicknesses. Common oscillatory WSS behavior was detected in all stent models. Comparing stent and vessel wall surfaces, maximum WSS values (in the order of 1Pa) were located on the stent surface area. WSS spatial distribution on the vascular wall surface showed decreasing values from the center of the vessel wall portion delimited by the stent struts to the wall regions close to the struts. The hemodynamic effects induced by two different thickness values for the same stent design were investigated, too, and a reduced extension of low WSS region (<0.5Pa) was observed for the model with a thicker strut.     

Computer modeling of deployment and mechanical expansion of neurovascular flow diverter in patient-specific intracranial aneurysms

Flow diverter (FD) is an emerging neurovascular device based on self-expandable braided stent for treating intracranial aneurysms. Variability in FD outcome has underscored a need for investigating the hemodynamic effect of fully deployed FD in patient-specific aneurysms. Image-based computational fluid dynamics, which can provide important hemodynamic insight, requires accurate representation of FD in deployed states. We developed a finite element analysis (FEA) based workflow for simulating mechanical deployment of FD in patient-specific aneurysms. We constructed FD models of interlaced wires emulating the Pipeline Embolization Device, using 3D finite beam elements to account for interactions between stent strands, and between the stent and other components. The FEA analysis encompasses all steps that affect the final deployed configuration including stent crimping, delivery and expansion. Besides the stent, modeling also includes key components of the FD delivery system such as microcatheter, pusher, and distal coil. Coordinated maneuver of these components allowed the workflow to mimic clinical operation of FD deployment and to explore clinical strategies. The workflow was applied to two patient-specific aneurysms. Parametric study indicated consistency of the deployment result against different friction conditions, but excessive intra-stent friction should be avoided. This study demonstrates for the first time mechanical modeling of braided FD stent deployment in cerebral vasculature to produce realistic deployed configuration, thus paving the way for accurate CFD analysis of flow diverters for reliable prediction and optimization of treatment outcome.     

Hemodynamics in coronary arteries with overlapping stents

Coronary artery stenosis is commonly treated by stent placement via percutaneous intervention, at times requiring multiple stents that may overlap. Stent overlap is associated with increased risk of adverse clinical outcome. While changes in local blood flow are suspected to play a role therein, hemodynamics in arteries with overlapping stents remain poorly understood. In this study we analyzed six cases of partially overlapping stents, placed ex vivo in porcine left coronary arteries and compared them to five cases with two non-overlapping stents. The stented vessel geometries were obtained by micro-computed tomography of corrosion casts. Flow and shear stress distribution were calculated using computational fluid dynamics. We observed a significant increase in the relative area exposed to low wall shear stress (WSS<0.5 Pa) in the overlapping stent segments compared both to areas without overlap in the same samples, as well as to non-overlapping stents. We further observed that the configuration of the overlapping stent struts relative to each other influenced the size of the low WSS area: positioning of the struts in the same axial location led to larger areas of low WSS compared to alternating struts. Our results indicate that the overlap geometry is by itself sufficient to cause unfavorable flow conditions that may worsen clinical outcome. While stent overlap cannot always be avoided, improved deployment strategies or stent designs could reduce the low WSS burden.     

Hemodynamics in stented vertebral artery ostial stenosis based on computational fluid dynamics simulations

Hemodynamic factors may affect the potential occurrence of in-stent restenosis (ISR) after intervention procedure of vertebral artery ostial stenosis (VAOS). The purpose of the present study is to investigate the influence of stent protrusion length in implantation strategy on the local hemodynamics of the VAOS. CTA images of a 58-year-old female patient with posterior circulation transient ischemic attack were used to perform a 3D reconstruction of the vertebral artery. Five models of the vertebral artery before and after the stent implantation were established. Model 1 was without stent implantation, Model 2–5 was with stent protruding into the subclavian artery for 0, 1, 2, 3 mm, respectively. Computational fluid dynamics simulations based on finite element analysis were employed to mimic the blood flow in arteries and to assess hemodynamic conditions, particularly the blood flow velocity and wall shear stress (WSS). The WSS and the blood flow velocity at the vertebral artery ostium were reduced by 85.33 and 35.36% respectively after stents implantation. The phenomenon of helical flow disappeared. Hemodynamics comparison showed that stent struts that protruded 1 mm into the subclavian artery induced the least decrease in blood speed and WSS. The results suggest that stent implantation can improve the hemodynamics of VAOS, while stent struts that had protruded 1 mm into the subclavian artery would result in less thrombogenesis and neointimal hyperplasia and most likely decrease the risk of ISR.     

Simulation of a balloon expandable stent in a realistic coronary artery—Determination of the optimum modelling strategy

Computational models of stent deployment in arteries have been widely used to shed light on various aspects of stent design and optimisation. In this context, modelling of balloon expandable stents has proved challenging due to the complex mechanics of balloon–stent interaction and the difficulties involved in creating folded balloon geometries. In this study, a method to create a folded balloon model is presented and utilised to numerically model the accurate deployment of a stent in a realistic geometry of an atherosclerotic human coronary artery. Stent deployment is, however, commonly modelled by applying an increasing pressure to the stent, thereby neglecting the balloon. This method is compared to the realistic balloon expansion simulation to fully elucidate the limitations of this procedure. The results illustrate that inclusion of a realistic balloon model is essential for accurate modelling of stent deformation and stent stresses. An alternative balloon simulation procedure is presented however, which overcomes many of the limitations of the applied pressure approach by using elements which restrain the stent as the desired diameter is achieved. This study shows that direct application of pressure to the stent inner surface may be used as an optimal modelling strategy to estimate the stresses in the vessel wall using these restraining elements and hence offer a very efficient alternative approach to numerically modelling stent deployment within complex arterial geometries. The method is limited however, in that it can only predict final stresses in the stented vessel and not those occurring during stent expansion, in which case the balloon expansion model is required.     

Numerical simulation of aneurysmal haemodynamics with calibrated porous-medium models of flow-diverting stents

Modelling flow-diverting (FD) stents as porous media (PM) markedly improves the efficiency of computational fluid dynamics (CFD) simulations in the study of intracranial aneurysm treatment. Nonetheless, the parameters of PM models adopted for simulations up until now were rarely calibrated to match the represented FD structure. We therefore sought to evaluate the PM parameters for a representative variety of commercially available stents, so characterising the flow-diversion behaviours of different FD devices on the market.We generated fully-resolved geometries for treatments using PED, Silk+, FRED, and dual PED stents. We then correspondingly derived the calibrated PM parameters—permeability (k) and inertial resistance factor (C₂)—for each stent design from CFD simulations, to ensure the calibrated PM model has identical flow resistance to the FD stent it represents. With each of the calibrated PM models respectively deployed in two aneurysms, we studied the flow-diversion effects of these stent configurations.This work for the first time reported several sets of parameters for PM models, which is vital to address the current knowledge gap and rectify the errors in PM model simulations, thereby setting right the modelling protocol for future studies using PM models. The flow resistance parameters were strongly affected by porosity and effective thickness of the commercial stents, and thus accounted for in the PM models. Flow simulations using the PM stent models revealed differences in aneurysmal mass flowrate (MFR) and energy loss (EL) between various stent designs.This study improves the practicability of FD simulation by using calibrated PM models, providing an individualised method with improved simulation efficiency and accuracy.     

Particle Imaging Velocimetry Evaluation of Intracranial Stents in Sidewall Aneurysm: Hemodynamic Transition Related to the Stent Design

We investigated the flow modifications induced by a large panel of commercial-off-the-shelf (COTS) intracranial stents in an idealized sidewall intracranial aneurysm (IA). Flow velocities in IA silicone model were assessed with and without stent implantation using particle imaging velocimetry (PIV). The use of the recently developed multi-time-lag method has allowed for uniform and precise measurements of both high and low velocities at IA neck and dome, respectively. Flow modification analysis of both regular (RSs) and flow diverter stents (FDSs) was subsequently correlated with relevant geometrical stent parameters. Flow reduction was found to be highly sensitive to stent porosity variations for regular stents RSs and moderately sensitive for FDSs. Consequently, two distinct IA flow change trends, with velocity reductions up to 50% and 90%, were identified for high-porosity RS and low-porosity FDS, respectively. The intermediate porosity (88%) regular braided stent provided the limit at which the transition in flow change trend occurred with a flow reduction of 84%. This transition occurred with decreasing stent porosity, as the driving force in IA neck changed from shear stress to differential pressure. Therefore, these results suggest that stents with intermediate porosities could possibly provide similar flow change patterns to FDS, favourable to curative thrombogenesis in IAs.     

Finite element comparison of performance related characteristics of balloon expandable stents

Cardiovascular stents are commonly made from 316L stainless steel and are deployed within stenosed arterial lesions using balloon expansion. Deployment involves inflating the balloon and plastically deforming the stent until the required diameter is obtained. This plastic deformation induces static stresses in the stent, which will remain for the lifetime of the device. In order to determine these stresses, finite element models of the unit cells of geometrically different, commercially available balloon expandable stents have been created, and deployment and elastic recoil have been simulated. In this work the residual stresses associated with deployment and recoil are compared for the various stent geometries, with a view to establishing appropriate initial stress states for fatigue loading for the stents. The maximum, minimum, and mean stresses induced in the stent due to systolic/diastolic pressure are evaluated, as are performance measures such as radial and longitudinal recoil.     

Fluid structure interaction with contact surface methodology for evaluation of endovascular carotid implants for drug-resistant hypertension treatment

Drug-resistant hypertensive patients may be treated by mechanical stimulation of stretch-sensitive baroreceptors located in the sinus of carotid arteries. To evaluate the efficacy of endovascular devices to stretch the carotid sinus such that the induced strain might trigger baroreceptors to increase action potential firing rate and thereby reduce systemic blood pressure, numerical simulations were conducted of devices deployed in subject-specific carotid models. Two models were chosen--a typical physiologic carotid and a diminutive atypical physiologic model representing a clinically worst case scenario--to evaluate the effects of device deployment in normal and extreme cases, respectively. Based on the anatomical dimensions of the carotids, two different device sizes were chosen out of five total device sizes available. A fluid structure interaction (FSI) simulation methodology with contact surface between the device and the arterial wall was implemented for resolving the stresses and strains induced by device deployment. Results indicate that device deployment in the carotid sinus of the physiologic model induces an increase of 2.5% and 7.5% in circumferential and longitudinal wall stretch, respectively, and a maximum of 54% increase in von Mises arterial stress at the sinus wall baroreceptor region. The second device, deployed in the diminutive carotid model, induces an increase of 6% in both circumferential and longitudinal stretch and a 50% maximum increase in von Mises stress at the sinus wall baroreceptor region. Device deployment has a minimal effect on blood-flow patterns, indicating that it does not adversely affect carotid bifurcation hemodynamics in the physiologic model. In the smaller carotid model, deployment of the device lowers wall shear stress at sinus by 16% while accelerating flow entering the external carotid artery branch. Our FSI simulations of carotid arteries with deployed device show that the device induces localized increase in wall stretch at the sinus, suggesting that this will activate baroreceptors and subsequently may control hypertension in drug-resistant hypertensive patients, with no consequential deleterious effects on the carotid sinus hemodynamics.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were 'treated' with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were ‘treated’ with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Local hemodynamic analysis after coronary stent implantation based on Euler-Lagrange method

Coronary stents are deployed to treat the coronary artery disease (CAD) by reopening stenotic regions in arteries to restore blood flow, but the risk of the in-stent restenosis (ISR) is high after stent implantation. One of the reasons is that stent implantation induces changes in local hemodynamic environment, so it is of vital importance to study the blood flow in stented arteries. Based on regarding the red blood cell (RBC) as a rigid solid particle and regarding the blood (including RBCs and plasma) as particle suspensions, a non-Newtonian particle suspensions model is proposed to simulate the realistic blood flow in this work. It considers the blood’s flow pattern and non-Newtonian characteristic, the blood cell-cell interactions, and the additional effects owing to the bi-concave shape and rotation of the RBC. Then, it is compared with other four common hemodynamic models (Newtonian single-phase flow model, Newtonian Eulerian two-phase flow model, non-Newtonian single-phase flow model, non-Newtonian Eulerian two-phase flow model), and the comparison results indicate that the models with the non-Newtonian characteristic are more suitable to describe the realistic blood flow. Afterwards, based on the non-Newtonian particle suspensions model, the local hemodynamic environment in stented arteries is investigated. The result shows that the stent strut protrusion into the flow stream would be likely to produce the flow stagnation zone. And the stent implantation can make the pressure gradient distribution uneven. Besides, the wall shear stress (WSS) of the region adjacent to every stent strut is lower than 0.5 Pa, and along the flow direction, the low-WSS zone near the strut behind is larger than that near the front strut. What’s more, in the regions near the struts in the proximal of the stent, the RBC particle stagnation zone is easy to be formed, and the erosion and deposition of RBCs are prone to occur. These hemodynamic analyses illustrate that the risk of ISR is high in the regions adjacent to the struts in the proximal and the distal ends of the stent when compared with struts in other positions of the stent. So the research can provide a suggestion on the stent design, which indicates that the strut structure in these positions of a stent should be optimized further.

     

Mis-sizing of stent promotes intimal hyperplasia: impact of endothelial shear and intramural stress

Stent can cause flow disturbances on the endothelium and compliance mismatch and increased stress on the vessel wall. These effects can cause low wall shear stress (WSS), high wall shear stress gradient (WSSG), oscillatory shear index (OSI), and circumferential wall stress (CWS), which may promote neointimal hyperplasia (IH). The hypothesis is that stent-induced abnormal fluid and solid mechanics contribute to IH. To vary the range of WSS, WSSG, OSI, and CWS, we intentionally mismatched the size of stents to that of the vessel lumen. Stents were implanted in coronary arteries of 10 swine. Intravascular ultrasound (IVUS) was used to size the coronary arteries and stents. After 4 wk of stent implantation, IVUS was performed again to determine the extent of IH. In conjunction, computational models of actual stents, the artery, and non-Newtonian blood were created in a computer simulation to yield the distribution of WSS, WSSG, OSI, and CWS in the stented vessel wall. An inverse relation (R(2) = 0.59, P < 0.005) between WSS and IH was found based on a linear regression analysis. Linear relations between WSSG, OSI, and IH were observed (R(2) = 0.48 and 0.50, respectively, P < 0.005). A linear relation (R(2) = 0.58, P < 0.005) between CWS and IH was also found. More statistically significant linear relations between the ratio of CWS to WSS (CWS/WSS), the products CWS × WSSG and CWS × OSI, and IH were observed (R(2) = 0.67, 0.54, and 0.56, respectively, P < 0.005), suggesting that both fluid and solid mechanics influence the extent of IH. Stents create endothelial flow disturbances and intramural wall stress concentrations, which correlate with the extent of IH formation, and these effects were exaggerated with mismatch of stent/vessel size. These findings reveal the importance of reliable vessel and stent sizing to improve the mechanics on the vessel wall and minimize IH.