In this paper, we derive the effective model describing a thin-domain flow with permeable boundary through which the fluid is injected into the domain. We start with incompressible Stokes system and perform the rigorous asymptotic analysis. Choosing the appropriate scaling for the injection leads to a compressible effective model. In this paper, we derive the effective model describing a thin-domain flow with permeable boundary through which the fluid is injected into the domain. We start with incompressible Stokes system and perform the rigorous asymptotic analysis. Choosing the appropriate scaling for the injection leads to a compressible effective model.

The present study tests the efficacy of the well-known viscous drag reduction strategy of imposing spanwise wall oscillations to reduce pressure drag contributions in transitional and fully rough turbulent wall flow. This is achieved by conducting a series of direct numerical simulations of a turbulent flow over two-dimensional (spanwise-aligned) semi-cylindrical rods, placed periodically along the streamwise direction with varying streamwise spacing. Surface oscillations, imposed at fixed viscous-scaled actuation parameters optimum for smooth wall drag reduction, are found to yield substantial drag reduction ($\gtrsim$25 %) for all the rough wall cases, maintained at matched roughness Reynolds numbers. While the total drag reduction is due to a drop in both viscous and pressure drag in the case of transitionally rough flow (i.e. with large inter-rod spacing), it is associated solely with pressure drag reduction for the fully rough cases (i.e. with small inter-rod spacing), with the latter being reported for the first time. The study finds that pressure drag reduction in all cases is caused by the attenuation of the vortex shedding activity in the roughness wake, in response to wall oscillation frequencies that are of the same order as the vortex shedding frequencies. Contrary to speculations in the literature, this study confirms that the mechanism behind pressure drag reduction, achieved via imposition of spanwise oscillations, is independent of the viscous drag reduction. This mechanism is responsible for weakening of the Reynolds stresses and increase in base pressure in the roughness wake, explaining the pressure drag reduction observed by past studies, across varying roughness heights and geometries.

We investigate the role of inter-scale interactions in the high-Reynolds-number skin-friction drag reduction strategy reported by Marusic et al. (Nat. Commun., vol. 12, 2021). The strategy involves imposing relatively low-frequency streamwise travelling waves of spanwise velocity at the wall to actuate the drag generating outer scales. This approach has proven to be more energy efficient than the conventional method of directly targeting the drag producing inner scales, which typically requires actuation at higher frequencies. Notably, it is observed that actuating the outer scales at low frequencies leads to a substantial attenuation of the major drag producing inner scales, suggesting that the actuations affect the nonlinear inner–outer coupling inherently existing in wall-bounded flows. In the present study, we find that increased drag reduction, through imposition of spanwise wall oscillations, is always associated with an increased coupling between the inner and outer scales. This enhanced coupling emerges through manipulation of the phase relationships between these triadically linked scales, with the actuation forcing the entire range of energy-containing scales, from the inner (viscous) to the outer (inertial) scales, to be more in phase. We also find that a similar enhancement of this nonlinear coupling, via manipulation of the inter-scale phase relationships, occurs with increasing Reynolds number for canonical turbulent boundary layers. This indicates improved efficacy of the energy-efficient drag reduction strategy at very high Reynolds numbers, where the energised outer scales are known to more strongly superimpose and modulate the inner scales. Leveraging the inter-scale interactions, therefore, offers a plausible mechanism for achieving energy-efficient drag reduction at high Reynolds numbers.

We present experimental evidence that the superstructures in turbulent boundary layers comprise smaller, geometrically self-similar coherent motions. The evidence comes from identifying and analysing instantaneous superstructures from large-scale particle image velocimetry datasets acquired at high Reynolds numbers, capable of capturing streamwise elongated motions extending up to 12 times the boundary layer thickness. Given the challenge in identifying the constituent motions of the superstructures based on streamwise velocity signatures, a new approach is adopted that analyses the wall-normal velocity fluctuations within these very long motions, which reveals the constituent motions unambiguously. The conditional streamwise energy spectra of the Reynolds shear stress and the wall-normal fluctuations, corresponding exclusively to the superstructure region, are found to exhibit the well-known distance-from-the-wall scaling in the intermediate-scale range. It suggests that geometrically self-similar motions are the constituent motions of these very-large-scale structures. Investigation of the spatial organization of the wall-normal momentum-carrying eddies, within the superstructures, also lends empirical support to the concatenation hypothesis for the formation of these structures. The association between the superstructures and self-similar motions is reaffirmed on comparing the vertical coherence of the Reynolds-shear-stress-carrying motions, by computing conditionally averaged two-point correlations, which are found to match with the mean correlations. The mean vertical coherence of these motions, investigated for the log region across three decades of Reynolds numbers, exhibits a unique distance-from-the-wall scaling invariant with Reynolds number. The findings support modelling of these dynamically significant motions via data-driven coherent structure-based models.

Turbulent drag reduction (DR) through streamwise travelling waves of the spanwise wall oscillation is investigated over a wide range of Reynolds numbers. Here, in Part 1, wall-resolved large-eddy simulations in a channel flow are conducted to examine how the frequency and wavenumber of the travelling wave influence the DR at friction Reynolds numbers $Re_\tau = 951$ and $4000$. The actuation parameter space is restricted to the inner-scaled actuation (ISA) pathway, where DR is achieved through direct attenuation of the near-wall scales. The level of turbulence attenuation, hence DR, is found to change with the near-wall Stokes layer protrusion height $\ell _{0.01}$. A range of frequencies is identified where the Stokes layer attenuates turbulence, lifting up the cycle of turbulence generation and thickening the viscous sublayer; in this range, the DR increases as $\ell _{0.01}$ increases up to $30$ viscous units. Outside this range, the strong Stokes shear strain enhances near-wall turbulence generation leading to a drop in DR with increasing $\ell _{0.01}$. We further find that, within our parameter and Reynolds number space, the ISA pathway has a power cost that always exceeds any DR savings. This motivates the study of the outer-scaled actuation pathway in Part 2, where DR is achieved through actuating the outer-scaled motions.

We present measurements of turbulent drag reduction (DR) in boundary layers at high friction Reynolds numbers in the range of $4500 \le Re_\tau \le 15\ 000$. The efficacy of the approach, using streamwise travelling waves of spanwise wall oscillations, is studied for two actuation regimes: (i) inner-scaled actuation (ISA), as investigated in Part 1 of this study, which targets the relatively high-frequency structures of the near-wall cycle, and (ii) outer-scaled actuation (OSA), which was recently presented by Marusic et al. (Nat. Commun., vol. 12, 2021) for high-$Re_\tau$ flows, targeting the lower-frequency, outer-scale motions. Multiple experimental techniques were used, including a floating-element balance to directly measure the skin-friction drag force, hot-wire anemometry to acquire long-time fluctuating velocity and wall-shear stress, and stereoscopic particle image velocimetry to measure the turbulence statistics of all three velocity components across the boundary layer. Under the ISA pathway, DR of up to 25 % was achieved, but mostly with net power saving (NPS) losses due to the high-input power cost associated with the high-frequency actuation. The low-frequency OSA pathway, however, with its lower input power requirements, was found to consistently result in positive NPS of 5–10 % for moderate DRs of 5–15 %. The results suggest that OSA is an attractive pathway for energy-efficient DR in high-Reynolds-number applications.

We study the asymptotic behaviour of the periodically mixed Zaremba problem. We cover the part of the boundary by a chess board with a small period (square size) $\varepsilon$ and impose the Dirichlet condition on black and the Neumann condition on white squares. As $\varepsilon \to 0$, we get the effective boundary condition which is always of the Dirichlet type. The Dirichlet data on the boundary, however, depend on the ratio between the magnitudes of the two boundary values.

In this study, we develop an analytical model to predict the turbulent boundary layer downstream of a step-change in the surface roughness where upstream flow conditions are given. We first revisit the classical model of Elliott (Trans. Am. Geophys. Union, vol. 39, 1958, pp. 1048–1054), who modelled the velocity distribution within and above the internal layer with a simple piecewise logarithmic profile, and evolved the velocity profile using the streamwise momentum equation. Elliott's model was originally developed for an atmospheric surface layer, and to make the model applicable to a spatially developing turbulent boundary layer with finite thickness, we propose a number of more physical refinements, including adding a wake function to the velocity profile, considering the growth of the entire boundary layer in the streamwise direction, and using a more realistic shear stress profile in the momentum equation. In particular, we implement the blending model (Li et al., J. Fluid Mech., vol. 923, 2021, p. A18) to account for the deviation of the mean flow within the internal layer from a canonical velocity profile based on the local wall condition. These refinements lead to improved agreement between the prediction and the measurement, especially in the vicinity of the rough-to-smooth change.

A large-scale spanwise and wall-normal array of sonic anemometers in the atmospheric surface layer is used to acquire all three components of instantaneous fluctuating velocity as well as temperature in a range of stability conditions. These data permit investigation of the three-dimensional statistical structure of turbulence structures. Based on a similar dataset, Krug et al. (Boundary-Layer Meteorol., vol. 172, 2019, pp. 199–214) reported a self-similar range of wall-attached turbulence structures under both unstable and near-neutral stability conditions. They considered only a wall-normal array and thus assessed statistical structure in the wall-normal direction, in relation to the streamwise wavelength. The present work extends the view of a self-similar range of turbulence structures, by including the statistical structure in the spanwise direction. Moreover, by analysing the phase shift between synchronized measurements in the spectral domain, it is inferred how a scale-dependent inclination angle in the streamwise/wall-normal plane varies with stability. Results suggest that the self-similar wall-attached structures have similar aspect ratios between streamwise/wall-normal scales and streamwise/spanwise scales such that $\lambda _x/\Delta z : \lambda _x/\Delta y \approx 1$ for both near-neutral and unstable conditions. Under the most unstable conditions, coherent structures with $\lambda _x/\delta = 1$ are inclined at angles as high as $65^{\circ }$ relative to the solid boundary, while larger scales $\lambda _x/\delta = 6$ exhibit inclination angles of approximately $35^{\circ }$. For near-neutral stability conditions, the angle tends towards $12^{\circ }$ for all scales. It is noted that in the near-neutral condition, the structure inclination angle and the aspect ratio – and thus the statistical modelling of coherent structures in the atmospheric surface layer – are highly sensitive to the value of the stability parameter.