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
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):
- 1706346
- NSF-PAR ID:
- 10317509
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 937
- ISSN:
- 0022-1120
- 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 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
-
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
-
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
-
more » « less
Abstract High fidelity near-wall velocity measurements in wall bounded fluid flows continue to pose a challenge and the resulting limitations on available experimental data cloud our understanding of the near-wall velocity behavior in turbulent boundary layers. One of the challenges is the spatial averaging and limited spatial resolution inherent to cross-correlation-based particle image velocimetry (PIV) methods. To circumvent this difficulty, we implement an explicit no-slip boundary condition in a wavelet-based optical flow velocimetry (wOFV) method. It is found that the no-slip boundary condition on the velocity field can be implemented in wOFV by transforming the constraint to the wavelet domain through a series of algebraic linear transformations, which are formulated in terms of the known wavelet filter matrices, and then satisfying the resulting constraint on the wavelet coefficients using constrained optimization for the optical flow functional minimization. The developed method is then used to study the classical problem of a turbulent channel flow using synthetic data from a direct numerical simulation (DNS) and experimental particle image data from a zero pressure gradient, high Reynolds number turbulent boundary layer. The results obtained by successfully implementing the no-slip boundary condition are compared to velocity measurements from wOFV without the no-slip condition and to a commercial PIV code, using the velocity from the DNS as ground truth. It is found that wOFV with the no-slip condition successfully resolves the near-wall profile with enhanced accuracy compared to the other velocimetry methods, as well as other derived quantities such as wall shear and turbulent intensity, without sacrificing accuracy away from the wall, leading to state of the art measurements in the
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2022.118
