Highlyresolved, direct numerical simulations of turbulent channel flows with sub
Kolmogorov grid resolution are performed to investigate the characteristics of wallbounded turbulent flows in the presence of sinusoidal wall waviness. The wall waviness serves as a simplified model to study the effects of welldefined geometric parameters of roughness on the characteristics of wallbounded turbulent flows. In this study, a twodimensional 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 wavywall simulations, the Reynolds number based on the friction velocity was kept constant. To study the effects of wave amplitude and wavelength on turbulence, twodimensional 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 smoothwall turbulent channel flow onedimensional pro files of these quantities are computed by averaging them over one wavelength of the wave profile. A strong effect of the wallwaviness 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
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 welldefined geometric features of roughnessamplitude, $a$ , and wavelength, $\lambda$ , on the flow. The flow in the nearwall 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 nearwall 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 nearwall region. The shearlayer 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 shearlayer velocity scale.
more »
« less
 Award ID(s):
 1706346
 NSFPAR ID:
 10317509
 Date Published:
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 937
 ISSN:
 00221120
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this


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 wallnormal location is obtained through the maximization of the power spectral density associated with the linearized Navier–Stokes equations with a turbulent mean profile and deltacorrelated 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 wallnormal 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 nearwall region is influenced by largescale structures further away from the wall. These structures resemble Townsend’s attached eddies in the crossplane, 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 nearwall 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

null (Ed.)Dimensional analysis suggests that the dissipation length scale ( $\ell _{{\it\epsilon}}=u_{\star }^{3}/{\it\epsilon}$ ) is the appropriate scale for the shearproduction range of the secondorder streamwise structure function in neutrally stratified turbulent shear flows near solid boundaries, including smooth and roughwall 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 loglaw dependence of the secondorder streamwise structure function on ( $r/\ell _{{\it\epsilon}}$ ), where $r$ is the streamwise spatial separation. This scaling law is confirmed by largeeddy 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 loglaw 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 shearproduction 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

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:xwiley:00359009:media:qj3839:qj3839math0001 (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

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 nearwall, asymptotic behaviour of turbulent fluxes using a large direct numerical simulation (DNS) database of turbulent channel flow at higher than usual wallnormal resolutions. We vary the Mach number at a constant friction Reynolds number to directly assess compressibility effects. We observe that the nearwall asymptotic behaviour for compressible turbulent flow is different from the corresponding incompressible flow even if the mean density variations are taken into account and semilocal scalings are used. For Mach numbers near the incompressible regimes, the nearwall asymptotic behaviour follows the wellknown theoretical behaviour. When the Mach number is increased, turbulent fluxes containing wallnormal components show a decrease in the slope owing to increased dilatation effects. We observe that $R_{vv}$ approaches its highMachnumber 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