We present a theoretical stability analysis for an expanding accretion shock that does not involve a rarefaction wave behind it. The dispersion equation that determines the eigenvalues of the problem and the explicit formulae for the corresponding eigenfunction profiles are presented for an arbitrary equation of state and finite-strength shocks. For spherically and cylindrically expanding steady shock waves, we demonstrate the possibility of instability in a literal sense, a power-law growth of shock-front perturbations with time, in the range of $h_c< h<1+2 {\mathcal {M}}_2$, where $h$ is the D'yakov-Kontorovich parameter, $h_c$ is its critical value corresponding to the onset of the instability and ${\mathcal {M}}_2$ is the downstream Mach number. Shock divergence is a stabilizing factor and, therefore, instability is found for high angular mode numbers. As the parameter $h$ increases from $h_c$ to $1+2 {\mathcal {M}}_2$, the instability power index grows from zero to infinity. This result contrasts with the classic theory applicable to planar isolated shocks, which predicts spontaneous acoustic emission associated with constant-amplitude oscillations of the perturbed shock in the range $h_c< h<1+2 {\mathcal {M}}_2$. Examples are given for three different equations of state: ideal gas, van der Waals gas and three-terms constitutive equation for simple metals.

We performed a linear stability analysis of a Newtonian ferrofluid cylinder surrounded by a Newtonian non-magnetic fluid in an azimuthal magnetic field. A wire is used at the centre of the ferrofluid cylinder to create this magnetic field. Isothermal conditions are considered and gravity is ignored. An axisymmetric perturbation is imposed at the interface between the two fluids and a dispersion relation is obtained allowing us to predict whether the flow is stable or unstable with respect to this perturbation. This relation is dependent on the Ohnesorge number of the ferrofluid, the dynamic viscosity ratio, the density ratio, the magnetic Bond number, the relative magnetic permeability and the dimensionless wire radius. Solutions to this dispersion relation are compared with experimental data from Arkhipenko et al. (Fluid Dyn., vol. 15, issue 4, 1981, pp. 477–481) and, more recently, Bourdin et al. (Phys. Rev. Lett., vol. 104, issue 9, 2010, 094502). A better agreement than the inviscid theory and the theory that only takes into account the viscosity of the ferrofluid is shown with the data of Arkhipenko et al. (Fluid Dyn., vol. 15, issue 4, 1981, pp. 477–481) and those of Bourdin et al. (Phys. Rev. Lett., vol. 104, issue 9, 2010, 094502) for small wavenumbers.

We study long surface and internal ring waves propagating in a stratified fluid over a parallel shear current. The far-field modal and amplitude equations for the ring waves are presented in dimensional form. We re-derive the modal equations from the formulation for plane waves tangent to the ring wave, which opens a way to obtaining important characteristics of the ring waves (group speed, wave action conservation law) and to constructing more general ‘hybrid solutions’ consisting of a part of a ring wave and two tangent plane waves. The modal equations constitute a new spectral problem, and are analysed for a number of examples of surface ring waves in a homogeneous fluid and internal ring waves in a stratified fluid. Detailed analysis is developed for the case of a two-layered fluid with a linear shear current where we study their wavefronts and two-dimensional modal structure. Comparisons are made between the modal functions (i.e. eigenfunctions of the relevant spectral problems) for the surface waves in homogeneous and two-layered fluids, as well as the interfacial waves described exactly and in the rigid-lid approximation. We also analyse the wavefronts of surface and interfacial waves for a large family of power-law upper-layer currents, which can be used to model wind generated currents, river inflows and exchange flows in straits. Global and local measures of the deformation of wavefronts are introduced and evaluated.

A coastal eddy is modelled as a barotropic vortex propagating along a coastal shelf. If the vortex speed matches the phase speed of any coastal trapped shelf wave modes, a shelf wave wake is generated leading to a flux of energy from the vortex into the wave field. Using a simple shelf geometry, we determine analytic expressions for the wave wake and the leading-order flux of wave energy. By considering the balance of energy between the vortex and wave field, this energy flux is then used to make analytic predictions for the evolution of the vortex speed and radius under the assumption that the vortex structure remains self-similar. These predictions are examined in the asymptotic limit of small rotation rate and shelf slope and tested against numerical simulations. If the vortex speed does not match the phase speed of any shelf wave, steady vortex solutions are expected to exist. We present a numerical approach for finding these nonlinear solutions and examine the parameter dependence of their structure.

It is well established that Lagrangian particle dispersion models, for inhomogeneous turbulent flows, must satisfy the ‘well-mixed condition’ of Thomson (J. Fluid Mech., vol. 180, 1987, pp. 529–556) in order to produce physically reasonable results. In more than one dimension, however, the well-mixed condition is not sufficient to define the dispersion model uniquely. The non-uniqueness, which is related to the rotational degrees of freedom of particle trajectories, permits models with trajectory curvatures and velocity autocorrelation functions which are clearly unphysical. A spin condition is therefore introduced to constrain the models. It requires an ensemble of particles with fixed initial position and velocity to have, at short times, expected angular momentum, measured relative to the mean position and velocity of an ensemble of fluid particles with initially random velocity, equal to the relative angular momentum of the mean flow at the ensemble mean location. The resulting unique model is found explicitly for the canonical example of inhomogeneous Gaussian turbulence and is characterised by accelerations which are exponential in the particle velocity. A simpler unique model with a quadratic acceleration is obtained using a weaker version of the spin condition. Unlike previous models, the unique models defined by the spin condition lead to particles having the correct (ensemble mean) angular speed in a turbulent flow in solid-body rotation. The properties of the new models are discussed in the settings of a turbulent channel flow and an idealised turbulent atmospheric boundary-layer flow.

The limit of power extraction by a device which makes use of constructive interference, i.e. local blockage, is investigated theoretically. The device is modelled using actuator disc theory in which we allow the device to be split into arrays and these then into sub-arrays an arbitrary number of times so as to construct an $n$-level multi-scale device in which the original device undergoes $n-1$ sub-divisions. The alternative physical interpretation of the problem is a planar system of arrayed turbines in which groups of turbines are homogeneously arrayed at the smallest $n\mathrm {th}$ scale, and then these groups are homogeneously spaced relative to each other at the next smallest $n-1\mathrm {th}$ scale, with this pattern repeating at all subsequent larger scales. The scale-separation idea of Nishino & Willden (J. Fluid. Mech., vol. 708, 2012b, pp. 596–606) is employed, which assumes mixing within a sub-array occurs faster than mixing of the by-pass flow around that sub-array, so that in the $n$-scale device mixing occurs from the inner scale to the outermost scale in that order. We investigate the behaviour of an arbitrary level multi-scale device, and determine the arrangement of actuator discs ($n\mathrm {th}$ level devices) which maximises the power coefficient (ratio of power extracted to undisturbed kinetic energy flux through the net disc frontal area). We find that this optimal arrangement is close to fractal, and fractal arrangements give similar results. With the device placed in an infinitely wide channel, i.e. zero global blockage, we find that the optimum power coefficient tends to unity as the number of device scales tends to infinity, a 27/16 increase over the Lanchester–Betz limit of $0.593$. For devices in finite width channels, i.e. non-zero global blockage, similar observations can be made with further uplift in the maximum power coefficient. We discuss the fluid mechanics of this energy extraction process and examine the scale distribution of thrust and wake velocity coefficients. Numerical demonstration of performance uplift due to multi-scale dynamics is also provided. We demonstrate that bypass flow remixing and ensuing energy losses increase the device power coefficient above the limits for single devices, so that although the power coefficient can be made to increase, this is at the expense of the overall efficiency of energy extraction which decreases as wake-scale remixing losses necessarily rise. For multi-scale devices in finite overall blockage two effects act to increase extractable power; an overall streamwise pressure gradient associated with finite blockage, and wake pressure recoveries associated with bypass-scale remixing.

The slow motion of a small buoyant sphere near a right dihedral corner made by tangentially sliding walls is investigated. Under creeping-flow conditions the force and torque on the sphere can be decomposed into eleven elementary types of motion involving simple particle translations, particle rotations and wall movements. Force and torque balances are employed to find the velocity and rotation of the particle as functions of its location. Depending on the ratio of the wall velocities and the gravitational settling velocity of the sphere, different dynamical regimes are identified. In particular, a non-trivial line attractor/repeller for the particle motion exists at a location detached from both the walls. The existence, location and stability of the corresponding two-dimensional fixed point are studied depending on the wall velocities and the buoyancy force. The impact of the line attractors/repellers on the motion of small particles in cavities and its relevance for corner cleaning applications are discussed.

Small viscous effects in high-Reynolds-number rotational flows always accumulate over time to have a leading-order effect. Therefore, the high-Reynolds-number limit for the Navier–Stokes equations is singular. It is important to investigate whether a solution of the Euler equations can approximate a real flow at large Reynolds number. These facts are often overlooked and, as a result, the Euler equations are used to simulate laminar rotational flows at large Reynolds number. Based on the Fredholm alternative, an asymptotic perturbation theory is described to establish secularity conditions determined by viscosity for an inviscid solution to approximate a real viscous fluid. Four important classical inviscid solutions are investigated using the theory with the following conclusions. The Stuart cats’ eyes and Mallier–Maslowe vortices are inconsistent with any real fluid at high Reynolds number; whereas Hill's spherical vortex is confirmed to be consistent with a steady state in the spherical core region and the Lamb–Chaplygin dipole is found to be consistent with a quasi-steady state in the circular core region. These solutions have been widely used for analysing the stability of vortex flows and wakes, and their interactions with shock waves or bubbles. Serendipitously, we have revealed an original exact solution of the Navier–Stokes equations which is time dependent, has non-zero nonlinear convective terms and is restricted to a finite domain with the decay rate depending on dipole radius.

Motivated by the industrial manufacture of organic light-emitting-diode displays, we formulate and analyse a mathematical model for the evolution of a thin droplet in a shallow axisymmetric well of rather general shape both before and after touchdown that accounts for the spatially non-uniform evaporation of the fluid, perform physical experiments using three cylindrical wells with different small aspect ratios, and validate the mathematical model by comparing the present experimental results with the corresponding theoretical predictions for a cylindrical well.

This study examines the impact of interfacial viscosity on the stability of an initially deformed droplet translating through an unbounded quiescent fluid. The boundary-integral formulation is employed to investigate the time evolution of a droplet in the Stokes flow limit. The droplet interface is modelled using the Boussinesq–Scriven constitutive relationship having surface shear viscosity $\eta _\mu$ and surface dilatational viscosity $\eta _\kappa$. We observe that, below a critical value of the capillary number, $Ca_C$, the initially perturbed droplet reverts to its spherical shape. Above $Ca_C$, the translating droplet deforms continuously, growing a tail at the rear end for initial prolate perturbations and a cavity for initial oblate perturbations. We find that surface shear viscosity inhibits the tail/cavity growth at the droplet's rear end and increases the $Ca_C$ compared with a clean droplet. In contrast, surface dilatational viscosity increases tail/cavity growth and lowers $Ca_C$ compared with a clean droplet. Surprisingly, both shear and dilatational surface viscosity appear to delay the time at which pinch off occurs, and hence satellite droplets form. Lastly, we explore the combined influence of surface viscosity and surfactant transport on droplet stability by assuming a linear dependence of surface tension on surfactant concentration and exponential dependence of interfacial viscosities on the surface pressure. We find that pressure-thinning/thickening effects significantly affect the droplet dynamics for surface shear viscosity but play a small role for surface dilatational viscosity. We lastly provide phase diagrams for the critical capillary number for different values of the droplet's viscosity ratio and initial Taylor deformation parameter.

Direct numerical simulations are carried out to study the effect of finite Weissenberg number up to $Wi=16$ on laminar and turbulent channel flows of an elastoviscoplastic (EVP) fluid, at a fixed bulk Reynolds number of $2800$. The incompressible flow equations are coupled with the evolution equation for the EVP stress tensor by a modified Saramito model that extends both the Bingham viscoplastic and the finite extensible nonlinear elastic-Peterlin (FENE-P) viscoelastic models. In turbulent flow, we find that drag decreases with both the Bingham and Weissenberg numbers, until the flow laminarises at high enough elastic and yield stresses. Hence, a higher drag reduction is achieved than in the viscoelastic flow at the same Weissenberg number. The drag reduction persists at Bingham numbers up to 20, in contrast to viscoplastic flow, where the drag increases in the laminar regime compared with a Newtonian flow. Moreover, elasticity affects the laminarisation of an EVP flow in a non-monotonic fashion, delaying it at lower and promoting it at higher Weissenberg numbers. A hibernation phenomenon is observed in the EVP flow, leading to large changes in the unyielded regions. Finally, plasticity is observed to affect both low- and high-speed streaks equally, attenuating the turbulent dissipation and the fragmentation of turbulent structures.

Chemically active colloids self-propel by catalysing the decomposition of molecular ‘fuel’ available in the surrounding solution. If the various molecular species involved in the reaction have distinct interactions with the colloid surface, and if the colloid has some intrinsic asymmetry in its surface chemistry or geometry, there will be phoretic flows in an interfacial layer surrounding the particle, leading to directed motion. Most studies of chemically active colloids have focused on spherical, axisymmetric ‘Janus’ particles, which (in the bulk, and in absence of fluctuations) simply move in a straight line. For particles with a complex (non-spherical and non-axisymmetric) geometry, the dynamics can be much richer. Here, we consider chemically active helices. Via numerical calculations and slender body theory, we study how the translational and rotational velocities of the particle depend on geometry and the distribution of catalytic activity over the particle surface. We confirm the recent finding of Katsamba et al. (J. Fluid Mech., vol. 898, 2020, p. A24) that both tangential and circumferential concentration gradients contribute to the particle velocity. The relative importance of these contributions has a strong impact on the motion of the particle. We show that, by a judicious choice of the particle design parameters, one can suppress components of angular velocity that are perpendicular to the screw axis, or even select for purely ‘sideways’ translation of the helix.