Surface bubbles in the ocean are critical in moderating several fluxes between the atmosphere and the ocean. In this paper, we experimentally investigate the drainage and lifetime of surface bubbles in solutions containing surfactants and salts, subjected to turbulence in the air surrounding them modelling the wind above the ocean. We carefully construct a set-up allowing us to repeatably measure the mean lifetime of a series of surface bubbles, while varying the solution and the wind speed or humidity of the air. To that end, we show that renewing the surface layer is critical to avoid a change of the physical properties of the interface. We show that the drainage of the bubbles is well modelled by taking into account the outwards viscous flow and convective evaporation. The mean lifetime of surface bubbles in solutions containing no salt is controlled by evaporation and independent on surfactant concentration. When salt is added, the same scaling is valid only at high surfactant concentrations. At low concentrations, the lifetime is always smaller and independent of wind speed, owing to the presence of impurities triggering a thick bursting event. When the mean lifetime is controlled by evaporation, the probability density of the lifetime is very narrow around its mean, while when impurities are present, a broad distribution is observed.

We present three-dimensional velocity gradient statistics from turbulent Rayleigh–Bénard convection experiments in a horizontally extended cell of aspect ratio 25, a paradigm for mesoscale convection with its organisation into large-scale patterns. The Rayleigh number ${\textit{Ra}}$ ranges from $3.7 \times 10^5$ to $4.8 \times 10^6$, the Prandtl number ${\textit{Pr}}$ from 5 to 7.1. Spatio-temporally resolved volumetric data are reconstructed from moderately dense Lagrangian particle tracking measurements. All nine components of the velocity gradient tensor from the experiments show good agreement with those from direct numerical simulations, both conducted at ${\textit{Ra}} = 1 \times 10^6$ and ${\textit{Pr}} = 6.6$. As expected, with increasing ${\textit{Ra}}$, the flow in the bulk approaches isotropic conditions in the horizontal plane. The focus of our analysis is on non-Gaussian velocity gradient statistics. We demonstrate that statistical convergence of derivative moments up to the sixth order is achieved. Specifically, we examine the probability density functions (PDFs) of components of the velocity gradient tensor, vorticity components, kinetic energy dissipation and local enstrophy at different heights in the bottom half of the cell. The probability of high-amplitude derivatives increases from the bulk to the bottom plate. A similar trend is observed with increasing ${\textit{Ra}}$ at fixed height. Both indicate enhanced small-scale intermittency of the velocity field. We also determine derivative skewness and flatness. The PDFs of the derivatives with respect to the horizontal coordinates are found to be more symmetric than the ones with respect to the vertical coordinate. The conditional statistical analysis of the velocity derivatives with respect to up-/down-welling regions and the rest did not display significant difference, most probably due to the moderate Rayleigh numbers. Furthermore, doubly logarithmic plots of the PDFs of normalised energy dissipation and local enstrophy at all heights show that the left tails follow slopes of 3 / 2 and 1 / 2, respectively, in agreement with numerical results. In general, the left tails of the dissipation and local enstrophy distributions show higher probability values with increasing proximity towards the plate, in comparison with those in the bulk.

We determine unsteady time-periodic flow perturbations that are optimal for enhancing the time-averaged rate of heat transfer between hot and cold walls (i.e. the Nusselt number Nu), under the constraint of fixed flow power (Pe$^2$, where Pe is the Péclet number). The unsteady flows are perturbations of previously computed optimal steady flows and are given by eigenmodes of the Hessian matrix of Nu, the matrix of second derivatives with respect to amplitudes of flow mode coefficients. Positive eigenvalues of the Hessian correspond to increases in Nu by unsteady flows, and occur at $Pe\geqslant 10^{3.5}$ and within a band of flow periods $\tau \sim Pe^{-1}$. For $\tau {\textit{Pe}}\leqslant 10^{0.5}$, the optimal flows are chains of vortices that move along the walls or along eddies enclosed by flow branches near the walls. At larger $\tau {\textit{Pe}}$, the vorticity distributions are often more complex and extend farther from the walls. The heat flux is enhanced at locations on the walls near the unsteady vorticity. We construct an iterative time-spectral solver for the unsteady temperature field, and find increases in Nu of up to 7 % at moderate-to-large perturbation amplitudes.

Shock-tube experiments are conducted to investigate the Atwood-number dependence of hydrodynamic instability induced by a strong shock with a Mach number exceeding 3.0. The compressible linear theory performs reliably under varying compressibility conditions. In contrast, the impulsive model significantly loses predictive accuracy at high shock intensities and Atwood numbers ($A_t$), particularly when specific heat ratio differences across the interface are pronounced. To address this limitation, we propose a modified impulsive model that offers favourable predictions over a wide range of compressibility conditions while retaining practical simplicity. In the nonlinear regime, increasing $A_t$ enhances both the shock-proximity and secondary-compression effects, which suppress bubble growth at early and late stages, respectively. Meanwhile, spike growth is promoted by the spike-acceleration and shock-proximity mechanisms. Several models reproduce spike growth across a wide range of $A_t$, whether physical or incidental. In contrast, no models reliably describe bubble evolution under all $A_t$ conditions, primarily due to neglecting compressibility effects that persist into the nonlinear regime. Building on these insights, we develop an empirical model that effectively captures bubble evolution over a wide $A_t$ range. Modal evolution is further shown to be strongly affected by compressibility-induced variations in interface morphology. The effect is particularly pronounced at moderate to high $A_t$, where it suppresses the fundamental mode growth while promoting higher-order harmonic generation.

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

Microswimmers and active colloids often move in confined systems, including those involving interfaces. Such interfaces, especially at the microscale, may deform in response to the stresses of the flow created by the active particle. We develop a theoretical framework to analyse the effect of a nearby membrane on the motion of an active particle whose flow fields are generated by force-free singularities. We demonstrate our results on a particle represented by a combination of a force dipole and a mass dipole, while the membrane resists deformation due to tension and bending rigidities. We find that the deformation either enhances or suppresses the motion of the active particle, depending on its orientation and the relative strengths between the fundamental singularities that describe its flow. Furthermore, the deformation can generate motion in new directions.

The interaction between a coherent vortex ring and an inertial particle is studied through a combination of experimental and numerical methods. The vortex ring is chosen as a model flow ubiquitous in various geophysical and industrial flows. A detailed description of the vortex properties together with the evolution of the particle kinematics during the interaction is addressed thanks to time-resolved particle image velocimetry and three-dimensional shadowgraphy visualisations. Complementary, direct numerical simulations are realised with a one-way coupling model for the particle, allowing for the identification of the elementary forces responsible for the interaction behaviours. The experimental and numerical results unequivocally demonstrate the existence of three distinct interaction regimes in the parameter range of the present study: simple deviation, strong deviation and capture. These regimes are delineated as functions of key controlled dimensionless parameters, namely, the Stokes number and the initial radial position of the particle relative to the vortex ring axis of propagation.

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

The flow in a rapidly rotating cylinder forced by the harmonic oscillations of a small sphere along the rotation axis is explored numerically. For oscillation frequencies less than twice the cylinder rotation frequency, the forced response flows feature conical shear layers emitted from the critical latitudes of the sphere. These latitudes are where the characteristics of the hyperbolic system, arrived at by ignoring nonlinear, viscous and forcing terms in the governing equations, are tangential to the sphere. These conical shear layers vary continuously with the forcing frequency so long as it remains inertial. At certain values of the forcing frequency, linear inviscid inertial modes of the cylinder are resonated. Of all possible inertial modes, only those whose symmetries are compatible with the symmetry of the forced system are resonated. This all occurs even in the linear limit of vanishingly small forcing amplitude. As the forcing amplitude is increased, nonlinearity leads to non-harmonic oscillations and a non-zero mean flow which features a Taylor columnar structure extending from the sphere to the two endwalls in an axially invariant fashion.

Pulsatile fluid flows through straight pipes undergo a sudden transition to turbulence that is extremely difficult to predict. The difficulty stems here from the linear Floquet stability of the laminar flow up to large Reynolds numbers, well above experimental observations of turbulent flow. This makes the instability problem fully nonlinear and thus dependent on the shape and amplitude of the flow perturbation, in addition to the Reynolds and Womersley numbers and the pulsation amplitude. This problem can be tackled by optimising over the space of all admissible perturbations to the laminar flow. In this paper, we present an adjoint optimisation code, based on a GPU implementation of the pseudo-spectral Navier–Stokes solver nspipe, which incorporates an automatic, optimal checkpointing strategy. We leverage this code to show that the flow is susceptible to two distinct instability routes: one in the deceleration phase, where the flow is prone to oblique instabilities, and another during the acceleration phase with similar mechanisms as in steady pipe flow. Instability is energetically more likely in the deceleration phase. Specifically, localised oblique perturbations can optimally exploit nonlinear effects to gain over nine orders of magnitude in energy at a peak Reynolds number of ${\textit{Re}}_{\textit{max}}\approx 4000$. These oblique perturbations saturate into regular flow patterns that decay in the acceleration phase or break down to turbulence depending on the flow parameters. In the acceleration phase, optimal perturbations are substantially less amplified, but generally trigger turbulence if their amplitude is sufficiently large.