Simulation of Ekman Boundary Layers by Large Eddy Model with Dynamic Mixed Subfilter Closure

Theoretical analysis of boundary layer turbulence has suggested a feasibility of sufficiently accurate turbulence resolving simulations at relatively coarse meshes. However, large eddy simulation (LES) codes, which employ traditional eddy-viscosity turbulence closures, fail to provide adequate turbulence statistics at coarse meshes especially within a surface layer. Manual tuning of parameters in these turbulence closures may correct low order turbulence statistics but severely harms spectra of turbulence kinetic energy (TKE). For more than decade, engineering LES codes successfully employ dynamic turbulence closures. A dynamic Smagorinsky turbulence closure (DSM) has been already tried in environmental LES. The DSM is able to provide adequate turbulence statistics at coarse meshes but it is not completely consistent with the LES equations. This paper investigates applicability of an advanced dynamic mixed turbulence closure (DMM) to simulations of Ekman boundary layers of high Reynolds number flows. The DMM differs from the DSM by explicit calculation of the Leonard term in the turbulence stress tensor. The Horizontal Array Turbulence Study (HATS) field program has revealed that the Leonard term is indeed an important component of the real turbulence stress tensor.