Simulation of particle-hemodynamics in a partially occluded artery segment with implications to the initiation of microemboli and secondary stenoses.

Computational results of laminar incompressible blood-particle flow analyses in an axisymmetric artery segment with a smooth local area constriction of 75 percent have been presented. The flow input waveform was sinusoidal with a nonzero average. The non-Newtonian behavior of blood was simulated with a modified Quemada model, platelet concentrations were calculated with a drift-flux model, and monocyte trajectories were described and compared for both Newtonian and Quemada rheologies. Indicators of "disturbed flow" included the time-averaged wall shear stress (WSS), the oscillatory shear index (OSI), and the wall shear stress gradient (WSSG). Implications of the vortical flow patterns behind the primary stenosis to the formation of microemboli and downstream stenoses are as follows. Elevated platelet concentrations due to accumulation in recirculation zones mixed with thrombin and ADP complexes assumed to be released upstream in high wall shear stress regions, could form microemboli, which are convected downstream. Distinct near-wall vortices causing a local increase in the WSSG and OSI as well as blood-particle entrainment with possible wall deposition, indicate sites susceptible to the onset of an additional stenosis proximal to the initial geometric disturbance.     

A numerical study on hemodynamics in the left coronary bifurcation with normal and hypertension conditions

In this study, a three-dimensional analysis of the non-Newtonian blood flow was carried out in the left coronary bifurcation. The Casson model and hyperelastic and rigid models were used as the constitutive equation for blood flow and vessel wall model, respectively. Physiological conditions were considered first normal and then compliant with hypertension disease with the aim of evaluating hemodynamic parameters and a better understanding of the onset and progression of atherosclerosis plaques in the coronary artery bifurcation. Two-way fluid–structure interaction method applying a fully implicit second-order backward Euler differencing scheme has been used which is performed in the commercial code ANSYS and ANSYS CFX (version 15.0). When artery deformations and blood pressure are associated, arbitrary Lagrangian–Eulerian formulation is employed to calculate the artery domain response using the temporal blood response. As a result of bifurcation, noticeable velocity reduction and backflow formation decrease shear stress and made it oscillatory at the starting point of the LCx branch which caused the shear stress to be less than 1 and 2 Pa in the LCx and the LAD branches, respectively. Oscillatory shear index (OSI) as a hemodynamic parameter represents the increase in residence time and oscillatory wall shear stress. Because of using the ideal 3D geometry and realistic physiological conditions, the values obtained for shear stress are more accurate than the previous studies. Comparing the results of this study with previous clinical investigations shows that the regions with low wall shear stress less than 1.20 Pa and with high OSI value more than 0.3 are in more potential risk to the atherosclerosis plaque development, especially in the posterior after the bifurcation.     

Numerical simulation of haemodynamics of the descending aorta in the non-diabetic and diabetic rabbits

Diabetes mellitus (DM) is a predisposing risk factor leading to macrovascular diseases. Changes in haemodynamics of the diabetic aortas remain largely unclear and relevant computational analyses are lacking in the literature. Ten adult rabbits (1.6–2.2 kg) were collected and the type I diabetic rabbit model was induced by injection of alloxan. A total of five control and five diabetic rabbit aortas were considered for subsequent numerical simulation. The CT scanning was performed to reconstruct three-dimensional model of the individual rabbit descending aorta. The flow velocity waveforms were measured by ultrasound machine and were set to be the inlet boundary conditions. The reconstructed aortas were then imported into ANSYS to perform mesh generation and computational analysis. Results showed that the distributions of haemodynamic indicators time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and transverse wall shear stress (transWSS) in the non-diabetic rabbit aortas were similar to those in the diabetic rabbit aortas. However, the mean values of TAWSS and transWSS in the non-diabetic rabbit aortas were significantly higher than those values in the diabetic rabbit aortas (TAWSS: p = 0.04; transWSS: p = 0.02). The back of right renal artery tended to have high OSI in both the non-diabetic and the diabetic rabbit aortas. Notably, the regions with high OSI tended to have intense disturbed flow and low TAWSS in the most diabetic rabbit aortas. The results suggest that diabetes leads to changes in haemodynamic parameters in the rabbit aortas. In particular, the lower TAWSS and the higher OSI within the diabetic aortas may further contribute to aortic wall remodeling.     

Flow patterns in three-dimensional porcine epicardial coronary arterial tree

The branching pattern of epicardial coronary arteries is clearly three-dimensional, with correspondingly complex flow patterns. The objective of the present study was to perform a detailed hemodynamic analysis using a three-dimensional finite element method in a left anterior descending (LAD) epicardial arterial tree, including main trunk and primary branches, based on computed tomography scans. The inlet LAD flow velocity was measured in an anesthetized pig, and the outlet pressure boundary condition was estimated based on scaling laws. The spatial and temporal wall shear stress (WSS), gradient of WSS (WSSG), and oscillatory shear index (OSI) were calculated and used to identify regions of flow disturbances in the vicinity of primary bifurcations. We found that low WSS and high OSI coincide with disturbed flows (stagnated, secondary, and reversed flows) opposite to the flow divider and lateral to the junction orifice of the main trunk and primary branches. High time-averaged WSSG occurs in regions of bifurcations, with the flow divider having maximum values. Low WSS and high OSI were found to be related through a power law relationship. Furthermore, zones of low time-averaged WSS and high OSI amplified for larger diameter ratio and high inlet flow rate. Hence, different focal atherosclerotic-prone regions may be explained by different physical mechanism associated with certain critical levels of low WSS, high OSI, and high WSSG, which are strongly affected by the diameter ratio. The implications of the flow patterns for atherogenesis are enumerated.     

Hemodynamic factors at the distal end-to-side anastomosis of a bypass graft with different POS:DOS flow ratios

A pulsatile flow in vitro model of the distal end-to-side anastomosis of an arterial bypass graft was used to examine the effects that different flow ratios between the proximal outlet segment (POS) and the distal outlet segment (DOS) have on the flow patterns and the distributions of hemodynamic factors in the anastomosis. Amberlite particles were tracked by flow visualization to determine overall flow patterns and velocity measurements were made with Laser Doppler anemometry (LDA) to obtain detailed hemodynamic factors along the artery floor and the graft hood regions. These factors included wall shear stress (WSS), spatial wall shear stress gradient (WSSG), and oscillatory index (OSI). Statistical analysis was used to compare these hemodynamic factors between cases having different POS:DOS flow ratios (Case 1-0:100, Case 2-25:75, Case 3-50:50). The results showed that changes in POS:DOS flow ratios had a great influence on the flow patterns in the anastomosis. With an increase in proximal outlet flow, the range of location of the stagnation point along the artery floor decreased, while the extent of flow separation along the graft hood increased. The statistical results showed that there were significant differences (p<0.05) for the mean WSS between cases along the graft hood, but no significant differences were detected along the artery floor. There were no significant differences for the spatial WSSG along both the artery floor and the graft hood. However, there were significant differences (p<0.05) in the mean OSI between Cases 1 and 2 and between Cases 1 and 3 both along the artery floor and along the graft hood. Comparing these mechanical factors with histological findings of intimal hyperplasia formation obtained by previous canine studies, the results of the statistical analysis suggest that regions exposed to a combination of low mean WSS and high OSI may be most prone to the formation of intimal hyperplasia.     

Large eddy simulation of LDL surface concentration in a subject specific human aorta

The development of atherosclerosis is correlated to the accumulation of lipids in the arterial wall, which, in turn, may be caused by the build-up of low-density lipoproteins (LDL) on the arterial surface. The goal of this study was to model blood flow within a subject specific human aorta, and to study how the LDL surface concentration changed during a cardiac cycle. With measured velocity profiles as boundary conditions, a scale-resolving technique (large eddy simulation, LES) was used to compute the pulsatile blood flow that was in the transitional regime. The relationship between wall shear stress (WSS) and LDL surface concentration was investigated, and it was found that the accumulation of LDL correlated well with WSS. In general, regions of low WSS corresponded to regions of increased LDL concentration and vice versa. The instantaneous LDL values changed significantly during a cardiac cycle; during systole the surface concentration was low due to increased convective fluid transport, while in diastole there was an increased accumulation of LDL on the surface. Therefore, the near-wall velocity was investigated at four representative locations, and it was concluded that in regions with disturbed flow the LDL concentration had significant temporal changes, indicating that LDL accumulation is sensitive to not only the WSS but also near-wall flow.     

Local hemodynamics affect monocytic cell adhesion to a three-dimensional flow model coated with E-selectin.

Monocyte adhesion to the endothelium depends on concentrations of receptors/ligands, local concentrations of chemoattractants, monocyte transport to the endothelial surface and hemodynamic forces. Monocyte adhesion to the inert surface of a three-dimensional perfusion model was shown to correlate inversely with wall shear stress, but was also affected by flow patterns which influenced the near-wall cell availability. We hypothesized that (a) under the same flow conditions, insolubilized E-selectin on the model's surface may mediate adhesive interactions at higher wall shear stresses, compared to an uncoated model, and (b) pulsatile flow may modify the adhesion profile obtained under steady flow. An axisymmetric flow model with a stenosis and a sudden expansion produced a range of wall shear stresses and a separated flow with recirculation and reattachment. Pre-activated U937 cells were perfused through the model under either steady (Re = 100, 140) or pulsatile (Remean = 107) flow. The velocity field was characterized through computational fluid dynamics and validated by inert particle tracking. Surface E-selectin greatly increased cell adhesion in all regions at Re = 100 and 140, compared to an uncoated model under the same flow conditions. In regions where the cells near the wall were abundant (taper and stenosis), adhesion to E-selectin correlated with the reciprocal of local wall shear stress when flow was steady. Pulsatile flow distributed the adherent cells more evenly throughout the coated model. Hence, characterizing both the local hemodynamics and the biological activity on the vessel wall is important in leukocyte adhesion.     

Computation in the rabbit aorta of a new metric – the transverse wall shear stress – to quantify the multidirectional character of disturbed blood flow

Spatial variation of the haemodynamic stresses acting on the arterial wall is commonly assumed to explain the focal development of atherosclerosis. Disturbed flow in particular is thought to play a key role. However, widely-used metrics developed to quantify its extent are unable to distinguish between uniaxial and multidirectional flows. We analysed pulsatile flow fields obtained in idealised and anatomically-realistic arterial geometries using computational fluid dynamics techniques, and in particular investigated the multidirectionality of the flow fields, capturing this aspect of near-wall blood flow with a new metric – the transverse wall shear stress (transWSS) – calculated as the time-average of wall shear stress components perpendicular to the mean flow direction. In the idealised branching geometry, multidirectional flow was observed downstream of the branch ostium, a region of flow stagnation, and to the sides of the ostium. The distribution of the transWSS was different from the pattern of traditional haemodynamic metrics and more dependent on the velocity waveform imposed at the branch outlet. In rabbit aortas, transWSS patterns were again different from patterns of traditional metrics. The near-branch pattern varied between intercostal ostia, as is the case for lesion distribution; for some branches there were striking resemblances to the age-dependent patterns of disease seen in rabbit and human aortas. The new metric may lead to improved understanding of atherogenesis.     

Hemodynamics and wall mechanics in human carotid bifurcation and its consequences for atherogenesis: investigation of inter-individual variation

Finite element simulations of fluid-solid interactions were used to investigate inter-individual variations in flow dynamics and wall mechanics at the carotid artery bifurcation, and its effects on atherogenesis, in three healthy humans (normal volunteers: NV1, NV2, NV4). Subject-specific calculations were based on MR images of structural anatomy and ultrasound measurements of flow at domain boundaries. For all subjects, the largest contiguous region of low wall shear stress (WSS) occurred at the carotid bulb, WSS was high (6-10 Pa) at the apex, and a small localized region of WSS > 10 Pa occurred close to the inner wall of the external carotid artery (ECA). NV2 and NV4 had a "spot" of low WSS distal to the bifurcation at the inner wall of the ECA. Low WSS patches in the common carotid artery (CCA) were contiguous with the carotid bulb low WSS region in NV1 and NV2, but not in NV4. In all three subjects, areas of high oscillatory shear index (OSI) were confined to regions of low WSS. Only NV4 exhibited high levels of OSI on the external adjoining wall of the ECA and CCA. For all subjects, the maximum wall shear stress temporal gradient (WSSTG) was highest at the flow divider (reaching 1,000 Pa/s), exceeding 300 Pa/s at the walls connecting the ECA and CCA, but remaining below 250 Pa/s outside of the ECA. In all subjects, (maximum principle) cyclic strain (CS) was greatest at the apex (NV1: 14%; NV2: 11%; NV4: 6%), and a second high CS region occurred at the ECA-CCA adjoining wall (NV1: 11%, NV2: 9%, NV4: 5%). Wall deformability was included in one simulation (NV2) to verify that it had little influence on the parameters studied. Location and magnitude of low WSS were similar, except for the apex (differences of up to 25%). Wall distensibility also influenced OSI, doubling it in most of the CCA, separating the single high OSI region of the carotid bulb into two smaller regions, and shrinking the ECA internal and external walls' high OSI regions. These observations provide further evidence that significant intra-subject variability exists in those factors thought to impact atherosclerosis.     

The influence of out-of-plane geometry on pulsatile flow within a distal end-to-side anastomosis

We present an experimental and computational investigation of time-varying flow in an idealized fully occluded 45 degrees distal end-to-side anastomosis. Two geometric configurations are assessed, one where the centerlines of host and bypass vessels lie within a plane, and one where the bypass vessel is deformed out of the plane of symmetry, respectively, termed planar and non-planar. Flow experiments were conducted by magnetic resonance imaging in rigid wall models and computations were performed using a high order spectral/hp algorithm. Results indicate a significant change in the spatial distribution of wall shear stress and a reduction of the time-averaged peak wall shear stress magnitude by 10% in the non-planar model as compared to the planar configuration. In the planar geometry the stagnation point follows a straight-line path along the host artery bed with a path length of 0.8 diameters. By contrast in the non-planar case the stagnation point oscillates about a center that is located off the symmetry plane intersection with the host artery bed wall, and follows a parabolic path with a 0.7 diameter longitudinal and 0.5 diameter transverse excursion. A definition of the oscillatory shear index (OSI) is introduced that varies between 0 and 0.5 and that accounts for a continuous range of wall shear stress vector angles. In both models, regions of elevated oscillatory shear were spatially associated with regions of separated or oscillating stagnation point flow. The mean oscillatory shear magnitude (considering sites where OSI>0.1) in the non-planar geometry was reduced by 22% as compared to the planar configuration. These changes in the dynamic behavior of the stagnation point and the oscillatory shear distribution introduced by out-of-plane graft curvature may influence the localization of vessel wall sites exposed to physiologically unfavorable flow conditions.     

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.     

Flow patterns and wall shear stress distribution in human internal carotid arteries: the geometric effect on the risk for stenoses

It has been widely observed that atherosclerotic stenosis occurs at sites with complex hemodynamics, such as arteries with high curvature or bifurcations. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present study, 3D sinusoidally pulsatile blood flow through the models of internal carotid artery (ICA) with different geometries was investigated with computational simulation. Three preferred sites of stenoses were found along the carotid siphon with low and highly oscillatory WSS. The risk for stenoses at these sites was scaled with the values of time-averaged WSS and oscillating shear index (OSI). The local risk for stenoses at every preferred site of stenoses was found different between 3 types of ICA, indicating that the geometry of the blood vessel plays significant roles in the atherogenesis. Specifically, the large curvature and planarity of the vessel were found to increase the risk for stenoses, because they tend to lower WSS and elevate OSI. Therefore, the geometric study makes it possible to estimate the stenosis location in the ICA siphon as long as the shape of ICA was measured.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

Experimental and numerical study of pulsatile flows through stenosis: wall shear stress analysis.

Different shapes of pulsatile flows through a model of stenosis are experimentally and numerically modeled to validate both methods and to determine the wall shear stress temporal evolution downstream from the stenosis. Two-dimensional velocity measurements are performed in a 75% severity stenosis using a pulsed Doppler ultrasonic velocimeter. Finite element package is employed for the transient numerical simulations. Polynomial method, based on the experimental velocity values, is proposed to determine the wall shear stress temporal evolution. There is a good agreement between the numerical and experimental results. The wall shear stress temporal analysis shows oscillating wall shear stress values during the cycle with high wall shear stress values at the throat of about 120 dyn/cm2, and low values downstream from the stenosis of about - 2.5 dyn/cm2. The key result of the study is that the presence of the stenosis leads the artery to work in a direction which is opposite to the direction of a healthy artery.     

Non-Newtonian blood flow in human right coronary arteries: steady state simulations

This study looks at blood flow through four different right coronary arteries, which have been reconstructed from bi-plane angiograms. Five non-Newtonian blood models, as well as the usual Newtonian model of blood viscosity, are used to study the wall shear stress in each of these arteries at a particular point in the cardiac cycle. It was found that in the case of steady flow in a given artery, the pattern of wall shear stress is consistent across all models. The magnitude of wall shear stress, however, is influenced by the model used and correlates with graphs of shear stress versus strain for each model. For mid-range velocities of around 0.2 m s(-1) the models are virtually indistinguishable. Local and global non-Newtonian importance factors are introduced, in an attempt to quantify the types of flows where non-Newtonian behaviour is significant. It is concluded that, while the Newtonian model of blood viscosity is a good approximation in regions of mid-range to high shear, it is advisable to use the Generalised Power Law model (which tends to the Newtonian model in those shear ranges in any case) in order to achieve better approximation of wall shear stress at low shear.     

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.     

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.