Analysis of Temporal Shear Stress Gradients During the Onset Phase of Flow Over a Backward-Facing Step

[+] Author and Article Information
Mark A. Haidekker, Charles R. White, John A. Frangos

University of California, San Diego, Department of Bioengineering, La Jolla, CA 92092-0412

J Biomech Eng 123(5), 455-463 (Apr 17, 2001) (9 pages) doi:10.1115/1.1389460 History: Received July 06, 2000; Revised April 17, 2001
Copyright © 2001 by ASME
Your Session has timed out. Please sign back in to continue.


Glagov,  S., Zarins,  C., Giddens,  D. P., and Ku,  D. N., 1988, “Hemodynamics and Atherosclerosis. Insights and Perspectives Gained From Studies of Human Arteries,” Arch. Pathol. Lab. Med., 112, pp. 1018–1031.
Davies,  P. F., Barbee,  K. A., Volin,  M. V., Robotewskyj,  A., Chen,  J., Joseph,  L., Griem,  M. L., Wernick,  M. N., Jacobs,  E., Polacek,  D. C., Depaola,  N., and Barakat,  A. I., 1999, “Spatial Relationships in Early Signaling Events of Flow-Mediated Endothelial Mechanotransduction,” Annu. Rev. Physiol., 59, pp. 527–549.
Davies,  P. F., Polacek,  D. C., Handen,  J. S., Helmke,  B. P., and Depaola,  N., 1999, “A Spatial Approach to Transcriptional Profiling: Mechanotransduction and the Focal Origin of Atherosclerosis,” Trends Biotechnol., 17, pp. 347–351.
Kuchan,  M. J., and Frangos,  J. A., 1994, “Role of Calcium and Calmodulin in Flow-Induced Nitric Oxide Production in Endothelial Cells,” Am. J. Physiol., 266, pp. C628–C636.
Frangos,  J. A., Eskin,  S. G., Mcintire,  L. V., and Ives,  C. L., 1985, “Flow Effects on Prostacyclin Production by Cultured Human Endothelial Cells,” Science, 227, pp. 1477–1479.
Levesque,  M. J., Nerem,  R. M., and Sprague,  E. A., 1990, “Vascular Endothelial Cell Proliferation in Culture and the Influence of Flow,” Biomaterials, 11, pp. 702–707.
Chiu,  J. J., Wang,  D. L., Chien,  S., Skalak,  R., and Usami,  S., 1998, “Effects of Disturbed Flow on Endothelial Cells,” ASME J. Biomech. Eng., 120, pp. 2–8.
Depaola,  N., Gimbrone,  M. A., Davies,  P. F., and Dewey,  C. F., 1992, “Vascular Endothelium Responds to Fluid Shear Stress Gradients,” Arterioscler. Thromb., 12, pp. 1254–1257.
Ojha,  M., 1994, “Wall Shear Stress Temporal Gradient and Anastomotic Intimal Hyperplasia,” Circ. Res., 74, pp. 1227–1231.
Ku,  D. N., Giddens,  D. P., Zarins,  C. K., and Glagov,  S., 1985, “Pulsatile Flow and Atherosclerosis in the Human Carotid Bifurcation. Positive Correlation Between Plaque Location and Low Oscillating Shear Stress,” Arteriosclerosis 5, 293–302.
Bao,  X., Lu,  C., and Frangos,  J. A., 1999, “Temporal Gradient in Shear but Not Steady Shear Stress Induces PDGF-A and MCP-1 Expression in Endothelial Cells: Role of NO, NF Kappa B, and Egr-1,” Arterioscler., Thromb., Vasc. Biol., 19, pp. 996–1003.
Bao,  X., Clark,  C. B., and Frangos,  J. A., 2000, “Temporal Gradient in Shear-Induced Signaling Pathway: Involvement of MAP Kinase, C-Fos, and Connexin43,” Am. J. Physiol., 278, pp. H1598–H1605.
Moore,  J. E., Xu,  C., Glagov,  S., Zarins,  C. K., and Ku,  D. N., 1994, “Fluid Wall Shear Stress Measurements in a Model of the Human Abdominal Aorta: Oscillatory Behavior and Relationship to Atherosclerosis,” Atherosclerosis, 110, pp. 225–240.
Truskey,  G. A., Barber,  K. M., Robey,  T. C., Olivier,  L. A., and Combs,  M. P., 1995, “Characterization of a Sudden Expansion Flow Chamber to Study the Response of Endothelium to Flow Recirculation,” ASME J. Biomech. Eng., 117, pp. 203–210.
Durst,  F., and Pereira,  J. C. F., 1988, “Time-Dependent Laminar Backward-Facing Step Flow in a Two-Dimensional Duct,” ASME J. Fluids Eng., 110, pp. 289–296.
Armaly,  B. F., Durst,  F., Pereira,  J. C. F., and Schönung,  B., 1982, “Experimental and Theoretical Investigation of Backward-Facing Step-Flow,” J. Fluid Mech., 127, pp. 473–496.
Kaiktsis,  L., Karniadakis,  G. E., and Orszag,  S. A., 1991, “Onset of Three-Dimensionality, Equilibria, and Early Transition in Flow Over a Backward-Facing Step,” J. Fluid Mech., 231, pp. 501–528.
Frangos,  J. A., Mcintire,  L. V., and Eskin,  S. G., 1988, “Shear Stress Induced Stimulation of Mammlian Cell Metabolism,” Biotechnol. Bioeng., 32, pp. 1053–1060.
White,  C. W., Haidekker,  M. A., and Frangos,  J. A., 2001, “Temporal Gradients in Shear, but Not Spatial Gradients, Stimulate Endothelial Cell Proliferation,” Circulation, 103, pp. 2508–2513.
Wright,  H. P., 1972, “Mitosis Patterns in Aortic Endothelium,” Atherosclerosis, 15, pp. 93–100.
Depaola,  N., Davies,  P. F., Pritchard,  W. F., Florez,  L., Harbeck,  N., and Polacek,  D. C., 1999, “Spatial and Temporal Regulation of Gap Junction Connexin43 in Vascular Endothelial Cells Exposed to Controlled Disturbed Flows in Vitro,” Proc. Natl. Acad. Sci. U.S.A., 96, pp. 3154–3159.
Nagel,  T., Resnick,  N., Dewey,  C. F., and Gimbrone,  M. A., 1999, “Vascular Endothelial Cells Respond to Spatial Gradients in Fluid Shear Stress by Enhanced Activation of Transcription Factors,” Arterioscler., Thromb., Vasc. Biol., 19, pp. 1825–1834.
Long,  Q., Xu,  X. Y., Ariff,  B., Thom,  S. A., Hughes,  A. D., and Stanton,  A. V., 2000, “Reconstruction of Blood Flow Patterns in a Human Carotid Bifurcation: a Combined CFD and MRI Study,” J. Magn. Reson Imaging, 11, pp. 299–311.
Pedersen,  E. M., Agerbaek,  M., Kristensen,  I. B., and Yoganathan,  A. P., 1997, “Wall Shear Stress and Early Atherosclerotic Lesions in the Abdominal Aorta in Young Adults,” Eur. J. Vasc. Endovasc Surg., 13, pp. 443–451.
Ziegler,  T., and Nerem,  R. M., 1994, “Effect of Flow on the Process of Endothelial Cell Division,” Arterioscler. Thromb., 14, pp. 636–643.
Lutz,  R. J., Cannon,  J. N., Bischoff,  K. B., Dedrick,  R. L., Stiles,  R. K., and Fry,  D. L., 1977, “Wall Shear Stress Distribution in a Model Canine Artery During Steady Flow,” Circ. Res., 41, pp. 391–399.
Tardy,  Y., Resnick,  N., Nagel,  T., Gimbrone,  M. A., and Dewey,  C. F. , 1997, “Shear Stress Gradients Remodel Endothelial Monolayers in Vitro Via a Cell Proliferation-Migration-Loss Cycle,” Arterioscler., Thromb., Vasc. Biol., 17, pp. 3102–3106.
Griebel, M., Dornseifer, T., and Neunhoeffer, T., 1998, “Numerical Simulations in Fluid Dynamics,” SIAM, Philadelphia.
Young,  D. F., 1979, “Fluid Mechanics of Arterial Stenoses,” ASME J. Biomech. Eng., 101, pp. 157–175.
Calvo,  W. J., Hajduczok,  G., Russell,  J. A., and Diamond,  S. L., 1999, “Inhibition of Nitric Oxide but Not Prostacyclin Prevents Poststenotic Dilatation in Rabbit Femoral Artery,” Circulation, 99, pp. 1069–1076.
Yongchareon,  W., and Young,  D. F., 1979, “Initiation of Turbulence in Models of Arterial Stenoses,” J. Biomech., 12, pp. 185–196.


Grahic Jump Location
Geometry of the backward-facing step flow chamber with the dimensions used in the computer simulation. The origin of the coordinate system is the lower right corner of the expansion step. In order to keep the grid size small, only the first 15 mm downstream of the expansion site were considered in the model, since the fully reattached parabolic-profile flow is established at that point. Also shown are the streamlines of the steady-state flow at Re=243 with the second stagnation point at 4.8 mm.
Grahic Jump Location
Upper panel: Calculated changes of the highest shear stress values with increasing grid size. Resolutions of 70×14,150×28,300×50,600×100, and 1200×200 grid nodes were chosen, and the steady-state flow at Re=243 computed. There is no significant change (less than 3 percent) between the two largest grid sizes. Lower panel: Different time step sizes (parameter κ in Eq. (4)) were used to verify the convergence of the discretization of the momentum equation. No significant changes were observed for κ≤0.6.
Grahic Jump Location
Streamlines of the developing flow in the flow chamber during a simulated 200 ms linear onset phase of the inflow velocity. After 200 ms, the inflow velocity reaches its final values of 3.5 ml/s. A steady state is reached approximately 300 ms after the initial onset of flow. During the final 100 ms, the recirculation eddy keeps developing, although the flow rate is constant.
Grahic Jump Location
Development of the recirculation eddy in dependency of time. The Reynolds number is ramped linearly over 200 ms (dashed line). After 200 ms, when the inflow velocity is kept constant, the recirculation zone keeps increasing in size. At 300 ms, the position of the second stagnation point lies within 2 percent of its steady-state value, X/S=8.6, which corresponds to 4.8 mm downstream of the expansion site.
Grahic Jump Location
Graphs of the bottom plate shear stress as a function of location and time, τ(x,t). The same function is shown as a three-dimensional graph (A) and a false-color contour plot (B) with yellow and red colors indicating the forward shear stress after the second stagnation point and blue colors indicating the recirculation zone. The final position of the second stagnation point is indicated by the white arrow. The onset phase of flow can easily be recognized by the rise of shear stress and the movement of the second stagnation point: Moving forward in time is equivalent to moving to the right. The sudden-expansion site is represented by the lower edge of the right graph or the right edge moving towards the right background in the left graph. The flow direction is indicated by black arrows. [Subscribers to the JOURNAL OF B IOMECHANICAL E NGINEERING may access Figs. 5(B), 6(B), and 7(B) in color on-line at http://ojps.aip.org/ASMEJournals/Biomechanical.]
Grahic Jump Location
Temporal partial derivative of the function shown in Fig. 5. As in Fig. 5, the same data are presented as a three-dimensional graph and a two-dimensional false-color contour plot. The graphs show that the highest temporal gradient of the shear stress occurs at the time the ramping ends. The system then reaches a steady state after approximately 300 to 350 ms, after which the temporal gradient remains zero. The zone exposed to the highest temporal gradients is located slightly upstream of the second stagnation point.
Grahic Jump Location
Spatial partial derivative of the function shown in Fig. 5. As in Fig. 5, the same data are presented as a three-dimensional graph and a two-dimensional false-color contour plot. It can be seen that the zone with the highest spatial gradients matches the second stagnation point. At the end of the ramping period, a maximum spatial gradient can be observed that is slightly higher than the steady-state gradient. As opposed to the temporal gradients, the spatial gradients do not vanish when the system reaches equilibrium.
Grahic Jump Location
The maximum temporal gradient of shear stress as a function of the ramp time. The data points were computed for the discrete ramp times 50, 100, 200, 500 ms and 1, 2, and 15 s. Due to scaling reasons, the 15 s data point is not shown. A nonlinear curve fit shows that the data points are well represented by the function τ̇max=A⋅tramp−B with A=221.5 and B=1.15(R2=0.999). With tramp→∞ the temporal gradient τ̇ vanishes.
Grahic Jump Location
The maximum spatial gradient of shear stress as a function of the ramp time. The data points were computed for the same discrete ramp times as in Fig. 8. A good fit can be achieved with an exponential decay function, τmax=s⋅exp(−tramp/Kt)+P with the constants s=42.7 kPa/m,Kt=8.6 ms, and P=40.11 kPa/m and a regression coefficient of R=0.996. As opposed to temporal gradients of shear stress, spatial gradients do not vanish with slow ramping times, but rather approach a lower boundary. Also, the increase of spatial gradients with shorter ramp times is by orders of magnitude less pronounced than in the case of temporal gradients.
Grahic Jump Location
Dependency of the maximum shear stress in the recirculation zone (□) and the forward reattachment zone (♦) on the expansion ratio. A high expansion ratio leads to an amplification of the maxima relative to the outflow shear stress.
Grahic Jump Location
Immunofluorescent micrographs of BrdU positive nuclei after 4 h exposure to recirculating flow in the backward-facing step chamber. Shown is a representative example from a paired set of HUVEC slides with cells exposed to a sudden onset of flow in the upper panel, and cells exposed to ramped flow (15 s) below. Each micrograph is a photo-construct that comprises three contiguous HPF. Micrographs span the region of flow recirculation and reattachment (2.2 mm to 9.0 mm downstream from the expansion site). The white dots represent BrdU positive proliferating nuclei. The scaling bar equals 1.1 mm, and the arrow indicates the position of the second stagnation point.
Grahic Jump Location
Cell exposure to temporal gradients compared to cell proliferation rate. Shown is the cell proliferation of cells exposed to shear stress at sudden onset of flow (200 ms, –▴–) relative to the proliferation rate for ramped (15 s) onset of flow (–▪–) 27, lower panel. The difference between cell proliferation rates with sudden and ramped onset of flow are statistically significant within the fields of view from 2.25 to 7.88 mm (shaded area). The total exposure of the cells to temporal shear stress gradients is shown in the middle panel for sudden onset of flow (solid line) versus ramped onset of flow (dashed line). Similarly, the maximum exposure of the cells to spatial gradients is shown in the upper panel for both flow conditions.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In