The helicity is a topological conserved quantity of the Euler equations which imposes significant constraints on the dynamics of vortex lines. In the compressible setting, the conservation law holds only under the assumption that the pressure is barotropic. Let us consider a volume $V$ containing a compressible fluid with density $\rho$, velocity field $\textbf{u}$ and vorticity $\boldsymbol{\omega}$. We show that by introducing a new definition of helicity density $h_{\rho }=(\rho {\boldsymbol {u}})\cdot \mbox {curl}\,(\rho {\boldsymbol {u}})$ the barotropic assumption on the pressure can be removed, although ${\int _{V}} h_{\rho }{\rm d}V$ is no longer conserved. However, we show for the non-barotropic compressible Euler equations that the new helicity density $h_{\rho }$ obeys an entropy-type relation (in the sense of hyperbolic conservation laws) whose flux ${\boldsymbol {J}}_{\rho }$ contains all the pressure terms and whose source involves the potential vorticity $q = \boldsymbol{\omega} \cdot \nabla \rho$. Therefore, the rate of change of ${\int _{V}} h_{\rho }{\rm d}V$ no longer depends on the pressure and is easier to analyse, as it depends only on the potential vorticity and kinetic energy as well as $\mbox {div}\,{\boldsymbol {u}}$. This result also carries over to the inhomogeneous incompressible Euler equations for which the potential vorticity $q$ is a material constant. Therefore, $q$ is bounded by its initial value $q_{0}=q({\boldsymbol {x}},\,0)$, which enables us to define an inverse resolution length scale $\lambda _{H}^{-1}$ whose upper bound is found to be proportional to $\|q_{0}\|_{\infty }^{2/7}$. In a similar manner, we also introduce a new cross-helicity density for the ideal non-barotropic magnetohydrodynamic (MHD) equations.

Gravity currents are a ubiquitous density-driven flow occurring in both the natural environment and in industry. They include: seafloor turbidity currents, primary vectors of sediment, nutrient and pollutant transport; cold fronts; and hazardous gas spills. However, while the energetics are critical for their evolution and particle suspension, they are included in system-scale models only crudely, so we cannot yet predict and explain the dynamics and run-out of such real-world flows. Herein, a novel depth-averaged framework is developed to capture the evolution of volume, concentration, momentum and turbulent kinetic energy from direct integrals of the full governing equations. For the first time, we show the connection between the vertical profiles, the evolution of the depth-averaged flow and the energetics. The viscous dissipation of mean-flow energy near the bed makes a leading-order contribution, and an energetic approach to entrainment captures detrainment of fluid through particle settling. These observations allow a reconsideration of particle suspension, advancing over 50 years of research. We find that the new formulation can describe the full evolution of a shallow dilute current, with the accuracy depending primarily on closures for the profiles and source terms. Critically, this enables accurate and computationally efficient hazard risk analysis and earth surface modelling.

We investigate the momentum fluxes between a turbulent air boundary layer and a growing–breaking wave field by solving the air–water two-phase Navier–Stokes equations through direct numerical simulations. A fully developed turbulent airflow drives the growth of a narrowbanded wave field, whose amplitude increases until reaching breaking conditions. The breaking events result in a loss of wave energy, transferred to the water column, followed by renewed growth under wind forcing. We revisit the momentum flux analysis in a high-wind-speed regime, characterized by the ratio of the friction velocity to wave speed $u_\ast /c$ in the range $[0.3\,{-}\,0.9]$, through the lens of growing–breaking cycles. The total momentum flux across the interface is dominated by pressure, which increases with $u_\ast /c$ during growth and reduces sharply during breaking. Drag reduction during breaking is linked to airflow separation, a sudden acceleration of the flow, an upward shift of the mean streamwise velocity profile and a reduction in Reynolds shear stress. We characterize the reduction of pressure stress and flow acceleration through an aerodynamic drag coefficient by splitting the analysis between growing and breaking stages, treating them as separate subprocesses. While drag increases with $u_\ast /c$ during growth, it decreases during breaking. Averaging over both stages leads to a saturation of the drag coefficient at high $u_\ast /c$, comparable to what is observed at high wind speeds in laboratory and field conditions. Our analysis suggests that this saturation is controlled by breaking dynamics.

We investigate the dynamics and the stability of the incompressible flow past a corrugated dragonfly-inspired airfoil in the two-dimensional (2-D) $\alpha {-}Re$ parameter space, where $\alpha$ is the angle of attack and $Re$ is the Reynolds number. The angle of attack is varied in the range of $-5^{\circ } \leqslant \alpha \leqslant 10^{\circ }$, and $Re$ (based on the free stream velocity and the airfoil chord) is increased up to $Re=6000$. The study relies on linear stability analyses and three-dimensional (3-D) nonlinear direct numerical simulations. For all $\alpha$, the primary instability consists of a Hopf bifurcation towards a periodic regime. The linear stability analysis reveals that two distinct modes drive the flow bifurcation for positive and negative $\alpha$, being characterised by a different frequency and a distinct triggering mechanism. The critical $Re$ decreases as $|\alpha |$ increases, and scales as a power law for large positive/negative $\alpha$. At intermediate $Re$, different limit cycles arise depending on $\alpha$, each one characterised by a distinctive vortex interaction, leading thus to secondary instabilities of different nature. For intermediate positive/negative $\alpha$, vortices are shed from both the top/bottom leading- and trailing-edge shear layers, and the two phenomena are frequency locked. By means of Floquet stability analysis, we show that the secondary instability consists of a 2-D subharmonic bifurcation for large negative $\alpha$, of a 2-D Neimark–Sacker bifurcation for small negative $\alpha$, of a 3-D pitchfork bifurcation for small positive $\alpha$ and of a 3-D subharmonic bifurcation for large positive $\alpha$. The aerodynamic performance of the dragonfly-inspired airfoil is discussed in relation to the different flow regimes emerging in the $\alpha {-}Re$ space of parameters.

Dispersion in spatio-temporal random flows is dominated by the competition between spatial and temporal velocity resets along particle paths. This competition admits a range of normal and anomalous dispersion behaviours characterised by the Kubo number, which compares the relative strength of spatial and temporal velocity resets. To shed light on these behaviours, we develop a Lagrangian stochastic approach for particle motion in spatio-temporally fluctuating flow fields. For space–time separable flows, particle motion is mapped onto a continuous time random walk (CTRW) for steady flow in warped time, which enables the upscaling and prediction of the large-scale dispersion behaviour. For non-separable flows, we measure Lagrangian velocities in terms of a new sampling variable, the average number of velocity transitions (both temporal and spatial) along pathlines, which renders the velocity series Markovian. Based on this, we derive a Lagrangian stochastic model that represents particle motion as a coupled space–time random walk, that is, a CTRW for which the space and time increments are intrinsically coupled. This approach sheds light on the fundamental mechanisms of particle motion in space–time variable flows, and allows for its systematic quantification. Furthermore, these results indicate that alternative strategies for the analysis of Lagrangian velocity data using new sampling variables may facilitate the identification of (hidden) Markov models, and enable the development of reduced-order models for otherwise complex particle dynamics.

Surfactant transport is central to a diverse range of natural phenomena with numerous practical applications in physics and engineering. Surprisingly, this process remains relatively poorly understood at the molecular scale. Here, we use non-equilibrium molecular dynamics (NEMD) simulations to study the spreading of sodium dodecyl sulphate on a thin film of liquid water. The molecular form of the control volume is extended to a coordinate system moving with the liquid–vapour interface to track surfactant spreading. We use this to compare the NEMD results to the continuum description of surfactant transport on an interface. By including the molecular details in the continuum model, we establish that the transport equation preserves substantial accuracy in capturing the underlying physics. Moreover, the relative importance of the different mechanisms involved in the transport process is identified. Consequently, we derive a novel exact molecular equation for surfactant transport along a deforming surface. Close agreement between the two conceptually different approaches, i.e. NEMD simulations and the numerical solution of the continuum equation, is found as measured by the surfactant concentration profiles, and the time dependence of the so-called spreading length. The current study focuses on a relatively simple specific solvent–surfactant system, and the observed agreement with the continuum model may not arise for more complicated industrially relevant surfactants and anti-foaming agents. In such cases, the continuum approach may fail to predict accompanying phase transitions, which can still be captured through the NEMD framework.

We investigate turbulent flow between two concentric cylinders, oriented either axially or azimuthally. The axial configuration corresponds to a concentric annulus, where curvature is transverse to the flow, while the azimuthal configuration represents a curved channel with longitudinal curvature. Using direct numerical simulations, we examine the effects of both types of curvature on turbulence, varying the inner radius from $r_i=0.025\delta$ to $r_i=95.5\delta$, where $\delta$ is the gap width. The bulk Reynolds number, based on bulk velocity and $\delta$, is set at $R_b\approx 5000$, ensuring fully turbulent conditions. Our results show that transverse curvature, although breaking the symmetry of axial flows, induces limited changes in the flow structure, leading to an increase in friction at the inner wall. In contrast, longitudinal curvature has a significant impact on the structure and statistics of azimuthal flows. For mild to moderate longitudinal curvatures ($r_i\gt 1.5\delta$), the convex wall stabilises the flow, reducing turbulence intensity, wall friction and turbulent kinetic energy (TKE) production. For extreme longitudinal curvatures ($r_i\leqslant 0.25\delta$), spanwise-coherent flow structures develop near the inner wall, leading to a complete redistribution of the TKE budget: production becomes negligible near the inner wall, while pressure–velocity correlations increase substantially. As a result, the mean TKE peaks near the inner wall, thereby weakening the stabilising effect of convex curvature.

We investigate the emergent three-dimensional (3-D) dynamics of a rapidly yawing spheroidal swimmer interacting with a viscous shear flow. We show that the rapid yawing generates non-axisymmetric emergent effects, with the active swimmer behaving as an effective passive particle with two orthogonal planes of symmetry. We also demonstrate that this effective asymmetry generated by the rapid yawing can cause chaotic behaviour in the emergent dynamics, in stark contrast to the emergent dynamics generated by rapidly rotating spheroids, which are equivalent to those of effective passive spheroids. In general, we find that the shape of the equivalent effective particle under rapid yawing is different to the average shape of the active particle. Moreover, despite having two planes of symmetry, the equivalent passive particle is not an ellipsoid in general, except for specific scenarios in which the effective shape is a spheroid. In these scenarios, we calculate analytically the equivalent aspect ratio of the effective spheroid. We use a multiple scales analysis for systems to derive the emergent swimmer behaviour, which requires solving a non-autonomous nonlinear 3-D dynamical system, and we validate our analysis via comparison to numerical simulations.

The interaction between porous structures and flows with mean and oscillatory components has many applications in fluid dynamics. One such application is the hydrodynamic forces on offshore jacket structures from waves and current, which have been shown to give a significant blockage effect, leading to a reduction in drag forces. To better understand this, we derived analytical expressions that describe the effect of current on drag forces from large waves, and conducted experiments that measured forces on a model jacket in collinear waves and currents. We utilised symmetry and phase-inversion techniques, relying on the underlying physics of wave structure interaction, to separate Morison drag and inertia-type forces and to decompose these forces into their respective frequency harmonics. We find that the odd harmonics of the drag force mostly contain the loads from waves, while even harmonics vary much more rapidly with the current speed flowing through the jacket. At the time of peak force, these current speeds were estimated to be 40 % of the undisturbed current and 50 % of the industry-standard estimates, a result that has significant implications for design and re-assessment of jackets. At times away from the peak force, when there are no waves and only current, the blockage effects are reduced. Hence, the variation in blocked current speeds appears to occur on a relatively fast time scale similar to the compact wave envelope. These findings may be generalisable to any jacket-type structure in flows with mean and high Keulegan–Carpenter number oscillatory components.

We investigate the effects of bottom roughness on bottom boundary-layer (BBL) instability beneath internal solitary waves (ISWs) of depression. Applying both two-dimensional (2-D) numerical simulations and linear stability theory, an extensive parametric study explores the effect of the Reynolds number, pressure gradient, roughness (periodic bump) height $h_b$ and roughness wavelength $\lambda _b$ on BBL instability. The simulations show that small $h_b$, comparable to that of laboratory-flume materials ($\sim$100 times less than the thickness of the viscous sublayer $\delta _v$), can destabilize the BBL and trigger vortex shedding at critical Reynolds numbers much lower than what occurs for numerically smooth surfaces. We identify two mechanisms of vortex shedding, depending on $h_b/\delta _v$. For $h_b/\delta _v \gtrapprox 1$, vortices are forced directly by local flow separation in the lee of each bump. Conversely, for $h_b/\delta _v \lessapprox 10^{-1}$ the roughness seeds perturbations in the BBL, which are amplified by the BBL flow. Roughness wavelengths close to those associated with the most unstable BBL mode, as predicted by linear instability theory, are preferentially amplified. This resonant amplification nature of the BBL flow, beneath ISWs, is consistent with what occurs in a BBL driven by surface solitary waves and by periodic monochromatic waves. Using the $N$-factor method for Tollmien–Schlichting waves, we propose an analogy between the roughness height and seed noise required to trigger instability. Including surface roughness, or more generally an appropriate level of seed noise, reconciles the discrepancies between the vortex-shedding threshold observed in the laboratory versus that predicted by otherwise smooth-bottomed 2-D spectral simulations.

Turbidity currents, which are stratified, sediment-laden bottom flows in the ocean or lakes, can run out for hundreds or thousands of kilometres in submarine channels without losing their stratified structure. Here, we derive a layer-averaged, two-layer model for turbidity currents, specifically designed to capture long-runout. A number of previous models have captured runout of only tens of kilometres, beyond which thickening of the flows becomes excessive, and the models without a lateral overspill mechanism fail. In our framework, a lower layer containing nearly all the sediment is a faster, gravity-driven flow that propels an upper layer, where sediment concentration is nearly zero. The thickness of the lower layer is controlled by competition between interfacial water entrainment due to turbulent mixing and water detrainment due to sediment settling at the interface. The detrainment mechanism, first identified in experiments, is the key feature that prevents excessive thickening of the lower layer and allows long-runout. Under normal flow conditions, we obtain an exact solution to the two-layer formulation, revealing a constant velocity and a constant thickening rate in each of the two layers. Numerical simulations applied

to gradually varied flows on both constant and exponentially declining bed slopes, with boundary conditions mimicking field observations, show that the predicted lower layer thickness after 200 km flow propagation compares with observed submarine channel depths, whereas previous models overestimate this thickness three- to fourfold. This formulation opens new avenues for modelling the fluid mechanics and morphodynamics of long-runout turbidity currents in the submarine setting.

Double-diffusive convection can arise when the fluid density is set by multiple species which diffuse at different rates. Different flow regimes are possible depending on the distribution of the diffusing species, including salt fingering and diffusive convection. Flows arising from diffusive convection commonly exhibit step-like density profiles with sharp density interfaces separated by well-mixed layers. The formation of density layers is also seen in stratified turbulence, where a framework based on sorted density coordinates (Winters & D’Asaro 1996 J. Fluid Mech. 317, 179–193) has been used to diagnose layer formation (Zhou et al. 2017 J. Fluid Mech. 823, 198–229; Taylor & Zhou 2017 J. Fluid Mech. 823, R5). In this framework, the evolution of the sorted density profile can be expressed solely in terms of the eddy diffusivity, $\kappa _e$. Here, we use the same framework to diagnose layer formation in two-dimensional simulations of double-diffusive convection. We show that $\kappa _e$ is negative everywhere, with the antidiffusive effect strongest in ‘well-mixed’ layers where a positive diffusion coefficient may be expected. By considering a decomposition of the budget of the square of the Brunt-Väisälä frequency $\partial N^2_*/\partial t$, we demonstrate that the density layers are maintained by fundamentally different processes than in single-component stratified turbulence. In more complicated flows where stratified turbulence and double-diffusive convection can coexist, this framework could provide a method to distinguish between the mechanisms responsible for generating density layers.

The flow over a cambered NACA 65(1)–412 airfoil at $Re\,=\,2\times 10^4$ is described based on a high-order direct numerical simulation. Simulations are run over a range of angles of attack, $\alpha$, where a number of instabilities in the unsteady, three-dimensional flow field are identified. The balance and competing effects of these instabilities are responsible for significant and abrupt (with respect to $\alpha$) changes in flow regime, with measurable consequences in time-averaged, integrated force coefficients, and in the far-wake footprint. At low $\alpha$, the flow is strongly influenced by vortex roll-up from the pressure side at the trailing edge. The interaction of this large-scale structure with shear and three-dimensional modal instabilities in the separated shear layer and associated wake region on the suction side, explains the transitions and bifurcations of the the flow states as $\alpha$ increases. The transition from a separation at low $\alpha$ to reattachment and establishment of a laminar separation bubble at the trailing edge at critical $\alpha$ is driven by instabilities within the separated shear layer that are absent at lower angles. Instabilities of different wavelengths are then shown to pave the path to turbulence in the near wake.

In soft porous media, deformation drives solute transport via the intrinsic coupling between flow of the fluid and rearrangement of the pore structure. Solute transport driven by periodic loading, in particular, can be of great relevance in applications including the geomechanics of contaminants in the subsurface and the biomechanics of nutrient transport in living tissues and scaffolds for tissue engineering. However, the basic features of this process have not previously been systematically investigated. Here, we fill this hole in the context of a one-dimensional model problem. We do so by expanding the results from a companion study, in which we explored the poromechanics of periodic deformations, by introducing and analysing the impact of the resulting fluid and solid motion on solute transport. We first characterise the independent roles of the three main mechanisms of solute transport in porous media – advection, molecular diffusion and hydrodynamic dispersion – by examining their impacts on the solute concentration profile during one loading cycle. We next explore the impact of the transport parameters, showing how these alter the relative importance of diffusion and dispersion. We then explore the loading parameters by considering a range of loading periods – from slow to fast, relative to the poroelastic time scale – and amplitudes – from infinitesimal to large. We show that solute spreading over several loading cycles increases monotonically with amplitude, but is maximised for intermediate periods because of the increasing poromechanical localisation of the flow and deformation near the permeable boundary as the period decreases.

The study examines supersonic square jets in a twin nozzle configuration with the aim of identifying and characterising emergent instability modes during overexpanded operation. Unlike screeching rectangular jets that undergo strong fluctuations normal to the wider jet dimension, the equilateral nature of the exit geometry in square nozzles leads to multiple instability states dictated by shock–turbulence interactions and nozzle operating conditions. Furthermore, strong coupling modes between the jets were identified that led to either phase locked or out of phase interactions of the inner shear layers. Results from experimental studies were examined using spatial and temporal decomposition techniques based on spectral methods to identify the resultants from triadic shock–turbulence interactions. The primary instability mode across both operating conditions were driven by optimal interactions while the harmonics were found to be associated with the suboptimal shock–turbulence interactions.

A millimetric droplet may bounce and self-propel across the surface of a vertically vibrating liquid bath, guided by the slope of its accompanying Faraday wave field. The ‘walker’, consisting of a droplet dressed in a quasi-monochromatic wave form, is a spatially extended object that exhibits many phenomena previously thought exclusive to the quantum realm. While the walker dynamics can be remarkably complex, steady and periodic states arise in which the energy added by the bath vibration necessarily balances that dissipated by viscous effects. The system energetics may then be characterised in terms of the exchange between the bouncing droplet and its guiding or ‘pilot’ wave. We here characterise this energy exchange by means of a theoretical investigation into the dynamics of the pilot-wave system when time-averaged over one bouncing period. Specifically, we derive simple formulae characterising the dependence of the droplet’s gravitational potential energy and wave energy on the droplet speed. Doing so makes clear the partitioning between the gravitational, wave and kinetic energies of walking droplets in a number of steady, periodic and statistically steady dynamical states. We demonstrate that this partitioning depends exclusively on the ratio of the droplet speed to its speed limit.

Analytical expressions are derived for the velocity field, and effective slip lengths, associated with pressure-driven longitudinal flow in a circular superhydrophobic pipe whose boundary is patterned with a general arrangement of longitudinal no-shear stripes not necessarily possessing any rotational symmetry. First, the flow in a superhydrophobic pipe with $M$ different no-shear stripes in general position is found for $M=1, 2, 3$. The method, which is based on use of so-called prime functions, is such that with these cases covered, generalisation to any $M \geqslant 1$ follows in a straightforward manner. It is shown how any one of these solutions can be generalised to solve for flow along superhydrophobic pipes where that pattern of $M$ menisci is repeated $N \geqslant 1$ times around the boundary in a rotational symmetric arrangement. The work provides an extension of the canonical pipe flow solution for an $N$-fold rotationally symmetric pattern of no-shear stripes due to Philip (Angew. Math. Phys., vol. 23, 1972, pp. 353–372). The novel solution method, and the solutions that it produces, have significance for a wide range of mixed boundary value problems involving Poisson’s equation arising in other applications.

We conducted a series of pore-scale numerical simulations on convective flow in porous media, with a fixed Schmidt number of 400 and a wide range of Rayleigh numbers. The porous media are modeled using regularly arranged square obstacles in a Rayleigh–Bénard (RB) system. As the Rayleigh number increases, the flow transitions from a Darcy-type regime to an RB-type regime, with the corresponding $Sh$–$Ra_D$ relationship shifting from sublinear scaling to the classical 0.3 scaling of RB convection. Here, $Sh$ and $Ra_D$ represent the Sherwood number and the Rayleigh–Darcy number, respectively. For different porosities, the transition begins at approximately $Ra_D = 4000$, at which point the characteristic horizontal scale of the flow field is comparable to the size of a single obstacle unit. When the thickness of the concentration boundary layer is less than approximately one-sixth of the pore spacing, the flow finally enters the RB regime. In the Darcy regime, the scaling exponent of $Sh$ and $Ra_D$ decreases as porosity increases. Based on the Grossman–Lohse theory (J. Fluid Mech. vol. 407, 2000, pp. 27–56; Phys. Rev. Lett. vol. 86, 2001, p. 3316), we provide an explanation for the scaling laws in each regime and highlight the significant impact of mechanical dispersion effects during the development of the plumes. Our findings provide some new insights into the validity range of the Darcy model.

We analyse the steady viscoelastic fluid flow in slowly varying contracting channels of arbitrary shape and present a theory based on the lubrication approximation for calculating the flow rate–pressure drop relation at low and high Deborah ($De$) numbers. Unlike most prior theoretical studies leveraging the Oldroyd-B model, we describe the fluid viscoelasticity using a FENE-CR model and examine how the polymer chains’ finite extensibility impacts the pressure drop. We employ the low-Deborah-number lubrication analysis to provide analytical expressions for the pressure drop up to $O(De^4)$. We further consider the ultra-dilute limit and exploit a one-way coupling between the parabolic velocity and elastic stresses to calculate the pressure drop of the FENE-CR fluid for arbitrary values of the Deborah number. Such an approach allows us to elucidate elastic stress contributions governing the pressure drop variations and the effect of finite extensibility for all $De$. We validate our theoretical predictions with two-dimensional numerical simulations and find excellent agreement. We show that, at low Deborah numbers, the pressure drop of the FENE-CR fluid monotonically decreases with $De$, similar to the previous results for the Oldroyd-B and FENE-P fluids. However, at high Deborah numbers, in contrast to a linear decrease for the Oldroyd-B fluid, the pressure drop of the FENE-CR fluid exhibits a non-monotonic variation due to finite extensibility, first decreasing and then increasing with $De$. Nevertheless, even at sufficiently high Deborah numbers, the pressure drop of the FENE-CR fluid in the ultra-dilute and lubrication limits is lower than the corresponding Newtonian pressure drop.