skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: Amplitude and wavelength scaling of sinusoidal roughness effects in turbulent channel flow at fixed
Direct numerical simulations are performed for incompressible, turbulent channel flow over a smooth wall and different sinusoidal wall roughness configurations at a constant $Re_\tau = 720$ . Sinusoidal walls are used to study the effects of well-defined geometric features of roughness-amplitude, $a$ , and wavelength, $\lambda$ , on the flow. The flow in the near-wall region is strongly influenced by both $a$ and $\lambda$ . Establishing appropriate scaling laws will aid in understanding the effects of roughness and identifying the relevant physical mechanisms. Using inner variables and the roughness function to scale the flow quantities provides support for Townsend's hypothesis, but inner scaling is unable to capture the flow physics in the near-wall region. We provide modified scaling relations considering the dynamics of the shear layer and its interaction with the roughness. Although not a particularly surprising observation, this study provides clear evidence of the dependence of flow features on both $a$ and $\lambda$ . With these relations, we are able to collapse and/or align peaks for some flow quantities and, thus, capture the effects of surface roughness on turbulent flows even in the near-wall region. The shear-layer scaling supports the hypothesis that the physical mechanisms responsible for turbulent kinetic energy production in turbulent flows over rough walls are greatly influenced by the shear layer and its interaction with the roughness elements. Finally, a semiempirical model is developed to predict the contribution of pressure and skin friction drag on the roughness element based purely on its geometric parameters and the corresponding shear-layer velocity scale.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Highly-resolved, direct numerical simulations of turbulent channel flows with sub- Kolmogorov grid resolution are performed to investigate the characteristics of wall-bounded turbulent flows in the presence of sinusoidal wall waviness. The wall waviness serves as a simplified model to study the effects of well-defined geometric parameters of roughness on the characteristics of wall-bounded turbulent flows. In this study, a two-dimensional wave profile with steepness ranging from 0.06 to 0.25 and wave amplitudes ranging from 9 to 36 wall units were considered. For the smooth and wavy-wall simulations, the Reynolds number based on the friction velocity was kept constant. To study the effects of wave amplitude and wavelength on turbulence, two-dimensional time and spanwise averaged distributions of the mean flow, turbulent kinetic energy, and Reynolds stresses as well as turbulent kinetic energy production and dissipation are examined. Furthermore, in order to provide a more direct comparison with the smooth-wall turbulent channel flow one-dimensional pro- files of these quantities are computed by averaging them over one wavelength of the wave profile. A strong effect of the wall-waviness and, in particular, the wave amplitude and wavelength on the characteristics of the turbulence was obtained. Wall waviness mainly affected the inner flow region while all recorded turbulent statistics collapsed in the outer flow region. Significant reductions in turbulent kinetic energy, production and dissipation were obtained with increasing wave amplitudes when reported in inner scale. While production is lower for all wavy wall cases considered here in comparison to the smooth wall, reducing the wavelength caused an increase in production and a decrease in dissipation. 
    more » « less
  2. This paper employs an input–output based approach to analyse the convective velocities and transport of fluctuations in turbulent channel flows. The convective velocity for a fluctuating quantity associated with streamwise–spanwise wavelength pairs at each wall-normal location is obtained through the maximization of the power spectral density associated with the linearized Navier–Stokes equations with a turbulent mean profile and delta-correlated Gaussian forcing. We first demonstrate that the mean convective velocities computed in this manner agree well with those reported previously in the literature. We then exploit the analytical framework to probe the underlying mechanisms contributing to the local convective velocity at different wall-normal locations by isolating the contributions of each streamwise–spanwise wavelength pair (flow scale). The resulting analysis suggests that the behaviour of the convective velocity in the near-wall region is influenced by large-scale structures further away from the wall. These structures resemble Townsend’s attached eddies in the cross-plane, yet show incomplete similarity in the streamwise direction. We then investigate the role of each linear term in the momentum equation to isolate the contribution of the pressure, mean shear, and viscous effects to the deviation of the convective velocity from the mean at each flow scale. Our analysis highlights the role of the viscous effects, particularly in regards to large channel spanning structures whose influence extends to the near-wall region. The results of this work suggest the promise of an input–output approach for analysing convective velocity across a range of flow scales using only the mean velocity profile. 
    more » « less
  3. null (Ed.)
    Dimensional analysis suggests that the dissipation length scale ( $\ell _{{\it\epsilon}}=u_{\star }^{3}/{\it\epsilon}$ ) is the appropriate scale for the shear-production range of the second-order streamwise structure function in neutrally stratified turbulent shear flows near solid boundaries, including smooth- and rough-wall boundary layers and shear layers above canopies (e.g. crops, forests and cities). These flows have two major characteristics in common: (i) a single velocity scale, i.e. the friction velocity ( $u_{\star }$ ) and (ii) the presence of large eddies that scale with an external length scale much larger than the local integral length scale. No assumptions are made about the local integral scale, which is shown to be proportional to $\ell _{{\it\epsilon}}$ for the scaling analysis to be consistent with Kolmogorov’s result for the inertial subrange. Here ${\it\epsilon}$ is the rate of dissipation of turbulent kinetic energy (TKE) that represents the rate of energy cascade in the inertial subrange. The scaling yields a log-law dependence of the second-order streamwise structure function on ( $r/\ell _{{\it\epsilon}}$ ), where $r$ is the streamwise spatial separation. This scaling law is confirmed by large-eddy simulation (LES) results in the roughness sublayer above a model canopy, where the imbalance between local production and dissipation of TKE is much greater than in the inertial layer of wall turbulence and the local integral scale is affected by two external length scales. Parameters estimated for the log-law dependence on ( $r/\ell _{{\it\epsilon}}$ ) are in reasonable agreement with those reported for the inertial layer of wall turbulence. This leads to two important conclusions. Firstly, the validity of the $\ell _{{\it\epsilon}}$ -scaling is extended to shear flows with a much greater imbalance between production and dissipation, indicating possible universality of the shear-production range in flows near solid boundaries. Secondly, from a modelling perspective, $\ell _{{\it\epsilon}}$ is the appropriate scale to characterize turbulence in shear flows with multiple externally imposed length scales. 
    more » « less
  4. Large Eddy Simulations (LES) of neutral flow over regular arrays of cuboids are conducted to explore connections between momentum (z 0m ) and scalar (z 0s ) roughness lengths in urban environments, and how they are influenced by surface geometry. As LES resolves the obstacles but not the micro‐scale boundary layers attached to them, the aforementioned roughness lengths are analyzed at two distinct spatial scales. At the micro‐scale (roughness of individual facets, e.g. roofs), it is assumed that both momentum and scalar transfer are governed by accepted arguments for smooth walls that form the basis for the LES wall model. At the macro‐scale, the roughness lengths are representative of the aggregate effects of momentum and scalar transfer over the resolved roughness elements of the whole surface, and hence they are directly computed from the LES. The results indicate that morphologically‐based parameterizations for macro‐scale z 0m are adequate overall. The relation between the momentum and scalar macro‐roughness values, as conventionally represented by log(z 0m /z 0s ) and assumed to scale with urn:x-wiley:00359009:media:qj3839:qj3839-math-0001 (where Re * is a roughness Reynolds number), is then interpreted using surface renewal theory (SRT). SRT predicts n = 1/4 when only Kolmogorov‐scale eddies dominate the scalar exchange, whereas n = 1/2 is predicted when large eddies limit the renewal dynamics. The latter is found to better capture the LES results. However, both scaling relations indicate that z 0s decreases when z 0m increases for typical urban geometries and scales. This is opposite to how their relation is usually modeled for urban canopies (i.e. z 0s /z 0m is a fixed value smaller than unity). 
    more » « less
  5. The asymptotic behaviour of Reynolds stresses close to walls is well established in incompressible flows owing to the constraint imposed by the solenoidal nature of the velocity field. For compressible flows, thus, one may expect a different asymptotic behaviour, which has indeed been noted in the literature. However, the transition from incompressible to compressible scaling, as well as the limiting behaviour for the latter, is largely unknown. Thus, we investigate the effects of compressibility on the near-wall, asymptotic behaviour of turbulent fluxes using a large direct numerical simulation (DNS) database of turbulent channel flow at higher than usual wall-normal resolutions. We vary the Mach number at a constant friction Reynolds number to directly assess compressibility effects. We observe that the near-wall asymptotic behaviour for compressible turbulent flow is different from the corresponding incompressible flow even if the mean density variations are taken into account and semi-local scalings are used. For Mach numbers near the incompressible regimes, the near-wall asymptotic behaviour follows the well-known theoretical behaviour. When the Mach number is increased, turbulent fluxes containing wall-normal components show a decrease in the slope owing to increased dilatation effects. We observe that $R_{vv}$ approaches its high-Mach-number asymptote at a lower Mach number than that required for the other fluxes. We also introduce a transition distance from the wall at which turbulent fluxes exhibit a change in scaling exponents. Implications for wall models are briefly presented. 
    more » « less