Test of large-eddy simulation models for homogeneous-shear turbulence in stratification

Numerical tests of various subgrid-scale (SGS) models were conducted for turbulence in thermally stratified homogeneous-shear flow at a relatively low Reynolds number. Compared with a direct numerical simulation (DNS), we found that nondynamic isotropic SGS models are not able to represent the energy spectrum very well because the energy decays considerably during the transition between an initial random stage and a stage of coherent turbulent structures. Dynamic models performed well for simulating the energy spectrum and the change of GS properties with time; anisotropy is not a necessary feature under the present simulation conditions, although one of the special features of stratified turbulence is anisotropy. This may be because the present grids for large-eddy simulation were fine enough to resolve the patches of counter-gradient heat fluxes, which play an important role in the evolution of turbulent energy in stratified turbulence. With respect to the domain-averaged values of SGS stresses, only the dynamic two-parameter mixed (DTM) model produced results of the same order of magnitude as those of filtered DNS. This is because of the terms arising from re-decomposing of the SGS stresses in the DTM model. It was also found that this incompetence in simulating the SGS stress is not necessary to simulate GS energy evolution, as is known for wall turbulence.