Application of Large-Eddy Simulation to the Study of Pulsatile Flow in a Modeled Arterial Stenosis

[+] Author and Article Information
R. Mittal, S. P. Simmons

Department of Mechanical Engineering, University of Florida, Gainesville, FL 32611

H. S. Udaykumar

Department of Mechanical Engineering, University of Iowa, Iowa City, IA 52242

J Biomech Eng 123(4), 325-332 (Mar 26, 2001) (8 pages) doi:10.1115/1.1385840 History: Received October 22, 2000; Revised March 26, 2001
Copyright © 2001 by ASME
Your Session has timed out. Please sign back in to continue.


Fung, Y. C., 1997, Biomechanics: Circulation, 2nd ed., Springer-Verlag, New York.
Pedley, T. J., 1980, The Fluid Mechanics of Large Blood Vessels, Cambridge University Press, Cambridge.
White, F. M., 1991, Viscous Fluid Flow, McGraw-Hill, New York.
Winter,  D. C., and Nerem,  R. M., 1984, “Turbulence in Pulsatile Flows,” Ann. Biomed. Eng., 12, pp. 357–369.
Davies,  C., and Carpenter,  P. W., 1997, “Numerical Simulation of the Evolution of Tollmien–Schlichting Waves Over Compliant Panels,” J. Fluid Mech., 335 pp. 361–392.
Ku,  D. N., 1997, “Blood Flow in Arteries,” Annu. Rev. Fluid Mech., 29, pp. 399-343; Arteriosclerosis, 5, pp. 289–302.
Nichols, W. W., and O’Rourke, M. F., 1998, McDonald’s Blood Flow in Arteries Theoretical, Experimental and Clinical Principles 4th ed., Oxford University Press, New York.
Roshko, A., 1954, “On the Development of Turbulent Wakes From Vortex Streets,” NACA Report 1191, pp. 124–132.
Bloor,  M. S., 1964, “The Transition to Turbulence in the Wake of a Circular Cylinder,” J. Fluid Mech., 21, p. 290.
Giddens,  D. P., Zarins,  C. K., and Glagov,  S., 1993, “The Role of Fluid Mechanics in the Localization and Detection of Atherosclerosis,” ASME J. Biomech. Eng., 115, pp. 588–594.
Downing,  J. M., and Ku,  D. N., 1997, “Effects of Frictional Losses and Pulsatile Flow on the Collapse of Stenotic Arteries,” ASME J. Biomech. Eng., 119, pp. 317–324.
Lees,  X. X., and Dewey,  Y. Y, 1970, “Phonoangiography: A New Noninvasive Diagnostic Method for Studying Arterial Disease,” Physica A, 67, No. 2, pp. 935–942.
Cloutier,  G, Allard,  L., and Durand,  L.-G., 1996, “Characterization of Blood Flow Turbulence With Pulsed-Wave and Power Doppler Ultrasound Imaging,” ASME J. Biomech. Eng., 118, pp. 318–325.
Rogallo,  R. S., and Moin,  P, 1984, “Numerical Simulation of Turbulent Flows,” Annu. Rev. Fluid Mech., 16, pp. 99–137.
Wilcox, D. C., 1998, Turbulence Modeling for CFD, DCW Industries.
Tutty,  O. R., 1992, “Pulsatile Flow in a Constricted Channel,” ASME J. Biomech. Eng., 114, pp. 50–54.
Germano,  M., Piomelli,  U., Moin,  P., and Cabot,  W., 1991, “A Dynamic Subgrid-Scale Eddy Viscosity Model,” Phys. Fluids, A3, pp. 1760–1765.
Lesieur,  M., and Metais,  O., 1994, “New Trends in Large-Eddy Simulation,” Annu. Rev. Fluid Mech., 28, pp. 45–82.
Mittal,  R., and Moin,  P., 1997, “Suitability of Upwind-Biased Schemes for Large-Eddy Simulation of Turbulent Flows,” AIAA J., 35, No. 8, pp. 1415–1417.
Kaltenbach,  H.-J., Fatica,  M., Mittal,  R., Lund,  T. S., and Moin,  P., 1999, “Study of Flow in Planar Asymmetric Diffuser Using Large-Eddy Simulation,” J. Fluid Mech., 390, pp. 151–186.
Smagorinsky,  J., 1963, “General Circulation Experiments With the Primitive Equations,” Mon. Weather Rev., 91, No. 3, pp. 99–164.
Lilly,  D. K., 1992, “A Proposed Modification of the Germano Subgrid-Scale Closure,” Method. Phys. Fluids, A4, pp. 633–635.
Piomelli,  U., 1999, “Large-Eddy Simulation: Achievements and Challenges,” Prog. Aerosp. Sci., 35, No. 4, pp. 335–362.
Chorin,  A. J., 1968, “Numerical Solution of the Navier–Stokes Equations,” Math. Comput., 22, pp. 745–757.
Holman,  E., 1955, “The Development of Arterial Aneurysms,” J. Surgical Gynecology Obstetrics, 100, pp. 599–611.
Ku, D. N, Giddens, D. P., Zarins, C. K., and Glagov, S., 1985, Pulsatile Flow and Atherosclerosis in the Human Carotid Bifurcation.
Friedman,  M. H., Hutchins,  G. M., Bargeron,  C. B., Deters,  O. J., and Mark,  F. F., 1981, “Correlation Between Intimal Thickness and Fluid Shear in Human Arteries,” Atherosclerosis, 39, pp. 425–436.
Stein,  P. D., Walburn,  F. J., and Sabbah,  H. N., 1982, “Turbulent Stresses in the Region of Aortic and Pulmonary Valves,” ASME J. Biomech. Eng., 104, pp. 238–244.
Nerem,  R. M, 1992, “Vascular Fluid Mechanics, the Arterial Wall and Atherosclerosis,” ASME J. Biomech. Eng., 114, pp. 274–282.
Zarins, C. K., and Glagov, S., 1994, “Pathophysiology of Human Atherosclerosis,” in: F. J. Vieth, R. W. Hobson II, R. A. Williams, et al., eds., Vascular Surgery—Principles and Practice, 2nd ed., McGraw-Hill, New York, pp. 21–39.
Khalifa,  A. M. A., and Giddens,  D. P., 1981, “Characterization and Evolution of Post-Stenotic Flow Disturbances,” J. Biomech., 14, pp. 279–296.
Lieber,  B. B., and Giddens,  D. P., 1990, “Post-stenotic Core Flow Behavior in Pulsatile Flow and Its Effects on Wall Shear Stress,” J. Biomech., 23, No. 6, pp. 597–605.
Reynolds,  W. C., and Hussain,  A. K. M. F., 1972, “The Mechanics of an Organized Wave in Turbulent Shear Flow. Part 3. Theoretical Models and Comparisons With Experiments,” J. Fluid Mech., 54, pp. 263–288.
Khalifa,  A. M. A., and Giddens,  D. P., 1978, “Analysis of Disorder in Pulsatile Flows With Applications to Poststenotic Blood Velocity Measurements in Dogs,” J. Biomech., 11, pp. 129–141.
Tennekes, H., and Lumley, J. L., 1972, A First Course in Turbulence, MIT Press, Cambridge, MA.
Prasad,  A., and Williamson,  C. H. K., 1996, “The Instability of the Separated Shear Layer From a Bluff Body,” Phys. Fluids, 8, No. 6pp. 1347–1349.


Grahic Jump Location
Schematic of the stenosis model employed in the current simulations
Grahic Jump Location
Two-dimensional view of the mesh in the vicinity of the stenosis. Only every second grid point in the x1 and every fourth point in the x2 direction is shown.
Grahic Jump Location
Comparison of mean velocity profiles obtained from two different simulations. The solid and dotted lines correspond to 240×64×16 and 240×96×16 meshes, respectively. The solid circles indicate locations where temporal variation of flow variable is extracted for frequency analysis.
Grahic Jump Location
Sequence of four spanwise-averaged, spanwise vorticity plots over one flow cycle for the Re=2000,St0=0.024 case. Dark and light shades represent clockwise and counterclockwise vorticity, respectively. (a) t/T=0, (b) t/T=0.25, (c) t/T=0.5, (d) t/T=0.75.
Grahic Jump Location
Three-dimensional isosurface plot of spanwise vorticity corresponding to Fig. 4(c)
Grahic Jump Location
Mean flow characteristics obtained by averaging in time and along the spanwise direction: (a) streamline plot; (b) contour plot of pressure
Grahic Jump Location
Nondimensionalized shear stress and pressure on the lower and upper walls of the channel: (a) skin friction coefficient Cf; (b) pressure coefficient Cp
Grahic Jump Location
Root mean square (rms) shear stress and pressure gradient on the lower and upper walls of the channel: (a) root mean square skin friction coefficient (Cf)rms; (b) root-mean-square pressure gradient (pτ)rms
Grahic Jump Location
(a) Temporal variation of streamwise velocity (u1) over one cycle at various locations on the channel centerline. The locations have been indicated by solid circles in Fig. 3. The plots have been offset in the vertical direction. (i) x/H=4.5;offset=0; (ii) x/H=6.1;offset=+2; (iii) x/H=9.1;offset=+6; (iv) x/H=12.3;offset=+8. (b) plot of velocity variation, its phase average (ũ1) and deviation from phase average (u1) at x/H=6.1. (c) Frequency spectra of (u1) corresponding to the variations in Fig. 9. The spectra have been offset by a factor in the vertical direction. (i) x/H=4.5;offset=×1; (ii) x/H=6.1;offset=×10; (iii) x/H=9.1;offset=×100; (iv) x/H=12.3;offset=×2000. Vertical line in the spectra indicates the Strouhal number(=0.94) of the high-frequency vortex shedding.
Grahic Jump Location
(a) Temporal variation of pressure coefficient over one cycle at various streamwise locations on the lower wall at streamline locations corresponding to previous plot. The plots have been offset in the vertical direction. (i) x/H=4.5;offset=0; (ii) x/H=6.1;offset=+4; (iii) x/H=9.1;offset=+8; (iv) x/H=12.3;offset=+12. (b) Frequency spectra of p corresponding to the pressure variations in Fig. 10(a). The spectra have been offset by a factor in the vertical direction. (i) x/H=4.5;offset=×1; (ii) x/H=6.1;offset=×10; (iii) x/H=9.1;offset=×103; (iv) x/H=12.3;offset=×〈5×105). Vertical line in the spectra indicates the Strouhal number (=0.94) of the high-frequency vortex shedding.
Grahic Jump Location
(a) Contour plot of fluctuation kinetic energy; (b) contour plot of averaged subgrid scale (SGS) viscosity. In both plots, dark shades represents higher magnitudes.




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