By comparing the budget of a data-driven quasi-linear approximation (DQLA) (Holford, Lee & Hwang, J. Fluid Mech., vol. 980, 2024, A12) and direct numerical simulation (DNS) (Lee & Moser, J. Fluid Mech., vol. 860, 2019, pp. 886–938), the energetics of linear models for wall-bounded turbulence are assessed. The DQLA is implemented with the linearised Navier–Stokes equations with a stochastic forcing term and an eddy viscosity diffusion model. The self-consistent nature of the DQLA allows for a global comparison across all wavenumbers to assess the role of the various terms in the linear model in replicating the features present in DNS. Starting from the steady-state second-order statistics of a Fourier mode, a spectral budget equation is derived, connecting Lyapunov-like equations to the transport budget equations obtained from DNS. It is found that the DQLA and DNS are in good qualitative agreement for the streamwise-elongated structures present in DNS, comparing well for production, viscous transport and wall-normal turbulent transport. However, the DQLA does not have an energy-conservative nonlinear term. This results in no dissipation under molecular viscosity, with energy instead being dissipated locally through the eddy viscosity model, which models the energy removal by the nonlinear term at integral length scales. Comparison of the pressure–strain statistics also highlights the absence of the streak instability, with production and forcing mainly being retained in the streamwise and wall-normal components or shifted to the spanwise component. It is demonstrated that the eddy viscosity diffusion term locally enforces a self-similar budget, making the model for the nonlinear term self-consistent with a logarithmic mean profile. Implications and recommendations to improve the current eddy viscosity enhanced linear models are also discussed concerning the comparison with DNS, as well as considerations with regard to pressure statistics to mimic the role of the streak instability through colour of turbulence models.

This DNS study considers transition to turbulence in plane Couette flow (pCf) with a rough stationary wall and a smooth moving wall. The roughness elements are square ribs of height $k=0.2h$ (where $h$ is the half-channel height). Two different pitch separations, $\lambda =2k$ and $10k$, are considered, i.e. d-type and k-type roughness, respectively. The transition in both rough pCf cases takes place through a stage of alternate laminar–turbulent bands aligned in an oblique fashion. However, roughness causes a shift in the transitional Reynolds number ($Re$) range. In the k-type roughness, stable bands are observed in the range $Re \in [300, 325]$, which is a downwards shift from the transitional $Re$ range for the smooth pCf ($Re \in [325,400]$). The d-type roughness, on the other hand, surprisingly shifts the transitional $Re$ range upwards to $Re \in [350,425]$. This peculiar behaviour is associated with the ability of the ribs to shed and regenerate vorticity. Large-scale components extracted using a filtering process relate to the transition bands and flow parallel to the oblique laminar–turbulent boundaries. Counter-rotating vortices are present in the turbulent regions of the flow field, which exist in tandem with the high- and low-velocity streaks. Another interesting observation is the secondary Reynolds shear stresses, $-\overline {v^{\prime }u^{\prime }}$ and $-\overline {w^{\prime }v^{\prime }}$, which are non-zero in the transitional regime, in contrast to the turbulent regime where they are negligible.

We study the effect of surface texture on an overlying turbulent flow for the case of textures made of an alternating slip/no-slip pattern, a common model for superhydrophobic surfaces, but also a particularly simple form of texture. For texture sizes $L^+ \gtrsim 25$, we have previously reported that, even though the texture effectively imposes homogeneous slip boundary conditions on the overlying, background turbulence, this is not its sole effect. The effective conditions only produce an origin offset on the background turbulence, which remains otherwise smooth-wall-like. For actual textures, however, as their size increases from $L^+ \gtrsim 25$ the flow progressively departs from this smooth-wall-like regime, resulting in additional shear Reynolds stress and increased drag, in a non-homogeneous fashion that could not be reproduced by the effective boundary conditions. This paper focuses on the underlying physical mechanism of this phenomenon. We argue that it is caused by the nonlinear interaction of the texture-coherent flow, directly induced by the surface topology, and the background turbulence, as it acts directly on the latter and alters it. This does not occur at the boundary where effective conditions are imposed, but within the overlying flow itself, where the interaction acts as a forcing on the governing equations of the background turbulence, and takes the form of cross-advective terms between the latter and the texture-coherent flow. We show this by conducting simulations where we remove the texture and introduce additional, forcing terms in the Navier–Stokes equations, in addition to the equivalent homogeneous slip boundary conditions. The forcing terms capture the effect of the nonlinear interaction on the background turbulence without the need to resolve the surface texture. We show that, when the forcing terms are derived accounting for the amplitude modulation of the texture-coherent flow by the background turbulence, they quantitatively capture the changes in the flow up to texture sizes $L^+ \approx 70{-}100$. This includes not just the roughness function but also the changes in the flow statistics and structure.

Modular floating solar farms exhibit periodic open surface coverages due to the strip configuration of floating modules that support the photovoltaic (PV) panels on top. The associated modulations in the surface boundary layer and its turbulence characteristics are investigated in the present study under fully developed open channel flows. Different coverage percentages of 100 % (i.e. full cover), 60 %, 30 % and 0 % (i.e. open surface) were tested and measurements were obtained using particle image velocimetry. The results showed that the turbulence statistics are similar when the coverage decreases from 100 % to 60 %. However, with 30 %, both the turbulence intensities and Reynolds stresses increase substantially, reaching up to 50 % higher compared with the 100 % coverage, and the boundary layer thickness increases by more than 25 %. The local skin friction beneath the openings increases by 50 %. Analysis of spanwise vortices and premultiplied spectra indicates that the periodic coverage elongates the hairpin vortex packets and reduces their inclination angle, imposing limitations on sustainable coherent structures. At 30 %, flow detachment and smaller-scale vortices become dominant, reducing the mean velocities and increasing the turbulence intensities. Decreasing coverage percentage with flow detachment also shifts the energy transfer to higher wavenumbers, increasing energy dissipation and decreasing the bulk flow velocity. The kinetic energy and Reynolds stress carried by very large-scale motions decreases from 40 %–50 % with the 100 % and 60 % coverage to around 30 %–40 % with the 30 % coverage. Further research studies involving spanwise heterogeneity, higher Reynolds number and varying submergence of PV modules are needed for environmental considerations.

The results of an experimental investigation of smooth-body adverse pressure gradient (APG) turbulent boundary layer flow separation and reattachment over a two-dimensional ramp are presented. These results are part of a larger archival smooth-body flow separation data set acquired in partnership with NASA Langley Research Center and archived on the NASA Turbulence Modeling Resource website. The experimental geometry provides initial canonical turbulent boundary layer growth under nominally zero pressure gradient conditions prior to encountering a smooth, two-dimensional, backward facing ramp geometry onto which a streamwise APG that is fully adjustable is imposed. Detailed surface and off-surface flow field measurements are used to fully characterize the smooth-body APG turbulent boundary layer separation and reattachment at multiple spanwise locations over the ramp geometry. Unsteady aspects of the flow separation are characterized. It is shown that the first and second spatial derivatives of the streamwise static surface pressure profile are sufficient to determine key detachment and reattachment locations. The imposed streamwise APG gives rise to inflectional mean velocity profiles and the associated formation of an embedded shear layer, which is shown to play a dominant role in the subsequent flow development. Similarity scaling is developed for both the mean velocity and turbulent stresses that is found to provide self-similar collapse of profiles for different regions of the ramp flow. Despite the highly non-equilibrium flow environment, a new similarity scaling proved capable of providing self-similar turbulent stress profiles over the full streamwise extent of flow separation and downstream reattachment.

The causal relevance of local flow conditions in open-channel turbulence is analysed using ensembles of interventional experiments in which the effect of perturbing the flow within a small cell is monitored at some future time. When this is done using the relative amplification of the perturbation energy, causality depends on the flow conditions within the cell before it is perturbed, and can be used as a probe of the flow dynamics. The key scaling parameter is the ambient shear, which is also the dominant diagnostic variable for wall-attached perturbations. Away from the wall, the relevant variables are the streamwise and wall-normal velocities. Causally significant cells are associated with sweeps that carry the perturbation towards the stronger shear near the wall, whereas irrelevant ones are associated with ejections that carry it towards the weaker shear in the outer layers. Causally significant and irrelevant cells are themselves organised into structures that share many characteristics with classical sweeps and ejections, such as forming spanwise pairs whose dimensions and geometry are similar to those of classical quadrants. At the wall, this is consistent with causally significant configurations in which a high-speed streak overtakes a low-speed one, and causally irrelevant ones in which the two streaks pull apart from each other. It is argued that this is probably associated with streak meandering.

Wind tunnel measurements of the incident turbulent velocity fields and axial forces on a horizontal axis turbine and porous disc analogues are reported. The models were tested in both a simulated atmospheric boundary layer (ABL) and in grid turbulence, allowing for a range of turbulence length scale to rotor diameter ratios to be considered. A theoretical framework to account for the combined effect of distortion and potential flow blocking in the induction zone is presented. In the case of very large length-scale turbulence to diameter ratios, where distortion effects are minimal, a quasi-steady approach is adopted for the effect of blocking. For the small length-scale ratio limit, the method is developed from the classical analyses for rapid distortion of turbulence and blockage from flow through a porous sheet of resistance. For general length-scale ratios, an efficient prediction method based on interpolation between the two length-scale ratio extremes is established. For very large length-scale ratios, a quasi-steady theory without distortion is appropriate for a rotor or disc in a simulated ABL. The small length-scale theory is applicable for tests conducted in grid turbulence. The results of the study can inform the prediction and interpretation of typical measurements of turbulence within the induction zone and the fluctuating loads on a rotor, at both prototype and full scale. This is of particular importance to fatigue load assessments.

Resolvent analysis provides a framework to predict coherent spatio-temporal structures of the largest linear energy amplification, through a singular value decomposition (SVD) of the resolvent operator, obtained by linearising the Navier–Stokes equations about a known turbulent mean velocity profile. Resolvent analysis utilizes a Fourier decomposition in time, which has thus far limited its application to statistically stationary or time-periodic flows. This work develops a variant of resolvent analysis applicable to time-evolving flows, and proposes a variant that identifies spatio-temporally sparse structures, applicable to either stationary or time-varying mean velocity profiles. Spatio-temporal resolvent analysis is formulated through the incorporation of the temporal dimension to the numerical domain via a discrete time-differentiation operator. Sparsity (which manifests in localisation) is achieved through the addition of an $l_1$-norm penalisation term to the optimisation associated with the SVD. This modified optimisation problem can be formulated as a nonlinear eigenproblem and solved via an inverse power method. We first showcase the implementation of the sparse analysis on a statistically stationary turbulent channel flow, and demonstrate that the sparse variant can identify aspects of the physics not directly evident from standard resolvent analysis. This is followed by applying the sparse space–time formulation on systems that are time varying: a time-periodic turbulent Stokes boundary layer and then a turbulent channel flow with a sudden imposition of a lateral pressure gradient, with the original streamwise pressure gradient unchanged. We present results demonstrating how the sparsity-promoting variant can either change the quantitative structure of the leading space–time modes to increase their sparsity, or identify entirely different linear amplification mechanisms compared with non-sparse resolvent analysis.

We report direct numerical simulations results of the rough-wall channel, focusing on roughness with high $k_{rms}/k_a$ statistics but small to negative $Sk$ statistics, and we study the implications of this new dataset on rough-wall modelling. Here, $k_{rms}$ is the root mean square, $k_a$ is the first-order moment of roughness height, and $Sk$ is the skewness. The effects of packing density, skewness and arrangement of roughness elements on mean streamwise velocity, equivalent roughness height ($z_0$) and Reynolds and dispersive stresses have been studied. We demonstrate that two-point correlation lengths of roughness height statistics play an important role in characterizing rough surfaces with identical moments of roughness height but different arrangements of roughness elements. Analysis of the present as well as historical data suggests that the task of rough-wall modelling is to identify geometric parameters that distinguish the rough surfaces within the calibration dataset. We demonstrate a novel feature selection procedure to determine these parameters. Further, since there is no finite set of roughness statistics that distinguish between all rough surfaces, we argue that obtaining a universal rough-wall model for making equivalent sand-grain roughness ($k_s$) predictions would be challenging, and that each rough-wall model would have its applicable range. This motivates the development of group-based rough-wall models. The applicability of multi-variate polynomial regression and feedforward neural networks for building such group-based rough-wall models using the selected features has been shown.

In hypersonic flight the shock wave and turbulent boundary layer interaction (STBLI) sharply increases wall heat transfer that intensifies the aerodynamic heating problems. In this work the STBLI is modelled by compression ramp flow with a Mach number of 5, a Reynolds number based on momentum thickness of 4652 and a wall to recovery temperature ratio of 0.5. The aerodynamic heat generation and transport mechanisms are investigated in the interaction based on theoretical analysis and direct numerical simulation (DNS) that agrees with previous studies. A prediction correlation of wall heat flux in STBLI is deduced theoretically and validated by some representative data including the present DNS, which improves the prediction accuracy and can be applied to a wider $Ma$ range compared with the canonical Q-P theory. The correlation indicates that the sharp increase of wall heat transfer in the STBLI can be explained by the boundary layer compression and the convection transport enhancement. Based on the DNS results, the aerodynamic heat generation and transport mechanisms are revealed in the separation, recirculation and reattachment zones in the STBLI. From this perspective, the peak heat flux can be further explained by the enhancement of near-wall turbulent energy dissipation, compression aerodynamic heat generation and the near-wall turbulent transport. The generation and transport of compression aerodynamic heat reveal the underlying mechanism of the strong correlation between the peak heat flux ratios and the pressure ratios in STBLIs.

Accurate prediction of the global properties of wall-bounded turbulence holds significant importance for both fundamental research and engineering applications. In atmospheric boundary layers, the relationship between friction drag and geostrophic wind is described by the geostrophic drag law (GDL). We use carefully designed large-eddy simulations to study nocturnal stable atmospheric boundary layers (NSBLs), which are characterized by a negative potential temperature flux at the surface and neutral stratification higher up. Our simulations explore a wider range of the Kazanski–Monin parameter, $\mu = L_f / L_s = [16.7, 193.3]$, with $L_f$ the Ekman length scale and $L_s$ the surface Obukhov length. We show collapse of the GDL coefficients onto single curves as functions of $\mu$, thereby validating the GDL's applicability to NSBLs over a very wide $\mu$ range. We show that the boundary-layer height $h$ scales with $\sqrt {L_f L_s}$, while both the streamwise and spanwise wind gradients scale with $u_*^2 / (h^2 f)$, where $u_*$ represents the friction velocity and $f$ the Coriolis parameter. Leveraging these insights, we developed new analytical expressions for the GDL coefficients, significantly enhancing our understanding of the GDL for turbulent boundary layers. These formulations facilitate the analytical prediction of the geostrophic drag coefficient and cross-isobaric angle.

This work introduces a formulation of resolvent analysis that uses wavelet transforms rather than Fourier transforms in time. Under this formulation, resolvent analysis may extend to turbulent flows with non-stationary mean states. The optimal resolvent modes are augmented with a temporal dimension and are able to encode the time-transient trajectories that are most amplified by the linearised Navier–Stokes equations. We first show that the wavelet- and Fourier-based resolvent analyses give equivalent results for statistically stationary flow by applying them to turbulent channel flow. We then use wavelet-based resolvent analysis to study the transient growth mechanism in the near-wall region of a turbulent channel flow by windowing the resolvent operator in time and frequency. The computed principal resolvent response mode, i.e. the velocity field optimally amplified by the linearised dynamics of the flow, exhibits characteristics of the Orr mechanism, which supports the claim that this mechanism is key to linear transient energy growth. We also apply this method to non-stationary parallel shear flows such as an oscillating boundary layer, and three-dimensional channel flow in which a sudden spanwise pressure gradient perturbs a fully developed turbulent channel flow. In both cases, wavelet-based resolvent analysis yields modes that are sensitive to the changing mean profile of the flow. For the oscillating boundary layer, wavelet-based resolvent analysis produces oscillating principal forcing and response modes that peak at times and wall-normal locations associated with high turbulent activity. For the turbulent channel flow under a sudden spanwise pressure gradient, the resolvent modes gradually realign themselves with the mean flow as the latter deviates. Wavelet-based resolvent analysis thus captures the changes in the transient linear growth mechanisms caused by a time-varying turbulent mean profile.

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the $-u^{\prime }w^{\prime }$ co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

This study investigates the correlation between the fluctuating wall heat flux, and the distribution and transport of Reynolds shear stress and turbulent heat flux in compressible boundary layers at Mach number 5.86 and friction Reynolds number 420, with a relatively weaker and a stronger wall cooling imposed. As illustrated from the probability density functions of the wall-heat-flux perturbations, with increasing wall cooling, the extreme wall heat flux is intensified and tends to be more negatively skewed. To examine the role of the extreme events in the transport of the momentum and heat, conditional analysis of the extreme positive and negative wall-heat-flux-perturbation events is conducted. In most regions of the boundary layer, the positive events are predominantly associated with an increase in Reynolds shear stress and a decrease in turbulent heat flux. Joint probability density functions of velocity and wall-heat-flux perturbations in the near-wall region indicate that the extreme positive events tend to be more correlated with ejections, which is particularly evident in the stronger wall-cooling case. To further shed light on the underlying mechanisms of the connections between wall heat flux and transport budgets, a transport equation for turbulent heat flux is derived, in a similar manner to that for Reynolds shear stress. The energy balance is inspected, with conditional analysis applied to budget terms and mean flow properties so as to quantify the correlation between wall-heat-flux fluctuations and energy evolution.

Turbulent flows over rough surfaces can be encountered in a wide range of engineering applications. Despite the progress made after several decades of studies, the prediction of drag and roughness function from the surface geometrical parameters remains an open question. Several methods have shown encouraging results. However, they lack generality and present some scatter in the data. In this paper we propose a new parameter, the effective distribution ($ED$), which lays foundation on the effective slope with some changes to take into account the sheltering effect of large roughness elements and the drag induced by pinnacles higher than the average roughness elements. To develop this new correlation between geometrical features of the wall and the drag, we performed a set of simulations of the turbulent flow over a rough surface made of triangular elements varying their height and spatial distribution. The $ED$ correlates quite well both with the drag and the roughness function for a wide range of cases having different mean roughness height, skewness and kurtosis. To further validate the $ED$, and assessing how it can be generalized to real rough wall, an irregular wall made from the superposition of random sinusoidal function was considered. Results were consistent with the correlation here presented.

Superhydrophobic (SHPo) surfaces can capture a thin layer of air called a plastron under water to reduce skin friction. Although a ~30 % drag reduction has been recently reported with longitudinal micro-trench SHPo surfaces under a boat and in a towing tank, the results lacked the consistency to establish a clear trend. Designed based on Yu et al. (J. Fluid Mech, vol. 962, 2023, A9), this work develops and tests a series of high-performance SHPo surface coupons that can sustain a pinned plastron underneath a passenger motorboat revamped to reach 14 knots. Importantly, plastrons in a pinned state, not just their existence, are confirmed during flow experiments for the first time. All the drag-reduction data measured on different longitudinal micro-trenches are found to collapse if plotted against slip length in wall units. In comparison, aligned posts and transverse trenches show less and little drag reduction, respectively, confirming the adverse effect of the spanwise slip in turbulent flows. This report not only verifies SHPo surfaces can provide a consistent drag reduction at high speeds in open sea but also shows that one may predict the amount of drag reduction in turbulent flows using the physical slip length obtained for Stokes flows.

In the present study, we propose a Reynolds analogy model for compressible wall turbulence. This model is demonstrated to be able to alleviate the defects of the generalized Reynolds analogy model (GRA) (Zhang et al., J. Fluid Mech., vol. 739, 2014, pp. 392–420), and maintains its success in describing the mean velocity–temperature relation. Furthermore, the present model is superior to the GRA in depicting the relationship between their fluctuating fields and also bridges the gap between the phenomenological model and the mathematical representation of the Reynolds analogy. The key points of the present model are validated by analysing the data of compressible wall-bounded turbulence with different Mach numbers, Reynolds numbers and wall thermal conditions.

Thwaites (Aeronaut. Q., vol. 1, 1949, pp. 245–280) developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations (Alber, 9th Aerospace Sciences Meeting, 1971). The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

Direct numerical simulations of oscillatory flow over a bed made of ripples have been performed. Two oscillatory flow forcing mechanisms have been compared: (i) a sinusoidal external pressure gradient (pressure-driven flow); and (ii) a sinusoidal velocity boundary conditions on the rippled bed (shear-driven flow). In the second case, the oscillations of the bed are such that when observed from a reference frame fixed with the bed, the free stream follows the same harmonic oscillation as in the pressure-driven case. While the outer layers have the same dynamics in the two cases, close to the bed differences are observed during the cycle, mostly because the large form drag across the ripples cannot be reproduced in the shear-driven case. A comparison against experimental data from an oscillating tray apparatus provides a relatively good agreement for the phase-averaged flow when the same forcing is considered (i.e. a shear-driven flow). The pressure-driven case has a comparable error to the shear-driven numerical results over the crest of the ripples, whereas the discrepancy is larger at the troughs. The discrepancies between the two cases are more limited for time-averaged flow quantities, such as the mean flow pattern and the time-averaged Reynolds stress distribution. This suggests that numerical or experimental shear-driven configurations may capture well the net effects of coastal transport processes (which occur in pressure-driven oscillatory flow), but care should be exercised in interpreting phase-dependent dynamics near the troughs. More work is needed to fully assess the sensitivity to the forcing mechanisms in different flow regimes.

This study aims to quantify how turbulence in a channel flow mixes momentum in the mean sense. We applied the macroscopic forcing method (Mani & Park, Phys. Rev. Fluids, 2021, 054607) to direct numerical simulation (DNS) of a turbulent channel flow at $Re_\tau =180$ using two different forcing strategies that are designed to separately assess the anisotropy and non-locality of momentum mixing. In the first strategy, the leading term of the Kramers–Moyal expansion of the eddy viscosity is quantified, revealing all 81 tensorial coefficients that essentially characterise the local-limit eddy viscosity. The results indicate the following: (1) the eddy viscosity has significant anisotropy, (2) Reynolds stresses are generated by both the mean strain rate and mean rotation rate tensors associated with the momentum field and (3) the local-limit eddy viscosity generates asymmetric Reynolds stress tensors. In the second strategy, the eddy viscosity is quantified as an integration kernel revealing the non-local influence of the mean momentum gradient at each wall-normal coordinate on all nine components of the Reynolds stresses over the channel width. Our results indicate that while the shear component of the Reynolds stress is reasonably reproduced by the local mean gradients, other components of the Reynolds stress are highly non-local. These results provide an understanding of anisotropy and non-locality requirements for closure modelling of momentum transport in attached wall-bounded turbulent flows.