 Award ID(s):
 1005363
 Publication Date:
 NSFPAR ID:
 10218845
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 801
 Page Range or eLocationID:
 459 to 474
 ISSN:
 00221120
 Sponsoring Org:
 National Science Foundation
More Like this

We discuss the Onsager theory of wallbounded turbulence, analysing the momentum dissipation anomaly hypothesized by Taylor. Turbulent drag laws observed with both smooth and rough walls imply ultraviolet divergences of velocity gradients. These are eliminated by a coarsegraining operation, filtering out smallscale eddies and windowing out nearwall eddies, thus introducing two arbitrary regularization lengthscales. The regularized equations for resolved eddies correspond to the weak formulation of the Navier–Stokes equation and contain, in addition to the usual turbulent stress, also an inertial drag force modelling momentum exchange with unresolved nearwall eddies. Using an Onsagertype argument based on the principle of renormalization group invariance, we derive an upper bound on wall friction by a function of Reynolds number determined by the modulus of continuity of the velocity at the wall. Our main result is a deterministic version of Prandtl’s relation between the Blasius − 1 / 4 drag law and the 1/7 powerlaw profile of the mean streamwise velocity. At higher Reynolds, the von Kármán–Prandtl drag law requires instead a slow logarithmic approach of velocity to zero at the wall. We discuss briefly also the largeeddy simulation of wallbounded flows and use of iterative renormalization group methods to establish universal statistics inmore »

We employ numerically implicit subgridscale modeling provided by the wellknown streamlined upwind/Petrov–Galerkin stabilization for the finite element discretization of advection–diffusion problems in a Large Eddy Simulation (LES) approach. Whereas its original purpose was to provide sufficient algorithmic dissipation for a stable and convergent numerical method, more recently, it has been utilized as a subgridscale (SGS) model to account for the effect of small scales, unresolvable by the discretization. The freestream Mach number is 2.5, and direct comparison with a DNS database from our research group, as well as with experiments from the literature of adiabatic supersonic spatially turbulent boundary layers, is performed. Turbulent inflow conditions are generated via our dynamic rescaling–recycling approach, recently extended to highspeed flows. Focus is given to the assessment of the resolved Reynolds stresses. In addition, flow visualization is performed to obtain a much better insight into the physics of the flow. A weak compressibility effect is observed on thermal turbulent structures based on twopoint correlations (IC vs. supersonic). The Reynolds analogy (u′ vs. t′) approximately holds for the supersonic regime, but to a lesser extent than previously observed in incompressible (IC) turbulent boundary layers, where temperature was assumed as a passive scalar. A much longermore »

The inertial subrange of turbulent scales is commonly reflected by a power law signature in ensemble statistics such as the energy spectrum and structure functions – both in theory and from observations. Despite promising findings on the topic of fractal geometries in turbulence, there is no accepted image for the physical flow features corresponding to this statistical signature in the inertial subrange. The present study uses boundary layer turbulence measurements to evaluate the selfsimilar geometric properties of velocity isosurfaces and investigate their influence on statistics for the velocity signal. The fractal dimension of streamwise velocity isosurfaces, indicating statistical selfsimilarity in the size of ‘wrinkles’ along each isosurface, is shown to be constant only within the inertial subrange of scales. For the transition between the inertial subrange and production range, it is inferred that the largest wrinkles become increasingly confined by the overall size of largescale coherent velocity regions such as uniform momentum zones. The selfsimilarity of isosurfaces yields powerlaw trends in subsequent onedimensional statistics. For instance, the theoretical 2/3 powerlaw exponent for the structure function can be recovered by considering the collective behaviour of numerous isosurface level sets. The results suggest that the physical presence of inertial subrange eddies ismore »

Turbulent processes in the ocean surface boundary layer (OSBL) play a key role in weather and climate systems. This study explores a Lagrangian analysis of wavedriven OSBL turbulence, based on a largeeddy simulation (LES) model coupled to a Lagrangian stochastic model (LSM). Langmuir turbulence (LT) is captured by Craik–Leibovich wave forcing that generates LT through the Craik–Leibovich type 2 (CL2) mechanism. Breaking wave (BW) effects are modeled by a surface turbulent kinetic energy flux that is constrained by wind energy input to surface waves. Unresolved LES subgridscale (SGS) motions are simulated with the LSM to be energetically consistent with the SGS model of the LES. With LT, Lagrangian autocorrelations of velocities reveal three distinct turbulent time scales: an integral, a dispersive mixing, and a coherent structure time. Coherent structures due to LT result in relatively narrow peaks of Lagrangian frequency velocity spectra. With and without waves, the highfrequency spectral tail is consistent with expectations for the inertial subrange, but BWs substantially increase spectral levels at high frequencies. Consistently, over short times, particlepair dispersion results agree with the Richardson–Obukhov law, and nearsurface dispersion is significantly enhanced because of BWs. Over longer times, our dispersion results are consistent with Taylor dispersion. Inmore »

Abstract Direct numerical simulations are performed to compare the evolution of turbulent stratified shear layers with different density gradient profiles at a high Reynolds number. The density profiles include uniform stratification, twolayer hyperbolic tangent profile and a composite of these two profiles. All profiles have the same initial bulk Richardson number ( $$Ri_{b,0}$$ R i b , 0 ); however, the minimum gradient Richardson number and the distribution of density gradient across the shear layer are varied among the cases. The objective of the study is to provide a comparative analysis of the evolution of the shear layers in term of shear layer growth, turbulent kinetic energy as well as the mixing efficiency and its parameterization. The evolution of the shear layers in all cases shows the development of Kelvin–Helmholtz billows, the transition to turbulence by secondary instabilities followed by the decay of turbulence. Comparison among the cases reveals that the amount of turbulent mixing varies with the density gradient distribution inside the shear layer. The minimum gradient Richardson number and the initial bulk Richardson number do not correlate well with the integrated TKE production, dissipation and buoyancy flux. The bulk mixing efficiency for fixed $$Ri_{b,0}$$ R i b ,more »