We present direct numerical simulation (DNS) and modelling of incompressible, turbulent, generalized Couette–Poiseuille flow. A particular example is specified by spherical coordinates $(Re,\theta,\phi )$, where $Re = 6000$ is a global Reynolds number, $\phi$ denotes the angle between the moving plate, velocity-difference vector and the volume-flow vector and $\tan \theta$ specifies the ratio of the mean volume-flow speed to the plate speed. The limits $\phi \to 0^\circ$ and $\phi \to 90^\circ$ give alignment and orthogonality, respectively, while $\theta \to 0^\circ,\ \theta \to 90^\circ$ correspond respectively to pure Couette flow in the $x$ direction and pure Poiseuille flow at angle $\phi$ to the $x$ axis. Competition between the Couette-flow shear and the forced volume flow produces a mean-velocity profile with directional twist between the confining walls. Resultant mean-speed profiles relative to each wall generally show a log-like region. An empirical flow model is constructed based on component log and log-wake velocity profiles relative to the two walls. This gives predictions of four independent components of shear stress and also mean-velocity profiles as functions of $(Re,\theta,\phi )$. The model captures DNS results including the mean-flow twist. Premultiplied energy spectra are obtained for symmetric flows with $\phi =90^\circ$. With increasing $\theta$, the energy peak gradually moves in the direction of increasing $k_x$ and decreasing $k_z$. Rotation of the energy spectrum produced by the faster moving velocity near the wall is also observed. Rapid weakening of a spike maxima in the Couette-type flow regime indicates attenuation of large-scale roll structures, which is also shown in the $Q$-criterion visualization of a three-dimensional time-averaged flow.

High-fidelity simulations are conducted to investigate the turbulent boundary layers around a finite-span NACA0012 wing with a rounded wing-tip geometry at a chord-based Reynolds number of $Re_c=200\,000$ and at various angles of attack up to $10^\circ$. The study aims to discern the differences between the boundary layers on the finite-span wing and those on infinite-span wings at equivalent angles of attack. The finite-span boundary layers exhibit: (i) an altered streamwise and a non-zero spanwise pressure gradient as a result of the variable downwash induced by the wing-tip vortices (an inviscid effect typical of finite-span wings); (ii) differences in the flow history at different wall-normal distances, caused by the variable flow angle in the wall-normal direction (due to constant pressure gradients and variable momentum normal to the wall); (iii) laminar flow entrainment into the turbulent boundary layers near the wing tip (due to a laminar–turbulent interface); and (iv) variations in boundary layer thickness across the span, attributed to the variable wall-normal velocity in that direction (a primarily inviscid effect). These physical effects are then used to explain the differences in the Reynolds stress profiles and other boundary layer quantities, including the reduced near-wall peak of the streamwise Reynolds stress and the elevated Reynolds stress levels near the boundary layer edge, both observed in the finite-span wings. Other aspects of the flow, such as the downstream development of wing-tip vortices and their interactions with the surrounding flow, are reserved for future investigations.

The first two parts of this study showed by experiment and by large-eddy simulation that, in mixing layers formed of processions of spanwise-oriented vortex structures, there is a fundamental change in the dynamics of the large-structure growth at what has been called the ‘mixing’ transition. This third part examines the relationship of the post-transition ‘organised’ flow to the mixing layer of classical turbulence theory. Further large-eddy simulations are presented which, like some of the experiments reported in Part 1, have captured the mixing layer in both its classical and organised turbulent states, thereby allowing them to be characterised separately. The simulation results are then used to calibrate model-free integral analyses of the dynamics of the mixing layer's spatial growth. It is shown that the organised and classical flows are alternative self-preserving turbulent flow states involving fundamentally different exchanges of mass, momentum and energy, one or other of which emerges naturally in particular realisations depending on the initial conditions.

Two-dimensional (2-D) quadrant analysis is generally used for investigating flow and sediment dynamics around a rigid structure in open channel flows. Given that particle distribution around rigid obstacles is not spatially uniform and changes in time, while vortices evolve to become three-dimensional (3-D) structures, 2-D quadrant analysis might be unsuitable to completely determine the sediment transport. Hence, 3-D quadrant and 3-D octant analyses should be considered, using the 3-D instantaneous velocity data and relative 3-D bursting process to define sediment transport surrounding the submerged square and circular cylinders. The turbulent kinetic energy (TKE), transition probabilities, occurrence probabilities, stress fraction and angles of inclination of 3-D bursting events are considered to quantify the coherent structures surrounding the cylinders and their interaction with bed particles. Experiments were conducted at the Hydraulics and Water Resources Engineering Laboratory, School of Infrastructure, Indian Institute of Technology, Bhubaneswar, Odisha, India, and velocity data were recorded at different cross-sections around the submerged cylinders using an acoustic Doppler velocimeter. Results show that the TKE is greater for internal ejection, external ejection, and internal sweep, external sweep in the upstream of the circular and square cylindrical structures. On the other part, the TKE is significant for internal ejection, external ejection, and internal sweep, external sweep in the downstream of the aligned square cylindrical structures, which justifies the highest scour depth that occurred upstream of the circular and square cylindrical structures and downstream of the aligned square cylindrical structure. The transition probability of the bursting events was determined using the Markov process from the measured velocity data to investigate the consecutive occurrence of bursting events. Further, the importance of sweeps and ejections on sediment erosion surrounding the cylinders within the scour hole at various stages of its development was investigated via 3-D quadrant analysis of the bursting occurrences. The outcomes show that external sweep and internal ejection events are active mechanisms for bed particle transport surrounding the cylinder. The maximum transition probability values are found around aligned cylindrical structures in comparison with the circular and square cylindrical structures in the transverse direction. This depicts the formation of a trailing vortex on both sides of the aligned square cylindrical object. The results reveal that the effect of inclination angles with respect to the water flow is greater for internal ejection and external sweep from upstream to downstream within the scour hole surrounding the cylindrical structures at various phases of development as horseshoe vortices and downflow develop upstream of the cylindrical structures while trailing vortices and wake vortices form at the top and downstream of the cylindrical structures. Internal and external ejection have a higher stress fraction than an internal and external sweep for square cylinders with alignment angles of 0°, 20° and circular cylinders over underdeveloped and developed scoured beds, respectively. With the higher percentage of fractional contributions for internal sweeps, the external sweep is predicted close to the cylindrical objects in comparison with the internal ejection and external ejection events because of the formation and warping of the horseshoe vortex close to the cylindrical objects, suggesting a significant probability of 3-D bursting occurrences with sediment movement near the cylindrical structures.

The Reynolds-averaged Navier–Stokes (RANS) models depend on empirical constants to close the Reynolds stress terms. The empirical constants were obtained using experiments conducted at low Reynolds numbers several decades ago. In this paper, we revisit the turbulent viscosity parameter $C_\mu$, based on the stress–intensity ratio $c^2 = {|\overline {uw}|}/{k}$. Here, $\overline {|uw|}$ and $k$ are the absolute values of the Reynolds stress and turbulent kinetic energy, respectively. Through a priori comparisons, we find that the currently accepted value of $C_\mu = 0.09$ does not agree with the latest direct numerical simulation (DNS) and experimental datasets of wall-bounded turbulent planar flows. Therefore, a new value is suggested by averaging $c^2$ in the equilibrium region, where the production ($\mathcal {P}$) of $k$ is within 10 % of the dissipation rate ($\epsilon$), and consequently, $c^4 \approx C_\mu$. We evaluate flows up to friction Reynolds number $Re_\tau \approx 10\,000$ and find that with increasing $Re_\tau$, $C_\mu$ approaches a value of 0.06, which is almost 50 % lower than the prevalent value of 0.09. Finally, we perform an a priori test with the new (proposed) value of $C_\mu = 0.06$ to show that the estimated turbulent viscosity $\nu _T$ for wall-bounded flows is in much closer agreement with the exact (DNS) values than when $\nu _T$ is estimated using $C_\mu = 0.09$.

Large-scale coherent structures in incompressible turbulent pipe flow are studied for a wide range of Reynolds numbers ($Re_\tau =180, 550, 1000, 2000$ and $5200$). Employing the Karhunen–Loève decomposition and a novel approach based on the Voronoi diagram, we identify and classify statistically coherent structures based on their location, dimensions and $Re_{\tau }$. With increasing $Re_{\tau }$, two distinct classes of structures become more energetic, namely wall-attached and detached eddies. The Voronoi methodology is shown to delineate these two classes without the need for specific criteria or thresholds. At the highest $Re_{\tau }$, the attached eddies scale linearly with the wall-normal distance with a slope of approximately $l_y\sim 1.2y/R$, while the detached eddies remain constant at the size of $l_y \approx 0.26R$, with a progressive shift towards the pipe centre. We extract these two classes of structures and describe their spatial characteristics, including radial size, helix angle and azimuthal self-similarity. The spatial distribution could help explain the differences in mean velocity between pipe and channel flows, as well as in modelling large and very-large-scale motions (LSM and VLSM). In addition, a comprehensive description is provided for both wall-attached and detached structures in terms of LSM and VLSM.

We use three-dimensional direct numerical simulations of homogeneous isotropic turbulence in a cubic domain to investigate the dynamics of heavy, chiral, finite-size inertial particles and their effects on the flow. Using an immersed-boundary method and a complex collision model, four-way coupled simulations have been performed, and the effects of particle-to-fluid density ratio, turbulence strength and particle volume fraction have been analysed. We find that freely falling particles on the one hand add energy to the turbulent flow but, on the other hand, they also enhance the flow dissipation: depending on the combination of flow parameters, the former or the latter mechanism prevails, thus yielding enhanced or weakened turbulence. Furthermore, particle chirality entails a preferential angular velocity which induces a net vorticity in the fluid phase. As turbulence strengthens, the energy introduced by the falling particles becomes less relevant and stronger velocity fluctuations alter the solid phase dynamics, making the effect of chirality irrelevant for the large-scale features of the flow. Moreover, comparing the time history of collision events for chiral particles and spheres (at the same volume fraction) suggests that the former tend to entangle, in contrast to the latter which rebound impulsively.

A high-resolution numerical simulation of an air–water turbulent upward bubbly flow in a pipe is performed to investigate the turbulence characteristics and bubble interaction with the wall. We consider three bubble equivalent diameters and three total bubble volume fractions. The bulk and bubble Reynolds numbers are $Re_{bulk}= u_{bulk} D/\nu _w = 5300$ and $Re_{bub}= (\langle u_{bub}\rangle - u_{bulk}) d_{eq}/\nu _w = 533\unicode{x2013}1000$, respectively, where $u_{bulk}$ is the water bulk velocity, $\langle u_{bub}\rangle$ is the overall bubble mean velocity, $D$ is the pipe diameter and $\nu _w$ is the water kinematic viscosity. The mean water velocity near the wall significantly increases due to bubble interaction with the wall, and the root-mean-square water velocity fluctuations are proportional to $\bar {\psi }(r)^{0.4}$, where $\bar {\psi } (r)$ is the mean bubble volume fraction. For the cases considered, the bubble-induced turbulence suppresses the shear-induced turbulence and becomes the dominant flow characteristic at all radial locations including near the wall. Rising bubbles near the wall mostly bounce against the wall rather than slide along the wall or hang around the wall without collision. Low-speed streaks observed in the near-wall region in the absence of bubbles nearly disappear due to the bouncing bubbles. These bouncing bubbles generate counter-rotating vortices in their wake, and increase the skin friction by sweeping high-speed water towards the wall. We also suggest an algebraic Reynolds-averaged Navier–Stokes model considering the interaction between shear-induced and bubble-induced turbulence. This model provides accurate predictions for a wide range of liquid bulk Reynolds numbers.

The influence of outer large-scale motions (LSMs) on near-wall structures in compressible turbulent channel flows is investigated. To separate the compressibility effects, velocity fluctuations are decomposed into solenoidal and dilatational components using the Helmholtz decomposition method. Solenoidal velocity fluctuations manifest as near-wall streaks and outer large-scale structures. The spanwise drifting of near-wall solenoidal streaks is found to be driven by the outer LSMs, while LSMs have a trivial influence on the spanwise density of solenoidal streaks, consistent with the outer LSM impacts found in incompressible flows (Zhou et al., J. Fluid Mech., vol. 940, 2022, p. A23). Dilatational motions are characterized by the near-wall small-scale travelling-wave packets and the large-scale parts in the outer region. The streamwise advection velocity of the near-wall structures remains at $16 \sim 18u_{\tau }$, hardly influenced by Mach numbers, Reynolds numbers and wall temperatures. The spanwise drifting of near-wall dilatational structures, quantified by the particle image velocimetry method, follows a mechanism distinct from solenoidal streaks. This drifting velocity is notably larger than those of the solenoidal streaks, and the influence of outer LSMs is not the primary trigger for this drifting.

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.

Comprehensive coherent structures around a surface-mounted low aspect ratio square cylinder in uniform flow with an oblique angle of $45^{\circ }$ were investigated for cylinder-width-based Reynolds numbers of 3000 and 10 000 by direct numerical simulation based on a topology-confined mesh refinement framework. High-resolution simulations and the critical-point concept were scrutinized to reveal for the first time the reasonable and compatible topologies of flow separation and complete near-wall structures, due to their extensive impact on various engineering applications. Large-scale horseshoe vortices are observed at two notable foci in the viscous sublayer. Within this layer, a wall-parallel jet is formed by downflow intruding into the bottom surface at the half-saddle point, then deflecting in the upstream direction and finally penetrating the bottom surface until the half-saddle point. A pair of conical vortices on the cylinder's top surface switch themselves on two sides of the square cylinder, where the switching frequency is identical with that of the sway of the side shear layer. The undulation of the Kelvin–Helmholtz instability is identified in the instantaneous development of a conical vortex and side shear layer, where the scaling of the ratio of the Kelvin–Helmholtz and von Kármán frequencies follows the power-law relation obtained by Lander et al. (J. Fluid Mech., vol. 849, 2018, pp. 1096–1119). Large-scale arch-shaped vortex is often detected in the intermediate wake region of a square cylinder, involving two interconnected portions, such as the leg portion separated from leeward surfaces and head portion rolled up from the top surface. The leg portion of the arch-shaped vortex was rooted by two foci near the bottom-surface plane.

We study the behaviour of the streamwise velocity variance in turbulent wall-bounded flows using a direct numerical simulation (DNS) database of pipe flow up to friction Reynolds number ${{Re}}_{\tau } \approx 12000$. The analysis of the spanwise spectra in the viscous near-wall region strongly hints to the presence of an overlap layer between the inner- and the outer-scaled spectral ranges, featuring a $k_{\theta }^{-1+\alpha }$ decay (with $k_{\theta }$ the wavenumber in the azimuthal direction, and $\alpha \approx 0.18$), hence shallower than suggested by the classical formulation of the attached-eddy model. The key implication is that the contribution to the streamwise velocity variance $(\langle{u}^2\rangle)$ from the largest scales of motion (superstructures) slowly declines as ${{Re}}_{\tau }^{-\alpha }$, and the integrated inner-scaled variance follows a defect power law of the type $\langle u^2 \rangle ^+ = A - B \, {{Re}}_{\tau }^{-\alpha }$, with constants $A$ and $B$ depending on $y^+$. The DNS data very well support this behaviour, which implies that strict wall scaling is restored in the infinite-Reynolds-number limit. The extrapolated limit distribution of the streamwise velocity variance features a buffer-layer peak value of $\langle u^2 \rangle ^+ \approx 12.1$, and an additional outer peak with larger magnitude. The analysis of the velocity spectra also suggests a similar behaviour of the dissipation rate of the streamwise velocity variance at the wall, which is expected to attain a limiting value of approximately $0.28$, hence slightly exceeding the value $0.25$ which was assumed in previous analyses (Chen & Sreenivasan, J. Fluid Mech., vol. 908, 2021, R3). We have found evidence suggesting that the reduced near-wall influence of wall-attached eddies is likely linked to the formation of underlying turbulent Stokes layers.

Flexible canopy flows are often encountered in natural scenarios, e.g. when crops sway in the wind or when submerged kelp forests are agitated by marine currents. Here, we provide a detailed characterisation of the turbulent flow developed above and between the flexible filaments of a fully submerged dense canopy and we describe their dynamical response to the turbulent forcing. We investigate a wide range of flexibilities, encompassing the case in which the filaments are completely rigid and standing upright as well as that where they are fully compliant to the flow and deflected in the streamwise direction. We are thus able to isolate the effect of the canopy flexibility on the drag and on the inner–outer flow interactions, as well as the two flapping regimes of the filaments already identified for a single fibre. Furthermore, we offer a detailed description of the Reynolds stresses throughout the wall-normal direction resorting to the Lumley triangle formalism, and we show the multi-layer nature of turbulence inside and above the canopy. The relevance of our investigation is thus twofold: the fundamental physical understanding developed here paves the way towards the investigation of more complex and realistic scenarios, while we also provide a thorough characterisation of the turbulent state that can prove useful in the development of accurate turbulence models for RANS and LES.

A turbulent two-vortex system (T-2VS) is obtained by inserting analytical model wake vortices into very weak homogeneous isotropic turbulence (HIT) and by evolving them in time using large-eddy simulation until a turbulent state at statistical equilibrium is reached. The T-2VS is characterised as follows: circulation distribution of the vortices; energy of the mean and fluctuating fields; energy dissipation rate. It is also verified that essentially the same T-2VS is obtained when varying the initial model or initial HIT perturbation. A wall-resolved simulation of the T-2VS further interacting with a smooth ground is then performed at $Re_\varGamma = 2 \times 10^5$; this is $10 \times$ higher than in previous works, which allows us to better capture the high Reynolds number behaviour. The high release height of the T-2VS also ensures a physically correct approach to the ground. The results are compared with the literature and also to what is obtained for the case of non-turbulent vortices interacting with the same ground at the same Reynolds number. The flow topologies are discussed, and significant differences are highlighted regarding the separation of the boundary layer generated at the ground, and the way this secondary vorticity interacts with the primary vortices and makes them decay. The vortex trajectories are also measured, together with their circulation distribution and global circulation evolution, and the differences are discussed.

Polymer chains in turbulent flows are generally modelled as dumbbells, i.e. two beads joined by a nonlinear spring. The dumbbell only maps a single spatial configuration, described by the polymer end-to-end vector, thus a multi-bead FENE (finitely extensible nonlinear elastic) chain seems a natural improvement for a more accurate characterisation of the polymer spatial conformation. At a large Weissenberg number, a comparison with the more accurate Kuhn chain reveals that the multi-bead FENE chain drastically overestimates the probability of folded configurations. Surprisingly, the dumbbell turns out to be the only meaningful bead-spring model to coarse-grain a polymer macromolecule in turbulent pipe flows.

Considerable efforts have been devoted to the understanding of the small-scale characteristics in turbulent flows. While the universality of small-scale quantities has been established for incompressible flows, their extension to high-pressure transcritical flows remains an open area of research. To address this question, we investigate the real-fluid thermodynamic effects on small-scale velocity statistics of high-pressure transcritical wall-bounded turbulence. We show that in the locally isotropic region for transcritical flows, low-order moments of small-scale statistics collapse for all cases and Kolmogorov's (1941) theory holds. However, real-fluid thermodynamic effects introduce deviations in the tails of the probability density function for the velocity derivative and, consequently, high-order moments of velocity gradients and dissipation rate in transcritical flows cannot collapse in the locally isotropic region. Analysis of the intermittency shows that the low-order structure functions in transcritical flows follow extended self-similarity, while the dependence of the intermittency factor on real-fluid effects is observed for high-order structure functions. The real-fluid effects on intermittency are explained by turbulent structures related to rare events. Additionally, the dissipation rate moments for transcritical flows follow a universal scaling with Reynolds number, and the scaling exponents are different from those of incompressible flows. These results extend the small-scale universality in incompressible flows (Schumacher et al., Proc. Natl Acad. Sci. USA, vol. 111, 2014, pp. 10961–10965) to realistic flows with significant changes in thermodynamic properties, and provide a physical underpinning of the scaling laws of small-scale statistics at transcritical conditions.

Direct numerical simulations are performed to investigate the characteristics of a turbulent core–annular flow with water-lubricated high viscosity oil in a horizontal pipe. Six different superficial velocity ratios ($\,j_w/j_o = 0.057\unicode{x2013}0.41$) are examined by changing the water superficial velocity $j_w$ for a fixed oil superficial velocity $j_o$. The pressure drops in the pipe and the shapes of the phase interface agree well with those from previous experiments. The oil core flow is almost a plug flow, and the gaps between the phase interface and pipe wall are narrow and wide near the upper and lower surfaces of the pipe, respectively, due to the buoyancy. Within a narrow gap, water is confined mostly in a valley region of the wavy phase interface and hardly goes through its crest. On the other hand, water near the phase interface at a wide gap convects downstream almost at the core speed and the flow near the wall is similar to that of single-phase wall-bounded turbulent flow. The annular flow is characterized by three different regimes depending on the clearance Reynolds number ($Re_c$) based on the core velocity and local gap size: laminar Couette flow driven by the core for $Re_c \le 600$, transitional flow for $600 < Re_c < 2500$ and turbulent flow for $Re_c \ge 2500$. The minimum pressure drop occurs at $j_w/j_o = 0.11$ in the early transition regime. For all $j_w/j_o$ considered, the negative lift force acting on the core comes from the pressure force which balances the buoyancy.

A previously developed method for large-eddy simulations (LES), based on spectral eddy-viscosity models, is generalised to the physical space representation. The method estimates the subgrid-scale (SGS) energy transfer using a similarity-type model expression for the SGS tensor obtained using Gaussian filtering of velocity fields advanced in the simulations. Following steps for the spectral space representation, the SGS transfer in the physical space is used to obtain a spatially varying eddy viscosity at each time step in LES. The computed eddy viscosity is employed to model the SGS stress tensor in the familiar Boussinesq form for use in LES. The method is tested in LES of isotropic turbulence at high Reynolds numbers where the inertial range dynamics is expected and for lower Reynolds number decaying turbulence under conditions of the classical Comte-Bellot and Corrsin experiments. In both cases the agreement with reference data is very good and the SGS transfer computed for the proposed eddy-viscosity model is highly correlated with the transfer computed for the similarity-type stress tensor.

This article extends the previous investigation of the spatial evolution of energy-containing motions in wall-bounded turbulent flows (Kannadasan et al., J. Fluid Mech., vol. 955, 2023, R1) by examining their scale-interactions through spectral analysis based on the spanwise scale decomposition of turbulent kinetic energy and the Reynolds stress transport equation. The energy-containing motions located at the inflow of a turbulent channel flow are artificially removed and the interscale transport mechanisms involved in their spatial evolution are studied. This scale interaction analysis reveals the presence of a significant inverse transfer of streamwise Reynolds stress from the near-wall streaks to larger scales in the spatial evolution of energy-containing motions. This transfer is due to the spanwise variation of streamwise velocity fluctuations, represented by $\partial u'/\partial z$, which is the primary mechanism of streak instability. The analysis presented in this study also shows that the inverse cascade of spanwise energy may correspond to the regeneration of streamwise vortices in the process of reactivating the self-sustaining mechanism in the spatial evolution of energy-containing motions.

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.