A phenomenological description is presented to explain the intermediate and low-frequency/large-scale contributions to the wall-shear-stress (${\tau }_w$) and wall-pressure ($\,{p}_w$) spectra of canonical turbulent boundary layers, both of which are well known to increase with Reynolds number, albeit in a distinct manner. The explanation is based on the concept of active and inactive motions (Townsend, J. Fluid Mech., vol. 11, issue 1, 1961, pp. 97–120) associated with the attached-eddy hypothesis. Unique data sets of simultaneously acquired ${\tau }_w$, ${p}_w$ and velocity-fluctuation time series in the log region are considered, across a friction-Reynolds-number ($Re_{\tau }$) range of $ {O}(10^3) \lesssim Re_{\tau } \lesssim {O}(10^6)$. A recently proposed energy-decomposition methodology (Deshpande et al., J. Fluid Mech., vol. 914, 2021, A5) is implemented to reveal the active and inactive contributions to the ${\tau }_w$- and $p_w$-spectra. Empirical evidence is provided in support of Bradshaw's (J. Fluid Mech., vol. 30, issue 2, 1967, pp. 241–258) hypothesis that the inactive motions are responsible for the non-local wall-ward transport of the large-scale inertia-dominated energy, which is produced in the log region by active motions. This explains the large-scale signatures in the ${\tau }_w$-spectrum, which grow with $Re_{\tau }$ despite the statistically weak signature of large-scale turbulence production, in the near-wall region. For wall pressure, active and inactive motions respectively contribute to the intermediate and large scales of the $p_w$-spectrum. Both these contributions are found to increase with increasing $Re_{\tau }$ owing to the broadening and energization of the wall-scaled (attached) eddy hierarchy. This potentially explains the rapid $Re_{\tau }$-growth of the $p_w$-spectra relative to ${\tau }_w$, given the dependence of the latter only on the inactive contributions.

Adverse pressure gradient (APG) turbulent boundary layers (TBL) require an understanding of the details of the pressure gradient, or history effect, to characterize the associated variation of spatiotemporal turbulent statistics. The streamwise-varying mean pressure gradient is reflected in the streamwise developing mean flow field and thus resolvent analysis, which captures the amplification of the Navier–Stokes equations linearized about the turbulent mean, can be used to understand linear amplification in APG TBLs. In particular, by using a biglobal approach in which the amplification is characterized by a temporal frequency and spanwise wavenumber, the streamwise and wall-normal inhomogeneities of the APG TBL can be resolved and related to the APG history. The linear response is able to identify multiscale phenomena, identifying a near-wall peak with $\lambda _{z}^+\approx 100$ for zero pressure gradient TBLs and mild to moderate APG TBLs as well as large-scale modes whose amplification increases with APG strength and Reynolds number. It is shown that the monotonic growth in the turbulent statistics with increasing APG is reflected in the linear growth in the associated resolvent amplification. Collapse in the Reynolds stresses is obtained through an augmented hybrid velocity scale, which replaces the local APG strength measure in the hybrid velocity scale presented in Romero et al. (Intl J. Heat Fluid Flow, vol. 93, 2022, 108885) with a velocity that encapsulates the pressure gradient history. While this resolvent approach is applicable to any APG TBL, it is shown from a scaling analysis of the linearized Navier–Stokes equations that the linear growth observed in the resolvent amplification with the history effect is limited to near-equilibrium APG TBLs.

The amplitude modulation coefficient, $R$, that is widely used to characterize nonlinear interactions between large- and small-scale motions in wall-bounded turbulence is not compatible with detecting the convective nonlinearity of the Navier–Stokes equations. Through a spectral decomposition of $R$ and a simplified model of triadic convective interactions, we show that $R$ suppresses the signature of convective scale interactions, but is strongly influenced by linear interactions between large-scale motions and the background mean flow. We propose an additional coefficient that is specifically designed for the detection of convective nonlinearities, and we show how this new coefficient, $R_T$, quantifies the turbulent kinetic energy transport involved in turbulent scale interactions and reveals a classical energy cascade across widely separated scales.

Prandtl's secondary flows of the second kind generated by laterally varying roughness are studied using the linearised Reynolds-averaged Navier–Stokes approach proposed by Zampino et al. (J. Fluid Mech., vol. 944, 2022, p. A4). The momentum equations are coupled to the Spalart–Allmaras model while the roughness is captured by adapting established strategies for homogeneous roughness to heterogeneous surfaces. Linearisation of the governing equations yields a framework that enables a rapid exploration of the parameter space associated with heterogeneous surfaces, in the limiting case of small spanwise variations of the roughness properties. Channel flow is considered, with longitudinal high- and low-roughness strips arranged symmetrically. By varying the strip width, it is found that linear mechanisms play a dominant role in determining the size and intensity of secondary flows. In this setting, secondary flows may be interpreted as the time-averaged output response of the turbulent mean flow subjected to a steady forcing produced by the wall heterogeneity. In fact, the linear model predicts that secondary flows are most intense when the strip width is about 0.7 times the half-channel height, in excellent agreement with available data. Furthermore, a unified framework to analyse combinations of heterogeneous roughness properties and laterally varying topographies, common in applications, is discussed. Noting that the framework assumes small spanwise variations of the surface properties, two separate secondary-flow-inducing source mechanisms are identified, i.e. the lateral variation of the virtual origin from which the turbulent structure develops and the lateral variation of the streamwise velocity slip, capturing the acceleration/deceleration perceived by the bulk flow over troughs and crests of non-planar topographies.

We explore the application of the reference map technique, originally developed for Eulerian simulation of solid mechanics, in Lagrangian kinematics of turbulent flows. Unlike traditional methods based on explicit particle tracking, the reference map facilitates the calculation of flow maps and gradients without the need for particles. This is achieved through an Eulerian update of the reference map, which records the take-off positions of fluid particles. This approach is found to be mathematically equivalent to the work of Leung (J. Comput. Phys., vol. 230, issue 9, 2011, pp. 3500–3524), who computed the flow map of simple two-dimensional flows using an Eulerian approach. We discuss important modifications necessary for its first application to complex three-dimensional turbulent flows, including the conservative, low-dissipation update of the flow map and the treatment of periodic boundary conditions. We first demonstrate the accuracy of finite-time Lyapunov exponent (FTLE) calculations based on the reference map against the standard particle-based approach in a two-dimensional Taylor–Green vortex. Then we apply it to turbulent channel flow at $Re_\tau =180$, where Lagrangian coherent structures identified as ridges of the backward-time FTLE are found to bound vortical regions of flow, consistent with Eulerian coherent structures from the $Q$-criterion. The reference map also proves suitable for material surface tracking despite not explicitly tracking particles. This capability can provide valuable insights into the Lagrangian landscape of turbulent momentum transport, complementing Eulerian velocity field analysis. The evolution of initially wall-normal material surfaces in the viscous sublayer, buffer layer and log layer sheds light on the Reynolds stress-generating events from a Lagrangian perspective.

Reynolds-averaged models for solving the Navier–Stokes equations are implicitly based on Kolmogorov's theory for describing energy transfers between the different turbulent scales, which means that all the energy produced at large scales is transferred at a constant rate to the smallest turbulent scales where it is dissipated. As a result, these models use a single scale to describe the turbulence spectrum, which in cases of non-equilibrium turbulence does not provide an adequate description of the transfers actually observed. This is particularly the case for wall-bounded flows at high Reynolds numbers, such as turbulent channel flows. Taking up an approach developed by Schiestel (2007 Modeling and Simulation of Turbulent Flows, ISTE Ltd and John Wiley & Sons), which aims to define a Reynolds-averaged Navier–Stokes model transporting several scales of turbulence, a two-scale Reynolds stress model (RSM) was developed in order to take into account the interactions between the inner and outer regions of wall-bounded flows. The results obtained with the model are compared with the direct numerical simulations (DNS) of Lee & Moser (J. Fluid Mech., vol. 860, 2019, pp. 886–938) in a turbulent channel for several friction Reynolds numbers up to $Re_{\tau }=5200$, for which partial integrations in spectral space were carried out, highlighting distinct behaviours between small and large scales of turbulence. The model developed provides an accurate description of the contributions at small and large scales and thus reproduces the high-Reynolds-number effects observed in DNS data. In addition, comparisons with the DNS data served to validate a large part of the closure relations used for the various terms in the two-scale RSM.

Rotating convection is considered on the tilted $f$-plane where gravity and rotation are not aligned. For sufficiently large rotation rates, $\Omega$, the Taylor–Proudman effect results in the gyroscopic alignment of anisotropic columnar structures with the rotation axis giving rise to rapidly varying radial length scales that vanishes as $\Omega ^{-1/3}$ for $\Omega \rightarrow \infty$. Compounding this phenomenon is the existence of viscous (Ekman) layers adjacent to the impenetrable bounding surfaces that scale as $\Omega ^{-1/2}$. In this investigation, these constraints are relaxed upon utilising a non-orthogonal coordinate representation of the fluid equations where the upright coordinate aligns with rotation axis. This exposes the problem to asymptotic perturbation methods that permit: (i) relaxation of the constraints of gyroscopic alignment; (ii) the filtering of Ekman layers through the uncovering of parameterised velocity pumping boundary conditions; and (iii) the development of reduced quasi-geostrophic systems valid in the limit $\Omega \rightarrow \infty$. Linear stability investigations reveal excellent quantitative agreement between results from parameterised or unapproximated mechanical boundary conditions. For no-slip boundaries, it is demonstrated that the associated Ekman pumping alters convective onset through an enhanced destabilisation of large spatial scales. The range of unstable modes at a fixed thermal forcing is thus significantly extended with a direct dependence on $\Omega$. This holds true even for geophysical and astrophysical regimes characterised by extreme values of the non-dimensional Ekman number $E$. The nonlinear regime is explored via the global heat and momentum transport of single-mode solutions to the quasi-geostrophic systems which indicate $O(1)$ changes which do not scale with the size of $E$.

The seminal Bolgiano–Obukhov (BO) theory established the fundamental framework for turbulent mixing and energy transfer in stably stratified fluids. However, the presence of BO scalings remains debatable despite their being observed in stably stratified atmospheric layers and convective turbulence. In this study, we performed precise temperature measurements with 51 high-resolution loggers above the seafloor for 46 h on the continental shelf of the northern South China Sea. The temperature observation exhibits three layers with increasing distance from the seafloor: the bottom mixed layer (BML), the mixing zone and the internal wave zone. A BO-like scaling $\alpha =-1.34\pm 0.10$ is observed in the temperature spectrum when the BML is in a weakly stable stratified ($N\sim 0.0018$ rad s$^{-1}$) and strongly sheared ($Ri\sim 0.0027$) condition, whereas in the unstably stratified convective turbulence of the BML, the scaling $\alpha =-1.76\pm 0.10$ clearly deviated from the BO theory but approached the classical $-$5/3 scaling in isotropic turbulence. This suggests that the convective turbulence is not the promise of BO scaling. In the mixing zone, where internal waves alternately interact with the BML, the scaling follows the Kolmogorov scaling. In the internal wave zone, the scaling $\alpha =-2.12 \pm 0.15$ is observed in the turbulence range and possible mechanisms are provided.

Boundary layers are present in many natural and industrial fluid flows. The concept of boundary layers can be traced back to Leonardo da Vinci's paintings of pipe flow, where he was aware of a higher velocity away from the walls. During the 19th century, the physics of boundary conditions had been extensively debated, and the well-known Maxwell–Navier slip length was proposed in 1823. In most cases, the no-slip boundary condition is valid at a fluid–solid interface. However, with the advancement of measurement techniques, slip lengths ranging from nanometre to micrometre scales were experimentally measured, raising questions regarding the applicability of the no-slip condition. In 2003, Lauga & Stone (J. Fluid Mech., vol. 489, 2003, pp. 55–77) proposed a simple model to elucidate the effect of surface heterogeneities on the slip length, elegantly bridging the microscopic structure of the wall-boundary conditions to the macroscopic effective slip length.

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.

A previous paper of the authors (Duck & Stephen, J. Fluid Mech., vol. 917, 2021, A56) considered the effect of three-dimensional, temporally periodic, linear and incompressible disturbances on a Blasius boundary layer, in particular when the disturbance wavelength is both comparable to and longer than the boundary-layer thickness. This previous study revealed that, unlike the two-dimensional counterpart, a mode exists that exhibits regimes of downstream spatial growth. In this paper we extend the analysis to the compressible regime, based on the boundary-region equations methodology. The aforementioned unstable mode is seen to persist into the compressible regime, and is studied using a combination of numerical and asymptotic methods. The paper adopts several approaches. First is a numerical approach in which the spatial development of the disturbances is assessed. This then leads to a consideration of the far-downstream behaviour, using (several) asymptotic limits. Of some note, in addition to unstable modes found in the incompressible case, is the existence of a further class of instability, not found in the incompressible case (which is also analysed asymptotically), corresponding to what amounts to an inviscid instability. The far-downstream analysis enables a (sub-)classification into entropy and non-entropy modes. The former, according to this analysis, are spatially damped, with one caveat, as revealed by our marching procedure, which highlights how spatial development of disturbances can be important.

The present study investigates streamwise ($\overline {u^2}$) energy-transfer mechanisms in the inner and outer regions of turbulent boundary layers (TBLs). Particular focus is placed on the $\overline {u^2}$ production, its inter-component and wall-normal transport as well as dissipation, all of which become statistically significant in the outer region with increasing friction Reynolds number ($Re_{\tau }$). These properties are analysed using published data sets of zero, weak and moderately strong adverse-pressure-gradient (APG) TBLs across a decade of $Re_{\tau }$, revealing similarity in energy-transfer pathways for all these TBLs. It is found that both the inner and outer peaks of $\overline {u^2}$ are always associated with local maxima in the $\overline {u^2}$ production and its inter-component transport, and the regions below/above each of these peaks are always dominated by wall-ward/away-from-wall transport of $\overline {u^2}$, thereby classifying the $\overline {u^2}$ profiles into four distinct regimes. This classification reveals existence of phenomenologically similar energy-transfer mechanisms in the ‘inner’ and ‘outer’ regions of moderately strong APG TBLs, which meet at an intermediate location coinciding with the minimum in $\overline {u^2}$ profiles. Conditional averaging suggests existence of similar phenomena even in low $Re_{\tau }$ canonical and/or weak APG TBLs, albeit with the outer-region mechanisms weaker than those in the inner region. This explains the absence of their $\overline {u^2}$ outer peak and the dominance of $\overline {u^2}$ wall-normal transport away from the wall, which potentially originates from the inner region. Given that the wall-ward/away-from-wall transport of $\overline {u^2}$ is governed by the $Q_4$(sweeps)/$Q_2$(ejections) quadrants of the Reynolds shear stress, it is argued that the emergence of the $\overline {u^2}$ outer peak corresponds with the statistical dominance of $Q_4$ events in the outer region. Besides unravelling the dynamical significance of $Q_2$ and $Q_4$ events in the outer region of TBLs, the present analysis also proposes new phenomenological arguments for testing on canonical wall-turbulence data at very high $Re_{\tau }$.

High-frequency observation data, including all three components of instantaneous fluctuating velocity, temperature, as well as particulate matter 10 ($PM_{10}$), collected from the unstable atmospheric surface layer at $z/L = -0.11$ and $-$0.12, $L$ being the Obukhov length, during sand and dust storms (SDS), were used to explore the scaling of vertical coherence and the logarithmic energy profile for wall-attached eddies. The present results demonstrate good agreement with the self-similar range of the wall-attached features for velocity and temperature components, as well as for $PM_{10}$ at lower heights ($z<15$ m) during SDS. Following the idea depicted by Davenport (Q. J. R. Meteorol., vol. 372, 1961, pp. 194–211), an empirically derived transfer kernel comprises implicit filtering via a scale-dependent gain and phase, parametrically defined as $|H_L^2(f)|=\exp (c_1-c_2\delta /\lambda _x)$, where $c_1$ and $c_2$ are parameters, $\delta$ is the boundary layer thickness and $\lambda _x$ is the streamwise wavelength. Linear coherence spectrum analysis is applied as a filter to separate the coherent and incoherent portions. After this separation procedure, the turbulence intensity decay for wall-attached eddies is described in a log–linear manner, which also identifies how the scaling parameter differs between the measured components. These findings present abundant features of wall-attached eddies during SDS which further are used to improve/enrich existing near-wall models.

In this experimental study, we explore the dynamics of the thermal boundary layer in liquid metal Rayleigh–Bénard convection, covering the parameter ranges of $0.026 \leq$ Prandtl numbers $(Pr) \leq 0.033$ and Rayleigh numbers ($Ra$) up to $2.9\times 10^9$. Our research focuses on characterising the thermal boundary layer near the top plate of a cylindrical convection cell with an aspect ratio of 0.5, distinguishing between two distinct regions: the shear-dominated region around the centre of the top plate and a location near the side wall where the boundary layer is expected to be affected by the impact or ejection of thermal plumes. The dependencies of the boundary layer thickness on $Ra$ at these positions reveal deviating scaling exponents with the difference diminishing as $Ra$ increases. We find stronger fluctuations in the boundary layer and increasing deviation from the Prandtl–Blasius–Pohlhausen profile with increasing $Ra$, as well as in the measurements outside the centre region. Our data illustrate the complex interplay between flow dynamics and thermal transport in low-$Pr$ convection.

The influence of symmetry-breaking effects of ridge-type roughness on secondary currents in turbulent channel flow is investigated using direct numerical simulations. The ridges have triangular cross-section, which is systematically varied from isosceles to right-angled triangle, introducing an imbalance to the slopes of the ridges’ lateral surfaces while the streamwise homogeneity of the surfaces is maintained. In all cases, secondary current vortices are produced, but asymmetric ridge cross-sections break the symmetry of these vortices. As a result of the asymmetry-induced misalignment and imbalance in the secondary current vortices, net spanwise flow emerges. The magnitude of the spanwise flow increases with the slope ratio of the ridge lateral surfaces and significantly modifies the mean flow topology, leading to the merging of critical points in the case of the right-angled triangular ridge shape. Within the cavities, the net spanwise flow is accompanied by a non-zero mean spanwise pressure gradient, while from the perspective of the outer flow, the scalene ridge surfaces have a similar effect as a wall that is slowly moving in the spanwise direction. Overall, the present results suggest the existence of a special type of Prandtl's secondary currents of the second kind, namely those that result in net spanwise flow.

This study proposes a novel super-resolution (or SR) framework for generating high-resolution turbulent boundary layer (TBL) flow from low-resolution inputs. The framework combines a super-resolution generative adversarial neural network (SRGAN) with down-sampling modules (DMs), integrating the residual of the continuity equation into the loss function. The DMs selectively filter out components with excessive energy dissipation in low-resolution fields prior to the super-resolution process. The framework iteratively applies the SRGAN and DM procedure to fully capture the energy cascade of multi-scale flow structures, collectively termed the SRGAN-based energy cascade reconstruction framework (EC-SRGAN). Despite being trained solely on turbulent channel flow data (via ‘zero-shot transfer’), EC-SRGAN exhibits remarkable generalization in predicting TBL small-scale velocity fields, accurately reproducing wavenumber spectra compared to direct numerical simulation (DNS) results. Furthermore, a super-resolution core is trained at a specific super-resolution ratio. By leveraging this pretrained super-resolution core, EC-SRGAN efficiently reconstructs TBL fields at multiple super-resolution ratios from various levels of low-resolution inputs, showcasing strong flexibility. By learning turbulent scale invariance, EC-SRGAN demonstrates robustness across different TBL datasets. These results underscore the potential of EC-SRGAN for generating and predicting wall turbulence with high flexibility, offering promising applications in addressing diverse TBL-related challenges.

The scale-dependent variability of convective velocities and structure inclination angles in wall-bounded turbulence was studied experimentally via space–time energy spectrum measurements. We found that the variability of convection velocities for large-scale motions (LSMs) decreased inversely with streamwise wavenumbers, and that the variability trend was not fully explained by earlier applications of Kraichnan's ‘random-sweeping’ model of turbulence that assumed perfect scale separation. By analytically extending the random-sweeping model to allow for a dominant large scale in the random-sweeping signal that can interact with other LSMs, we showed how scale interactions can explain the variability trend in convection velocities for LSMs. The variability in convection velocities was also shown to correlate with the scale-dependent inclination angles of coherent structures that were obtained via cross-spectral analysis. Large-scale motions tended to exhibit shallower inclination to the wall with increasing convection velocity, while small-scale motions and those far from the wall exhibited the reverse behaviour. We proposed that these two opposite relationships between inclination angle and convection velocity can be explained in terms of a balance between opposing effects of the mean shear and the coherent structure geometry. Descriptions and models of convection velocity variability effects are useful both for modelling turbulence spectra and explaining the geometry of coherent structures.

A new wall-wake law is proposed for the streamwise turbulence in the outer region of a turbulent boundary layer. The formulation pairs the logarithmic part of the profile (with a slope $A_1$ and additive constant $B_1$) to an outer linear part, and it accurately describes over 95 % of the boundary layer profile at high Reynolds numbers. Once the slope $A_1$ is fixed, $B_1$ is the only free parameter determining the fit. Most importantly, $B_1$ is shown to follow the same trend with Reynolds number as the wake factor in the wall-wake law for the mean velocity, which is tied to changes in scaling of the mean flow and the turbulence that occur at low Reynolds number.

Two sets of measurements utilizing hot-wire anemometry and oil-film interferometry for flat-plate turbulent boundary layers, exposed to various controlled adverse and favourable pressure gradients, are used to evaluate history effects of the imposed and varying free-stream gradients. The results are from the NDF wind tunnel at Illinois Tech (IIT) and the MTL wind tunnel at KTH, over the range $800 < Re_\tau < 22\,000$ (where $Re_{\tau }$ is the friction Reynolds number). The streamwise pressure-gradient parameter $\beta \equiv (-\ell /\tau _{w}) \cdot (\partial P_{e}/\partial x)$ varied between $-2 < \beta < 7$, where $\ell$ is an outer length scale for boundary layers equivalent to the half-height of channel flow and the radius of pipe flow, and is estimated for each boundary-layer profile; note that $\tau_w$ is the wall-shear stress and $P_e$ is the free-stream static pressure. Extracting from each profile the three parameters of the overlap region, following the recent work of Monkewitz & Nagib (J. Fluid Mech., vol. 967, 2023, p. A15) that led to an overlap region of combined logarithmic and linear parts, we find minimum history effects in the overlap region. Thus, the overlap region in this range of pressure-gradient boundary layers appears to be in ‘quasiequilibrium’.

The topology of turbulent coherent structures is known to regulate the transport of energy, mass and momentum in the atmospheric boundary layer (ABL). While previous research has primarily focused on characterizing the structure of turbulence in stationary ABL flows, real-world scenarios frequently deviate from stationarity, giving rise to nuanced and poorly understood changes in the turbulence geometry and associated transport mechanisms. This study sheds light on this problem by examining topological changes in ABL turbulence induced by non-stationarity and their effects on momentum transport. Results from a large-eddy simulation of pulsatile open channel flow over an array of surface-mounted cuboids are examined. The analysis reveals that the flow pulsation triggers a phase-dependent shear rate, and the ejection-sweep pattern varies with the shear rate during the pulsatile cycle. From a turbulence structure perspective, it is attributed to the changes in the geometry of hairpin vortices. An increase (decrease) in the shear rate intensifies (relaxes) these structures, leading to an increase (decrease) in the frequency of ejections and an amplification (reduction) of their percentage contribution to the total momentum flux. Furthermore, the size of the hairpin packets undergoes variations, which depend on the geometry of the constituting hairpin vortices, yet the packet inclination preserves its orientation throughout the pulsatile cycle. These observations reinforce the important role non-stationarity holds in shaping the structure of ABL turbulence and the momentum transport mechanisms it governs.