An eddy-viscosity model is proposed and applied in large-eddy simulation of turbulent shear flows with quite satisfactory results. The model is essentially not more complicated than the Smagorinsky model, but is constructed in such a way that its dissipation is relatively small in transitional and near-wall regions. The model is expressed in first-order derivatives, does not involve explicit filtering, averaging, or clipping procedures, and is rotationally invariant for isotropic filter widths. Because of these highly desirable properties the model seems to be well suited for engineering applications. In order to provide a foundation of the model, an algebraic framework for general three-dimensional flows is introduced. Within this framework several types of flows are proven to have zero energy transfer to subgrid scales. The eddy viscosity is zero in the same cases; the theoretical subgrid dissipation and the eddy viscosity have the same algebraic structure. In addition, the model is based on a fundamental realizability inequality for the theoretical subgrid dissipation. Results are shown for a transitional and turbulent mixing layer at high Reynolds number and a turbulent channel flow. In both cases the present model is found to be more accurate than the Smagorinsky model and as good as the standard dynamic model. Unlike the Smagorinsky model, the present model is able to adequately handle not only turbulent but also transitional flow.

Direct numerical simulation databases have been used to study the effect of compressibility on mixing layers. The simulations cover convective Mach numbers from 0.2 to 1.2 and all contain a fully resolved turbulent energy cascade to small spatial scales. Statistical information is extracted from the databases to determine reasons for the reduced growth rate that is observed as the convective Mach number is increased. It is found that the dilatational contribution to dissipation is negligible even when eddy shocklets are observed in the flow. Also pressure-dilatation is not found to be significant. Using an accurate relation between the momentum thickness growth rate and the production of turbulence kinetic energy together with integrated equations for the Reynolds stress tensor it is shown that reduced pressure fluctuations are responsible for the changes in growth rate via the pressure-strain term. A deterministic model for the required pressure fluctuations is given based on the structure of variable-density vortices and the assumption that the limiting eddies are sonic. Simple anisotropy considerations are used to close the averaged equations. Good agreement with turbulence statistics obtained from the simulations is found.

DOI to the publisher's website. • The final author version and the galley proof are versions of the publication after peer review. • The final published version features the final layout of the paper including the volume, issue and page numbers. Link to publication General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal. If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the "Taverne" license above, please follow below link for the End User Agreement:

No abstract

Direct numerical simulation (DNS) and large-eddy simulation (LES) of particle-laden turbulent channel flow, in which the particles experience a drag force, are performed. In this flow turbophoresis leads to an accumulation of particles near the walls. It is shown that the turbophoresis in LES is reduced, in case the subgrid effects in the particle equations of motion are ignored. To alleviate this problem an inverse filtering model is proposed and incorporated into the particle equations. The model is shown to enhance the turbophoresis in actual LES, such that a good agreement with the DNS prediction is obtained.