We investigated experimentally the settling behaviour of vertically aligned spherical particles within various quiescent media at different release frequencies. The particles had a diameter of $d = 4$ mm and density of $\rho _s = 2200$ kg m$^{-3}$, and were released near the free surface of water, ethanol, a G60 water–glycerine mixture (60 % glycerine by weight) and oil media at frequencies of $f_P = 4$, 6 and 8 Hz, thereby allowing study of Galileo numbers, $Ga \in [16, 976]$. Particle tracking velocimetry quantified the motion of nearly 800 particles in a 600 mm high tank, and particle image velocimetry examined flow patterns around the particles. Results revealed that the centre of mass of the particle trajectories exhibited preferential in-plane motions, with significant lateral dispersion and large $Ga$ in water and ethanol, and nearly vertical paths with low $Ga$ in the G60 mixture and oil media. Varying degrees of particle separation resulted in higher terminal velocities than for a single particle. Hence, particle drag decreased in all cases, with the oil medium showing the highest drag reduction under the closest particle separation, reaching up to nearly 70 % of that for the single particle. The vertical and lateral pair dispersions, $R^2_z$ and $R^2_L$, exhibited ballistic scaling, with dependences on the initial separation, $r_0$, and the type of medium. With large $Ga$, $R^2_z$ displayed a ballistic regime followed by a slower rate, whereas with small $Ga$, $R^2_z$ maintained a consistent ballistic regime throughout settling. Finally, normalized $R^2_z$ demonstrated distinct scaling (exponent 2/3 and 1) dependent on the normalized initial separation and $Ga$.

Metachronal rowing is a biological propulsion mechanism employed by many swimming invertebrates (e.g. copepods, ctenophores, krill and shrimp). Animals that swim using this mechanism feature rows of appendages that oscillate in a coordinated wave. In this study, we used observations of a swimming ctenophore (comb jelly) to examine the hydrodynamic performance and vortex dynamics associated with metachronal rowing. We first reconstructed the beating kinematics of ctenophore appendages based on a high-speed video of a metachronally coordinated row. Following the reconstruction, two numerical models were developed and simulated using an in-house immersed-boundary-method-based computational fluid dynamics solver. The two models included the original geometry (16 appendages in a row) and a sparse geometry (8 appendages, formed by removing every other appendage along the row). We found that appendage tip vortex interactions contribute to hydrodynamic performance via a vortex-weakening mechanism. Through this mechanism, appendage tip vortices are significantly weakened during the drag-producing recovery stroke. As a result, the swimming ctenophore produces less overall drag, and its thrust-to-power ratio is significantly improved (up to 55.0 % compared with the sparse model). Our parametric study indicated that such a propulsion enhancement mechanism is less effective at higher Reynolds numbers. Simulations were also used to investigate the effects of substrate curvature on the unsteady hydrodynamics. Our results illustrated that, compared with a flat substrate, arranging appendages on a curved substrate can boost the overall thrust generation by up to 29.5 %. These findings provide new insights into the fluid dynamic principles of propulsion enhancement underlying metachronal rowing.

The gravity flow of a granular material between two vertical walls separated by a width $2W$ is simulated using the discrete element method (DEM). Periodic boundary conditions are applied in the flow (vertical) and the other horizontal directions. The mass flow rate is controlled by specifying the average solids fraction $\bar {\phi }$, the ratio of the volume of the particles to the volume of the channel. A steady fully developed state can be achieved for a narrow range of $\bar {\phi }$, $\bar {\phi }_{max} \geq \bar {\phi } \geq \bar {\phi }_{cr}$, and the material is in free fall for $\bar {\phi } < \bar {\phi }_{min}$. For an intermediate range of $\bar {\phi }$ ($\bar {\phi }_{cr} > \bar {\phi } \geq \bar {\phi }_{min}$), there are oscillations in the horizontal coordinate of the centre of mass, velocity components and stress. As $\bar {\phi }$ decreases in the range $\bar {\phi }_{cr} > \bar {\phi } \geq \bar {\phi }_{min}$, the amplitude of the oscillations increases proportional to $\sqrt {\bar {\phi }_{cr} - \bar {\phi }}$ and the frequency appears to tend to a non-zero value as $\bar {\phi } \rightarrow \bar {\phi }_{cr}$, indicating a supercritical Hopf bifurcation. The relation between the dominant frequency and the higher harmonics of the position, velocity and stress fluctuations are explained using the momentum balance. It is found that dissipation in the inter-particle and particle–wall contacts is critical for the presence of an oscillatory state.

The recent Lagrangian relaxation towards equilibrium (LaRTE) approach (Fowler et al., J. Fluid Mech., vol. 934, 2022, A44) is a wall model for large-eddy simulations (LES) that isolates quasi-equilibrium wall-stress dynamics from non-equilibrium responses to time-varying LES inputs. Non-equilibrium physics can then be modelled separately, such as the laminar Stokes layers that form in the viscous region and generate rapid wall-stress responses to fast changes in the pressure gradient. To capture additional wall-stress contributions due to near-wall turbulent eddies, a model term motivated by the attached eddy hypothesis is proposed. The total modelled wall stress thus includes contributions from various processes operating at different time scales (i.e. the LaRTE quasi-equilibrium plus laminar and turbulent non-equilibrium wall stresses) and is called the multi-time-scale (MTS) wall model. It is applied in LES of turbulent channel flow subject to a wide range of unsteady conditions from quasi-equilibrium to non-equilibrium. Flows considered include pulsating and linearly accelerating channel flow for several forcing frequencies and acceleration rates, respectively. We also revisit the sudden spanwise pressure gradient flow (considered in Fowler et al., J. Fluid Mech., vol. 934, 2022, A44) to review how the newly introduced model features affect this flow. Results obtained with the MTS wall model show good agreement with direct numerical simulation data over a vast range of conditions in these various non-equilibrium channel flows. To further understand the MTS model, we also describe and test the instantaneous-equilibrium limit of the MTS wall model. In this limit, good wall-stress predictions are obtained with reduced model complexity but providing less complete information about the wall-stress physics.

Algorithms for reconstructing and predicting nonlinear ocean wave fields from remote measurements are presented. Three types of synthetic observations are used to quantify the influence of remote measurement modulation mechanisms on the algorithms’ performance. First, the observations correspond to randomly distributed surface elevations. Then, they are related to a marine radar model – the second type takes the wave shadowing modulation into account whereas the third one also includes the tilt modulation. The observations are numerically generated based on unidirectional waves of various steepness values. Linear and weakly nonlinear prediction algorithms based on analytical models are considered, as well as a highly nonlinear algorithm relying on the high-order spectral (HOS) method. Reconstructing surfaces from shadowed observations is found to have an impact limited to the non-visible regions, while tilt modulation affects the reconstruction more generally due to the indirect, more complex extraction of wave information. It is shown that the accuracy of the surface reconstruction mainly depends on the correct modelling of the wave shape nonlinearities. Modelling the nonlinear correction of the dispersion relation, in particular the frequency-dependent wave phase effects in the case of irregular waves, substantially improves the prediction. The suitability of the algorithms for severe wave conditions in finite depth and using non-perfect observations is assessed through wave tank experiments. It shows that only the third-order HOS solution predicts the right amplitude and phase of an emerging extreme wave, emphasizing the relevance of the corresponding physical modelling.

We study the spectral energy transfer due to wave–triad interactions in the Garrett–Munk spectrum of internal gravity waves based on a numerical evaluation of the collision integral in the wave kinetic equation. Our numerical evaluation builds on the reduction of the collision integral on the resonant manifold for a horizontally isotropic spectrum. We evaluate directly the downscale energy flux available for ocean mixing, whose value is in close agreement with the finescale parameterization. We further decompose the energy transfer into contributions from different mechanisms, including local interactions and three types of non-local interactions, namely parametric subharmonic instability, elastic scattering (ES) and induced diffusion (ID). Through analysis on the role of each mechanism, we resolve two long-standing paradoxes regarding the mechanism for forward cascade in frequency and zero ID flux for the GM76 spectrum. In addition, our analysis estimates the contribution of each mechanism to the energy transfer in each spectral direction, and reveals new understanding of the importance of local interactions and ES in the energy transfer.

We propose a constitutive equation for flows of gases in high-rate regimes where the Navier–Stokes theory breaks down. The model generalizes the Navier–Stokes relation and agrees well with that model in all lower rate flows examined. Our proposed constitutive relation is calibrated with the method of objective molecular dynamics (OMD) using families of compressible and incompressible flows of Lennard-Jones argon. The constitutive relation makes use of the higher-order objective strain rates due to Rivlin and Ericksen (J. Rat. Mech. Anal., vol. 4, 1955, pp. 323–425). The constitutive relation is fully frame-indifferent, and the macroscopic flows corresponding to the OMD simulations are exact solutions for the proposed model. The model is shown to agree with atomistic results much better than the Navier–Stokes equations in the transition regime. The success of our model indicates that it is not higher gradients that become important in the high-rate regime, but rather higher rates of change of the strain rate tensor. While somewhat more complicated to implement than the Navier–Stokes relation, the proposed model is expected to be compatible with existing methods of computational fluid dynamics and may extend those methods to higher rate regimes, while preserving their ability to handle large spatial scales.

In this paper we report on an experimental study focusing on the manifestation and dynamics of the large-scale circulation (LSC) in turbulent liquid metal convection. The experiments are performed inside a cylinder of aspect ratio $\varGamma = 0.5$ filled with the ternary alloy GaInSn, which has a Prandtl number of $Pr = 0.03$. The large-scale flow structures are classified and characterized at Rayleigh numbers of ${Ra} = 9.33 \times 10^6, 5.31 \times 10^7$ and $6.02 \times 10^8$ by means of the contactless inductive flow tomography which enables the full reconstruction of the three-dimensional (3-D) flow structures in the entire convection cell. This is complemented with the multi-thermal-probe method for capturing the azimuthal temperature variation induced by the LSC at the sidewall. We use proper orthogonal decomposition (POD) to identify the dominating modes of the turbulent convection. The analysis reveals that a single-roll structure of the LSC alternates in short succession with double-roll structures or a three-roll structure. This is accompanied by dramatic fluctuations of the Reynolds number, whose instantaneous values can deviate by more than 50 % from the time-average value. No coherent oscillations are observed, whereas a correlation analysis indicates a residual contribution of the torsion and sloshing modes. Results of the POD analysis suggest a stabilization of the single-roll LSC with increasing $Ra$ at the expense of flow structures with multiple rolls. Moreover, the relative lifetime of all identified flow states, measured in units of free-fall times, increases with rising $Ra$.

We investigate steady and oscillatory flow through a hexagonal close-packed arrangement of spheres in the framework of the volume-averaged momentum equation. We quantify the friction and pressure drag based on a direct numerical simulation dataset. Using the pressure decomposition of Graham (J. Fluid Mech., vol. 881, 2019), the pressure drag can be further split up into an accelerative, a viscous and a convective contribution. For the accelerative pressure, a closed-form expression can be given in terms of the potential flow solution. We investigate the contributions of the different drag components to the volume-averaged momentum budget and their Reynolds number scaling. For steady flow, we find that the friction and viscous pressure drag are proportional to $Re$ at low Reynolds numbers and scale with $Re^{1.4}$ for high Reynolds numbers. This is close to the steady laminar boundary layer scaling. For the convective pressure drag, we find a cubic scaling at low and a quadratic scaling at high Reynolds numbers. The Reynolds stresses have a minor contribution to the momentum budget. For oscillatory flow at low and medium Womersley numbers, the amplitudes of the drag components are similar to the steady cases at the same Reynolds number. At high Womersley numbers, the drag components behave quite differently and the friction and viscous pressure drag are relatively insignificant. The drag components are not in phase with the forcing and the superficial velocity; the phase lag increases with the Womersley number. This suggests that new models beyond the current quasisteady approaches need to be developed.

The refraction of surface gravity waves by currents leads to spatial modulations in the wave field and, in particular, in the significant wave height. We examine this phenomenon in the case of waves scattered by a localised current feature, assuming (i) the smallness of the ratio between current velocity and wave group speed, and (ii) a swell-like, highly directional wave spectrum. We apply matched asymptotics to the equation governing the conservation of wave action in the four-dimensional position–wavenumber space. The resulting explicit formulas show that the modulations in wave action and significant wave height past the localised current are controlled by the vorticity of the current integrated along the primary direction of the swell. We assess the asymptotic predictions against numerical simulations using WAVEWATCH III for a Gaussian vortex. We also consider vortex dipoles to demonstrate the possibility of ‘vortex cloaking’ whereby certain currents have (asymptotically) no impact on the significant wave height. We discuss the role of the ratio of the two small parameters characterising assumptions (i) and (ii) above, and show that caustics are significant only for unrealistically large values of this ratio, corresponding to unrealistically narrow directional spectra.

We carry out direct numerical simulation (DNS) of turbulent flow in rough pipes. Two types of irregular roughness are investigated, namely a grit-blasted and a graphite surface. A wide range of Reynolds numbers is tested, from the laminar up to the fully rough regime, attempting to replicate Nikuradse's pioneering study. Despite the large relative roughness, outer-layer similarity is achieved at high Reynolds number as hypothesised by Townsend, with deviations from the smooth wall case of 4 % for the grit-blasted surface and 13 % for the graphite surface. This makes it possible to define a roughness function and the equivalent sand-grain roughness. The results are compared with those obtained in plane channels, with small differences pointing to the residual influence of the duct cross-sectional shape in the presence of relatively large roughness. The computed friction factors behave similar to those Nikuradse's chart, with differences in terms of the friction factor in the laminar region and of the critical Reynolds number, which are partly absorbed by using the hydraulic radius as reference length scale. The distributions of the velocity fluctuations intensities point to a isotropisation of turbulence in the near-wall region resulting from the roughness, with influence of the roughness geometry. Comparison of the computed equivalent sand-grain roughness height suggest that existing correlations suffer from poor predictive power, at least for surfaces with large relative roughness.

A turbulent circular disk wake database (Chongsiripinyo & Sarkar, J. Fluid Mech., vol. 885, 2020) at Reynolds number $Re = U_\infty D/\nu = 5 \times 10^{4}$ is interrogated to identify the presence of large-scale streaks – coherent elongated regions of streamwise velocity. The unprecedented streamwise length – until $x/D \approx 120$ – of the simulation enables investigation of the near and far wakes. The near wake is dominated by the vortex shedding (VS) mode residing at azimuthal wavenumber $m=1$ and Strouhal number $St = 0.135$. After filtering out the VS structure, conclusive evidence of large-scale streaks with frequency $St \rightarrow 0$, equivalently streamwise wavenumber $k_x \rightarrow 0$ in the wake, becomes apparent in visualizations and spectra. These streaky structures are found throughout the simulation domain beyond $x/D \approx 10$. Conditionally averaged streamwise vorticity fields reveal that the lift-up mechanism is active in the near as well as the far wake, and that ejections contribute more than sweeps to events of intense $-u'_xu'_r$. Spectral proper orthogonal decomposition is employed to extract the energy and the spatiotemporal features of the large-scale streaks. The streak energy is concentrated in the $m=2$ azimuthal mode over the entire domain. Finally, bispectral mode decomposition is conducted to reveal strong interaction between $m=1$ and $St = \pm 0.135$ modes to give the $m=2$, $St \rightarrow 0$ streak mode. Our results indicate that the self-interaction of the VS mode generates the $m=2$, $St \rightarrow 0$ streamwise vortices, which leads to streak formation through the lift-up process. To the authors’ knowledge, this is the first study that reports and characterizes large-scale low-frequency streaks and the associated lift-up mechanism in a turbulent wake.

Self-propulsion of chemically active droplets and phoretic disks has been studied widely; however, most research overlooks the influence of disk shape on swimming dynamics. Inspired by experimentally observed prolate composite droplets and elliptical camphor disks, we employ simulations to investigate the phoretic dynamics of an elliptical disk that emits solutes uniformly in the creeping flow regime. By varying the disk's eccentricity $e$ and the Péclet number $Pe$, we distinguish five disk behaviours: stationary, steady, orbiting, periodic and chaotic. We perform a linear stability analysis (LSA) to predict the onset of instability and the most unstable eigenmode when a stationary disk transitions spontaneously to steady self-propulsion. In addition to the LSA, we use an alternative approach to determine the perturbation growth rate, illustrating the competing roles of advection and diffusion. The steady motion features a transition from a puller-type to a neutral-type swimmer as $Pe$ increases, which occurs as a bimodal concentration profile at the disk surface shifts to a polarized solute distribution, driven by convective solute transport. An elliptical disk achieves an orbiting motion through a chiral symmetry-breaking instability, wherein it repeatedly follows a circular path while simultaneously rotating. The periodic swinging motion, emerging from a steady motion via a supercritical Hopf bifurcation, is characterized by a wave-like trajectory. We uncover a transition from normal diffusion to superdiffusion as eccentricity $e$ increases, corresponding to a random-walking circular disk and a ballistically swimming elliptical counterpart, respectively.

The objective of the present study is to provide a numerical database of thermal boundary layers and to contribute to the understanding of the dynamics of passive scalars at different Prandtl numbers. In this regard, a direct numerical simulation (DNS) of an incompressible zero-pressure-gradient turbulent boundary layer is performed with the Reynolds number based on momentum thickness ${\textit {Re}}_{\theta }$ ranging up to $1080$. Four passive scalars, characterized by the Prandtl numbers ${\textit {Pr}} = 1,2,4,6$ are simulated using the pseudo-spectral code SIMSON (Chevalier et al., SIMSON : a pseudo-spectral solver for incompressible boundary layer flows. Tech. Rep. TRITA-MEK 2007:07. KTH Mechanics, Stockholm, Sweden, 2007). To the best of our knowledge, the present DNS provides the thermal boundary layer with the highest Prandtl number available in the literature. It corresponds to that of water at ${\sim }24\,^{\circ }{\rm C}$, when the fluid temperature is considered as a passive scalar. Turbulence statistics for the flow and thermal fields are computed and compared with available numerical simulations at similar Reynolds numbers. The mean flow and scalar profiles, root-mean-squared velocity and scalar fluctuations, turbulent heat flux, turbulent Prandtl number and higher-order statistics agree well with the numerical data reported in the literature. Furthermore, the pre-multiplied two-dimensional spectra of the velocity and of the passive scalars are computed, providing a quantitative description of the energy distribution at the different length scales for various wall-normal locations. The energy distribution of the heat-flux fields at the wall is concentrated on longer temporal structures with increasing Prandtl number. This is due to the thinner thermal boundary layer as thermal diffusivity decreases and, thereby, the longer temporal structures exhibit a different footprint at the wall.

The Earth's landscape hosts a variety of patterns resulting from the interaction between a sediment-carrying fluid and an erodible boundary. Here, the morphodynamics of river bifurcations is interpreted as a second-order phase transition. A consolidated one-dimensional bifurcation model is re-examined in the light of classical Landau theory of critical phenomena. The transition from a balanced to an unbalanced flux partition is described in terms of an order parameter. The equilibrium states of the system are shown to be minima of a morphodynamic potential function. Finally, the role of a generic external forcing is investigated. A threshold value of the forcing is shown to set bounds separating two different morphodynamic responses to allogenic and autogenic dynamics.

The present paper finds that the coexistence of multiple primary instability waves may cause a non-trivial nonlinear interaction and breakdown process, which has not been reported before. In the considered Mach 6 flat-plate boundary layer, a global resolvent analysis reports three optimal disturbances (local maxima): a high-frequency planar wave, a low-frequency oblique wave and a stationary streak. For the dominant planar and oblique waves, a parabolised stability equation analysis identifies the initial non-modal transient growth and downstream modal growth. Initiated by these two optimal disturbances jointly, the complete linear and nonlinear instability processes until breakdown to turbulence are shown with direct numerical simulation. Owing to the transient growth, the oblique wave may be more significant than the planar wave in the breakdown. The oblique wave and scales of nonlinear interactions are pronounced in the outer layer, whose significance may not be comprehensively characterised by the wall pressure measurement. Fourier modes characterising the oblique-wave oblique breakdown, the planar-wave fundamental resonance, the planar-wave subharmonic resonance and the combination resonance related to a detuned mode are observed successively. The detuned mode seems to dominate the near-wall dynamics in the late nonlinear stage, characterised by $\varLambda$-like structures. Meanwhile, the existence of this detuned mode is independent of the initial amplitude ratio and the absolute amplitude of the oblique and planar waves. Weakly nonlinear stability analyses demonstrate that the detuned mode is mainly a consequence of the secondary instability under the combination of planar and oblique primary waves. Wave vector plots reveal the resonant state of multiple triads. Energy budget and amplitude-correction analyses provide a clear physical image of energy transfer.

This study examined the impact of coherent structures on the aerodynamic forces exerted on a NACA0012 aerofoil with angles of attack $7.5^{\circ }$ and $10^{\circ }$, and a chord-based Reynolds number $50\,000$. The study utilized the spectral proper orthogonal decomposition (SPOD) algorithm to identify the coherent structures, and vorticity force analysis to quantify their impact on lift and drag forces. Results showed that at $10^{\circ }$, the zeroth frequency of the first SPOD mode had a significant impact on drag and lift forces due to a large vortex structure that caused a strong flow along the suction side of the aerofoil. The first and second frequencies of the first SPOD mode represented asymmetric vortex pairs and a series of vortex pairs that determined the leading-edge separation, respectively. At $7.5^{\circ }$, the zeroth frequency of the first mode corresponded to an oscillating near-wall stream that followed the reattachment flow pattern, while the first frequency corresponded to a counter-rotating vortex pair that originated where the flow reattaches. Finally, the second frequency of the first mode corresponded to smaller counter-rotating vortex pairs at the shear layer originated near the reattachment point. These findings suggest that coherent structures have a significant impact on aerodynamic forces exerted on aerofoils, and can be identified and quantified using the SPOD algorithm and vorticity force analysis.

Microvortex generators are passive control devices smaller than the boundary layer thickness that energise the boundary layer to prevent flow separation with limited induced drag. In this work, we use direct numerical simulations (DNS) to investigate the effect of the Reynolds number in a supersonic turbulent boundary layer over a microramp vortex generator. Three friction Reynolds numbers are considered, up to $Re_\tau = 2000$, for fixed free stream Mach number $M_\infty =2$ and fixed relative height of the ramp with respect to the boundary layer thickness. The high-fidelity data set sheds light on the instantaneous and highly three-dimensional organisation of both the wake and the shock waves induced by the microramp. The full access to the flow field provided by DNS allows us to develop a qualitative model of the near wake, explaining the internal convolution of the Kelvin–Helmholtz vortices around the low-momentum region behind the ramp. The overall analysis shows that numerical results agree excellently with recent experimental measurements in similar operating conditions and confirms that microramps effectively induce a significantly fuller boundary layer even far downstream of the ramp. Moreover, results highlight significant Reynolds number effects, which in general do not scale with the ramp height. Increasing Reynolds number leads to enhanced coherence of the typical vortical structures in the field, faster and stronger development of the momentum deficit region, increased upwash between the primary vortices from the sides of the ramp – and thus increased lift-up of the wake – and faster transfer of momentum towards the wall.

When a blunt body impacts an air–water interface, large hydrodynamic forces often arise, a phenomenon many of us have unfortunately experienced in a failed dive or ‘belly flop’. Beyond assessing risk to biological divers, an understanding and methods for remediation of such slamming forces are critical to the design of numerous engineered naval and aerospace structures. Herein we systematically investigate the role of impactor elasticity on the resultant structural loads in perhaps the simplest possible scenario: the water entry of a simple harmonic oscillator. Contrary to conventional intuition, we find that ‘softening’ the impactor does not always reduce the peak impact force, but may also increase the force as compared with a fully rigid counterpart. Through our combined experimental and theoretical investigation, we demonstrate that the transition from force reduction to force amplification is delineated by a critical ‘hydroelastic’ factor that relates the hydrodynamic and elastic time scales of the problem.

Non-stationarity is the rule in the atmospheric boundary layer (ABL). Under such conditions, the flow may experience departures from equilibrium with the underlying surface stress, misalignment of shear stresses and strain rates, and three-dimensionality in turbulence statistics. Existing ABL flow theories are primarily established for statistically stationary flow conditions and cannot predict such behaviours. Motivated by this knowledge gap, this study analyses the impact of time-varying pressure gradients on mean flow and turbulence over urban-like surfaces. A series of large-eddy simulations of pulsatile flow over cuboid arrays is performed, programmatically varying the oscillation amplitude $\alpha$ and forcing frequency $\omega$. The analysis focuses on both longtime-averaged and phase-dependent flow dynamics. Inspection of longtime-averaged velocity profiles reveals that the aerodynamic roughness length $z_0$ increases with $\alpha$ and $\omega$, whereas the displacement height $d$ appears to be insensitive to these parameters. In terms of oscillatory flow statistics, it is found that $\alpha$ primarily controls the oscillation amplitude of the streamwise velocity and Reynolds stresses, but has a negligible impact on their wall-normal structure. On the other hand, $\omega$ determines the size of the region affected by the unsteady forcing, which identifies the so-called Stokes layer thickness $\delta _s$. Within the Stokes layer, phase-averaged resolved Reynolds stress profiles feature substantial variations during the pulsatile cycle, and the turbulence is out of equilibrium with the mean flow. Two phenomenological models have been proposed that capture the influence of flow unsteadiness on $z_0$ and $\delta _s$, respectively.