Vortex reconnection under Biot–Savart evolution is investigated geometrically and numerically using a tent model consisting of vortex filaments initially in the form of two tilted hyperbolic branches; the vortices are antiparallel at their points of nearest approach. It is shown that the tips of these vortices approach each other, accelerating as they do so to form a finite-time singularity at the apex of the tent. The minimum separation of the vortices and the maximum velocity and axial strain rate exhibit nearly self-similar Leray scaling, but the exponents of the velocity and strain rate deviate slightly from their respective self-similar values of $-1/2$ and $-1$ ; this deviation is associated with the appearance of distinct minima of curvature leading to cusp structures at the tips. The writhe and twist of each vortex are both zero at all times up to the instant of reconnection. By way of validation of the model, the structure of the eigenvalues and eigenvectors of the rate-of-strain tensor is investigated: it is shown that the second eigenvalue $\unicode[STIX]{x1D706}_{2}$ has dipole structure around the vortex filaments. At the tips, it is observed that $\unicode[STIX]{x1D706}_{2}$ is positive and the corresponding eigenvector is tangent to the filament, implying persistent stretching of the vortex.

Fundamental solutions (Green’s functions) are derived for the regularised 13-moment system (R13) of rarefied gas dynamics, for small departures from equilibrium; these solutions show the presence of Knudsen layers, associated with exponential decay terms, that do not feature in the solution of lower-order systems (e.g. the Navier–Stokes–Fourier equations). Incorporation of these new fundamental solutions into a numerical framework based on the method of fundamental solutions (MFS) allows for efficient computation of three-dimensional gas microflows at remarkably low computational cost. The R13-MFS approach accurately recovers analytic solutions for low-speed flow around a stationary sphere and heat transfer from a hot sphere (for which a new analytic solution has been derived), capturing non-equilibrium flow phenomena missing from lower-order solutions. To demonstrate the potential of the new approach, the influence of kinetic effects on the hydrodynamic interaction between approaching solid microparticles is calculated. Finally, a programme of future work based on the initial steps taken in this article is outlined.

We present the results of a combined experimental and theoretical study of drop coalescence in the presence of an initial temperature difference $\unicode[STIX]{x0394}T_{0}$ between a drop and a bath of the same liquid. We characterize experimentally the dependence of the residence time before coalescence on $\unicode[STIX]{x0394}T_{0}$ for silicone oils with different viscosities. Delayed coalescence arises above a critical temperature difference $\unicode[STIX]{x0394}T_{c}$ that depends on the fluid viscosity: for $\unicode[STIX]{x0394}T_{0}>\unicode[STIX]{x0394}T_{c}$ , the delay time increases as $\unicode[STIX]{x0394}T_{0}^{2/3}$ for all liquids examined. This observed dependence is rationalized theoretically through consideration of the thermocapillary flows generated within the drop, the bath and the intervening air layer.

We explore the dynamics of starting plumes by analysis of a series of new small-scale laboratory experiments combined with a theoretical model for mass, momentum, and buoyancy conservation. We find that the head of the plume ascends with a speed which is approximately 0.6 times the characteristic speed of the fluid in the following steady plume, in accord with Turner (J. Fluid Mech., vol. 13 (03), 1962, pp. 356–368), and so the fluid released from the source eventually catches the head of the flow. On reaching the top of the plume it recirculates and mixes in the plume head. We estimate that approximately $0.61\pm 0.04$ of the total buoyancy released from the source accumulates in the plume head, with the remainder in the following steady plume. Using measurements of the volume of the head, we estimate that a fraction $0.16\pm 0.08$ of the volume of the head is entrained directly from the ambient, with the remainder of the fluid in the head being supplied by the following steady plume. These results imply that the buoyancy force exerted on the plume head plus the momentum flux supplied by the following plume exceeds the rate of change of momentum of the plume head even including the added mass of the plume head. We propose that the difference is associated with a drag force resulting from the displacement of ambient fluid around the plume head. Using our experimental data, we estimate that the drag coefficient $C_{d}$ has a value $4.2\pm 1.4$ , with the range in values associated with the uncertainty in our estimate of entrainment of fluid directly into the plume head. As a test, the proposed model is shown to provide a reasonable description of a starting plume rising through a stratified environment in the region below the maximum height of rise of the associated steady plume, although, above this point, the shape of the plume head changes and the model breaks down.

The flow in a decelerating turbulent round jet is investigated using direct numerical simulation. The simulations are initialised with a flow field from a statistically stationary turbulent jet. Upon stopping the inflow, a deceleration wave passes through the jet, behind which the velocity field evolves towards a new statistically unsteady self-similar state. Assumption of unsteady self-similar behaviour leads to analytical relations concerning the evolution of the centreline mean axial velocity and the shapes of the radial profiles of the velocity statistics. Consistency between these predictions and the simulation data supports the use of the assumption of self-similarity. The mean radial velocity is predicted to reverse in direction near to the jet centreline as the deceleration wave passes, contributing to an approximately threefold increase in the normalised mass entrainment rate. The shape of the mean axial velocity profile undergoes a relatively small change across the deceleration transient, and this observation provides direct evidence in support of previous models that have assumed that the mean axial velocity profile, and in some cases also the jet spreading angle, remain approximately constant within unsteady jets.

We theoretically investigate the pattern formation observed when a fluid flows over a solid substrate that can dissolve or melt. We use a turbulent mixing description that includes the effect of the bed roughness. We show that the dissolution instability at the origin of the pattern is associated with an anomaly at the transition from a laminar to a turbulent hydrodynamic response with respect to the bed elevation. This anomaly, and therefore the instability, disappears when the bed becomes hydrodynamically rough, above a threshold roughness-based Reynolds number. This suggests a possible mechanism for the selection of the pattern amplitude. The most unstable wavelength scales as observed in nature on the thickness of the viscous sublayer, with a multiplicative factor that depends on the dimensionless parameters of the problem.

We consider turbulence driven by a large-scale horizontal shear in Kolmogorov flow (i.e. with sinusoidal body forcing) and a background linear stable stratification with buoyancy frequency $N_{B}^{2}$ imposed in the third, vertical direction in a fluid with kinematic viscosity $\unicode[STIX]{x1D708}$ . This flow is known to be organised into layers by nonlinear unstable steady states, which incline the background shear in the vertical and can be demonstrated to be the finite-amplitude saturation of a sequence of instabilities, originally from the laminar state. Here, we investigate the next order of motions in this system, i.e. the time-dependent mechanisms by which the density field is irreversibly mixed. This investigation is achieved using ‘recurrent flow analysis’. We identify (unstable) periodic orbits, which are embedded in the turbulent attractor, and use these orbits as proxies for the chaotic flow. We find that the time average of an appropriate measure of the ‘mixing efficiency’ of the flow $\mathscr{E}=\unicode[STIX]{x1D712}/(\unicode[STIX]{x1D712}+{\mathcal{D}})$ (where ${\mathcal{D}}$ is the volume-averaged kinetic energy dissipation rate and $\unicode[STIX]{x1D712}$ is the volume-averaged density variance dissipation rate) varies non-monotonically with the time-averaged buoyancy Reynolds numbers $\overline{Re}_{B}=\overline{{\mathcal{D}}}/(\unicode[STIX]{x1D708}N_{B}^{2})$ , and is bounded above by $1/6$ , consistently with the classical model of Osborn (J. Phys. Oceanogr., vol. 10 (1), 1980, pp. 83–89). There are qualitatively different physical properties between the unstable orbits that have lower irreversible mixing efficiency at low $\overline{Re}_{B}\sim O(1)$ and those with nearly optimal $\mathscr{E}\lesssim 1/6$ at intermediate $\overline{Re}_{B}\sim 10$ . The weaker orbits, inevitably embedded in more strongly stratified flow, are characterised by straining or ‘scouring’ motions, while the more efficient orbits have clear overturning dynamics in more weakly stratified, and apparently shear-unstable flow.

There is a long-standing debate in the literature of stratified flows over topography concerning the correct dimensionless number to refer to as a Froude number. Common definitions using external quantities of the flow include $U/(ND)$ , $U/(Nh_{0})$ , and $Uk/N$ , where $U$ and $N$ are, respectively, scales for the background velocity and buoyancy frequency, $D$ is the depth, and $h_{0}$ and $k^{-1}$ are, respectively, height and width scales of the topography. It is also possible to define an internal Froude number $Fr_{\unicode[STIX]{x1D6FF}}=u_{0}/\sqrt{g^{\prime }\unicode[STIX]{x1D6FF}}$ , where $u_{0}$ , $g^{\prime }$ , and $\unicode[STIX]{x1D6FF}$ are, respectively, the characteristic velocity, reduced gravity, and vertical length scale of the perturbation above the topography. For the case of hydrostatic lee waves in a deep ocean, both $U/(ND)$ and $Uk/N$ are insignificantly small, rendering the dimensionless number $Nh_{0}/U$ the only relevant dynamical parameter. However, although it appears to be an inverse Froude number, such an interpretation is incorrect. By non-dimensionalizing the stratified Euler equations describing the flow of an infinitely deep fluid over topography, we show that $Nh_{0}/U$ is in fact the square of the internal Froude number because it can identically be written in terms of the inner variables, $Fr_{\unicode[STIX]{x1D6FF}}^{2}=Nh_{0}/U=u_{0}^{2}/(g^{\prime }\unicode[STIX]{x1D6FF})$ . Our scaling also identifies $Nh_{0}/U$ as the ratio of the vertical velocity scale within the lee wave to the group velocity of the lee wave, which we term the vertical Froude number, $Fr_{vert}=Nh_{0}/U=w_{0}/c_{g}$ . To encapsulate such behaviour, we suggest referring to $Nh_{0}/U$ as the lee-wave Froude number, $Fr_{lee}$ .

Most commonly used models for turbulent mixing in the ocean rely on a background stratification against which turbulence must work to stir the fluid. While this background stratification is typically well defined in idealized numerical models, it is more difficult to capture in observations. Here, a potential discrepancy in ocean mixing estimates due to the chosen calculation of the background stratification is explored using direct numerical simulation data of breaking internal waves on slopes. Two different methods for computing the buoyancy frequency $N$ , one based on a three-dimensionally sorted density field (often used in numerical models) and the other based on locally sorted vertical density profiles (often used in the field), are used to quantify the effect of $N$ on turbulence quantities. It is shown that how $N$ is calculated changes not only the flux Richardson number $R_{f}$ , which is often used to parameterize turbulent mixing, but also the turbulence activity number or the Gibson number $Gi$ , leading to potential errors in estimates of the mixing efficiency using $Gi$ -based parameterizations.

The response of the free liquid surface in a Hele-Shaw cell subjected to a horizontal oscillation is investigated. We study the low-oscillation-amplitude regime and we show, by varying the fluid viscosity, $\unicode[STIX]{x1D708}$ , and the forcing frequency, $\unicode[STIX]{x1D714}$ , that the ratio between the Stokes viscous length, $\sqrt{2\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714}}$ , and the cell thickness greatly affects the amplitude and phase lag of the gravity waves. In particular, the sloshing system undergoes an underdamped/overdamped transition for sufficiently large viscosities. A consistent theoretical model, based on a modification of Darcy’s law to include unsteadiness, is then introduced to rationalize the experimental observations. Contrary to traditional sloshing wave theory, the viscous flow dissipation comes at leading order in the analysis, rather than as a higher-order asymptotic correction to the inviscid sloshing dynamics. The analytical expression for the resonance curves agrees well with experimental results without tunable parameters.

Accurate prediction of the spatial evolution of turbulent wake flows under pressure gradient conditions is required in some engineering applications such as the design of high-lift devices and wind farms over topography. In this paper, we aim to develop an analytical model to predict the evolution of a turbulent planar wake under an arbitrary pressure gradient condition. The model is based on the cross-stream integration of the streamwise momentum equation and uses the self-similarity of the mean flow. We have also made an experimentally supported assumption that the ratio of the maximum velocity deficit to the wake width is independent of the imposed pressure gradient. The asymptotic response of the wake to the pressure gradient is also investigated. After its derivation, the model is successfully validated against experimental data by comparing the evolution of the wake width and maximum velocity deficit. The inputs of the model are the imposed pressure gradient and the wake width under zero pressure gradient. The model does not require any parameter tuning and is deemed to be practical, computationally fast, accurate enough, and therefore useful for the scientific and engineering communities.

This experimental study focuses on the effect of horizontal boundaries with pyramid-shaped roughness elements on the heat transfer in rotating Rayleigh–Bénard convection. It is shown that the Ekman pumping mechanism, which is responsible for the heat transfer enhancement under rotation in the case of smooth top and bottom surfaces, is unaffected by the roughness as long as the Ekman layer thickness $\unicode[STIX]{x1D6FF}_{E}$ is significantly larger than the roughness height  $k$ . As the rotation rate increases, and thus $\unicode[STIX]{x1D6FF}_{E}$ decreases, the roughness elements penetrate the radially inward flow in the interior of the Ekman boundary layer that feeds the columnar Ekman vortices. This perturbation generates additional thermal disturbances which are found to increase the heat transfer efficiency even further. However, when $\unicode[STIX]{x1D6FF}_{E}\approx k$ , the Ekman boundary layer is strongly perturbed by the roughness elements and the Ekman pumping mechanism is suppressed. The results suggest that the Ekman pumping is re-established for $\unicode[STIX]{x1D6FF}_{E}\ll k$ as the faces of the pyramidal roughness elements then act locally as a sloping boundary on which an Ekman layer can be formed.

Wavepackets obtained as solutions of the flow equations linearised around the mean flow have been shown in recent work to yield good agreement, in terms of amplitude and phase, with those educed from turbulent jets. Compelling agreement has been demonstrated, for the axisymmetric and first helical mode, up to Strouhal numbers close to unity. We here extend the range of validity of wavepacket models to Strouhal number $St=4.0$ and azimuthal wavenumber $m=4$ by comparing solutions of the parabolised stability equations with a well-validated large-eddy simulation of a Mach 0.9 turbulent jet. The results show that the near-nozzle dynamics can be correctly described by the homogeneous linear model, the initial growth rates being accurately predicted for the entire range of frequencies and azimuthal wavenumbers considered. Similarly to the lower-frequency wavepackets reported prior to this work, the high-frequency linear waves deviate from the data downstream of their stabilisation locations, which move progressively upstream as the frequency increases.

Wetting of a planar solid substrate is investigated in the presence of a macroscopic particle in the complete wetting regime. A drop of silicone oil spreads on the substrate and its macroscopic edge is incident on the particle at the late stage of spreading. The drop–particle interaction is observed in detail by shadowgraph and interferometry. Although the spreading drop edge is pinned by the particle for a short time, a sharp acceleration occurs when the liquid starts wetting the extra surface area offered by the particle and forming a meniscus. This process yields a net gain in spreading speed. A theoretical model based on the classical wetting dynamics dictated by Cox’s law is developed. It predicts that the capillary energy of the meniscus gives rise to a rapid motion of the liquid edge, showing good agreement with the dynamics observed in the experiments.

We present detailed experiments on transient growth of turbulent spots induced by external forcing in plane Couette–Poiseuille flow, which are studied in the framework of linear transient growth. The experimental investigation is supplemented with full theoretical analysis. We compare quantitatively the experimental and theoretical results, including maximal gain and the time at which it occurs. We also present the limits of validity for the application of the linear theory at high amplitude perturbation and Reynolds number, showing experiments with self-sustained states.

A generalized reciprocal theorem is formulated for the motion and hydrodynamic force moments of an active particle in an arbitrary background flow of a (weakly nonlinear) complex fluid. This formalism includes as special cases a number of previous calculations of the motion of both passive and active particles in Newtonian and non-Newtonian fluids.

Turbulent wakes past bluff bodies commonly present asymmetric flow states reminiscent of bifurcations in the laminar regime. Understanding the sensitivity of these states to flow forcing is crucial to the modelling and control of flow symmetry properties. In this study, the near wake of a rectangular bluff body in proximity to a wall is disturbed by the use of passive devices located between the model and the wall, upstream of the massive flow separation occurring at the blunt trailing edges. Due to the proximity to the boundary, the wake initially presents wall-normal asymmetry and a negative wall-normal pressure gradient along the base. The application of disturbances with variable size, however, sets flow symmetry along the wall-normal plane, leading to the intermittent spanwise wake reversals reported recently in the literature. A further increase in the size of perturbation suppresses wake switching, and wall-normal asymmetry is recovered, but with a positive wall-normal pressure gradient. The dynamical features of this bifurcation scenario can be retrieved using two coupled symmetry-breaking models for spanwise and wall-normal pressure gradients. This confirms the high sensitivity of the separated flow to external perturbations. More importantly, the results unify observations of the bluff-body wake topologies covered in previous investigations.

The nonlinear stability of viscous, immiscible multilayer flows in plane channels driven both by a pressure gradient and gravity is studied. Three fluid phases are present with two interfaces. Weakly nonlinear models of coupled evolution equations for the interfacial positions are derived and studied for inertialess, stably stratified flows in channels at small inclination angles. Interfacial tension is demoted and high-wavenumber stabilisation enters due to density stratification through second-order dissipation terms rather than the fourth-order ones found for strong interfacial tension. An asymptotic analysis is carried out to demonstrate how these models arise. The governing equations are $2\times 2$ systems of second-order semi-linear parabolic partial differential equations (PDEs) that can exhibit inertialess instabilities due to interaction between the interfaces. Mathematically this takes place due to a transition of the nonlinear flux function from hyperbolic to elliptic behaviour. The concept of hyperbolic invariant regions, found in nonlinear parabolic systems, is used to analyse this inertialess mechanism and to derive a transition criterion to predict the large-time nonlinear state of the system. The criterion is shown to predict nonlinear stability or instability of flows that are stable initially, i.e. the initial nonlinear fluxes are hyperbolic. Stability requires the hyperbolicity to persist at large times, whereas instability sets in when ellipticity is encountered as the system evolves. In the former case the solution decays asymptotically to its uniform base state, while in the latter case nonlinear travelling waves can emerge that could not be predicted by a linear stability analysis. The nonlinear analysis predicts threshold initial disturbances above which instability emerges.

A one-dimensional, unsteady nozzle flow is modelled to identify the sources of indirect noise in multicomponent gases. First, from non-equilibrium thermodynamics relations, it is shown that a compositional inhomogeneity advected in an accelerating flow is a source of sound induced by inhomogeneities in the mixture (i) chemical potentials and (ii) specific heat capacities. Second, it is shown that the acoustic, entropy and compositional linear perturbations evolve independently from each other and they become coupled through mean-flow gradients and/or at the boundaries. Third, the equations are cast in invariant formulation and a mathematical solution is found by asymptotic expansion of path-ordered integrals with an infinite radius of convergence. Finally, the transfer functions are calculated for a supersonic nozzle with finite spatial extent perturbed by a methane–air compositional inhomogeneity. The proposed framework will help identify and quantify the sources of sound in nozzles with relevance, for example, to aeronautical gas turbines.

In low-density axisymmetric jets, the onset of global instability is known to depend on three control parameters, namely the jet-to-ambient density ratio $S$ , the initial momentum thickness $\unicode[STIX]{x1D703}_{0}$ and the Reynolds number $Re$ . For sufficiently low values of $S$ and $\unicode[STIX]{x1D703}_{0}$ , these jets bifurcate from a steady state (a fixed point) to a self-excited oscillatory state (a limit cycle) when $Re$ increases above a critical value corresponding to the Hopf point, $Re_{H}$ . In the literature, this Hopf bifurcation is often regarded as supercritical. In this experimental study, however, we find that under some conditions, there exists a hysteretic bistable region at $Re_{SN}<Re<Re_{H}$ , where $Re_{SN}$ denotes a saddle-node point. This shows that, contrary to expectations, the Hopf bifurcation can also be subcritical, which we explore by evaluating the coefficients of a truncated Landau model. The existence of subcritical bifurcations implies the potential for triggering and the need for weakly nonlinear analyses to be performed to at least fifth order if one is to be able to predict saturation and bistability. We conclude by proposing a universal scaling for $Re_{H}$ in terms of $S$ and $\unicode[STIX]{x1D703}_{0}$ . This scaling, which is insensitive to the super/subcritical nature of the bifurcations, can be used to predict the onset of self-excited oscillations, providing further evidence to support Hallberg & Strykowski’s concept (J. Fluid Mech., vol. 569, 2006, pp. 493–507) of universal global modes in low-density jets.