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.

Linear and nonlinear stability analyses are performed to determine critical Rayleigh numbers (${Ra}_{cr}$) for a Rayleigh–Bénard convection configuration with an imposed bottom boundary heat flux that varies harmonically in time with zero mean. The ${Ra}_{cr}$ value depends on the non-dimensional frequency $\omega$ of the boundary heat-flux modulation. Floquet theory is used to find ${Ra}_{cr}$ for linear stability, and the energy method is used to find ${Ra}_{cr}$ for two different types of nonlinear stability: strong and asymptotic. The most unstable linear mode alternates between synchronous and subharmonic frequencies at low $\omega$, with only the latter at large $\omega$. For a given frequency, the linear stability ${Ra}_{cr}$ is generally higher than the nonlinear stability ${Ra}_{cr}$, as expected. For large $\omega$, ${Ra}_{cr} \omega ^{-2}$ approaches an $O(10)$ constant for linear stability but zero for nonlinear stability. Hence the domain for subcritical instability becomes increasingly large with increasing $\omega$. The same conclusion is reached for decreasing Prandtl number. Changing temperature and/or velocity boundary conditions at the modulated or non-modulated plate leads to the same conclusions. These stability results are confirmed by selected direct numerical simulations of the initial value problem.

We study the inverse flux of waves in one of the simplest geophysical fluid dynamics models: one-dimensional rotating shallow water equations. Based on direct numerical integration of the governing equations, we find that waves injected at small scales get transferred upscale predominantly via resonant quartic interactions between wave modes. The waves’ upscale transfer is non-local and involves turbulent transfer between disparate scales of the flow. Our analysis reveals that the upscale transfer of waves is extremely intermittent and is a result of localized-in-time bursts in wave action flux. These intermittent events of flux bursts lead to shallower waves’ spectrum and relatively higher amplitude wave fields in physical space. On examining statistics of the flow fields, we find that low-energy high wavenumbers more or less comply with the assumptions used in wave turbulence theory, such as uniformly distributed wave phases and the Gaussian distribution of fields, while non-uniform distribution of wave phases and non-Gaussian statistics dominate at large scales or low wavenumbers that contain a major share of the flow energy. Our findings point out that the one-dimensional rotating shallow water equations, despite being a simple geophysical fluid dynamic model, harbour complex and intricate features associated with the upscale transfer of waves that have not been recognized in the past.

Three wave resonant triad interactions in two space and one time dimensions form a well-known system of first-order quadratically nonlinear evolution equations that arise in many areas of physics. In deep water waves, they were first derived by Simmons in 1969 and later shown to be exactly solvable by Ablowitz & Haberman in 1975. Specifically, integrability was established by introducing a system of six wave interactions whose symmetry reduction leads to the well-known three wave equations. Here, it is shown that the six wave interaction and classical three wave equations satisfying triad resonance conditions in finite-depth gravity waves can be derived from the non-local integro-differential formulation of the free surface gravity wave equation with surface tension. These quadratically nonlinear six wave interaction equations and their reductions to the classical and non-local complex as well as real reverse space–time three wave interaction equations are integrable. Limits to infinite and shallow water depth are also discussed.

Supraglacial lakes play a central role in storing melt water, enhancing surface melt and ultimately in driving ice flow and ice shelf melt through injecting water into the subglacial environment and through facilitating fracturing. Here, we develop a model for the drainage of supraglacial lakes through the dissipation-driven incision of a surface channel. The model consists of the St Venant equations for flow in the channel, fed by an upstream lake reservoir, coupled with an equation for the evolution of channel elevation due to advection, uplift and downward melting. After reduction to a ‘stream power’-type hyperbolic model, we show that lake drainage occurs above a critical rate of water supply to the lake due to the backward migration of a shock that incises the lake seal. The critical water supply rate depends on advection velocity and uplift (or more precisely, drawdown downstream of the lake) as well as model parameters such as channel wall roughness and the parameters defining the relationship between channel cross-section and wetted perimeter. Once lake drainage does occur, it can either continue until the lake is empty, or terminate early, leading to oscillatory cycles of lake filling and draining, with the latter favoured by large lake volumes and relatively small water supply rates.

The large-scale pulsations of shock-induced separation with length scale that significantly exceeds the incoming boundary layer thickness are investigated. The shock–boundary layer interaction (SBLI) unit is generated by an inward-turning axisymmetric compression ramp at an inflow Mach number of 2.5. A substantial region surrounding the centre azimuth exhibited mean and dynamic flow features that are consistent with two-dimensional separation. Two-dimensional highly resolved maps of surface pressure field are obtained using fast-response pressure-sensitive paint fluorescence imaging at 40 kHz repetition rate. The measurement domain covered significant regions of the incoming boundary layer through the relaxing boundary layer downstream of the reattachment as well as over 25 boundary layer thicknesses in the azimuthal direction. These measurements provide new insights into the spanwise coupling of the SBLI unit in addition to its inherent dynamics. The power spectral density (PSD) of the centreline pressure exhibits very good agreement with theoretical models and complementary measurements using fast-response pressure transducers, which served to validate the pressure field measurements. Detailed examination of the PSD reveals strong agreement with the literature, which includes the peak Strouhal number of the separation and reattachment shock motions as well as the downward frequency shift along the separation bubble. Furthermore, the pressure fluctuation maps reveal streamwise-elongated structures just downstream of the ramp leading edge that persist well downstream of the reattachment. A time sequence of conditional average pressure fluctuation maps is constructed surrounding isolated pressure excursions in the intermittent region. This sequence, along with two-point cross correlation analysis, provides critical information about the flow processes that drive the separation bubble pulsations in the SBLI units with large separation scales. Overall, the imbalance in the mass within the separation bubble appears to be the critical mechanism that drives the separation bubble pulsations. Furthermore, the pressure perturbations originating at azimuthally offset locations are also observed to influence the separation bubble dynamics.

We propose a compressible multi-scale model that (i) captures the dynamics of both large vapour cavities (resolved vapour) and micro-bubbles (unresolved vapour), and (ii) accounts for medium compressibility. The vapour mass, momentum and energy in the compressible homogeneous mixture equations are explicitly decomposed into constituent resolved and unresolved components that are independently treated. The homogeneous mixture of liquid and resolved vapour is tracked as a continuum in an Eulerian sense. The unresolved vapour terms are expressed in terms of subgrid bubble velocities and radii that are tracked in a Lagrangian sense using a novel ‘$kR$-$RP$ equation’ (k, constant multiple; R, bubble size; RP, Rayleigh-Plesset). The $kR$-$RP$ equation is formally derived in terms of the pressure at a finite distance ($kR$) from the bubble while accounting for the effects of neighbouring bubbles; $p(kR)$ may therefore be either a near-field or far-field pressure. The equation exactly recovers the classical Rayleigh–Plesset and Keller–Miksis equations in the limits that $k$ and $c$ (speed of sound) become very large. Also, the results are independent of $k$ for a single bubble for all $k$, and for multiple bubbles when $kR < d$ (where $d$ denotes separation distance). Numerical results show this robustness of the model to the choice of $k$, which can be different for each bubble. The multi-scale model is validated for the collapse of a single resolved/unresolved bubble. Its ability to capture inter-bubble interactions is demonstrated for multiple bubbles exposed to an acoustic pulse. The model is then applied to a problem where resolved and unresolved bubbles co-exist. Finally, it is validated using a cluster of $1200$ bubbles exposed to a strong acoustic pulse. The results show the impact of the bubble cluster on the transmitted and reflected waves and the shielding effect where bubbles at the edge of the cluster shield the interior bubbles by dampening the incident acoustic wave.

The near-field characteristics of highly buoyant plumes, commonly referred to as lazy plumes, remain relatively poorly understood across a range of flow conditions, particularly compared with our understanding of far-field characteristics. Here, we perform fully resolved three-dimensional numerical simulations of round helium plumes to characterize the effects of different inlet Richardson, ${Ri}_0$, and Reynolds, ${Re}_0$, numbers on first- and second-order statistical moments as well as average vertical fluxes in the near field. For sufficiently high ${Re}_0$ at a particular ${Ri}_0$, heavy air can penetrate the core of the plume, reminiscent of spikes in the classical Rayleigh–Taylor instability. In the most turbulent simulation, this penetration becomes so strong that a recirculation zone forms along the centreline of the plume. Vertical fluxes are found to scale linearly with vertical distance from the plume inlet, consistent with experimental and numerical observations (Jiang & Luo, Flow Turbul. Combust., vol. 64, 2000, pp. 43–69; Kaye & Hunt, Intl J. Heat Fluid Flow, vol. 30, 2009, pp. 1099–1105). We analytically derive this linear scaling from the governing equations by making a radial entrainment hypothesis whereby ambient fluid is entrained, on average, only in the radial direction at a finite distance from the inlet. Through this derivation, we identify physical mechanisms that can cause these relationships to remain only approximately valid for the present simulations. Lastly, we identify near-field power-law scaling relations for the flux magnitudes based on ${Ri}_0$, and also examine vertical profiles of the non-dimensional Richardson number flux. Ultimately, insights from the present simulations are used to define near-, intermediate- and far-field regions in buoyant plumes.

We study the effect of insoluble surfactants on the impact of surfactant-free droplets onto surfactant-laden thin liquid films via a fully three-dimensional direct numerical simulation approach that employs a hybrid interface-tracking/level-set method, and by taking into account surfactant-induced Marangoni stresses due to gradients in interfacial surfactant concentration. Our numerical predictions for the temporal evolution of the surfactant-free crown are validated against the experimental work by Che & Matar (Langmuir, vol. 33, 2017, pp. 12140–12148). We focus on the ‘crown-splash regime’, and we observe that the crown dynamics evolves through various stages: from the growth of linear modes (through a Rayleigh–Plateau instability) to the development of nonlinearities leading to primary and secondary breakup events (through droplet shedding modulated by an end-pinching mechanism). We show that the addition of surfactants does not affect the wave selection via the Rayleigh–Plateau instability. However, the presence of surfactants plays a key role in the late stages of the dynamics as soon as the ligaments are driven out from the rim. Surfactant-induced Marangoni stresses delay the end-pinching mechanisms to result in longer ligaments prior to their capillary singularity. Our results indicate that Marangoni stresses bridge the gap between adjacent protrusions promoting the adjacent protrusions' collision and the merging of ligaments. Finally, we demonstrate that the addition of surfactants leads to surface rigidification and consequently to the retardation of the flow dynamics.

Although natural and industrial flows often transport and mix multiple scalars, relatively few studies of turbulent multi-scalar mixing have been undertaken. In the present work, a novel three-wire thermal-anemometry-based probe – capable of simultaneously measuring velocity, helium concentration and temperature – is used to investigate the evolution of multiple scalars ($\phi _1$, $\phi _2$, $\phi _3$) and velocity in turbulent coaxial jets. The jets consist of (i) a centre jet containing a mixture of helium and air ($\phi _1$), (ii) an annular jet containing pure (unheated) air ($\phi _2$) and (iii) a coflow of (pure) heated air ($\phi _3$). Axial measurements are made at three different momentum flux ratios ($M=0.77, 2.1, 4.2$). Increasing $M$ was observed to result in complex, competing effects. Larger momentum flux ratios cause the potential core of the centre jet to decrease in size, greater scalar fluctuations and more rapid correlation of $\phi _1$ and $\phi _2$. However, at the same time, certain statistics, including those describing the velocity field, evolve more slowly. Moreover, the flow near the beginning of the fully merged region appears to be less mixed at higher values of $M$. The present work finally demonstrates that differences can be observed in the evolution and mixing of coaxial jets between those in which $M<1$ and those in which $M>1$, thereby presenting an opportunity by which the mixing process in coaxial jets may be controlled.

A linear stability analysis is employed to investigate thermal effects in shear flows. The cases analysed consist of unstably stratified horizontal boundary layers under a mixed convection regime, where forced and free convection mechanisms compete. Governing equations are given by the incompressible Navier–Stokes equations with the Oberbeck–Boussinesq approximation, where the base flow comes from their boundary layer approximation. Modal and non-modal analyses are used to investigate the behaviour of small-amplitude disturbances superposed to this base flow. An evaluation of the inertial, shearing and buoyancy mechanisms in the mixed convection boundary layer stability is performed through variations in the Reynolds, Prandtl and Richardson numbers. On the one hand, the spectra lead to the parametric conditions for the time-asymptotic onset of instability, which is still caused by Tollmien–Schlichting (TS) waves as in the traditional Blasius case. However, thermal effects have a destabilizing effect on them, more so for liquids than gases. On the other hand, the pseudospectra obtained from a resolvent analysis indicate the existence of transient growth at this same onset. However, contrary to the traditional Blasius case, thermal effects cause it to be dominated by the continuous frequency spectrum instead of the discrete TS modes. In order to elucidate this qualitative change, a componentwise input–output analysis is employed to quantify the receptivity to specific external disturbances. It shows that thermal effects directly impact the conversion of thermal to kinetic linear disturbance energy, causing a strong amplification of the flow response due to the non-normality of the linear operator. Results reveal that heating from below causes the forcing and response modes of the input–output analysis to have a free-stream spatial support due to non-modal excitation of the continuous spectrum. Such a behaviour suggests that the unstably stratified boundary layer is susceptible to free-stream thermal disturbances, which can potentially impact bypass transition.

The dynamics of thin, non-circular droplets evaporating in the diffusion-limited regime is examined. The challenging non-rectilinear mixed boundary problem this poses is solved using a novel asymptotic approach and an asymptotic expansion for the evaporative flux from the free surface of the droplet is found. While theoretically valid only for droplets that are close to circular, it is demonstrated that the methodology can successfully be applied to droplets with a wide variety of footprint shapes, including polygons and highly non-convex domains. As our solution for the flux fundamentally represents a novel result in potential theory, the applications are numerous, as the mixed boundary value problem arises in fields as diverse as electrostatics and contact mechanics. Here, we demonstrate the practicality of our result by considering the analytically tractable case of deposition of solute from large droplets in detail, including a matched asymptotic analysis to resolve the pressure, streamlines and deposition up to second order.

In a dilute suspension where sinking spheroids or motile gyrotactic micro-organisms are modelled as orientable and negatively buoyant particles, we have found analytical solutions to their steady distributions under any arbitrary continuous vertical shear flow. The two-way coupling between their distribution and the vertical flow is nonlinear, enabling the uniform base state to bifurcate into a structure reminiscent of the streamers in settling spheroid suspensions and gyrotactic plumes. This bifurcation depends on a single parameter that is proportional to the average number of particles on a horizontal cross-section. In a three-dimensional axisymmetric system, the plume structure blows up when the parameter is above a threshold. We discuss how this singularity is analogous to the chemotactic collapse of a Keller–Segel model, and the significance that this analogy entails.

We revisit the classical problem of calculating the pancake-like shape of a sessile drop at large Bond numbers. Starting from a formulation where drop volume and contact angle are prescribed, we develop an asymptotic scheme which systematically produces approximations to the two key pancake parameters, height and radius. The scheme is based on asymptotic matching of a ‘flat region’ where capillarity is negligible and an ‘edge region’ near the contact line. Major simplifications follow from the distinction between algebraically and exponentially small terms, together with the use of two exact integral relations. The first represents a force balance in the vertical direction. The second, which can be interpreted as a radial force balance on the drop edge (up to exponentially small terms), generalises an approximate force balance used in classical treatments. The resulting approximations for the geometric pancake parameters, which go beyond known leading-order results, are compared with numerical calculations tailored to the pancake limit. These, in turn, are facilitated by an asymptotic approximation for the exponentially small apex curvature, which we obtain using a Wentzel–Kramers–Brillouin method. We also consider the comparable two-dimensional problem, where similar integral balances explicitly determine the pancake parameters in closed form up to an exponentially small error.

Thin-film flow down a fibre exhibits rich dynamics and is relevant to applications such as desalination, fibre coating and fog harvesting. These flows are subject to instabilities that result in dynamic bead-on-fibre patterns. We perform an experimental study of shear-thinning flow down fibres using 20 different xanthan gum solutions as our working liquid. The bead-on-fibre morphology can be oriented either symmetrically or asymmetrically on the fibre, and this depends upon the surface tension, fibre diameter and liquid rheology, as defined by the Ostwald power-law index. For highly shear-thinning liquids, it is possible for the pattern to be complex and exhibit simultaneously both asymmetric large beads and symmetric small beads in the isolated and convective flow regimes. We quantify the transition between flow regimes and bead dynamics for the asymmetric morphology, and compare with Newtonian flow, as it depends upon the experimental parameters. Finally, the dimensionless bead frequency is shown to scale with the Bond number for all of our experimental data (symmetric and asymmetric).

The effect of surfactants on the flow characteristics during rapid drop formation in a microchannel is investigated using high-speed imaging, micro-particle image velocimetry and numerical simulations; the latter are performed using a three- dimensional multiphase solver that accounts for the transport of soluble surfactants in the bulk and at the interface. Drops are generated in a flow-focusing microchannel, using silicone oil ($4.6$ mPa s) as the continuous phase and a 52 % w/w glycerol solution as the dispersed phase. A non-ionic surfactant (Triton X-100) is dissolved in the dispersed phase at concentrations below and above the critical micelle concentration. Good agreement is found between experimental and numerical data for the drop size, drop formation time and circulation patterns. The results reveal strong circulation patterns in the forming drop in the absence of surfactants, whose intensity decreases with increasing surfactant concentration. The surfactant concentration profiles in the bulk and at the interface are shown for all stages of drop formation. The surfactant interfacial concentration is large at the front and the back of the forming drop, while the neck region is almost surfactant free. Marangoni stresses develop away from the neck, contributing to changes in the velocity profile inside the drop.

We study the propagation of water waves in a three-dimensional device alternating open canals and resonant canals with subwavelength resonances. The dispersion of water waves in such a medium is obtained by analysing the full three-dimensional problem and combining Bloch–Floquet analysis with an asymptotic technique. We obtain the closed forms of the dispersions for resonant canals containing one or two resonators, which depend on only two functions associated with symmetric and antisymmetric modes, and on a geometric parameter analogous to the hopping parameter in topological systems. The analysis of the complete band structure reveals frequency ranges alternating between elliptical and hyperbolic dispersions; in particular, the hyperbolic regime gives rise to a negative effective water depth with a consequent negative refraction. Throughout the course of our study, our theoretical results are validated by comparison with numerical calculations of the full three-dimensional problem.

This study investigates the dynamics of the simultaneous impact of two droplets on a dry substrate. We develop a new micro-controlled droplet generator that releases two equally sized water droplets simultaneously on-demand, with no trailing droplets. The impact Weber number, based on impact droplet size and velocity, and the inter-droplet spacing relative to the impact droplet size are varied in the ranges of 54 to 155 and 1.32 to 2.25, respectively, leading to the strong interaction of the spreading lamellae that form a central uprising sheet, which eventually deposits or breaks into tiny droplets. We analyse the impact processes for both deposition and splashing of the uprising sheets. Simultaneous high-speed imaging from two orthogonal views of the droplet impacts quantifies the three-dimensional structure of the sheet morphology, including the temporal evolution of the rim-bounded ‘semilunar’ shape, surface waves, rim corrugations and finger formation, and deposition or splashing of the liquid sheet. The characteristics of the sheet surface waves and the rim instabilities are quantified. Novel scaling is developed for the maximum sheet height, sheet width and thickness, which considers the geometrical constraints and mass balance of the interacting lamellae to describe the temporal evolution of a ‘semilunar’ uprising sheet and is in good agreement with the measurements. The uprising sheet splashing generates larger droplets than those from splashing of single-droplet impacts, and it occurs due to the end-pinching of sheet fingers and at conditions that single-droplet impacts lead only to liquid deposition.

The root systems of mangroves, a tree species found in intertidal tropical and subtropical coastal zones, provide a natural barrier that dissipates wave energy effectively and reduces sediment erosion. Here, we use a combination of experiments and numerical simulations to examine the wake and drag characteristics of porous arrays of cylinders, which serve as simplified models of mangrove root networks. Optimal arrangements of the arrays are obtained by coupling Navier–Stokes simulations with a multi-objective optimization algorithm, which seeks configurations that minimize wake enstrophy and maximize drag on the porous structure. These optimal configurations are investigated using particle image velocimetry, and the internal and external flows around the porous arrays are analysed using a combination of proper orthogonal decomposition and Lagrangian particle tracking. Large variations in drag and enstrophy are observed by varying the relative positions of the cylinders, which indicates that the geometrical arrangement of porous arrays plays a prominent role in determining wake and drag characteristics. A sensitivity analysis suggests that enstrophy is more sensitive than drag to specific cylinder placement, and depends on distinctive flow patterns that develop in the interior due to interactions among neighbouring cylinders. Arrays with higher drag involve a combination of larger projected frontal area and minimal flux through the interior, leading to increased wake enstrophy, which is unfavourable for particle deposition and erosion. Based on the analysis of characteristics associated with the optimal arrays, several manually designed arrays are tested, and they display the expected behaviour with regard to drag and enstrophy.