Libration driven elliptical instability Phys. Fluids 24, 061703 (2012) On integral length scales in anisotropic turbulence Phys. Fluids 24, 061702 (2012) An experimental study of small Atwood number Rayleigh-Taylor instability using the magnetic levitation of paramagnetic fluids Phys. Fluids 24, 052106 (2012) Direct measurement of the spectral transfer function of a laser based anemometer Rev. Sci. Instrum. 83, 033111 (2012) Additional information on Phys. Fluids A One major drawback of the eddy viscosity subgrid-scale stress models used in large-eddy simulations is their inability to represent correctly with a single universal constant different turbulent fields in rotating or sheared flows, near solid walls, or in transitional regimes. In the present work a new eddy viscosity model is presented which alleviates many of these drawbacks. The model coefficient is computed dynamically as the calculation progresses rather than input apriori. The model is based on an algebraic identity between the subgrid-scale stresses at two different filtered levels and the resolved turbulent stresses. The subgrid-scale stresses obtained using the proposed model vanish in laminar flow and at a solid boundary, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The results of large-eddy simulations of transitional and turbulent channel flow that use the proposed model are in good agreement with the direct simulation data.

The dynamic subgrid-scale (SGS) model of Germano et al. [Phys. Fluids A 3, 1760 (1991)] is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

No abstract

No abstract

Most subgrid-scale (SGS) models for large-eddy simulations (LES) are absolutely dissipative (that is, they remove energy from the large scales at each point in the physical space). The actual SGS stresses, however, may transfer energy to the large scales (backscatter) at a given location. Recent work on the LES of transitional flows [ Piomelli et al., Phys. Fluids A 2,257 (1990) ] has shown that failure to account for this phenomenon can cause inaccurate prediction of the growth of the perturbations. Direct numerical simulations of transitional and turbulent channel flow and compressible isotropic turbulence are used to study the backscatter phenomenon. In all flows considered roughly 50% of the grid points were experiencing backscatter when a Fourier cutoff filter was used. The backscatter fraction was less with a Gaussian filter, and intermediate with a box filter in physical space. Moreover, the backscatter and forward scatter contributions to the SGS dissipation were comparable, and each was often much larger than the total SGS dissipation. The SGS dissipation (normalized by total dissipation) increased with filter width almost independently of filter type. The amount of backscatter showed an increasing trend with Reynolds number. In the near-wall region of the channel, events characterized by strong Reynolds shear stress correlated fairly well with areas of high SGS dissipation (both forward and backward). In compressible isotropic turbulence similar results were obtained, independent of fl uctuation Mach number.

a b s t r a c tFlows in which shock waves and turbulence are present and interact dynamically occur in a wide range of applications, including inertial confinement fusion, supernovae explosion, and scramjet propulsion. Accurate simulations of such problems are challenging because of the contradictory requirements of numerical methods used to simulate turbulence, which must minimize any numerical dissipation that would otherwise overwhelm the small scales, and shock-capturing schemes, which introduce numerical dissipation to stabilize the solution. The objective of the present work is to evaluate the performance of several numerical methods capable of simultaneously handling turbulence and shock waves. A comprehensive range of high-resolution methods (WENO, hybrid WENO/central difference, artificial diffusivity, adaptive characteristic-based filter, and shock fitting) and suite of test cases (Taylor-Green vortex, Shu-Osher problem, shock-vorticity/entropy wave interaction, Noh problem, compressible isotropic turbulence) relevant to problems with shocks and turbulence are considered. The results indicate that the WENO methods provide sharp shock profiles, but overwhelm the physical dissipation. The hybrid method is minimally dissipative and leads to sharp shocks and well-resolved broadband turbulence, but relies on an appropriate shock sensor. Artificial diffusivity methods in which the artificial bulk viscosity is based on the magnitude of the strain-rate tensor resolve vortical structures well but damp dilatational modes in compressible turbulence; dilatation-based artificial bulk viscosity methods significantly improve this behavior. For well-defined shocks, the shock fitting approach yields good results.

A mixing transition has been observed in a high-resolution large-eddy simulation of Rayleigh-Taylor instability. During the transition, an inertial range forms in the velocity spectrum and the rate of growth of the mixing zone is reduced. By measuring growth of the layer in units of dominant initial wavelength, criteria are established for reaching the hypothetical self-similar state of the mixing layer. A model including mixing effects is derived for the growth rate. The rate of growth of the mixing layer is determined by the net mass flux through the plane associated with the initial location of the interface. All of the information necessary for predicting the growth rate is contained in this plane. The model incorporates an effective Atwood number and provides a good match to the simulation data. The model suggests reduced growth for miscible fluids, compared to the immiscible case.

No abstract