Search

Heat transfer characteristics of random suspensions of 0.25 aspect ratio (AR) cylinders are investigated for Reynolds numbers (Re) between 10 and 300 and solid fraction (φ) ranging from 0.1 to 0.3 using particle resolved simulations. The effect of particle inclination with respect to flow and particle clustering on heat transfer is investigated. The Nusselt number decreases with an increase in inclination angle and the dependence becomes stronger as φ and Re increase. On the other hand, while prolate ellipsoid suspensions of AR 2.5 follow the same trend, the Nusselt number increases with inclination angle as AR increases to 5 and 10 and as φ increases. Local particle clustering nominally decreases the Nusselt number because of the dominance of thermal wakes. At low φ⁠, this effect is felt only at low Re⁠, but as φ increases, the effect spreads to higher Re⁠. Similar but weaker trends are also found in suspensions of prolate ellipsoids of AR 2.5, 5, and 10. High AR, low Re prolate ellipsoids exhibit the greatest dependence of Nusselt number on local solid fraction. Implementation of two independent definitions of reference length, i.e., volume equivalent sphere diameter deq for ellipsoids and diameter dp of the cylindrical particle in the correlation of Tavassoli et al. (2015, "Direct Numerical Simulation of Fluid-Particle Heat Transfer in Fixed Random Arrays of Non-Spherical Particles," Chem. Eng. Sci., 129, pp. 42-48) provides good estimates of the respective suspension mean Nusselt numbers. [ABSTRACT FROM AUTHOR]

Large eddy simulation calculations are conducted for flow in a channel with dimples and protrusions on opposite walls with both surfaces heated at three Reynolds numbers, [Re.sub.H] = 220, 940, and 9300, ranging from laminar, weakly turbulent, to fully turbulent, respectively. Turbulence generated by the separated shear layer in the dimple and along the downstream rim of the dimple is primarily responsible for heat transfer augmentation on the dimple surface. On the other hand, augmentation on the protrusion surface is mostly driven by flow impingement and flow acceleration between protrusions, while the turbulence generated in the wake has a secondary effect. Heat transfer augmentation ratios of 0.99 at [Re.sub.H] = 220, 2.9 at [Re.sub.H] = 940, and 2.5 at [Re.sub.H] = 9300 are obtained. Both skin friction and form losses contribute to pressure drop in the channel. Form losses increase from 45% to 80% with increasing Reynolds number Friction coefficient augmentation ratios of 1.67, 4.82, and 6.37 are obtained at [Re.sub.H] = 220, 940, and 9300, respectively. Based on the geometry studied, it is found that dimples and protrusions may not be viable heat transfer augmentation surfaces when the flow is steady and laminar

Detailed investigation of film cooling for a cylindrical leading edge is carried out using large eddy simulation (LES). The paper focuses on the effects of coolant to mainstream blowing ratio on flow features and, consequently, on the adiabatic effectiveness and heat transfer coefficient. With the advantage of obtaining unique, accurate, and dynamic results from LES, the influential coherent structures in the flow are identified. Describing the mechanism of jet-mainstream interaction, it is shown that as the blowing ratio increases, a more turbulent shear layer and stronger mainstream entrainment occur. The combined effects lead to a lower adiabatic effectiveness and higher heat transfer coefficient. Surface distribution and span-averaged profiles are shown for both adiabatic effectiveness and heat transfer (presented by Frossling number). Results are in good agreement with the experimental data of Ekkad et al. [1998, 'Detailed Film Cooling Measurement on a Cylindrical Leading Edge Model: Effect of Free-Steam Turbulence and Coolant Density,' ASME J. Turbomach., 120, pp. 799-807].

The problem of accurately predicting the flow and heat transfer in the ribbed internal cooling duct of a rotating gas turbine blade is addressed with the use of large eddy simulations (LES). Four calculations of the developing flow region of a rotating duct with ribs on opposite walls are used to study changes in the buoyancy parameter at a constant rotation rate. The Reynolds number is 20,000, the rotation number is 0.3, and the buoyancy parameter is varied between 0.00, 0.25, 0.45, and 0.65. Previous experimental studies have noted that leading wall heat transfer augmentation decreases as the buoyancy parameter increases with low buoyancy, but heat transfer then increases with high buoyancy. However, no consistent physical explanation has been given in the literature. The LES results from this study show that the initial decrease in augmentation with buoyancy is a result of larger separated regions at the leading wall However, as the separated region spans the full pitch between ribs with an increase in buoyancy parameter, it leads to increased turbulence and increased entrainment of mainstream fluid, which is redirected toward the leading wall by the presence of a rib. The impinging mainstream fluid results in heat transfer augmentation in the region immediately upstream of a rib. The results obtained from this study are in very good agreement with previous experimental results. [DOI: 10.1115/1.2437779]

A ribbed square duct (P/e= 10, e/[D.sub.h]=0.10) subjected to sand ingestion is studied using large-eddy simulations (LES). Particle sizes of 10 [micro]m, 50 [micro]m, and 100 [micro]m with nondimensional response times (normalized by friction velocity and hydraulic diameter) of 0.06875, 1.71875, and 6.875, respectively are considered. The calculations are performed for a nominal bulk Reynolds number of 20,O00 under fully-developed conditions. Distributions of impingement density, impingement velocities and angles, together with fractional energy transfer are presented for each surface. It is found that about 40% of the total number of 10 micron particles are concentrated in the vicinity (within 0.05 Oh) of the duct surfaces, compared to 25-30% of the 50 and 100 micron particles. The 10 micron particles are more sensitive to the primary and secondary flow velocities than the larger particles. While the 10 micron particles exhibit high energy transfer to the surface near the rib side-wall junction and immediately upstream of the rib, the larger particles exhibit more uniform distributions. The largest fraction of incoming particulate energy is transferred to the front face of the rib and is between one to two orders of magnitude larger than the other surfaces. As particle size increases, substantial particle energy is also transferred to the back face of the rib by particles bouncing off the front face and carrying enough momentum to impinge on the back face of the preceding rib. [DOI: 10.1115/1.2720509]

A nonlocal pressure equation is derived from mean-field free energy theory for calculating liquid-vapor systems. The proposed equation is validated analytically by showing that it reduces to van der Waals' square-gradient approximation under the assumption of slow density variations. The proposed nonlocal pressure is implemented in the mean-field free energy lattice Boltzmann method (LBM). The LBM is applied to simulate equilibrium liquid-vapor interface properties and interface dynamics of capillary waves and oscillating droplets in vapor. Computed results are validated with Maxwell constructions of liquid-vapor coexistence densities, theoretical relationship of variation of surface tension with temperature, theoretical planar interface density profiles, Laplace's law of capillarity, dispersion relationship between frequency and wave number of capillary waves, and the relationship between radius and the oscillating frequency of droplets in vapor. It is shown that the nonlocal pressure formulation gives excellent agreement with theory.

The capabilities of the detached eddy simulation (DES) and the unsteady Reynolds averaged Navier-Stokes (URANS) versions of the 1988 k-[omega] model in predicting the turbulent flow field in a two-pass internal cooling duct with normal ribs is presented. The flow is dominated by the separation and reattachment of shear layers; unsteady vorticity induced secondary flows and strong streamline curvature. The techniques are evaluated in predicting the developing flow at the entrance to the duct and downstream of the 180 deg bend, fully developed regime in the first pass, and in the 180 deg bend. Results of mean flow quantities, secondary flows, and the average friction factor are compared to experiments and large-eddy simulations (LES). DES predicts a slower flow development than LES, whereas URANS predicts it much earlier than LES computations and experiments. However, it is observed that as fully developed conditions are established, the capability of the base model in predicting the flow is enhanced by the DES formulation. DES accurately predicts the flow both in the fully developed region as well as the 180 deg bend of the duct. URANS fails to predict the secondary flows in the fully developed region of the duct and is clearly inferior to DES in the 180 deg bend. [DOI: 10.1115/1.2353279]

Large eddy simulation of the 180 deg bend in a stationary ribbed duct is presented. The domain studied includes three ribs upstream of the bend region and three ribs downstream of the bend with an outflow extension added to the end, using a total of 8.4 million cells. Two cases are compared to each other: one includes a rib in the bend and the other does not. The friction factor, mean flow, turbulence, and heat transfer are compared in the two cases to help explain the benefits and disadvantages of the wide number of flow effects seen in the bend, including flow separation at the tip of the dividing wall, counter-rotating Dean vortices, high heat transfer at areas of flow impingement, and flow separation at the upstream and downstream corners of the bend. Mean flow results show a region of separated flow at the tip of the dividing region in the case with no rib in the bend, but no separation region is observed in the case with a rib. A pair of counter-rotating Dean vortices in the middle of the bend is observed in both cases. Turbulent kinetic energy profiles show a 30% increase in the midplane of the bend when the rib is added. High gradients of heat transfer augmentation are observed on the back wall and downstream outside wall, where mean flow impingement occurs. This heat transfer is increased with the presence of a rib. Including a rib in the bend increases the friction factor in the bend by 80%, and it increases the heat transfer augmentation by approximately 20%, resulting in a trade-off between pressure drop and heat transfer.

Detached Eddy Simulation (DES) of a hydrodynamic and thermally developed turbulent flow is presented for a stationary duct with square ribs aligned normal to the main flow direction. The rib height to channel hydraulic diameter (e/[D.sub.h]) is 0.1, the rib pitch to rib height (P/e) is 10 and the calculations have been carried out for a bulk Reynolds number of 20,000. DES calculations are carried out on a [96.sup.3] grid, a [64.sup.3] grid, and a [48.sup.3] grid to study the effect of grid resolution. Based on the agreement with earlier LES computations, the [64.sup.3] grid is observed to be suitable for the DES computation. DES and RANS calculations carried out on the [64.sup.3] grid are compared to LES calculations on [96.sup.3]/ [128.sup.3] grids and experimental measurements. The flow and heat transfer characteristics for the DES cases compare well with the LES results and the experiments. The average friction and the augmentation ratios are consistent with experimental results, predicting values within 10% of the measured quantities, at a cost lower than the LES calculations. RANS fails to capture some key features of the flow. [DOI: 10.1115/1.2033010]

Results from large eddy simulations (LES) of fully developed flow in a 90 deg ribbed duct are presented with rib pitch-to-height ratio P/e=10 and a rib height-to-hydraulic-diameter ratio e/[D.sub.h] = 0.1. Three rotation numbers Ro = 0.18, 0.36, and 0.68 are studied at a nominal Reynolds number based on bulk velocity of 20 000. Centrifugal buoyancy effects are included at two Richardson numbers of Ri=12, 28 (Buoyancy parameter, Bo =0.12 and 0.30) for each rotation case. Heat transfer augmentation on the trailing side of the duct due to the action of Coriolis forces alone asymptotes to a value of 3.7 [+ or -] 5% by Ro = 0.2. On the other hand, augmentation ratios on the leading surface keep decreasing with an increase in rotation number with values ranging from 1.7 at Ro=0.18 to 1.2 at Ro=0.67. Secondary flow cells augment the heat transfer coefficient on the smooth walls by 20% to 30% over a stationary duct. Centrifugal buoyancy further strengthens the secondary flow cells in the duct cross-section which leads to an additional increase of 10% to 15%. Buoyancy also accentuates the augmentation of turbulence near the trailing wall of the duct and increases the heat transfer augmentation ratio 10% to 20% over the action of Coriolis forces alone. However, it does not have any significant effect at the leading side of the duct. The overall effect of buoyancy on heat transfer augmentation for the ribbed duct is found to be less than 10% over the effect of Coriolis forces alone. Friction on the other hand is augmented 15% to 20% at the highest buoyancy number studied. Comparison with available experiments in the literature show excellent agreement.

The impact of the cylindrical pin fins with tip clearances on the flow and heat transfer at low Reynolds number is discussed. The tip gaps are shown to reduce the form losses and friction factor of the system, hence affecting the heat transfer in the system.

The characteristics of detached eddy simulation (DES) are studied by application of internal cooling of turbine blades. The detached eddy simulation predicts the flow and heat transfer in ribbed ducts, rotation induced coriolis forces and centrifugal buoyancy.

Sand transport and deposition is investigated in a two-pass internal cooling ribbed geometry at near engine conditions. Large-eddy simulation (LES) calculations are performed for bulk Reynolds number of 25,000 to calculate flow field and heat transfer. Constant wall temperature boundary condition is used to investigate the effect of temperature on particle deposition. Three different wall temperatures of 950°C, 1000°C, and 1050°C are considered. Particle sizes in range 5-25 βm are considered. A new deposition model which accounts for particle composition, temperature, impact velocity and angle and material properties of particle and suiface is developed and applied. Calculated impingement and deposition patterns are discussed for different exposed surfaces in the two pass geometry. Other than the leading rib faces, the highest particle impingement and deposition is observed in the bend region and first quarter of the second pass. Significant deposition is observed in the two pass geometry for all three wall temperatures considered. Particle impingement and hence deposition is dominated by larger particles except in the downstream half of the bend region. In total, approximately 38%, 59%, and 67% of the injected particles deposit in the two passes, for the three wall temperatures of 950 °C, 1000°C, and 1050°C, respectively. While particle impingement is highest for wall temperature o f950 ° C , higher deposition is observed for 1000 °C and 1050 °C cases. Deposition increases significantly with wall temperature. For 1000°C, roughly 12% of the impacting particles deposit. For 1050°C, approximately 23% of the particles deposit on impact. For alt the three cases, the second pass experiences higher deposition compared to the first pass due to higher turbulence and direct impingement. [ABSTRACT FROM AUTHOR]

The incidence of stent late restenosis is high (Zwart et al., 2010, "Coronary Stent Thrombosis in the Current Era: Challenges and Opportunities for Treatment," Curr. Treat. Options Cardiovasc. Med., 12(1), pp. 46-57) despite the extensive use of stents, and is most prevalent at the proximal and distal ends of the stent. Elastic modulus change in stented coronary arteries subject to the motion of the myocardium is not studied extensively. It is our objective to understand and reveal the mechanism by which changes in elastic modulus and geometry contribute to the generation of nonphysiological wall shear stress (WSS). Such adverse hemodynamic conditions could have an effect on the onset of restenosis. Three-dimensional (3D), spatiotemporally resolved computational fluid dynamics (CFD) simulations of pulsatile flow with moving wall boundaries and fluid structure interaction (FSI) were carried out for a helical artery with physiologically relevant flow parameters. To study the effect of coronary artery (CA) geometry change on stent elastic modulus mismatch, models where the curvature, torsion and both curvature and torsion change were examined. The elastic modulus is increased by a factor of two, five, and ten in the stented section for all three modes of motion. The changes in elastic modulus and arterial geometry cause critical variations in the local pressure and velocity gradients and secondary flow patterns. The pressure gradient change is 47%, with respect to the unstented baseline when the elastic modulus is increased to 10. The corresponding WSS change is 15.4%. We demonstrate that these changes are attributed to the production of vorticity (vorticity flux) caused by the wall movement and elastic modulus discontinuity. The changes in curvature dominate torsion changes in terms of the effects to local hemodynamics. The elastic modulus discontinuities along with the dynamic change in geometry affected the secondary flow patterns and vorticity flux, which in turn affects the WSS. [ABSTRACT FROM AUTHOR]

Understanding flow through real porous media is of considerable importance given their significance in a wide range of applications. Direct numerical simulations of such flows are very useful in their fundamental understanding. Past works have focused mainly on ordered and disordered arrays of regular shaped structures such as cylinders o1 spheres to emulate porous media. More recently, extension of these studies to more realistic pore spaces are available in the literature highlighting the enormous potential of such studies in helping the fundamental understanding of pore-level flow physics. In an effort to advance the simulation of realistic porous media flows further, an immersed boundary method (IBM) framework capable of simulating flows through arbitrary surface contours is used in conjunction with a stochastic reconstruction procedure based on simulated annealing. The developed framework is tested in a two-dimensional channel with two types of" porous sections--one created using a random assembly of square blocks and another using the stochastic' reconstruction procedure. Numerous simulations are performed to demonstrate the capability of the developed framework. The computed pressure drops across the porous section are compared with predictions from the Darcy-Forchheimer equation for media composed of different structure sizes. Finally, the developed methodology is applied to study C02 diffusion in porous spherical particles of varying porosities. [ABSTRACT FROM AUTHOR]

A numerical study is performed to investigate the deposition of Syngas ash in the leading edge region of a turbine vane. The leading edge of the vane is modeled as a symmetric semicyUnder with a flat afterbody. Three rows of coolant holes located at stagnation and at ±21.3 deg from stagnation are simulated at blowing ratios of 0.5, 1.0, 1.5, and 2.0. Large eddy simulation (LES) is used to model the flow field of the coolant jet-mainstream interaction and Syngas ash particles are modeled using a discrete particle method. The capture efficiency for eight different ash compositions of particle sizes 5 and 10 microns are investigated. Under the conditions of the current simulations, both ash particles have Stokes numbers less than unity and hence are strongly affected by the flow and thermal field generated by the coolant interaction with the mainstream. Because of this, the cool-ant jets at stagnation are quite successful in pushing the particles away from the surface and minimizing deposition in the stagnation region. Among all of the ash samples, the ND ash sample shows the highest capture efficiency due to its low softening temperature. For the 5 micron particles, when the blowing ratio increases from 1.5 to 2.0, the percent-age of the capture efficiency increases as more numbers of particles are transported to the surface by strong mainstream entrainment by the coolant jets. The deposition results are also estimated using the discrete random walk (DRW) model and are compared to that obtained from the LES calculations. For both particle sizes, the DRW model under-predicts the capture efficiency when compared to the LES calculations and the difference increases with the increasing blowing ratio and decreases with increasing particle size. [ABSTRACT FROM AUTHOR]

Cardiovascular diseases are the number one cause of death in the world, making the understanding of hemodynamics and the development of treatment options imperative. The effect of motion of the coronary artery due to the motion of the myocardium is not extensively studied. In this work, we focus our investigation on the localized hemodynamic effects of dynamic changes in curvature and torsion. It is our objective to understand and reveal the mechanism by which changes in curvature and torsion contribute towards the observed wall shear stress distribution. Such adverse hemodynamic conditions could have an effect on circumferential intimal thickening. Three-dimensional spatiotemporally resolved computational fluid dynamics (CFD) simulations of pulsatile flow with moving wall boundaries were carried out for a simplified coronary artery with physiologically relevant flow parameters. A model with stationary walls is used as the baseline control case. In order to study the effect of curvature and torsion variation on local hemodynamics, this baseline model is compared to models where the curvature, torsion, and both curvature and torsion change. The simulations provided detailed information regarding the secondary flow dynamics. The results suggest that changes in curvature and torsion cause critical changes in local hemodynamics, namely, altering the local pressure and velocity gradients and secondary flow patterns. The wall shear stress (WSS) varies by a maximum of 22% when the curvature changes, by 3% when the torsion changes, and by 26% when both the curvature and torsion change. The oscillatory shear stress (OSI) varies by a maximum of 24% when the curvature changes, by 4% when the torsion changes, and by 28% when both the curvature and torsion change. We demonstrate that these changes are attributed to the physical mechanism associating the secondary flow patterns to the production of vorticity (vorticity flux) due to the wall movement. The secondary flow patterns and augmented vorticity flux affect the wall shear stresses. As a result, this work reveals how changes in curvature and torsion act to modify the near wall hemodynamics of arteries. [ABSTRACT FROM AUTHOR]

Large-eddy simulations are used to investigate Coriolis forces effect on flow structure and heat transfer in a rotating dimpled channel. Two geometries with two dimple depths are considered, δ=0.2 and 0.3 of channel height, for a wide range of rotation number Rob =0.0-0.70, based on mean bulk velocity and channel height. It is found that the turbulent flow is destabilized near the trailing side and stabilized near the leading side, with secondary flow structures generated in the channel under the effect of Coriolis forces. Higher heat transfer levels are obtained at the trailing surface of the channel, especially in regions offlow reattachment and boundary layer regeneration at the dimple surface. Coriolis forces showed a stronger effect on the flow structure for the shallow dimple geometry (δ=0.2) compared with the deeper dimple where the growth and shrinkage of the flow recirculation zone in the dimple cavity with rotation were more pronounced than the deep dimple geometry (δ = 0.3). Under the action of rotation, heat transfer augmentation increased by 57% for δ = 0.2 and by 70% for δ = 0.3 on the trailing side and dropped by 50% for δ = 0.2 and by 45% for δ = 0.3 on the leading side from that of the stationary case. [ABSTRACT FROM AUTHOR]

The paper investigates the deposition and erosion caused by Syngas ash particles in a film cooled leading edge region of a representative turbine vane. The carrier phase is predicted using large eddy simulation for three blowing ratios of 0.4, 0.8, and 1.2. Ash particle sizes of 1 μm, 3 μm, 5 μm, 7 μm, and 10 μm are investigated using Lagrangian dynamics. The 1 μm particles with momentum Stokes number, Stp=0.03 (based on approach velocity and leading edge diameter), follow the flow streamlines around the leading edge and few particles reach the blade surface. The 10 μm particles, on the other hand with a high momentum Stokes number, Stp=0.03, directly impinge on the surface, with blowing ratio having a minimal effect. The 3 μm, 5 μm, and 7 μm particles with Stp=0.03, 0.8 and 1.4, respectively, show some receptivity to coolant flow and blowing ratio. On a number basis, 85-90% of the 10 μm particles, 70-65% of 7 μm particles, 40-50% of 5 μm particles, 15% of 3 μm particles, and less than 1% of 1 μm particles deposit on the surface. Overall there is a slight decrease in percentage of particles deposited with increase in blowing ratio. On the other hand, the potential for erosive wear is highest in the coolant hole and is mostly attributed to 5 μm and 7 μm particles. It is only at BR=1.2 that 10 μm particles contribute to erosive wear in the coolant hole. [ABSTRACT FROM AUTHOR]

A nonlocal pressure equation is derived from mean-field free energy theory for calculating liquid-vapor systems. The proposed equation is validated analytically by showing that it reduces to van der Waals' square-gradient approximation under the assumption of slow density variations. The proposed nonlocal pressure is implemented in the mean-field free energy lattice Boltzmann method (LBM). The LBM is applied to simulate equilibrium liquid-vapor interface properties and interface dynamics of capillary waves and oscillating droplets in vapor. Computed results are validated with Maxwell constructions of liquid-vapor coexistence densities, theoretical relationship of variation of surface tension with temperature, theoretical planar interface density profiles, Laplace's law of capillarity, dispersion relationship between frequency and wave number of capillary waves, and the relationship between radius and the oscillating frequency of droplets in vapor. It is shown that the nonlocal pressure formulation gives excellent agreement with theory. [ABSTRACT FROM AUTHOR]

Previous investigations on the performance of straight pins, pins with tip clearance, and profiled fins showed that closely packed cylindrical pin fins are very competitive with the modified pins. Therefore, the objective of this paper is to investigate the effect of pin density on performance. Steady/time-dependent calculations are performed to investigate the effect of pin density on friction and heat transfer. Pins packed at distances of SD=1.1, 1.2, 1.3, 1.4, 1.5, 2.0, 2.5, and 3 pin diameters (D) are investigated for 10≤ ReD≤600. Two performance measures are used to compare the different pin fin densities. The first measure is to maximize heat transfer capacity for a given pumping power compared with a plane channel. The second measure used is based on entropy generation minimization (EGM), where the objective is to reduce the total irreversibility of the pin fin array to obtain an optimal spacing. Based on the performance measure of maximizing heat capacity, it is shown that for plain channels operating in the laminar range using denser pin packing has distinct advantages with SD=1.1 providing the best augmentation. However, the augmentation in heat capacity becomes relatively independent of the pin density for a channel operating in the turbulent regime. Based on the EGM method, at ReD>200, SD=1.3, 1.4, and 1.5 are the most suitable, with the least entropy generation observed at SD=1.4. At ReD<200, SD=1.1, 1.2, and 1.3 are also suitable for keeping entropy generation low. [ABSTRACT FROM AUTHOR]

Cardiovascular diseases are the number one cause of death in the world, making the understanding of hemodynamics and development of treatment options imperative. The most common modality for treatment of occlusive coronary artery diseases is the use of stents. Stent design profoundly influences the postprocedural hemodynamic and solid mechanical environment of the stented artery. However, despite their wide acceptance, the incidence of stent late restenosis is still high (Zwart et al., 2010, "Coronary Stent Thrombosis in the Current Era: Challenges and Opportunities for Treatment," Current Treatment Options in Cardiovascular Medicine, 12(1), pp. 46-57), and it is most prevailing at the proximal and distal ends of the stent. In this work, we focus our investigation on the localized hemodynamic effects of compliance mismatch due to the presence of a stent in an artery. The compliance mismatch in a stented artery is maximized at the proximal and distal ends of the stent. Hence, it is our objective to understand and reveal the mechanism by which changes in compliance contribute to the generation of nonphysiological wall shear stress (WSS). Such adverse hemodynamic conditions could have an effect on the onset of restenosis. Three-dimensional, spatiotemporally resolved computational fluid dynamics simulations of pulsatile flow with fluid-structure interaction were carried out for a simplified coronary artery with physiologically relevant flow parameters. A model with uniform elastic modulus is used as the baseline control case. In order to study the effect of compliance variation on local hemodynamics, this baseline model is compared with models where the elastic modulus was increased by two-, five-, and tenfold in the middle of the vessel. The simulations provided detailed information regarding the recirculation zone dynamics formed during flow reversals. The results suggest that discontinuities in compliance cause critical changes in local hemodynamics, namely, altering the local pressure and velocity gradients. The change in pressure gradient at the discontinuity was as high as 90%. The corresponding changes in WSS and oscillatory shear index calculated were 9% and 15%, respectively. We demonstrate that these changes are attributed to the physical mechanism associating the pressure gradient discontinuities to the production of vorticity (vorticity flux) due to the presence of the stent. The pressure gradient discontinuities and augmented vorticity flux are affecting the wall shear stresses. As a result, this work reveals how compliance variations act to modify the near wall hemodynamics of stented arteries.