In this experimental work, a conus impacts a deep liquid pool at a speed varying from 1.3 to $19.0\ {\rm cm}\ {\rm s}^{-1}$. Two liquids (2.5 % butanol–water solution or distilled water) and four coni made from duralumin with a diameter of 180 mm and different deadrise angles $\beta$ ($2^{\circ }$, $3^{\circ }$, $4^{\circ }$ and 5$^{\circ }$) are tested. An air cushion is trapped between the conus solid surface and the liquid. Several types of bubble patterns after the collapse of the air cushion are observed: one or multiple bubbles near the conus centre (vertex), irregular trails of bubbles on the conus surface and a ring of bubbles in a ‘necklace’-shaped arrangement. With a total internal reflection set-up and appropriate image post-processing, the external and internal radii of the ring-shaped wetted area are estimated for each frame. The external (internal) radius increases (decreases) in time following a linear (exponential) law. The speed of the outer border of the wetted area is in agreement with the Wagner theory for a body impacting onto a liquid. The initial radius of the annular touchdown region is estimated as the intersection of the relevant fitting curves. In the studied range of parameters, the initial radius obeys a universal scaling law, which follows from the air–water lubrication–inertia balance.

The swinging motion of the eigenmodes of a free inviscid drop has been known for nearly a century. Yet, as the drop sits on a solid substrate, getting flattened by gravity, analytical solutions waver due to the non-spherical base state and the dynamics of the three-phase contact line. The recent paper by Zhang et al. (J. Fluid Mech., vol. 962, 2023, A10) investigated the effect of gravity on the harmonic modes of sessile droplets for free and pinned contact line conditions. An effective boundary element method has been used to solve both axisymmetric and non-axisymmetric modes for a variety of Bond numbers and static contact angles, also revising on the way a debated capillary instability.

Hydrodynamic instabilities in the interstage of turbomachinery and their relation to fan broadband noise are investigated over wide operating conditions. By applying compressible linear stability analysis to the velocity profiles of NASA's high-speed fan rig, source diagnostic test (SDT), three distinct unstable regimes are found: the first mode amplifies disturbances over a wide frequency range including the rotor speed. The second mode rotates much slower and develops inside the outer-wall boundary layer. The third mode spins in the opposite direction. An approach equivalent to dynamic mode decomposition (DMD) successfully extracts these modes from a database of improved delayed detached-eddy simulation (IDDES) solving the SDT geometry at approach and cut-back conditions: the extracted DMD modes generally capture the characteristics of eigenfunctions predicted by the linear stability analysis. At the approach condition, the first two unstable modes seem to interact inside the outer-wall boundary layer, while at the cut-back condition, disturbances associated with the first regime migrate away from the wall. Moreover, outer-wall pressure fluctuations of the IDDES database are Fourier decomposed in the interstage: coherent structures extracted using proper orthogonal decomposition are found to be dominated by duct acoustic modes, and they can be mapped not only to interaction tones but also to broadband noise likely associated with the interaction between disturbances in the rotor wakes and quasi-stationary structures near the outer wall. This interaction noise appears to be radiated to both inlet and exhaust with similar azimuthal-mode distributions; however, the importance relative to other noise sources decreases with the engine speed.

Studies of Kelvin–Helmholtz (KH) instability have typically modelled the initial flow as an isolated shear layer. In geophysical cases, however, the instability often occurs near boundaries and may therefore be influenced by boundary proximity effects. Ensembles of direct numerical simulations are conducted to understand the effect of boundary proximity on the evolution of the instability and the resulting turbulence. Ensemble averages are used to reduce sensitivity to small variations in initial conditions. Both the transition to turbulence and the resulting turbulent mixing are modified when the shear layer is near a boundary: the time scales for the onset of instability and turbulence are longer, and the height of the KH billow is reduced. Subharmonic instability is suppressed by the boundary because phase lock is prevented due to the diverging phase speeds of the KH and subharmonic modes. In addition, the disruptive influence of three-dimensional secondary instabilities on pairing is more profound as the two events coincide more closely. When the shear layer is far from the boundary, the shear-aligned convective instability is dominant; however, secondary central-core instability takes over when the shear layer is close to the boundary, providing an alternate route for the transition to turbulence. Both the efficiency of the resulting mixing and the turbulent diffusivity are dramatically reduced by boundary proximity effects.

We derive a mathematical model of an electromagnetic (EM) field generated by tsunamigenic seabed deformation over an ocean of constant depth. We solve the governing Maxwell equations for the EM field, coupled with a potential flow model of Cauchy–Poisson type for the transient fluid motion forced by seabed deformation. Our new model advances previous studies, where simplified formulae without direct forcing were assumed for the wave field. Using complex integration and large-time asymptotics, we obtain a novel analytical solution for the magnetic field propagating at large distance from the seabed deformation in two dimensions. We show that this magnetic field is made of two terms, one proportional to an Airy function, and thus propagating similarly to the surface gravity wave, and one proportional to a Scorer function, which exhibits a phase lag with respect to the surface gravity wave. Such a phase lag explains the time difference between the arrival of the EM field and the surface gravity wave generated by seabed deformation, which were observed in recent measurements and numerical results. Finally, we discuss the opportunity to detect EM fields as precursors of surface gravity waves in tsunami early warning systems. We introduce a novel non-dimensional parameter to identify the propagation regime of the magnetic field, i.e. self-induction versus diffusion-dominated. We show that tsunami early warning via EM field is possible for diffusion-dominated regimes when the water depth is less than 2 km. Our findings provide a rigorous analytical explanation of existing observations and numerical results.

We present numerical analyses of two-dimensional electrohydrodynamic (EHD) flows of a dielectric liquid between a wire electrode and two plate electrodes with a Poiseuille flow, using direct numerical simulation and global stability analysis. Both conduction and injection mechanisms for charge generation are considered. In this work we focused on the intensity of the cross-flow and studied the EHD flows without a cross-flow, with a weak cross-flow and with a strong cross-flow. (1) In the case without a cross-flow, we investigated its nonlinear flow structures and linear dynamics. We found that the flow in the conduction regime is steady, whereas the flow in the injection regime is oscillatory, which can be explained by a global stability analysis. (2) The EHD flow with a weak cross-flow is closely related to the flow phenomena in an electrostatic precipitator (ESP). Our analyses indicate that increasing the cross-flow intensity or the electric Reynolds number leads to a less stable flow. Based on these results, we infer that one should adopt a relatively low voltage and weak cross-flow in the wire-plate EHD flow to avoid flow instability, which may hold practical implications for ESP. (3) The case of strong cross-flow is examined to study the EHD effect on the wake flow. By comparing the conventional cylindrical wake with the EHD wake in linear and nonlinear regimes, we found that the EHD effect brings forward the vortex shedding in wake flows. Besides, the EHD effect reduces the drag coefficient when the cross-flow is weak, but increases it when it is strong.

This paper discusses how the submergence ratio, defined as the ratio between the flow depth, D, and the height, h, of the solid rigid cylinders forming the array affects flow and turbulence structure inside and around a rectangular array of cylinders placed adjacent to one of the channel sidewalls. As the array becomes submerged, a vertical shear layer develops in between the top face of the array and the free surface, which strongly increases flow three-dimensionality and modifies how the momentum exchange between the array and the surrounding open water regions occurs with respect to the case of an emerged array where only a horizontal shear layer forms as part of the incoming flow approaching the array is deflected laterally. Eddy-resolving simulations are conducted for several values of the solid volume fraction, ϕ, and of the submergence ratio, 1.0 ≤ D/h ≤ 4.0. Similar to the limiting case of an emerged array (D/h = 1.0), the width- and depth-averaged streamwise velocity inside the array reaches a constant value after an initial adjustment region in the submerged-array cases. For D/h ≥ 1.33, the mean normal velocities through the top and side faces of the array do not become equal to zero downstream of the initial adjustment region. The flow inside the array reaches an equilibrium regime where the local flux of fluid leaving the array through its side face is balanced by the local flux of fluid entering the array through its top face. This regime is observed until close to the end of the array. For D/h ≥ 2.0, the horizontal shear layer vortices do not generate successive regions of high and low streamwise velocity and bed friction velocity inside the array, as is observed in the emerged cases. With respect to the emerged case, the size of the shear layer vortices and the shear layer width increase for low submergence ratios before decreasing rapidly for D/h ≥ 1.33. Significant three-dimensional effects are present inside the array and the horizontal shear layer for cases with both emerged and submerged arrays. In the ϕ = 0.08 cases, strong upwelling and downwelling motions are observed inside the array for D/h ≥ 1.33, while the circulation of the streamwise-oriented cell of secondary flow forming close to the lateral face of the array peaks when the submergence ratio is close to 1.33. For constant ϕ, the total streamwise drag force normalized with the height and width of the array increases with increasing submergence ratio. As the array submergence increases, the cylinders near the front of the array contribute less to the total force acting on the cylinders forming the array. For constant D/h, the total streamwise drag force acting on the array increases with ϕ for the submerged and emerged cases.

We experimentally evaluate whether active wall-normal surface deformations are suitable for the targeted control of very-large-scale motions (VLSMs) in a turbulent boundary layer at a friction Reynolds number of $Re_\tau =2600$. Circular surface deformations with a diameter $D$ roughly equal to the boundary layer thickness $\delta$ are generated periodically at a constant amplitude of $0.03\delta$ and at actuation frequencies of $St=0.05$ to 0.20, where $St$ is the Strouhal number based on $D$ and the free stream velocity $U_\infty$. The resulting impact on the flow was captured using high-speed particle image velocimetry and analysed using a triple decomposition. We find that the active surface deformations produce high- and low-speed streamwise velocity fluctuations that are concentrated along the centreline of the actuator. These motions have a negligible impact on the mean velocity profile downstream, i.e. they are truly high and low speed with respect to the unactuated base flow. The motions produced at $St\lesssim 0.1$ are comparable to synthetic VLSMs in terms of their lengths and widths but with a reduced wall-normal extent and rapidly decaying strength. These synthetic motions produce a strong modulation of the turbulence similar to that of the naturally occurring VLSMs. Most notably, we observe that synthetic high-speed motions with magnitudes of the order of $0.05U_\infty$ cause up to a 30 % reduction in turbulence production within the logarithmic layer. The strength and turbulence-modulating characteristics of the synthetic motions appear well suited for targeting the naturally occurring VLSMs locally using a control scheme.

A novel data-driven modal analysis method, reduced-order variational mode decomposition (RVMD), is proposed, inspired by the Hilbert–Huang transform and variational mode decomposition (VMD), to resolve transient or statistically non-stationary flow dynamics. First, the form of RVMD modes (referred to as an ‘elementary low-order dynamic process’, ELD) is constructed by combining low-order representation and the idea of intrinsic mode function, which enables the computed modes to characterize the non-stationary properties of space–time fluid flows. Then, the RVMD algorithm is designed based on VMD to achieve a low-redundant adaptive extraction of ELDs in flow data, with the modes computed by solving an elaborate optimization problem. Further, a combination of RVMD and Hilbert spectral analysis leads to a modal-based time-frequency analysis framework in the Hilbert view, providing a potentially powerful tool to discover, quantify and analyse the transient and non-stationary dynamics in complex flow problems. To provide a comprehensive evaluation, the computational cost and parameter dependence of RVMD are discussed, as well as the relations between RVMD and some classic modal decomposition methods. Finally, the virtues and utility of RVMD and the modal-based time-frequency analysis framework are well demonstrated via two canonical problems: the transient cylinder wake and the planar supersonic screeching jet.

We present the results of direct numerical simulations of a NACA 0012 airfoil, with Mach number 0.3 and angle of attack of $3^\circ$, examining the dynamics of the flow with increasing Reynolds numbers. Two-dimensional simulation results are obtained with chord-based Reynolds numbers in the range $3.2 \times 10^3 \leq Re \leq 2.70 \times 10^4$, where each simulation uses the last time step of the previous one as a starting point, to capture the evolution of dynamics as a function of $Re$. The development of the pressure fluctuations with time shows a transition from periodic to quasi-periodic attractor for $2.38 \times 10^4 \leq Re \leq 2.42 \times 10^4$, leading to the emergence of secondary tones in the wall and acoustic field pressure spectra, different from peaks related to the fundamental frequency $f_1$ and the respective harmonics; a second, incommensurate frequency $f_2$ appears, leading to several secondary tones with frequency $af_1 + bf_2$, with $a$ and $b$ integers. Further increase of the Reynolds number leads to the emergence of a tertiary frequency, $f_3$, indicating a route to chaos of the Ruelle–Takens–Newhouse type. Such a mechanism is related to the ladder-type characteristic structure of the tones, indicating that dynamic systems theory is an important tool for understanding airfoil tonal noise.

Thermoacoustic systems are complex systems where the interactions between the hydrodynamic, acoustic and heat release rate fluctuations lead to diverse dynamics such as chaos, intermittency and ordered dynamics. Such complex interactions cause catastrophically high-amplitude acoustic pressure oscillations and the emergence of order in the spatiotemporal dynamics, referred to as thermoacoustic instability. In this work, we use multilayer networks to study the spatial pattern of inter-subsystem interactions between the vorticity dynamics and thermoacoustic power generated due to acoustically coupled combustion in a bluff-body-stabilised turbulent dump combustor. We construct a two-layered network where the layers represent the thermoacoustic power and vorticity fields. The inter-layer links are determined using cross-variable short-window correlations between vorticity and thermoacoustic power fluctuations at any two locations in the flow field. Analysing the topology of inter-layer networks, using network properties such as degree correlations and link-rank distributions, helps us infer the spatial inhomogeneities in inter-subsystem interactions and unravel the fluid mechanical processes involved during different dynamical states. We show that, during chaotic dynamics, interactions between subsystems are non-localised and spread throughout the flow field of the combustor. During the state of thermoacoustic instability (order), we find that intense interactions occur in between regions of coherent vortex shedding and thermoacoustic power generation and we understand that these processes are strongly and locally coupled. Moreover, we discover that such dense inter-layer connections emerge in spatial pockets in the dump plane of the combustor during the state of intermittency much prior to the onset of order.

Utilizing a fan-stirred chamber and two-dimensional particle image velocimetry, we analyse the modification of homogeneous and isotropic turbulence ($50 \leq Re_\lambda \leq 140$, with supplementary data out to $Re_\lambda = 310$, where $Re_\lambda$ is the longitudinal Taylor Reynolds number) induced by both a non-volatile (water) and a volatile (ethanol) isolated and anchored droplet in the range $(0.3 \leq d/\eta \leq 5.1)$, where $d/\eta$ is the ratio of droplet diameter to the Kolmogorov length scale. The dissipation rate, $\varepsilon$, is calculated via the corrected spatial gradient method, and the resultant fields of both turbulent kinetic energy, $k$, and $\varepsilon$ are presented as spatial heat maps and as shell averages, ${\overline {k_{\Delta r}}}$ and ${\overline {\varepsilon _{\Delta r}}}$, vs the radial coordinate normalized by the droplet radius, $r/R$. The dissipation rate near the water droplet surface may exceed the corresponding unladen flow value by a factor of twenty or more. The normalized radius of recovery, $r^*$, which designates the radial location where ${\overline {k_{\Delta r}}}$ or ${\overline {\varepsilon _{\Delta r}}}$ has returned to within 10 % of the unladen value, is reasonably expressed as $r^* \propto (d/\lambda )^{-C_2}$ in either case, where $\lambda$ is the longitudinal Taylor microscale and $C_2$ is a positive empirical fitting parameter. Recovery of ${\overline {k_{\Delta r}}}$ and ${\overline {\varepsilon _{\Delta r}}}$ may take up to 14 normalized radii when $d/\lambda$ is small. Trend line extrapolation suggests that the attenuation region becomes negligible as $d/\lambda \to 1$. Ethanol, which evaporates up to five times faster than water, induces a much smaller dissipation spike near the surface. The mass ejection phenomenon appears to reduce the strong near-surface damping of the radial root-mean-square component. However, the radius of recovery trend for fields surrounding a volatile ethanol droplet falls directly in line with the non-volatile water droplet data for both $k$ and $\varepsilon$, indicating that droplet vaporization has little effect on the far-field return to isotropy.

A flat-plate turbulent boundary layer (TBL) is experimentally subjected to a family of 22 favourable–adverse pressure gradients (FAPGs) using a ceiling panel of variable convex curvature. We define FAPGs as a sequence of streamwise pressure gradients in the order of favourable followed by adverse, similar in sequence to the pressure gradients over the suction side of an airfoil. For the strongest pressure gradient case, the acceleration parameter, $K$, varied spatially from $6 \times 10^{-6}$ to $- 4.8 \times 10^{-6}$. The adverse pressure gradient (APG) region of this configuration is studied using particle image velocimetry in the streamwise–wall-normal plane. The statistics of the APG TBL show that the upstream favourable pressure gradient (FPG) exerts a strong and lasting downstream influence, and that the rapid spatial changes in the pressure gradients imposed cause an internal boundary layer to grow within the TBL for 15 of the cases studied. The internal layer typically occupies 20 $\%$ of the boundary layer thickness and dominates the boundary layer response, containing the peak turbulent production, peak strength and population of vortices, and most of the spectral energy content of the flow. The outer layer, on the other hand, develops in the APG region without considerable changes to the state dictated by the upstream FPG. These trends are in contrast to APG TBLs that originate from a zero pressure gradient region, where the outer/wake region is known to dominate TBL response. The observed changes are quantified across the family of FAPGs imposed.

We show that the energy required by a turbulent flow to displace a given amount of fluid through a straight duct in a given time interval can be reduced by modulating in time the pumping power. The control strategy is hybrid: it is passive, as it requires neither a control system nor control energy, but it manipulates how pumping energy is delivered to the system (as in active techniques) to increase the pumping efficiency. Our control employs a temporally periodic pumping pattern, where a short and intense acceleration (in which the pumping system is on) followed by a longer deceleration (in which the pumping system is off) makes the flow alternately visit a quasi-laminar and a turbulent state. The computational study is for a plane channel flow, and employs direct numerical simulations, which present specific computational challenges, for example the highly varying instantaneous value of the Reynolds number, and the importance of discretisation effects. Particular care is devoted to a meaningful definition of drag reduction in the present context. The ability of the forcing to yield significant savings is demonstrated. Since only a small portion of the parameter space is investigated, the best performance of the control technique remains to be assessed.

We numerically investigated the global linear instability and bifurcations in electro-thermo-convection (ETC) of a dielectric liquid confined in a two-dimensional (2-D) concentric annulus subjected to a strong unipolar injection. Seven kinds of solutions exist in this ETC system due to the complex bifurcations, i.e. saddle-node, subcritical and supercritical Hopf bifurcations. These bifurcation routes constitute at most four solution branches. Global linear instability analysis and energy analysis were conducted to explain the instability mechanism and transition of different solutions and to predict the local instability regions. The linearized lattice Boltzmann method (LLBM) for global linear instability analysis, first proposed by Pérez et al. (Theor. Comput. Fluid Dyn., vol. 31, 2017, pp. 643–664) to analyse incompressible flows, was extended here to solve the whole set of coupled linear equations, including the linear Navier–Stokes equations, the linear energy equation, Poisson's equation and the linear charge conservation equation. A multiscale analysis was also performed to recover the macroscopic linearized Navier–Stokes equations from the four different discrete lattice Boltzmann equations (LBEs). The LLBM was validated by calculating the linear critical value of 2-D natural convection; it has an error of 1.39% compared with the spectral method. Instability with global travelling wave behaviour is a unique behaviour in the annulus configuration electrothermohydrodynamic system, which may be caused by the baroclinity. Finally, the chaotic behaviour was quantitatively analysed through calculation of the fractal dimension and Lyapunov exponent.

When an immiscible oil drop is immersed in a stably stratified ethanol–water mixture, the Marangoni flow on the surface of the drop can experience an oscillatory instability, so that the drop undergoes a transition from levitating to bouncing. The onset of the instability and its mechanisms have been studied previously (Li et al., Phys. Rev. Lett., vol. 126, 2021, 124502; Li et al., J. Fluid Mech., vol. 932, 2022, A11), yet the bouncing motion of the drop itself, which is a completely different problem, has not yet been investigated. Here we study how the bouncing characteristics (jumping height, rising and sinking time) depend on the control parameters (drop radius, stratification strength, drop viscosity). We first record experimentally the bouncing trajectories of drops of different viscosities in different stratifications. Then a simplified dynamical analysis is performed to get the scaling relations of the jumping height and the rising and sinking times. The rising and sinking time scales are found to depend on the drag coefficient $C_D^S$ of the drop in the stratified liquid, which is determined empirically for the current parameter space (Zhang et al., J. Fluid Mech., vol. 875, 2019, 622–656). For low-viscosity (5 cSt) oil drops, the results on the drag coefficient match those from the literature (Yick et al., J. Fluid Mech., vol. 632, 2009, pp. 49–68; Candelier et al., J. Fluid Mech., vol. 749, 2014, pp. 184–200). For high-viscosity (100 cSt) oil drops, the parameter space had not been explored and the drag coefficients are not readily available. Numerical simulations are therefore performed to provide external verification for the drag coefficients, which well match with the experimental results.

We examine the structure of the continuous circular hydraulic jump and recirculation for a jet impinging on a disk. We use a composite mean-field thin-film approach consisting of subdividing the flow domain into three regions of increasing gravity strength: a developing boundary layer near impact, an intermediate supercritical viscous layer and a region comprising the jump and subcritical flow. Unlike existing models, the approach does not require any empirically or numerically adjusted boundary conditions. We demonstrate that the stress or corner singularity for a film draining at the edge is equivalent to an infinite slope of the film surface, which we impose as the downstream boundary condition. The model is validated against existing experiment and numerical simulation of the boundary-layer and Navier–Stokes equations. We find that the flow in the supercritical region remains insensitive to the change in gravity level but is greatly affected by viscosity. The existence of the jump is not necessarily commensurate with the presence of recirculation, which is strongly dependent on the upstream curvature and steepness of the jump.

The dynamics of suspensions of particles has been an active area of research since Einstein first calculated the leading-order correction to the viscosity of a suspension of spherical particles (Einstein, Proc. R. Soc., vol. A102, 1906, pp. 161–179). Since then, researchers have strived to develop an accurate description of the behaviours of suspensions that goes beyond just leading order in the particle volume fraction. Here, we consider the low-Reynolds-number behaviour of a suspension of spherical particles. Working from the Green's functions for the flow due to a single particle, we derive a continuum-level description of the dynamics of suspensions. Our analysis corrects an error in the derivation of these equations in the work of Jackson (Chem. Engng Sci., vol. 52, 1997, pp. 2457–2469) and leads to stable equations of motion for the particles and fluid. In addition, our resulting equations naturally give the sedimentation speed for suspended particles and correct a separate error in the calculation by Batchelor (J. Fluid Mech., vol. 52, 1972, pp. 245–268). Using the pair-correlation function for hard spheres, we are able to compute the sedimentation speed out to seventh order in the volume fraction, which agrees with experimental data up to 30 %–35 %, and also get higher-order corrections to the suspension viscosity, which agree with experiments up to $\sim$15 %. Then, using the pair distribution for spheres in shear flow, we find alterations to both the first and second normal stresses.

Direct numerical simulations of the turbulent flow in a channel are conducted up to $Re_{\tau }=1000$ to examine the influence of the friction Reynolds number on the translational and angular velocities of inertial, prolate ellipsoids. The quadrant distribution of the turbulent events seen by the particles is not significantly affected by the value of $Re_{\tau }$, but subtle modifications take place, depending on the position in the channel and on the particle relaxation time. Overall, the influence of $Re_{\tau }$ on the first and second statistical moments of the ellipsoids translational velocity is the same as that observed for the fluid velocity. The weak dependence of these statistics to the particle shape previously observed at low Reynolds number remains at higher values of $Re_{\tau }$. Similarly, the mean and root mean square (r.m.s.) of the angular velocity of the fluid seen by the particles weakly depend on particle shape and they have the same dependence to $Re_{\tau }$ as the angular velocity statistics of the carrier fluid. Particle angular velocity statistics are more strongly affected by the flow Reynolds number due to the evolution of the complex shape and inertia dependent rotation orbits with $Re_{\tau }$. In the near-wall region the average angular velocity of weakly inertial ellipsoids increases with $Re_{\tau }$ due to their stronger alignment with the mean fluid vorticity. Furthermore, the r.m.s. of the wall-normal component of the angular velocity of more inertial ellipsoids increases with $Re_{\tau }$ owing to the larger fluctuations of the angle between the particle major axis and the velocity-gradient plane.

The water surface response to subcritical turbulent flow over a backward-facing step (BFS) is studied via high-resolution large-eddy simulation (LES). The LES method is validated first using data of previously reported experiments. The LES-predicted water surface is decomposed into different types of gravity waves as well as turbulence-driven forced waves. Analysis of the LES data reveals the interplay between low-frequency large-scale turbulence structures, which are the result of flow separation from the step and reattachment behind the step, and the dynamics of the water surface. The water surface deformation is mainly the result of freely propagating gravity waves and forced waves, owing to turbulence in the form of rollers and/or hairpin vortices. Gravity waves with zero group velocity define the characteristic spatial and temporal scales of the surface deformations at higher frequencies, while large eddies determine their low-frequency modulation. These deformations are mainly confined in lateral bands that propagate downstream following the advection of the near-surface streamwise vortices (rollers) that are shed from the step. Steeper surface waves are observed in regions of negative perturbation velocity gradient and down-welling, downstream of the larger rollers, and are associated with thin isolated regions of high vorticity near the surface. The investigation of such a complex flow has shown that the decomposition of the water surface fluctuations into its different physical components may be used to identify the dynamics of the underlying flow structure.