On the basis of the generalized Saint-Venant equations, hydraulic jumps in laminar open channel flow are obtained as continuous shock structures. They arise as stable stationary solutions of the governing equations and are particularly well suited for a dynamical systems analysis, where they manifest themselves as near-parabolic trajectories in phase space. We derive an analytic approximate expression for the jump length as a function of the Froude and effective Reynolds numbers, reflecting the fact that – along with inertia – gravity and viscous diffusion both contribute to the balance of forces that shape the jump. The stability of the jumps is confirmed by means of a numerical experiment and an analysis of the propagation and attenuation of small perturbations.

We calculate analytically the self-similar Stokes flow driven by an externally applied electric field in a multiple coaxial Taylor cone consisting of an arbitrary number of immiscible leaky-dielectric or dielectric fluids. The proposed conical solutions open new avenues for innovative technological applications, some of which are briefly discussed.

The net erosion-deposition rate of an avalanche is fundamental to its dynamics and in determining its growth or decay. Small-scale experiments are performed by releasing a given volume of yellow sand onto a stationary erodible red sand layer on a rough inclined plane. Depending on the erodible layer depth and the slope angle, the avalanche is found to either decay, grow, propagate steadily or rapidly shed grains to produce secondary avalanches. The use of different coloured sand with identical properties shows that a particle exchange occurs, which eventually results in a flow that is comprised entirely of particles from the stationary layer rather than the initial release. It is notoriously difficult to model the erosion and deposition processes in granular flows, but it is shown that a two-dimensional depth-averaged avalanche model, with a hysteretic basal friction law, can reproduce all of the observed behaviours. The results illustrate how a continuous exchange of particles with the substrate layer is fundamentally important to the propagation of such avalanches. An investigation into long distance propagation behaviour reveals that avalanches can reach a steady state, the size and speed of which are independent of the initially released volume. In certain conditions avalanches can grow to steady states that are significantly more massive than the flows from which they are originally formed. This paper demonstrates the importance of correctly including erosion-deposition in operational forecast models of snow avalanches and other geophysical mass flows.

We carry out high-resolution particle image velocimetry experiments to characterize the flow field of fractal grids located at the exit section of a turbulent round jet. We comment on the mean flow organization and on the turbulence properties of such jets by comparing the results with those obtained with square grids, a regular grid (having the same effective mesh length) and a jet without turbulator. We find that, different from the case of decaying grid turbulence, a correction must be accounted for to properly scale the turbulence intensity profiles with a length scale based on grid parameters. We perform a low-order reconstruction of the velocity field based on the most energetic proper orthogonal decomposition modes and we compare the flow-field structure produced in the lee of fractal grids with a single square object and the jet without turbulator. The typical turbulence intensity profile detailed in Cafiero et al. (Phys. Fluids, vol. 27, 2015, 115103) for jets with fractal grids is produced by the interaction of small eddies shed by the central grid item. In the single square grid case, the turbulence is built upon the interaction between larger structures. Conversely, the interaction of the outward spreading wake with the external shear layer produces pairs of vortical structures, which we relate to the higher entrainment rate featured by jets with fractal turbulators. The secondary grid iterations have a disruptive effect on the turbulence transport, with a corresponding large correlation between the velocity fluctuations at the jet core with those at the jet shear layer.

We report laboratory results on bubble plumes released from a point source into a homogeneous liquid environment of depth $H$ and rotation rate $\varOmega$. The gas phase is characterised by the non-dimensional slip velocity $u_N=u_s/(BH^{-1})^{1/3}$, where $u_s$ is the slip velocity of the dominant bubble species and $B$ is the source buoyancy flux. The effects of the background rotation are characterised by the Rossby number $Ro=(B\varOmega ^{-3})^{1/4}/H$. We study the regimes $0.06\lessapprox u_N\lessapprox 0.36$ and $0.03\lessapprox Ro\lessapprox 0.3$. We establish that, from $\varOmega t\approx 2.3$, the growth of the maximum plume width $b_{max}$ is slowed down compared to the non-rotating $t^{3/4}$ power law, where $t$ is the time. At $\varOmega t\approx {\rm \pi}$, the plume axis starts to tilt laterally and causes a slowdown of the rise of the height $h_c$ of the plume silhouette centroid from the $t^{3/4}$ power law. These critical times do not depend on $u_N$. After $\varOmega t \approx {\rm \pi}$, the slip velocity counteracts the effects of rotation such that $h_c$ rises faster for larger $u_N$. The subsequent onset of the anticyclonic plume precession causes the disintegration of the rising plume front into vertical columnar structures for $Ro\lessapprox 0.15$. Once the plume reaches the free surface, the subsurface lateral dispersion of bubbles is increased compared to the non-rotating case. However, background rotation suppresses the lateral dispersion of bubbles at the free surface. We find that, asymptotically, the surface area $A$ affected by the bubbles scales as $(B\varOmega ^{-1})^{1/2}t$, and the proportionality factor reduces with an increasing $u_N$.

We found a multi-scale steady solution of the Boussinesq equations for Rayleigh–Bénard convection in a three-dimensional periodic domain between horizontal plates with a constant temperature difference. This was realised using a homotopy from the wall-to-wall optimal transport solution provided by Motoki et al. (J. Fluid Mech., vol. 851, 2018, R4). A connected steady solution, which is a consequence of bifurcation from a thermal conduction state at Rayleigh number $Ra\sim 10^{3}$, is tracked up to $Ra\sim 10^{7}$ using a Newton–Krylov iteration. The three-dimensional exact coherent thermal convection exhibits a scaling of $Nu\sim Ra^{0.31}$ (where $Nu$ is the Nusselt number) as well as multi-scale thermal plume and vortex structures, which are quite similar to those in turbulent Rayleigh–Bénard convection. The mean temperature profiles and the root-mean-square of the temperature and velocity fluctuations are in good agreement with those of the turbulent states. Furthermore, the energy spectrum follows Kolmogorov's $-5/3$ scaling law with a consistent prefactor, and the energy transfer to small scales with a nearly constant flux in the wavenumber space is in accordance with the turbulent energy transfer.

We analyse the model problem of rigid, spherical pendula impacting on rubber membranes at different Reynolds numbers to understand the contact dynamics in deformable bodies in a viscous fluid. We have investigated the problem both by laboratory and numerical experiments and a new contact model has been developed to perform the simulations. We have found that the collision dynamics depends on many parameters, the most important ones being the impact Stokes number ($St$) and the ratio of the membrane thickness ($e$) to sphere diameter ($2R$). For $e/(2R)=0.1$ the pendulum rebounds for $St \geqslant 18$ and this threshold increases as the membrane thickness decreases. Also, the membrane surface affects the dynamics of the impact by delaying the rebound for smaller membranes.

We present a numerical study of the dynamics of an elastic fibre in a shear flow at low Reynolds number, and seek to understand several aspects of the fibre's motion using the equations for slender-body theory coupled to the elastica. The numerical simulations are performed in the bead-spring framework including hydrodynamic interactions in two theoretical schemes: the generalized Rotne–Prager–Yamakawa model and a multipole expansion corrected for lubrication forces. In general, the two schemes yield similar results, including for the dominant scaling features of the shape that we identify. In particular, we focus on the evolution of an initially straight fibre oriented in the flow direction and show that the time scales of fibre bending, curling and rotation, which depend on its length and stiffness, determine the overall motion and evolution of the shapes. We document several characteristic time scales and curvatures representative of the shape that vary as power laws of the bending stiffness and fibre length. The numerical results are further supported by an interpretation using an elastica model.

Celebrating the life and work of Professor George K. Batchelor FRS Founding Editor of the Journal of Fluid Mechanics.

Graphical Abstract © Mark Hallworth. Reproduced with Permission.

Flowing granular materials exhibit several features that distinguish them from molecular fluids. A prominent feature is dilatancy, or volume deformation caused by shear deformation. Its significance in sustained flow has not been much appreciated, as its effect was thought to be confined to thin shear layers. However, it was demonstrated recently by Krishnaraj & Nott (Nat. Commun., vol. 7, 2016, 10630) that dilatancy drives a large-scale secondary flow in a cylindrical Couette device. They hypothesized that the combination of shear and gravity, when their directions are non-collinear, is necessary for the occurrence of the secondary flow. In this paper we investigate the phenomenon by considering a more complex primary flow generated in a split-bottom Couette device, wherein a part of the base of the container filled with a granular material is translated or rotated. It is known from previous studies that the height to which the material is filled determines the shape and extent of the shearing region in the primary flow. We show that the fill height also determines the form of the secondary flow, and argue that the two are intricately coupled and evolve together. Though the secondary flow is more complex than in a cylindrical Couette device, the mechanism is indeed what Krishnaraj & Nott hypothesized: the combined effect of dilatancy driven expansion, and gravity-driven flow down regions of low density. Unlike fluid instabilities that are typically driven by inertia, the secondary flow occurs at arbitrarily low shear rate and appears to be an integral part of the kinematic response.

Capsules, composed of a liquid core protected by a thin deformable membrane, offer high-potential applications in many fields of industry such as bioengineering. One of their limitations comes from the absence of models of capsule damage and/or rupture when they are subjected to an external flow. To assess when rupture is initiated, we develop a fluid–structure interaction (FSI) numerical model of a capsule in Stokes flow that accounts for potential damage of the capsule membrane. We consider the framework of continuum damage mechanics and model the membrane with an isotropic brittle damage model, in which the membrane damage state depends on the history of loading. The FSI problem is solved by coupling the finite element method, to solve for the membrane deformation, with the boundary integral method, to solve for the inner and outer fluid flows. The model is applied to an initially spherical capsule subjected to a simple shear flow. Damage initiates at a critical value of the capillary number, ratio of the fluid viscous forces to the membrane elastic forces and rupture occurs at a higher capillary number, when it reaches a threshold value. The material parameters introduced in the damage model do not influence the mode of damage but only the values of the critical and threshold capillary numbers. When the capillary number is larger than the critical value, damage develops in the two symmetric central regions containing the vorticity axis. It is indeed in these regions that the internal tensions are the highest on the membrane.

Turbulent mixing exerts a significant influence on many physical processes in the ocean. In a stably stratified Boussinesq fluid, this irreversible mixing describes the conversion of available potential energy (APE) to background potential energy (BPE). In some settings the APE framework is difficult to apply and approximate measures are used to estimate irreversible mixing. For example, numerical simulations of stratified turbulence often use triply periodic domains to increase computational efficiency. In this set-up, however, BPE is not uniquely defined and the method of Winters et al. (J. Fluid Mech., vol. 289, 1995, pp. 115–128) cannot be directly applied to calculate the APE. We propose a new technique to calculate APE in periodic domains with a mean stratification. By defining a control volume bounded by surfaces of constant buoyancy, we can construct an appropriate background buoyancy profile $b_\ast (z,t)$ and accurately quantify diapycnal mixing in such systems. This technique also permits the accurate calculation of a finite-amplitude local APE density in periodic domains. The evolution of APE is analysed in various turbulent stratified flow simulations. We show that the mean dissipation rate of buoyancy variance $\chi$ provides a good approximation to the mean diapycnal mixing rate, even in flows with significant variations in local stratification. When quantifying measures of mixing efficiency in transient flows, we find significant variation depending on whether laminar diffusion of a mean flow is included in the kinetic energy dissipation rate. We discuss how best to interpret these results in the context of quantifying diapycnal diffusivity in real oceanographic flows.

A concentrated, vertical monolayer of identical spherical squirmers, which may be bottom heavy, and which are subjected to a linear shear flow, is modelled computationally by two different methods: Stokesian dynamics, and a lubrication-theory-based method. Inertia is negligible. The aim is to compute the effective shear viscosity and, where possible, the normal stress differences as functions of the areal fraction of spheres $\phi$, the squirming parameter $\beta$ (proportional to the ratio of a squirmer's active stresslet to its swimming speed), the ratio $Sq$ of swimming speed to a typical speed of the shear flow, the bottom-heaviness parameter $G_{bh}$, the angle $\alpha$ that the shear flow makes with the horizontal and two parameters that define the repulsive force that is required computationally to prevent the squirmers from overlapping when their distance apart is less than a critical value. The Stokesian dynamics method allows the rheological quantities to be computed for values of $\phi$ up to $0.75$; the lubrication-theory method can be used for $\phi > 0.5$. For non-bottom-heavy squirmers, which are unaffected by gravity, the effective shear viscosity is found to increase more rapidly with $\phi$ than for inert spheres, whether the squirmers are pullers ($\beta > 0$) or pushers ($\beta < 0$); it also varies with $\beta$, although not by very much. However, for bottom-heavy squirmers the behaviour for pullers and pushers as $G_{bh}$ and $\alpha$ are varied is very different, since the viscosity can fall even below that of the suspending fluid for pushers at high $G_{bh}$. The normal stress differences, which are small for inert spheres, can become very large for bottom-heavy squirmers, increasing with $\beta$, and varying dramatically as the orientation $\alpha$ of the flow is varied from 0 to ${\rm \pi} /2$. A major finding is that, despite very different assumptions, the two methods of computation give overlapping results for viscosity as a function of $\phi$ in the range $0.5 < \phi < 0.75$. This suggests that lubrication theory, based on near-field interactions alone, contains most of the relevant physics, and that taking account of interactions with more distant particles than the nearest is not essential to describe the dominant physics.

Coating of filtration membranes with a polymer brush holds great promise for efficiently preventing the deposition of fouling particles. The polymer chains are compressed by incoming particles, carried with the permeation flow towards the membrane, and consequently exert a repulsive force that acts to keep the particles away from the membrane surface. Here, we theoretically investigate the effect of a polymer brush coating on the permeation-induced hydrodynamic force, $F_{h}$, pulling a particle towards the membrane, and its balance with the steric repulsion exerted by the compressing brush, resisting the particle's approach. Lubrication theory yields an ordinary differential equation for the pressure, from which $F_{h}$ is calculated numerically. Further, an asymptotic analysis is performed for the limiting cases of a dilute or dense brush, providing analytic expressions that demonstrate how brush properties affect $F_{h}$. Finally, the equilibrium position of a particle is evaluated by considering a balance between the opposing forces. The results provide an upper boundary for the brush properties, beyond which the brush is barely compressed under conditions typical of membrane filtration processes. Further increasing the brush density or thickness only decreases the total system permeance, resulting in increased energy consumption. The results shed light on the mechanisms by which a polymer brush affects the forces acting on a foulant particle, providing quantitative measures for assessing the potential efficacy of brush coatings.

We consider the solidification of idealised two-component mixtures comprising a solvent or suspending fluid and dissolved solute molecules or suspended colloidal particles, each considered as hard spheres. We review some fundamental thermodynamic ideas regarding relative motion between species and phase equilibria in such mixtures to show how the related solid–liquid phase diagrams depend on the size of the spheres. Using similarity solutions, we first describe freezing of the solvent to form a pure solid (here referred to as ‘ice’), with the solute rejected from the solid forming a boundary layer or dense particle layer ahead of the freezing front. We extend ideas of constitutional supercooling to the case of colloidal suspensions and show that, for a given temperature difference driving solidification, constitutional supercooling occurs only for an intermediate range of particle sizes. Constitutional supercooling promotes the formation of a mushy layer in which segregated ice separates regions of concentrated solute or particles on the microscale. We formulate a continuum model of the mushy layer that relies on a key observation that the regelative motion of concentrated clusters of particles is independent of the size and geometry of the cluster. Our modelling begins with a description of relative motion as a Fickian diffusive process. However, at high particle concentrations, we show that it is more convenient and more computationally tractable to use an equivalent formulation in terms of Darcy flow of the solvent. Within a mushy layer these diffusive fluxes correspond directly to the regelative flux of particle clusters at a rate determined by the local temperature and temperature gradient.

For a small sessile or pendant droplet it is generally assumed that gravity does not play any role once the Bond number is small. This is even assumed for evaporating binary sessile or pendant droplets, in which convective flows can be driven due to selective evaporation of one component and the resulting concentration and thus surface tension differences at the air–liquid interface. However, recent studies have shown that in such droplets gravity indeed can play a role and that natural convection can be the dominant driving mechanism for the flow inside evaporating binary droplets (Edwards et al., Phys. Rev. Lett., vol. 121, 2018, 184501; Li et al., Phys. Rev. Lett., vol. 122, 2019, 114501). In this study, we derive and validate a quasi-stationary model for the flow inside evaporating binary sessile and pendant droplets, which successfully allows one to predict the prevalence and the intriguing interaction of Rayleigh and/or Marangoni convection on the basis of a phase diagram for the flow field expressed in terms of the Rayleigh and Marangoni numbers.

Motivated by the desire to understand complex transient behaviour in fluid flows, we study the dynamics of an air bubble driven by the steady motion of a suspending viscous fluid within a Hele-Shaw channel with a centred depth perturbation. Using both experiments and numerical simulations of a depth-averaged model, we investigate the evolution of an initially centred bubble of prescribed volume as a function of flow rate and initial shape. The experiments exhibit a rich variety of organised transient dynamics, involving bubble breakup as well as aggregation and coalescence of interacting neighbouring bubbles. The long-term outcome is either a single bubble or multiple separating bubbles, positioned along the channel in order of increasing velocity. Up to moderate flow rates, the life and fate of the bubble are reproducible and can be categorised by a small number of characteristic behaviours that occur in simply connected regions of the parameter plane. Increasing the flow rate leads to less reproducible time evolutions with increasing sensitivity to initial conditions and perturbations in the channel. Time-dependent numerical simulations that allow for breakup and coalescence are found to reproduce most of the dynamical behaviour observed experimentally, including enhanced sensitivity at high flow rate. An unusual feature of this system is that the set of steady and periodic solutions can change during temporal evolution because both the number of bubbles and their size distribution evolve due to breakup and coalescence events. Calculation of stable and unstable solutions in the single- and two-bubble cases reveals that the transient dynamics is orchestrated by weakly unstable solutions of the system that can appear and disappear as the number of bubbles changes.

In this work, model closures of the multiphase Reynolds-averaged Navier–Stokes (RANS) equations are developed for homogeneous, fully developed gas–particle flows. To date, the majority of RANS closures are based on extensions of single-phase turbulence models, which fail to capture complex two-phase flow dynamics across dilute and dense regimes, especially when two-way coupling between the phases is important. In the present study, particles settle under gravity in an unbounded viscous fluid. At sufficient mass loadings, interphase momentum exchange between the phases results in the spontaneous generation of particle clusters that sustain velocity fluctuations in the fluid. Data generated from Eulerian–Lagrangian simulations are used in a sparse regression method for model closure that ensures form invariance. Particular attention is paid to modelling the unclosed terms unique to the multiphase RANS equations (drag production, drag exchange, pressure strain and viscous dissipation). A minimal set of tensors is presented that serve as the basis for modelling. It is found that sparse regression identifies compact, algebraic models that are accurate across flow conditions and robust to sparse training data.

Electrohydrodynamics of drops is a classic fluid mechanical problem where deformations and microscale flows are generated by application of an external electric field. In weak fields, electric stresses acting on the drop surface drive quadrupolar flows inside and outside and cause the drop to adopt a steady axisymmetric shape. This phenomenon is best explained by the leaky-dielectric model under the premise that a net surface charge is present at the interface while the bulk fluids are electroneutral. In the case of dielectric drops, increasing the electric field beyond a critical value can cause the drop to start rotating spontaneously and assume a steady tilted shape. This symmetry-breaking phenomenon, called Quincke rotation, arises due to the action of the interfacial electric torque countering the viscous torque on the drop, giving rise to steady rotation in sufficiently strong fields. Here, we present a small-deformation theory for the electrohydrodynamics of dielectric drops for the complete Melcher–Taylor leaky-dielectric model in three dimensions. Our theory is valid in the limits of strong capillary forces and highly viscous drops and is able to capture the transition to Quincke rotation. A coupled set of nonlinear ordinary differential equations for the induced dipole moments and shape functions are derived whose solution matches well with experimental results in the appropriate small-deformation regime. Retention of both the straining and rotational components of the flow in the governing equation for charge transport enables us to perform a linear stability analysis and derive a criterion for the applied electric field strength that must be overcome for the onset of Quincke rotation of a viscous drop.

Circular milling, a stunning manifestation of collective motion, is found across the natural world, from fish shoals to army ants. It has been observed recently that the plant-animal worm Symsagittifera roscoffensis exhibits circular milling behaviour, both in shallow pools at the beach and in Petri dishes in the laboratory. Here we investigate this phenomenon experimentally and theoretically, from a fluid dynamical viewpoint, focusing on the effect that an established circular mill has on the surrounding fluid. Unlike systems such as confined bacterial suspensions and collections of molecular motors and filaments that exhibit spontaneous circulatory behaviour, and which are modelled as force dipoles, the front–back symmetry of individual worms precludes a stresslet contribution. Instead, singularities such as source dipoles and Stokes quadrupoles are expected to dominate. We analyse a series of models to understand the contributions of these singularities to the azimuthal flow fields generated by a mill, in light of the particular boundary conditions that hold for flow in a Petri dish. A model that treats a circular mill as a rigid rotating disc that generates a Stokes flow is shown to capture basic experimental results well, and gives insights into the emergence and stability of multiple mill systems.