On the Subgrid-Scale Modeling of Compressible Turbulence

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented

On the large-eddy simulation of transitional wall-bounded flows

The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed

A framework for the evaluation of turbulence closures used in mesoscale ocean large-eddy simulations

We present a methodology to determine the best turbulence closure for an eddy-permitting ocean model through measurement of the error-landscape of the closure's subgrid spectral transfers and flux. We apply this method to 6 different closures for forced-dissipative simulations of the barotropic vorticity equation on a f-plane (2D Navier-Stokes equation). Using a high-resolution benchmark, we compare each closure's model of energy and enstrophy transfer to the actual transfer observed in the benchmark run. The error-landscape norms enable us to both make objective comparisons between the closures and to optimize each closure's free parameter for a fair comparison. The hyper-viscous closure most closely reproduces the enstrophy cascade, especially at larger scales due to the concentration of its dissipative effects to the very smallest scales. The viscous and Leith closures perform nearly as well, especially at smaller scales where all three models were dissipative. The Smagorinsky closure dissipates enstrophy at the wrong scales. The anticipated potential vorticity closure was the only model to reproduce the upscale transfer of kinetic energy from the unresolved scales, but would require high-order Laplacian corrections in order to concentrate dissipation at the smallest scales. The Lagrangian-averaged alpha-model closure did not perform successfully for forced 2D isotropic Navier-Stokes: small-scale filamentation is only slightly reduced by the model while small-scale roll-up is prevented. Together, this reduces the effects of diffusion.Comment: 44 pages, 21 figures, 1 Appendix, submitted to Ocean Modelin