0
Special Section Articles

# Large Eddy Simulations of a Stratified Shear Layer

[+] Author and Article Information
Hieu T. Pham

Mechanical and Aerospace Engineering,
University of California, San Diego,
La Jolla, CA 92092
e-mail: h8pham@ucsd.edu

Sutanu Sarkar

Fellow ASME
Mechanical and Aerospace Engineering,
University of California,
San Diego,
La Jolla, CA 92092
e-mail: sarkar@ucsd.edu

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received May 2, 2013; final manuscript received December 15, 2013; published online April 28, 2014. Assoc. Editor: Ye Zhou.

J. Fluids Eng 136(6), 060913 (Apr 28, 2014) (11 pages) Paper No: FE-13-1287; doi: 10.1115/1.4026416 History: Received May 02, 2013; Revised December 15, 2013

## Abstract

The performance of the large eddy simulation (LES) approach in predicting the evolution of a shear layer in the presence of stratification is evaluated. The LES uses a dynamic procedure to compute subgrid model coefficients based on filtered velocity and density fields. Two simulations at different Reynolds numbers are simulated on the same computational grid. The fine LES simulated at a low Reynolds number produces excellent agreement with direct numerical simulations (DNS): the linear evolution of momentum thickness and bulk Richardson number followed by an asymptotic approach to constant values is correctly represented and the evolution of the integrated turbulent kinetic energy budget is well captured. The model coefficients computed from the velocity and the density fields are similar and have a value in range of $0.01-0.02$. The coarse LES simulated at a higher Reynolds number Re = 50,000 shows acceptable results in terms of the bulk characteristics of the shear layer, such as momentum thickness and bulk Richardson number. Analysis of the turbulent budgets shows that, while the subgrid stress is able to remove sufficient energy from the resolved velocity fields, the subgrid scalar flux and thereby the subgrid scalar dissipation are underestimated by the model.

<>

## References

Smyth, W. D., and Moum, J. N., 2000, “Length Scales of Turbulence in Stably Stratified Mixing Layers,” Phys. Fluids, 12(6), pp. 1327–1342.
Smyth, W. D., and Moum, J. N., 2000, “Anisotropy of Turbulence in Stably Stratified Mixing Layers,” Phys. Fluids, 12(6), pp. 1343–1362.
Caulfield, C. P., and Peltier, W. R., 2000, “The Anatomy of the Mixing Transition in Homogeneous and Stratified Free Shear Layers,” J. Fluid Mech., 413, pp. 1–47.
Brucker, K., and Sarkar, S., 2007, “Evolution of an Initially Turbulent Stratified Shear Layer,” Phys. Fluids, 19, p. 101105.
Pham, H. T., Sarkar, S., and Brucker, K. A., 2009, “Dynamics of a Stratified Shear Layer Above a Region of Uniform Stratification,” J. Fluid Mech., 630, pp. 191–223.
Smyth, W. D., Moum, J. N., and Caldwell, D. R., 2001, “The Efficiency of Mixing in Turbulent Patches: Inferences From Direct Simulations and Microstructure Observations,” J. Phys. Oceanogr., 31, pp. 1969–1992.
Pham, H. T., and Sarkar, S., 2010, “Transport and Mixing of Density in a Continuously Stratified Shear Layer,” J. Turbulence, 11, pp. 1–23.
Moeng, C., and Sullivan, P., 1994, “A Comparison of Shear- and Buoyancy-Driven Planetary Boundary Layer Flows,” J. Atmos. Sci., 51(7), pp. 999–1022.
Chamecki, M., Meneveau, C., and Parlange, M., 2007, “The Local Structure of Atmospheric Turbulence and Its Effect on the Smagorinsky Model for Large Eddy Simulation,” J. Atmos. Sci., 64, pp. 1941–1958.
Scotti, A., 2010, “Large Eddy Simulation in the Ocean,” Intl. J. Comput. Fluid Dyn., 24, pp. 393–406.
Vreman, B., Geurts, B., and Kuerten, H., 1997, “Large-Eddy Simulation of the Turbulent Mixing Layer,” J. Fluid Mech., 339, pp. 357–390.
Rogers, M. M., and Moser, R. D., 1994, “Direct Simulation of a Self-Similar Turbulent Mixing Layer,” Phys. Fluids, 6, pp. 903–923.
Pantano, C., and Sarkar, S., 2002, “A Study of Compressibility Effects in the High-Speed Turbulent Shear Layer Using Direct Simulation,” J. Fluid Mech., 451, pp. 329–371.
Thorpe, S. A., 1973, “Experiments on Instability and Turbulence in a Stratified Shear Flow,” J. Fluid Mech., 61, pp. 731–751.
Koop, C. G., and Browand, F. K., 1979, “Instability and Turbulence in a Stratified Fluid With Shear,” J. Fluid Mech., 93, pp. 135–159.
Geyer, W. R., Lavery, A. C., Scully, M. E., and Trowbridge, J. H., 2010, “Mixing by Shear Instability at Higher Reynolds Number,” Geophys. Res. Lett., 37(22), pp. 1–5.
Smagorinsky, J., 1963, “General Circulation Experiments With the Primitive Equations I. The Basic Experiment,” Mon. Weather Rev., 91, pp. 99–164.
Germano, M., Piomelli, U., Moin, P., and Cabot, W., 1991, “A Dynamic Subgrid-Scale Eddy Viscosity Model,” Phys. Fluids, 3, pp. 1760–1765.
Bardina, J., Ferziger, J., and Reynolds, W., 1980, “Improved Subgrid Models for Large Eddy Simulation,” AIAA Paper No. 80-1357.
Tejada-Martinez, A. E., and Grosch, C. E., 2007, “Langmuir Turbulence in Shallow Water. Part 2. Large-Eddy Simulation,” J. Fluid Mech., 576, pp. 63–108.
Scotti, A., and Piomelli, U., 2001, “Numerical Simulation of Pulsating Turbulent Channel Flow,” Phys. Fluids, 13, pp. 1367–1384.
Salon, S., Armenio, V., and Crise, A., 2007, “A Numerical Investigation of the Stokes Boundary Layer in the Turbulent Regime,” J. Fluid Mech., 570, pp. 253–296.
Radhakrishnan, S., and Piomelli, U., 2008, “Large-Eddy Simulation of Oscillating Boundary Layers: Model Comparison and Validation,” J. Geophys. Res., 113, p. C02022.
Armenio, V., and Sarkar, S., 2002, “An Investigation of Stably-Stratified Turbulent Channel Flow Using Large Eddy Simulation,” J. Fluid Mech., 459, pp. 1–42.
Taylor, J. R., Sarkar, S., and Armenio, V., 2005, “Large Eddy Simulation of Stably Stratified Open Channel Flow,” Phys. Fluids, 17(116602), pp. 1–18.
Taylor, J. R., and Sarkar, S., 2007, “Internal Gravity Waves Generated by a Turbulent Bottom Ekman Layer,” J. Fluid Mech., 590, pp. 331–354.
Gayen, B., Sarkar, S., and Taylor, J. R., 2010, “Large Eddy Simulation of a Stratified Boundary Layer Under an Oscillatory Current,” J. Fluid Mech., 643, pp. 233–266.
Lilly, D. K., 1992, “A Proposed Modification of the Germano Subgrid-Scale Closure Method,” Phys. Fluids A, 4, pp. 633–635.
Basak, S., and Sarkar, S., 2006, “Dynamics of a Stratified Shear Layer With Horizontal Shear,” J. Fluid Mech., 568, pp. 19–54.
Brucker, K. A., and Sarkar, S., 2010, “A Comparative Study of Self-Propelled and Towed Wakes in a Stratified Fluid,” J. Fluid Mech., 652, pp. 373–404.
Slinn, D. N., and Riley, J. J., 1998, “A Model for the Simulation of Turbulent Boundary Layers in an Incompressible Stratified Flow,” J. Comput. Phys., 144, pp. 550–602.
Le Ribault, C., Sarkar, S., and Stanley, S. A., 1999, “Large Eddy Simulation of a Plane Jet,” Phys. Fluids, 11(10), pp. 3069–3083.

## Figures

Fig. 1

Snapshots of the density field at various times in the DNS of a stratified shear layer [7]. Development of Kelvin–Helmholtz billows in (a) breakdown of the billows into turbulence in (b) and the late-time layer of decaying turbulence in (c) are important characteristics of a turbulent stratified shear layer.

Fig. 2

Evolution of length scales at center of the shear layer in the LES models: (a) Kolmogorov scale η and Ozmidov scale LO and (b) energy-containing scale LEN and Ellison scale LE. Solid lines and dash lines denote the fine (low-Re) and coarse (high-Re) LES models, respectively. The scales are normalized by the horizontal grid spacing Δx=0.12.

Fig. 3

Comparison of wavenumber spectra between the DNS and the LES cases at t=180 and z=0: (a) velocity fields and (b) density fields. The energy spectrum are plotted against the streamwise wavenumber kx.

Fig. 4

Snapshots of the density field at various times in the fine (low-Re) LES. Similar to results of DNS, the LES captures the Kelvin–Helmholtz billows, their breakdown into turbulence, and the decay of turbulence.

Fig. 5

Temporal evolution of (a) the momentum thickness δθ and (b) the bulk Richardson number Rib is similar between the DNS and the fine (low-Re) LES

Fig. 6

Comparison of integrated turbulent kinetic energy budget over the computational domain in (a) the DNS and (b) the fine (low-Re) LES

Fig. 7

Budgets of (a) the turbulent kinetic energy and (b) the density variance 〈ρ'2〉 at t=120 in the fine (low-Re) LES

Fig. 8

(a) Dynamic subgrid coefficient Cd and (b) horizontally averaged eddy viscosity at various times in the fine (low-Re) LES

Fig. 9

(a) Dynamic subgrid coefficient Cθ and (b) horizontally averaged eddy diffusivity at various times in the fine (low-Re) LES

Fig. 10

Snapshots of the density field at various times in the coarse (high-Re) LES. The development of Kelvin–Helmholtz billows and the evolution of turbulence are well captured.

Fig. 11

Comparison of (a) the momentum thickness δθ and (b) the bulk Richardson number Rib between the fine (low-Re) and the coarse (high-Re) LES cases

Fig. 12

Budgets of (a) the turbulent kinetic energy K and (b) the density variance 〈ρ'2〉 at t=120 in the coarse (high-Re) LES

Fig. 13

(a) Difference in the eddy viscosity and eddy diffusivity at z=0 between the fine (low-Re) and coarse (high-Re) LES cases and (b) evolution of squared shear rate S2 and squared buoyancy frequency N2 at the same location

## Discussions

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 Proceedings Articles
Related eBook Content
Topic Collections