The local structure of atmospheric turbulence and its effect on the Smagorinsky model for large eddy simulation

M. Chamecki 1,3, C. Meneveau 2,3 & M.B. Parlange 4

1Department of Geography and Environmental Engineering 2Department of Mechanical Engineering 3Center for Environmental and Applied Fluid Mechanics The Johns Hopkins University, Baltimore MD 21218, 4 School of Architecture, Civil and Environmental Engineering, Ecole Polytechnique Federale de Lausanne, Lausanne, Switzerland.

ABSTRACT: Phenomena such as large-scale shear, buoyancy and the proximity to the ground surface significantly affect interactions among scales in atmospheric boundary layer turbulent flows. Hence, these phenomena affect parameters used in subgrid-scale (SGS) parameterizations used in large eddy simulation (LES) of the atmospheric boundary layer. The effects of these phenomena upon SGS parameters have, to date, been studied mostly as function of the global state of the flow. For instance, the Smagorinsky coefficient has been measured as function of the mean shear and stability condition of the atmosphere as determined from the average surface heat and momentum fluxes. However, in LES the global average field values are often difficult to determine a priori and the SGS parameters ideally must be expressed as function of local flow variables that characterize the instantaneous flow phenomena. With the goal of improving the Smagorinsky closure, in this study several dimensionless parameters characterizing the local structure and important dynamical characteristics of the flow are defined. These local parameters include enstrophy, vortex stretching, self-amplification of strain rate, and normalized temperature gradient and all are defined in such a way that they remain bounded under all circumstances. The dependence of the Smagorinsky coefficient on these local parameters is studied \textit{a priori} from field data measured in the atmospheric surface layer and, as a reference point, from direct numerical simulation of neutrally buoyant, isotropic turbulence. To capture the local effects in a statistically meaningful fashion, conditional averaging is used. Results show various important and inter related trends, such as significant increases of the coefficient in regions of large strain-rate self-amplification and vortex stretching. Results also show that the joint dependences on the parameters are rather complicated and can not be described by products of functions that depend on single parameters. Dependences on locally defined parameters are expected to improve the SGS model by sensitizing it to local flow conditions, and by enabling possible generalizations of the dynamic model based on conditional averaging.

(2007), J. Atmos. Sciences 64, 1941-1958.

full pdf article -- (© AMS, see http://www.ametsoc.org/AMS/)