In turbulence near a free surface, strong vortices attach to the surface, creating surface imprints visible as nearly circular ‘dimples’. By studying these imprints in direct numerical simulation (DNS) data we make two observations. First, the imprints of surface-attached vortices can be very effectively distinguished from other turbulent surface features using two physical features: they are nearly circular in shape, and persist for a long time compared with other pertinent time scales. Secondly, the instantaneous number of surface dimples from surface-attached vortices in an area, $N(t)$, is intimately related to its mean-square surface divergence, $\beta ^2(t)$. We develop a simple and physically transparent computer vision procedure which, using the properties of low eccentricity and longevity, detects and tracks vortices from their surface features only, with sensitivity and accuracy of $90\,\%$ or better. We compare $N(t)$ and $\beta ^2(t)$, finding a normalised cross-correlation of $0.90$, with changes in $N$ lagging around $0.8T_\infty$ behind those in $\beta ^2$ ($T_\infty$ is an integral time scale), confirming the common observation that vortices are spawned by strong upwelling events where $\beta ^2$ is large. These findings suggest that the rate of mass flux across the surface, being closely related to surface divergence, can be estimated remotely in some natural flows using visible free-surface dimples as proxy.

The transition to turbulence in conduits is among the longest-standing problems in fluid mechanics. Challenges in producing or saving energy hinge on understanding promotion or suppression of turbulence. While a global picture based on an intrinsically 3-D subcritical mechanism is emerging for 3-D turbulence, subcritical turbulence is yet to even be observed when flows approach two dimensions, e.g. under intense rotation or magnetic fields. Here, stability analysis and direct numerical simulations demonstrate a subcritical quasi-two-dimensional (quasi-2-D) transition from laminar flow to turbulence, via a radically different 2-D mechanism to the 3-D case, driven by nonlinear Tollmien–Schlichting waves. This alternative scenario calls for a new line of thought on the transition to turbulence and should inspire new strategies to control transition in rotating devices and nuclear fusion reactor blankets.

In this note we provide an analytical proof of the zero helicity condition for systems governed by the Gross–Pitaevskii equation (GPE). The proof is based on the hydrodynamic interpretation of the GPE, and the direct use of Noether's theorem by applying Kleinert's multi-valued gauge theory. As a by-product we also demonstrate the conservation and quantization of the circulation for the GPE.

In homogeneous and isotropic turbulence, the relative contributions of different physical mechanisms to the energy cascade can be quantified by an exact decomposition of the energy flux (Johnson, Phys. Rev. Lett., vol. 124, 2020, 104501; J. Fluid Mech., vol. 922, 2021, A3). We extend the formalism to the transfer of kinetic helicity across scales, important in the presence of large-scale mirror‐breaking mechanisms, to identify physical processes resulting in helicity transfer and quantify their contributions to the mean flux in the inertial range. All subfluxes transfer helicity from large to small scales. Approximately 50 % of the mean flux is due to the scale-local vortex flattening and vortex twisting. We derive a new exact relation between these effects, similar to the Betchov relation for the energy flux, revealing that the mean contribution of the former is three times larger than that of the latter. Multi-scale effects account for the remaining 50 % of the mean flux, with approximate equipartition between multi-scale vortex flattening, twisting and entangling.

When surface waves interact with ambient turbulence, the two affect each other mutually. Turbulent eddies get redirected, intensified and periodically stretched and compressed, while the waves suffer directional scattering. We study these mutual interactions experimentally in the water channel laboratory at the Norwegian University of Science and Technology (NTNU) Trondheim. Long groups of waves were propagated upstream on currents with identical mean flow but different turbulence properties, created by an active grid at the current inlet. The subsurface flow in the spanwise–vertical plane was measured with stereo particle-image velocimetry. Comparing the subsurface velocity fields before and after the passage of a wave group, a strong enhancement of streamwise vorticity is observed which increases rapidly towards the surface for $k_0z\gtrsim -0.3$ ($z$, vertical distance from still surface; $k_0$, carrier wavenumber) in qualitative agreement with theory. Next, we measure the broadening of the directional wave spectrum at increasing propagation distance. The rate of directional diffusion is greatest for the turbulent case with the highest energy at the longest length scales whereas the highest total turbulent kinetic energy overall did not produce the most scattering. The variance of directional spectra is found to increase linearly as a function of propagation time.

The interaction of supercritical turbulent flows with granular sediment beds is challenging to study both experimentally and numerically. This challenging task has hampered advances in understanding antidunes, the most characteristic bedform of supercritical flows. This article presents the first numerical attempt to simulate upstream-migrating antidunes with geometrically resolved particles and a liquid–gas interface. Our simulations provide data at a resolution higher than laboratory experiments, and they can therefore provide new insights into the mechanisms of antidune migration and contribute to a deeper understanding of the underlying physics. To manage the simulations’ computational costs and physical complexity, we employ the cumulant lattice Boltzmann method in conjunction with a discrete element method for particle interactions, as well as a volume-of-fluid scheme to track the deformable free surface of the fluid. By reproducing two flow configurations of previous experiments (Pascal et al., Earth Surf. Process. Landf., vol. 46, issue 9, 2021, pp. 1750–1765), we demonstrate that our approach is robust and accurately predicts the antidunes’ amplitude, wavelength and celerity. Furthermore, the simulated wall shear stress, a key parameter governing sediment transport, is in excellent agreement with the experimental measurements. The highly resolved data of fluid and particle motion from our simulation approach open new perspectives for detailed studies of morphodynamics in shallow supercritical flows.

We derive equations relating the pressure at a flat seabed and the free-surface profile for steady gravity waves with constant vorticity. The resulting set of nonlinear equations enables the recovery of the free surface from pressure measurements at the bed. Furthermore, the flow vorticity (unknown a priori) is determined solely from the bottom pressure as part of the recovery method. This approach is applicable even in the presence of stagnation points and its efficiency is illustrated via numerical examples.

Vortex rings are critical for thrust production underwater. In the ocean, self-propelled mesozooplankton generate vortices while swimming within a weakly stratified fluid. While large-scale biogenic transport has been observed during vertical migration in the wild and lab experiments, little focus has been given to the evolution of induced vortex rings as a function of their propagation direction relative to the density gradient. In this study, the evolution of an isolated vortex ring crossing the interface of a stable two-layer system is examined as a function of its translation direction with respect to gravity. The vortex ring size and position are visualized using planar laser-induced fluorescence (PLIF) and the induced vorticity field derived from particle image velocimetry (PIV) is examined. It is found that the production of baroclinic vorticity significantly affects the propagation of vortex rings crossing the density interface. As a result, any expected symmetry between vortex rings travelling from dense to light fluids and from light to dense fluids breaks down. In turn, the maximum penetration depth of the vortex ring occurs in the case in which the vortex propagates against the density gradient due to the misalignment of the pressure and density gradients. Our results have far-reaching implications for the characterization of local ecosystems in marine environments.

To explore the direction of inter-scale transfer of scalar variance between subgrid scale (SGS) and resolved scalar fields, direct numerical simulation data obtained earlier from two complex-chemistry lean hydrogen–air flames are analysed by applying Helmholtz–Hodge decomposition (HHD) to the simulated velocity fields. Computed results show backscatter of scalar (combustion progress variable $c$) variance, i.e. its transfer from SGS to resolved scales, even in a highly turbulent flame characterized by a unity-order Damköhler number and a ratio of Kolmogorov length scale to thermal laminar flame thickness as low as 0.05. Analysis of scalar fluxes associated with the solenoidal and potential velocity fields yielded by HHD shows that the documented backscatter stems primarily from the potential velocity perturbations generated due to dilatation in instantaneous local flames, with the backscatter being substantially promoted by a close alignment of the spatial gradient of mean scalar progress variable and the potential-velocity contribution to the local SGS scalar flux. The alignment is associated with the fact that combustion-induced thermal expansion increases local velocity in the direction of $\boldsymbol {\nabla } c$. These results call for development of SGS models capable of predicting backscatter of scalar variance in turbulent flames in large eddy simulations.

A closed expression for the average pressure difference (often called the macroscopic dynamic capillary pressure in the literature) is proposed for two-phase, Newtonian, incompressible, isothermal and creeping flow in homogeneous porous media. This upscaled equation complements the average equations for mass and momentum transport derived in a previous article. Consistently with this work, the expression is derived employing a simplified version of the volume-averaging method that makes use of elements of the adjoint method and Green's formula. The resulting equation for the average pressure difference is novel, as it shows that this quantity is controlled by the pressure gradient (and body forces) in each phase, as well as interfacial effects, and is applicable to situations in which the fluid–fluid interface is not necessarily at its steady position. The effective-medium quantities associated with the sources are all obtained from the solution of a single adjoint (or closure) problem to be solved on a (periodic) unit cell representative of the process. The average pressure difference predicted by the derived expression is validated through excellent comparisons with direct numerical simulations performed in a model porous structure.

Recent studies of pressure-driven flows of dilute polymer solutions in straight channels demonstrated the existence of two-dimensional coherent structures that are disconnected from the laminar state and appear through a subcritical bifurcation from infinity. These travelling-wave solutions were suggested to organise the phase-space dynamics of purely elastic and elasto-inertial chaotic channel flows. Here, we consider a wide range of parameters, covering the purely elastic and elasto-inertial cases, and demonstrate that the two-dimensional travelling-wave solutions are unstable when embedded in sufficiently wide three-dimensional domains. Our work demonstrates that studies of purely elastic and elasto-inertial turbulence in straight channels require three-dimensional simulations, and no reliable conclusions can be drawn from studying strictly two-dimensional channel flows.

Analysis of satellite altimetry and Argo float data leads Ni et al. (J. Geophys. Res., 125, 2020, e2020JC016479) to argue that mesoscale dipoles are widespread features of the global ocean having a relatively uniform three-structure that can lead to strong vertical exchanges. Almost all the features of the composite dipole they construct can be derived from a model for multipoles in the surface quasi-geostrophic equations for which we present a straightforward novel solution in terms of an explicit linear algebraic eigenvalue problem, allowing simple evaluation of the higher radial modes that appear to be present in the observations and suggesting that mass conservation may explain the observed frontogenetic velocities.

Fast microjets can emerge out of liquid pools from the rebounding of drop-impact craters, or when a bubble bursts at its surface. The fastest jets are the narrowest and are a source of aerosols both from the ocean and from a glass of champagne, of importance to climate and the olfactory senses. The most singular jets, which we observe experimentally at a maximum velocity of $137\pm 4\ {\rm m}\ {\rm s}^{-1}$ and a diameter of $12\ \mathrm {\mu }{\rm m}$, under reduced ambient pressure, are produced when a small dimple forms at the crater bottom and rebounds without pinching off a small bubble. The radial collapse and rebounding of this dimple is purely inertial, but highly sensitive to initial conditions. High-resolution numerical simulations reveal a new focusing mechanism, which drives the fastest jet within a converging conical channel, where an entrained air sheet provides effective slip at the outer boundary of the conically converging flow into the jet. This configuration bypasses any viscous cutoff of the jetting speed and explains the extreme sensitivity to initial conditions observed in detailed experiments of the phenomenon.

We consider transverse shear flow over a $2L$-periodic array of air-filled grooves, which are separated by solid slats. We address the scenario where the liquid partially invades the grooves. Following Crowdy's solution of the longitudinal problem (J. Fluid Mech. vol. 925, 2021, R2), we focus upon the simplified geometry of infinitely narrow slats. For flat menisci, the ratio $\lambda$ of the slip length to $L$ depends upon a single geometric parameter, namely the ratio $H$ of the invasion depth to $L$. We analyse the singular limit $H\ll 1$ using matched asymptotic expansions, with an outer region on the scale of the period and an inner region on the scale of a single slat. In the outer region, we employ the complex-variable formulation of Stokes flows, obtaining the complex velocity as a divergent series. This formal representation is transformed to an admissible closed-form expression. In the inner region, the problem is equivalent to that of Stokes flow past a finite line segment. Asymptotic matching yields $\lambda = [\ln (2/{\rm \pi} H)-1]/{\rm \pi}$.

A bio-inspired, passively deployable flap attached to an airfoil by a torsional spring of fixed stiffness can provide significant lift improvements at post-stall angles of attack. In this work, we describe a hybrid active–passive variant to this purely passive flow control paradigm, where the stiffness of the hinge is actively varied in time to yield passive fluid–structure interaction of greater aerodynamic benefit than the fixed-stiffness case. This hybrid active–passive flow control strategy could potentially be implemented using variable-stiffness actuators with less expense compared with actively prescribing the flap motion. The hinge stiffness is varied via a reinforcement-learning-trained closed-loop feedback controller. A physics-based penalty and a long–short-term training strategy for enabling fast training of the hybrid controller are introduced. The hybrid controller is shown to provide lift improvements as high as 136 % and 85 % with respect to the flapless airfoil and the best fixed-stiffness case, respectively. These lift improvements are achieved due to large-amplitude flap oscillations as the stiffness varies over four orders of magnitude, whose interplay with the flow is analysed in detail.

In the context of transition analysis, linear input–output analysis determines the worst-case disturbances to a laminar base flow based on a generic right-hand-side volumetric/boundary forcing term. The worst-case forcing is not physically realizable, and, to our knowledge, a generic framework for posing physically realizable worst-case disturbance problems is lacking. In natural receptivity analysis, disturbances are forced by matching (typically local) solutions within the boundary layer to outer solutions consisting of free-stream vortical, entropic and acoustic disturbances. We pose a scattering formalism to restrict the input forcing to a set of realizable disturbances associated with plane-wave solutions of the outer problem. The formulation is validated by comparing with direct numerical simulations of a Mach 4.5 flat-plate boundary layer. We show that the method provides insight into transition mechanisms by identifying those linear combinations of plane-wave disturbances that maximize energy amplification over a range of frequencies. We also discuss how the framework can be extended to accommodate scattering from shocks and in shock layers for supersonic flow.

‘Freeze-out’ of amplitude growth, i.e. the amplitude growth stagnation of a shocked helium–air interface, is realized through a reflected shock, which produces baroclinic vorticity of the opposite sign to that deposited by the first shock. Theoretically, a model is constructed to calculate the relations among the initial parameters for achieving freeze-out. In particular, if the amplitude growth is within the linear regime at the arrival of the reflected shock, the time interval between the impacts of two shock waves is linearly related to the initial perturbation wavelength, and is independent of the initial perturbation amplitude. Experimentally, an air–SF$_6$ (or air–argon) plane interface is adopted to produce a weak reflected shock. Seven experimental runs with specific initial conditions are examined. For all cases, freeze-out is achieved after the reflected shock impact under the designed conditions.

At high surface Péclet numbers, it is common to associate the presence of surfactants with surface immobilization, where a free surface becomes indistinguishable from a no-slip surface. A different mechanism has recently been proposed for longitudinal shear flow along a unidirectional trench (Baier & Hardt, J. Fluid Mech., vol. 949, 2022, A34) wherein, at high Marangoni numbers, the meniscus spanning the finite-length trench becomes a constant-shear-stress surface due to contamination by incompressible surfactant. That model predicts recirculating interfacial flows on the meniscus, a phenomenon that has been observed experimentally (Song et al., Phys. Rev. Fluids, vol. 3, issue 3, 2018, 033303). By finding an explicit solution to the constant-shear-stress model at all protrusion angles and calculating the effective slip length for a dilute mattress of such surfactant-laden trenches, we show that those effective slip lengths are almost indistinguishable from those for a surface whose menisci have the same deflection but have been completely immobilized (i.e. they are no-slip surfaces). This means that, despite the presence of non-trivial recirculating vortical flows on the menisci, the aggregate slip characteristics of such surfaces are that they have been effectively immobilized. This surprising result underscores the need for caution in comparing theory with experiments based on effective slip properties alone.

We derive a scaling law for the characteristic frequencies of wall pressure fluctuations in swept shock wave/turbulent boundary layer interactions in the presence of cylindrical symmetry, based on analysis of a direct numerical simulations database. Direct numerical simulations in large domains show evidence of spanwise rippling of the separation line, with typical wavelength proportional to separation bubble size. Pressure disturbances around the separation line are shown to be convected at a phase speed proportional to the cross-flow velocity. This information is leveraged to derive a simple model for low-frequency unsteadiness, which extends previous two-dimensional models (Piponniau et al., J. Fluid Mech., vol. 629, 2009, pp. 87–108), and which correctly predicts growth of the typical frequency with the sweep angle. Inferences regarding the typical frequencies in more general swept shock wave/turbulent boundary layer interactions are also discussed.

When a cylinder is mounted on an elastic support within a current, vortex-induced vibrations (VIV) may occur down to a Reynolds number (Re) close to $20$, based on the body diameter ($D$) and inflow velocity ($U$), i.e. below the critical value of $47$ reported for the onset of flow unsteadiness when the body is fixed. The impact of a forced rotation of the elastically mounted cylinder on the system behaviour is explored numerically for $Re \leqslant 30$, over wide ranges of values of the rotation rate (ratio between body surface velocity and $U$, $\alpha \in [0,5]$) and reduced velocity (inverse of the oscillator natural frequency non-dimensionalized by $D$ and $U$, $U^\star \in [2,30]$). The influence of the rotation is not monotonic, but the most prominent effect uncovered in this work is a substantial enhancement of the subcritical-Re, flow-induced vibrations beyond $\alpha =2$. This enhancement is twofold. First, the rotation results in a considerable expansion of the vibration/flow unsteadiness region in the $({Re},U^\star )$ domain, down to $Re=4$. Second, the elliptical orbits described by the rotating body are subjected to a major amplification, with a transition from VIV to responses whose magnitude tends to increase unboundedly with $U^\star$, even though still synchronized with flow unsteadiness. The emergence of such galloping-like oscillations close to the onset of vibrations disrupts the scenario of gradual vibration growth with Re, as amplitudes larger than $10$ body diameters may be observed at $Re=10$.