Secondary flows induced by spanwise heterogeneous surface roughness play a crucial role in determining engineering-relevant metrics such as surface drag, convective heat transfer and the transport of airborne scalars. While much of the existing literature has focused on idealized configurations with regularly spaced roughness elements, real-world surfaces often feature irregularities, clustering and topographic complexity for which the secondary flow response remains poorly understood. Motivated by this gap, we investigate multicolumn roughness configurations that serve as a regularized analogue of roughness clustering. Using large-eddy simulations, we systematically examine secondary flows across a controlled set of configurations in which cluster density and local arrangement are varied in an idealized manner, and observe that these variations give rise to distinct secondary flow polarities. Through a focused parameter study, we identify the spanwise gap between the edge-most roughness elements of adjacent columns, normalized by the channel half-height ($s_a/H$), as a key geometric factor governing this polarity. In addition to analysing the time-averaged structure, we investigate how variations in polarity affect the instantaneous dynamics of secondary flows. Here, we find that the regions of high- and low-momentum fluid created by the secondary flows alternate in a chaotic, non-periodic manner over time. Further analysis of the vertical velocity signal shows that variability in vertical momentum transport is a persistent and intrinsic feature of secondary flow dynamics. Taken together, these findings provide a comprehensive picture of how the geometric arrangement of roughness elements governs both the mean structure and temporal behaviour of secondary flows.

Jet vortex generators (JVGs) are a promising technique for controlling laminar separation in low-Reynolds-number aerofoils, such as those used in micro air vehicles (MAVs). While previous studies have demonstrated their aerodynamic benefits, the three-dimensional structure of the vortices they generate and their interaction with the boundary layer remain poorly characterised experimentally. In this study, volumetric velocity measurements are performed using the double-pulse Shake-the-Box (STB) technique on an SD7003 aerofoil equipped with skewed and pitched JVGs. Experiments are conducted at Reynolds numbers of 30 000 and 80 000, for angles of attack of 8$^{\circ}$, 10$^{\circ}$ and 14$^{\circ}$. The results provide the first experimental visualisation of the full three-dimensional vortex topology induced by JVGs, revealing asymmetric streamwise vortices that penetrate the separated shear layer and re-energise the near-wall region. In pre-stall conditions, the JVGs reshape the laminar separation bubble into a thinner and more stable structure, reducing its sensitivity to angle of attack. In stall conditions, they induce partial or full flow reattachment, delaying large-scale separation. The evolution of characteristic bubble parameters and the chordwise distribution of the shape factor $H = \delta ^{\ast }/\theta$, where $\delta ^{\ast }$ is the displacement thickness and $\theta$ is the momentum thickness, show a consistent trend of enhanced boundary-layer recovery. These findings offer new insight into the physical mechanisms underlying active separation control at low Reynolds numbers and establish a framework for evaluating vortex-based control strategies using volumetric diagnostics.

High Reynolds number effects of wall-bounded flows, involving interscale energy transfers between small and large scales of turbulence within and between the inner and outer regions, challenge the classical description of the structure of these flows and the ensuing turbulence models. The two-scale Reynolds stress model recently proposed by Chedevergne et al. (2024, J. Fluid Mech. vol. 1000), was able to reproduce the small- and large-scale contributions in turbulent channel flows that follow the scale separation performed by Lee & Moser (2019, J. Fluid Mech. vol. 860, pp. 886–938), by partitioning energy spectra at a given wavelength. However, the interscale interactions within the inner region were modelled in an ad hoc manner, but without physical relevance, making the two-scale Reynolds stress model less and less accurate for boundary layer applications as the Reynolds number was increased. In this study, by re-analysing direct numerical simulations data from Lee & Moser (2019), with the objective of modelling these scale interactions, crucial observations on energy transfers between large and small scales could be made. In particular, the analysis reveals the important role played by the spanwise component of the Reynolds stress in the logarithmic region. From the analysis undertaken, a revisited version of the two-scale model was thus proposed, focusing efforts on interscale transfer modelling. The resulting model is then successfully tested on high Reynolds number boundary layer configurations without pressure gradient, up to $\textit{Re}_{\tau }=20\,000$. The excellent agreement reflects the good prediction capabilities of the proposed model, and above all, the relevance of the modelling of the energy transfers within and between the inner and outer regions of wall-bounded flows.

Uniform momentum zones (UMZs) are widely used to describe and model the coherent structure of wall-bounded turbulent flows, but their detection has traditionally relied on relatively narrow fields of view which preclude fully resolving features at the scale of large-scale motions (LSMs). We refine and extend recent proposals to detect UMZs with moving-window fields of view by including physically motivated coherency criteria. Using synthetic data, we show how this updated moving-window approach can eliminate noise contamination that is likely responsible for the previously reported, high fractal dimension of UMZ interfaces. By applying the approach to channel flow direct numerical simulation (DNS), we identify a significant number of previously undetected, large-scale UMZ interfaces, including a small fraction of highly linear interfaces with well-defined streamwise inclination angles. We show that the inclination angles vary inversely with the size of the UMZ interfaces and that this relationship can be modelled by the opposing effects of shear-induced inclination and vortex-induced lift-up on hairpin packets. These geometric properties of large-scale UMZ interfaces play an important role in the development of improved stochastic models of wall-bounded turbulence.

Fully resolving turbulent flows remains challenging due to a turbulent systems’ multiscale complexity. Existing data-driven approaches typically demand expensive retraining for each flow scenario and struggle to generalize beyond their training conditions. Leveraging the universality of small-scale turbulent motions (Kolmogorov’s K41 theory), we propose a scale-oriented zonal generative adversarial network (SoZoGAN) framework for high-fidelity, zero-shot turbulence generation across diverse domains. Unlike conventional methods, SoZoGAN is trained exclusively on a single dataset of moderate-Reynolds-number homogeneous isotropic turbulence (HIT). The framework employs a zonal decomposition strategy, partitioning turbulent snapshots into subdomains based on scale-sensitive physical quantities. Within each subdomain, turbulence is synthesized using scale-indexed models pretrained solely on the HIT database. A SoZoGAN demonstrates high accuracy, cross-domain generalizability and robustness in zero-shot super-resolution of unsteady flows, as validated on untrained HIT, turbulent boundary layer and channel flow. Its strong generalization, demonstrated for homogeneous and inhomogeneous turbulence cases, suggests potential applicability to a wider range of industrial and natural turbulent flows. The scale-oriented zonal framework is architecture-agnostic, readily extending beyond generative adversarial networks to other deep learning models.

We investigate the dynamics of an oscillatory boundary layer developing over a bed of collisional and freely evolving sediment grains. We perform Euler–Lagrange simulations at Reynolds numbers ${\textit{Re}}_\delta = 200$, 400 and 800, density ratio $\rho _{\!p}/\rho _{\!f} = 2.65$, Galileo number ${\textit{Ga}} = 51.9$, maximum Shields numbers from $5.60 \times 10^{-2}$ to $2.43 \times 10^{-1}$, based on smooth wall configuration, and Keulegan–Carpenter number from $134.5$ to $538.0$. We show that the dynamics of the oscillatory boundary layer and sediment bed are strongly coupled due to two mechanisms: (i) bed permeability, which leads to flow penetration deep inside the sediment layer, a slip velocity at the bed–fluid interface, and the expansion of the boundary layer, and (ii) particle motion, which leads to rolling-grain ripples at ${\textit{Re}}_\delta = 400$ and ${\textit{Re}}_\delta = 800$. While at ${\textit{Re}}_\delta = 200$ the sediment bed remains static during the entire cycle, the permeability of the bed–fluid interface causes a thickening of the boundary layer. With increasing ${\textit{Re}}_\delta$, the particles become mobile, which leads to rolling-grain ripples at ${\textit{Re}}_\delta = 400$ and suspended sediment at ${\textit{Re}}_\delta = 800$. Due to their feedback force on the fluid, the mobile sediment particles cause greater velocity fluctuations in the fluid. Flow penetration causes a progressive alteration of the fluid velocity gradient near the bed interface, which reduces the Shields number based upon bed shear stress.

The spatio-temporal evolution of very large-scale coherent structures, also known as superstructures, is investigated in both smooth- and rough-wall boundary layers by means of direct numerical simulations up to a frictional Reynolds number of ${\textit{Re}}_\tau = 3\,150$. One smooth-wall and four rough-wall cases are considered, all developing over a region as long as $\sim$60 times the incoming boundary-layer thickness in the streamwise direction. Bio-inspired, biofouling-type topographies are employed for the rough-wall cases, following the previous work of Womack et al. (2022 J. Fluid Mech. vol. 933, p. A38) and Kaminaris et al. (2023 J. Fluid Mech. vol. 961, p. A23). We utilise three-dimensional time series, as well as multiple two-point correlation functions along the boundary layer to capture the detailed length- and time-scale evolution of the superstructures. The results suggest that the presence of roughness significantly amplifies both the strength and the streamwise growth rate of superstructures. Interestingly, however, their ratios relative to the local boundary-layer thickness, $\mathscr{L}_{\!x}/\delta$ and $\mathscr{L}_z/\delta$, remain constant and independent of the streamwise coordinate, indicating that such scaled length scales might constitute a possible flow invariant. Volumetric correlations revealed that all cases induce structures inclined with respect to the mean-flow direction, with those over the rough-wall topographies exhibiting steeper inclination angles. Finally, via proper orthogonal decomposition, pairs of counter-rotating roll modes were detected and found to flank the high- and low-speed superstructures, supporting the conjecture in the literature regarding the mechanisms responsible for the lateral momentum redistribution. The latter also suggests that the way momentum organises itself in high Reynolds number wall-bounded flows might be independent of the roughness terrain underneath.

This study applies the scaling patch approach to investigate the influence of pressure gradients on the mean-momentum balance in turbulent boundary layers (TBLs). Under strong pressure gradients, the force balance in the outer region is dominated by advective and pressure forces, with gradients of Reynolds stresses playing a minimal role. To retain the relevance of Reynolds stress gradients within the scaling patch framework, we propose a redistribution of the component $U_e \textrm {d}U_e/\textrm {d}x$ from the advective term to the pressure-gradient term. Here, $U_e$ is the mean streamwise velocity at the boundary layer edge. This reformulation enhances the outer-scaling framework of Wei & Knopp (2023 J. Fluid Mech. 958, 1–21), ensuring consistency across a wide range of pressure gradients, including those involving flow separation. Remarkably, the new outer-scaled gradient of Reynolds shear stress in TBLs under a pressure gradient closely resembles that observed in zero-pressure-gradient TBLs. In the inner region, the impact of pressure gradient is well captured by the Stratford–Mellor parameter $\beta _{\textit{in}}$. For weak pressure gradients ($|\beta _{\textit{in}}| \ll 0.07$), traditional inner scaling remains valid. However, for stronger pressure gradients $|\beta _{\textit{in}}| \gtrsim 0.07$, the near-wall dynamics is governed by a balance between pressure gradient and viscous force, as described by Stratford (1959 J. Fluid Mech. 5, 1–16) and Mellor (1966 J. Fluid Mech. 24, 255–274). In this sub-layer, viscosity and the imposed wall pressure gradient dictate the relevant velocity and length scales. Moreover, when $|\beta _{\textit{in}}| \gtrsim 0.7$ and the wall pressure $P_{w\textit{all}}$ gradient $\textrm { d}P_{w\textit{all}}/\textrm {d}x \gt 0$, a distinct sub-layer emerges outside the pressure–viscous balance region, characterised by a dominant balance between the imposed pressure gradient and the gradient of the Reynolds shear stress. In this region, the Reynolds shear stress increases linearly with distance from the wall. These findings provide new insights into the structure of TBLs under pressure gradients and establish a refined framework for modelling their dynamics.

A new statistical definition for the mean turbulent boundary layer (TBL) thickness is introduced, based on identification of the wall-normal location where the streamwise velocity skewness changes sign, from negative to positive, in the outermost region of the boundary layer. Importantly, this definition is independent of arbitrary thresholds, and broadly applicable, including to past single-point measurements. Furthermore, this definition is motivated by the phenomenology of streamwise velocity fluctuations near the turbulent/non-turbulent interface (TNTI), whose local characteristics are shown to be universal for TBLs under low free-stream turbulence conditions (i.e. with or without pressure gradients, surface roughness, etc.) through large-scale experiments, simulations and coherent structure-based modelling. The new approach yields a TBL thickness that is consistent with previous definitions, such as those based on Reynolds shear stress or ‘composite’ mean velocity profiles, and which can be used practically, e.g. to calculate integral thicknesses. Two methods are proposed for estimating the TBL thickness using this definition: one based on simple linear interpolation and the other on fitting a generalised Fourier model to the outer skewness profile. The robustness and limitations of these methods are demonstrated through analysis of several published experimental and numerical datasets, which cover a range of canonical and non-canonical TBLs. These datasets also vary in key characteristics such as wall-normal resolution and measurement noise, particularly in the critical TNTI region.

In this work, the correlations between streamwise velocity and temperature fluctuations are investigated in compressible turbulent channel flows from the perspective of coherent structures. The intense fluctuation structures and quadrant-event structures of both velocity and temperature have been identified, extracted separately and compared. Analyses show that although their structure sizes are similar in the whole channel, high correlation only exists in the near-wall region with a high overlapping rate of the instantaneous structures. The hierarchy of the temperature structures are passively formed following the dynamic process of the velocity such as ejections, which contributes to the remaining correlation in the outer layer. However, this passive scalar property cannot provide the production mechanism in the outer layer according to the budget analysis after scale decomposition, and the interscale energy transfer progress is also different from the velocity fluctuation field. Therefore, the temperature structures deviate from the velocity structures in the outer layer and cannot be carried by the following dynamic process of the velocity such as sweeps, passively, which can be found from the conditional averaged structures. All of these findings provide a new perspective for understanding the velocity–temperature relationship in compressible channel flows.

Interactions of turbulent boundary layers with a compliant surface are investigated experimentally at Reτ = 3300–8900. Integrating tomographic particle tracking with Mach–Zehnder interferometry enables simultaneous mapping of the compliant wall deformation and the three-dimensional velocity and pressure fields. Our initial study (J. Fluid. Mech. vol. 980, R2) shows that the flow–deformation correlations decrease with increasing Reτ, despite an order of magnitude increase in deformation amplitude. To elucidate the mechanisms involved, the same velocity, pressure and kinetic energy fields are decomposed to ‘wave-coherent’ and ‘stochastic’ parts using a Hilbert projection method. The phase dependent coherent variables, especially the pressure, are highly correlated with the wave, but decrease with increasing Reτ. While the coherent energy is 6 %–10 % of the stochastic level, the pressure root mean square is comparable near the wall. The energy flux between the coherent and stochastic parts and the pressure diffusion reverse sign at the critical layer. To explain the Reτ dependence, the characteristic deformation wavelength (three times the thickness) is compared with the scales of the energy-containing eddies in the boundary layer represented by the k−1 range in the energy spectrum. When the deformation wavelength is matched with the kxEuu peak at the present lowest Reτ, the flow–deformation correlations and coherent pressure become strong, even for submicron deformations. In this case, the flow and wall motion become phase locked, suggesting resonant behaviours. As Reτ increases, the wall wavelengths and spectral range of attached eddies are no longer matched, resulting in reduced correlations and lower coherent energy and pressure, despite larger deformation.

Steady flow at low Reynolds (Re) number through a planar channel with converging or diverging width is investigated in this study. Along the primary direction of flow, the small dimension of the channel cross-section remains constant while the sidewalls bounding the larger dimension are oriented at a constant angle. Due in part to ease of manufacturing, parallel-plate geometries such as this have found widespread use in microfluidic devices for mixing, heat exchange, flow control and flow patterning at small length scales. Previous analytical solutions for flows of this nature have required the converging or diverging aspect of the channel to be gradual. In this work, we derive a matched asymptotic solution, validated against numerical modelling results, that is valid for any sidewall angle, without requiring the channel width to vary gradually. To accomplish this, a cylindrical coordinate system defined by the angle of convergence between the channel sidewalls is considered. From the mathematical form of the composite expansion, a delineation between two secondary flow components emerges naturally. The results of this work show how one of these two components, originating from viscous shear near the channel sidewalls, corresponds to convective mixing, whereas the other component impresses the sidewall geometry on streamlines in the outer flow.

Sink flow boundary layers on smooth and rough walls were studied experimentally. In all cases a turbulent, zero-pressure-gradient boundary layer was subject to acceleration with K = 3.2 × 10–6, which suppressed the turbulence in the outer region and produced conditions similar to those in turbulent sink flow cases with lower K. In the smooth-wall case, after the momentum thickness Reynolds number had dropped to about 600, the near-wall turbulence then dropped, resulting in relaminarisation. In the rough-wall cases, the near-wall turbulence was sustained in spite of the strong favourable pressure gradient, and relaminarisation did not occur. A temporary equilibrium appears to occur that is similar to that seen with lower K, in spite of the ratio of the boundary-layer thickness to the roughness height dropping to less than 5. Mean velocity and Reynolds stress profiles, quadrant analysis and turbulence spectra are used to show the development of the boundary layer in response to the pressure gradient and the differences between the rough- and smooth-wall cases. This is believed to be the first study to consider the spatial evolution of constant-K rough-wall boundary layers with K large enough to cause relaminarisation in the smooth-wall case.

This paper considers the propagation, arrest and recession of a planar hydraulic fracture in a porous elastic medium whose footprint is constrained to a growing or shrinking rectangular region with a constant height. Hydraulic fractures with large aspect ratio rectangular footprints are frequently referred to as PKN fractures in recognition of the original researchers (Perkins & Kern 1961 J. Petrol. Tech. 13, 937–949) and (Nordgren 1972 J. Petrol Technol. 1972, 306–314) who first analyzed models of such fracture geometries. We investigate the one-dimensional non-local PKN approximation to a fully planar rectangular hydraulic fracture model in a three-dimensional elastic medium. By analysing the tip behaviour of the non-local PKN model, a transformation procedure is established to render the asymptotic equations for the dynamics of the steady semi-infinite PKN and plane strain models formally identical, which implies that all the existing multiscale plane strain asymptotes can be converted directly to the PKN case by making use of this transformation. Using this transformation, it is shown that the appropriate PKN asymptotes for the average aperture $\bar {w}$ with distance $\hat {x}$ to the fracture front are $\bar {w}\sim \hat {x}^{1/2},\ \hat {x}^{5/8}\ {\textrm{and}}\, \ \hat {x}^{2/3}$ in the toughness, leak-off and viscous modes of propagation, respectively; as well as the linear elastic fracture mechanics tip asymptote $\bar {w}\sim \hat {x}^{1/2}$ for arrest, which transitions to the linear asymptote tip $\bar {w}\sim \hat {x}$ for a fracture driven to recede due to fluid leak-off. Both the arrest and recession tip asymptotes share the intermediate leak-off asymptote $\bar {w}\sim \hat {x}^{3/4}$. A scaling analysis yields the arrest time, length and aperture as functions of a dimensionless injection-cessation time $\omega$. An asymptotic analysis of the non-local PKN model is used to establish the fundamental decoupling between dynamics and kinematics, which leads to the emergence of a similarity solution – termed the sunset solution – close to the time of collapse of the fracture. The multiscale PKN numerical solutions agree well with those for a fully planar multiscale rectangular hydraulic fracture model in a three-dimensional elastic medium. The scaling laws and the emergence of the sunset solution are confirmed by the PKN numerical model. The sunset solution also emerges in the fully planar numerical model and persists beyond the collapse time of the PKN model, by which time its footprints have separated from the upper and lower constraining sedimentary layer boundaries and have assumed self-similar elliptic shapes that shrink as they approach collapse.

In the rapidly rotating limit, we derive a balanced set of reduced equations governing the strongly nonlinear development of the convective wall-mode instability in the interior of a general container. The model illustrates that wall-mode convection is a multiscale phenomenon where the dynamics of the bulk interior diagnostically determine the small-scale dynamics within Stewartson boundary layers at the sidewalls. The sidewall boundary layers feedback on the interior via a nonlinear lateral heat-flux boundary condition, providing a closed system. Outside the asymptotically thin boundary layer, the convective modes connect to a dynamical interior that maintains scales set by the domain geometry. In many ways, the final system of equations resembles boundary-forced planetary geostrophic baroclinic dynamics coupled with barotropic quasi-geostrophic vorticity. The reduced system contains the results from previous linear instability theory but captured in an elementary fashion, providing a new avenue for investigating wall-mode convection in the strongly nonlinear regime. We also derive the dominant Ekman-flux correction to the onset Rayleigh number for large Taylor number, ${\textit {Ra}} \approx 31.8 \,{\textit{Ta}}^{1/2} - 4.43 \,{\textit{Ta}}^{5/12} + {\mathcal{O}}({\textit{Ta}}^{1/3})$ for no-slip boundaries. However, we find that the linear onset in a finite cylinder differs noticeably compared with a Cartesian channel. We demonstrate some of the reduced model’s nonlinear dynamics with numerical simulations in a cylindrical container.

A model is proposed for the one-dimensional spectrum and streamwise Reynolds stress in pipe flow for arbitrarily large Reynolds numbers. Constructed in wavenumber space, the model comprises four principal contributions to the spectrum: streaks, large-scale motions, very-large-scale motions and incoherent turbulence. It accounts for the broad and overlapping spectral content of these contributions from different eddy types. The model reproduces well the broad structure of the premultiplied one-dimensional spectrum of the streamwise velocity, although the bimodal shape that has been observed at certain wall-normal locations, and the $-5/3$ slope of the inertial subrange, are not captured effectively because of the simplifications made within the model. Regardless, the Reynolds stress distribution is well reproduced, even within the near-wall region, including key features of wall-bounded flows such as the Reynolds number dependence of the inner peak, the formation of a logarithmic region, and the formation of an outer peak. These findings suggest that many of these features arise from the overlap of energy content produced by both inner- and outer-scaled eddy structures combined with the viscous-scaled influence of the wall. The model is also used to compare with canonical turbulent boundary layer and channel flows, and despite some differences being apparent, we speculate that with only minor modifications to its coefficients, the model can be adapted to these flows as well.

Experiments were performed that (i) document the effect of the steady spanwise buffer layer blowing on the mean characteristics of the turbulent boundary layer for a range of momentum thickness Reynolds numbers from 4760 to 10 386, and (ii) document the effect of the buffer layer blowing on the unsteady characteristics and coherent vorticity in a boundary layer designed to provide sufficiently high spatial resolution. The spanwise buffer layer blowing of the order of $u_{\tau }$ is produced by a surface array of pulsating direct current (pulsed-DC) plasma actuators. This was found to substantially reduce the wall shear stress that was directly measured with a floating element coupled with a force sensor. The direct wall shear measurements agreed with values derived using the Clauser method to within $\pm 0.85$ %. The degree to which the buffer layer blowing affected $\tau _w$ was found to primarily depend on the inner variable spanwise spacing between the pulsed-DC actuator electrodes, i.e. ‘blowing sites’. Utilizing pairs of $[u,v]$ and $[u,w]$ hot-wire sensors, the latter experiments correlated significant reductions in the $\omega _y$ and $\omega _x$ vorticity components that resulted from the buffer layer blowing and translated into lower Reynolds stresses and turbulence production. The time scale to which these observed changes in the boundary layer characteristics would return to the baseline condition was subsequently documented. This revealed a recovery length of $x^+ \approx 86\,000$ that translated to a streamwise fetch of $x \approx 66\delta$. Finally, a comparison with the recent work by Cheng et al. (2021, J. Fluid Mech. vol. 918, A24) and Wei & Zhou (2024 in TSFP13, June 25–28, 2024) that followed our experimental approach to achieve comparable wall shear stress (drag) reductions has led to a new scaling based on the baseline boundary layer $\textit{Re}_{\tau }$ and buffer layer blowing velocity.

In this study, changes in the mean flow of a compressible turbulent boundary layer spatially evolving from low to ‘moderate’ Reynolds numbers are examined. All discussions are based on literature data and a direct numerical simulation (DNS) of a supersonic boundary layer specifically designed to be effectively free of spurious inflow effects in the range $4000 \lessapprox Re_\theta \lessapprox 5000$, which enables discussion of sensitive properties such as the turbulent wake. Most noticeably, the DNS data show the formation of a distinct ‘bend’ in the friction coefficient distribution reflected in sudden deviation from established low-Reynolds-number correlations. As will be shown, the bend is related to the surprisingly abrupt saturation of the turbulent wake, marking the change from low- to moderate-Reynolds-number behaviour; in previous studies, this trend was potentially obscured by data scatter in experiments and/or insufficient domain length in DNS. Moreover, the influence of the wake saturation on the formation of the early logarithmic overlap layer is assessed, which, if fully developed, leads to the onset of high-Reynolds-number behaviour further downstream.

Understanding the vertical coherence of the pressure structure and its interaction with velocity fields is critical for elucidating the mechanisms of acoustic generation and radiation in hypersonic turbulent boundary layers. This study employs linear coherence analysis to examine the self-similar coherent structures in the velocity and pressure fields within a Mach 6 hypersonic boundary layer, considering a range of wall-to-recovery temperature ratios. The influence of wall cooling on the geometric characteristics of these structures, such as inclination angles and three-dimensional aspect ratios, is evaluated. Specifically, the streamwise velocity exhibits self-similar coherent structures with the streamwise/wall-normal aspect ratio ranging from 16.5 to 38.7, showing a linear increases with decreasing wall temperatures. Similar linear dependence between the streamwise/wall-normal aspect ratio and the wall temperatures are observed for the Helmholtz-decomposed streamwise velocity and the pressure field. In terms of velocity–pressure coupling, the solenoidal component exhibits stronger interactions with the pressure fields in the near-wall region, while the dilatational component has stronger interactions with the pressure field at large scales with the increase of height. Such coupling generally follows the distance-from-the-wall scaling of the pressure field, except in cooled wall cases. Using the linear stochastic estimation, the pressure field across the boundary layer is predicted by inputting the near-wall pressure/velocity signal along with the transfer kernel. The result demonstrates that near-wall pressure signals provide the most accurate description of the pressure field in higher regions of the boundary layer. As wall-mounted sensors can measure near-wall pressure fluctuations, this study presents a potential approach to predict the off-wall pressure field correlated with the near-wall structures based on wall-pressure measurements.

We define ‘surface layer’ (SL) as an inertia-dominated turbulence region outside a viscous or roughness surface-adjacent sub-layer (SAS) that is characterised by linear scaling of specific coherence length scales on wall-normal distance, $z$. We generalise the mechanisms that underlie the formation of the classical inertial SL in the shear-dominated turbulent boundary layer (TBL) to wall-bounded turbulent flows with zero mean shear. Using particle image velocimetry data from two wind tunnel facilities, we contrast the classical TBL SL with a non-classical shear-free SL generated within grid turbulence advected over an impermeable plate using two grids with different turbulence length scales. Integral-scale variations with $z$ and other statistics are quantified. In both shear-dominated and shear-free SLs we observe well-defined linear increases in $z$ of the streamwise integral scale of vertical velocity fluctuations. In grid turbulence the shear-free SL initiates just above the SAS that confines friction-generated motions. By contrast, the TBL SL forms with non-zero mean shear rate that extends streamwise coherence lengths of streamwise fluctuations. In both flow classes only the integral scales of vertical fluctuating velocity increase linearly with $z$, indicating that the SL is generated by the blockage of vertical fluctuations in the vertical. Whereas the SAS in the TBL is much thinner than in the grid-turbulence flows, the generation of a shear-free SL by the interaction of turbulence eddies and a surface depends on the relative thinness of the SAS. We conclude that the common generalisable SL mechanism is direct blockage of vertical fluctuations by the impermeable surface.