Violent respiratory events such as coughs and sneezes play a key role in transferring respiratory diseases between infectious and susceptible individuals. We present the results of a combined experimental and theoretical investigation of the fluid dynamics of such violent expiratory events. Direct observation of sneezing and coughing events reveals that such flows are multiphase turbulent buoyant clouds with suspended droplets of various sizes. Our observations guide the development of an accompanying theoretical model of pathogen-bearing droplets interacting with a turbulent buoyant momentum puff. We develop in turn discrete and continuous models of droplet fallout from the cloud in order to predict the range of pathogens. According to the discrete fallout model droplets remain suspended in the cloud until their settling speed matches that of the decelerating cloud. A continuous fallout model is developed by adapting models of sedimentation from turbulent fluids. The predictions of our theoretical models are tested against data gathered from a series of analogue experiments in which a particle-laden cloud is ejected into a relatively dense ambient. Our study highlights the importance of the multiphase nature of respiratory clouds, specifically the suspension of the smallest drops by circulation within the cloud, in extending the range of respiratory pathogens.

Direct numerical simulation (DNS) of liquid film flow is used to compute fully developed solitary waves and to compare their characteristics with the predictions of low-dimensional models. Emphasis is placed on the regime of high inertia, where available models provide widely differing results. It is found that the parametric dependence of wave properties on inertia is highly non-trivial, and is satisfactorily approximated only by the four-equation model of Ruyer-Quil & Manneville (Eur. Phys. J. B, vol. 15, 2000, pp. 357–369). Detailed comparison of the asymptotic shapes of upstream and downstream tails is performed, and inherent limitations of all long-wave models are revealed. Local flow reversal in front of the main hump, which has been previously discussed in the literature, is shown to occur for an inertia range bounded from below and from above, and the boundaries are interpreted in terms of the capillary origin of the phenomenon. Computational results are reported for the entire range of Froude numbers, providing benchmark data for all wall inclinations.

We present a theoretical and experimental study of viscous flows injected into a porous medium that is confined vertically by horizontal impermeable boundaries and filled with an ambient fluid of different density and viscosity. General three-dimensional equations describing such flows are developed, showing that the dynamics can be affected by two separate contributions: spreading due to gradients in hydrostatic pressure, and that due to the pressure drop introduced by the injection. In the illustrative case of a two-dimensional injection of fluid at a constant volumetric rate, the injected fluid initially forms a viscous gravity current insensitive both to the depth of the medium and to the viscosity of the ambient fluid. Beyond a characteristic time scale, the dynamics transition to being dominated by the injection pressure, and the injected fluid eventually intersects the second boundary to form a second moving contact line. Three different late-time asymptotic regimes can emerge, depending on whether the viscosity of the injected fluid is less than, equal to or greater than that of the ambient fluid. With a less viscous injection, the flow undergoes a slow decay towards a similarity solution in which the two contact lines extend linearly in time with differing prefactors. Perturbations from this long-term state are shown to decay algebraically with time. Equal viscosities result in both contact lines approaching the same leading-order asymptotic position but with a first-order correction to the distance between them that expands as $t^{1/2}$ due to gravitational spreading. For a more viscous injection, the distance between the contact lines approaches a constant value, with perturbations decaying exponentially. Data from a new series of laboratory experiments confirm these theoretical predictions.

We present a fluid-mechanical model of the coalescence of a number of elastic objects due to surface tension. We consider an array of spring–block elements separated by thin liquid films, whose dynamics are modelled using lubrication theory. With this simplified model of elastocapillary coalescence, we present the results of numerical simulations for a large number of elements, $N=O(10^4)$. A linear stability analysis shows that pairwise coalescence is always the most unstable mode of deformation. However, the numerical simulations show that the cluster sizes actually produced by coalescence from a small white-noise perturbation have a distribution that depends on the relative strength of surface tension and elasticity, as measured by an elastocapillary number $K$. Both the maximum cluster size and the mean cluster size scale like $K^{-1/2}$ for small $K$. An analytical solution for the response of the system to a localized perturbation shows that such perturbations generate propagating disturbance fronts, which leave behind ‘frozen-in’ clusters of a predictable size that also depends on $K$. A good quantitative comparison between the cluster-size statistics from noisy perturbations and this ‘frozen-in’ cluster size suggests that propagating fronts may play a crucial role in the dynamics of coalescence.

Local linear stability analysis has been shown to provide valuable information about the response of jet diffusion flames to flow-field perturbations. However, this analysis commonly relies on several modelling assumptions about the mean flow prescription, the thermo-viscous-diffusive transport properties, and the complexity and representation of the chemical reaction mechanisms. In this work, the effects of these modelling assumptions on the stability behaviour of a jet diffusion flame are systematically investigated. A flamelet formulation is combined with linear stability theory to fully account for the effects of complex transport properties and the detailed reaction chemistry on the perturbation dynamics. The model is applied to a methane–air jet diffusion flame that was experimentally investigated by Füri et al. (Proc. Combust. Inst., vol. 29, 2002, pp. 1653–1661). Detailed simulations are performed to obtain mean flow quantities, about which the stability analysis is performed. Simulation results show that the growth rate of the inviscid instability mode is insensitive to the representation of the transport properties at low frequencies, and exhibits a stronger dependence on the mean flow representation. The effects of the complexity of the reaction chemistry on the stability behaviour are investigated in the context of an adiabatic jet flame configuration. Comparisons with a detailed chemical-kinetics model show that the use of a one-step chemistry representation in combination with a simplified viscous-diffusive transport model can affect the mean flow representation and heat release location, thereby modifying the instability behaviour. This is attributed to the shift in the flame structure predicted by the one-step chemistry model, and is further exacerbated by the representation of the transport properties. A pinch-point analysis is performed to investigate the stability behaviour; it is shown that the shear-layer instability is convectively unstable, while the outer buoyancy-driven instability mode transitions from absolutely to convectively unstable in the nozzle near field, and this transition point is dependent on the Froude number.

The gravity-driven flow of a viscous liquid film coating the inside of a tube is studied both theoretically and experimentally. As the film moves downward, small perturbations to the free surface grow due to surface tension effects and can form liquid plugs. A first-principles strongly nonlinear model based on long-wave asymptotics is developed to provide simplified governing equations for the motion of the film flow. Linear stability analysis on the basic solution of the model predicts the speed and wavelength of the most unstable mode, and whether the film is convectively or absolutely unstable. These results are found to be in remarkable agreement with the experiments. The model is also solved numerically to follow the time evolution of instabilities. For relatively thin films, these instabilities saturate as a series of small-amplitude travelling waves, while thicker films lead to solutions whose amplitude becomes large enough for the liquid surface to approach the centre of the tube in finite time, suggesting liquid plug formation. Next, the model’s periodic travelling wave solutions are determined by a continuation algorithm using the results from the time evolution code as initial seed. It is found that bifurcation branches for these solutions exist, and the critical turning points where branches merge determine film mean thicknesses beyond which no travelling wave solutions exist. These critical thickness values are in good agreement with those for liquid plug formations determined experimentally and numerically by the time-evolution code.

It is shown here that a subset of the implicit analytical shock solutions discovered by Becker and by Johnson can be inverted, yielding several exact closed-form solutions of the one-dimensional compressible Navier–Stokes equations for an ideal gas. For a constant dynamic viscosity and thermal conductivity, and at particular values of the shock Mach number, the velocity can be expressed in terms of a polynomial root. For a constant kinematic viscosity, independent of Mach number, the velocity can be expressed in terms of a hyperbolic tangent function. The remaining fluid variables are related to the velocity through simple algebraic expressions. The solutions derived here make excellent verification tests for numerical algorithms, since no source terms in the evolution equations are approximated, and the closed-form expressions are straightforward to implement. The solutions are also of some academic interest as they may provide insight into the nonlinear character of the Navier–Stokes equations and may stimulate further analytical developments.

The structure of a steady axisymmetric thermal plume rising through a very viscous fluid with strongly temperature-dependent viscosity of the form $\mu \propto \exp (-\gamma T)$ is investigated. An analytic asymptotic solution is derived for the fast-flowing core of the plume, which predicts that the excess centreline temperature decays exponentially as $\exp \{ - 12 \pi \kappa z/(\gamma A) \}$, where $\kappa $ is the thermal diffusivity, $z$ the height and $A$ the vertical heat flux. This rate of decay, which is found to be in good agreement with numerical simulations of the boundary-layer equations, is three times faster than that predicted by the oft-quoted model of Olson, Schubert and Anderson (J. Geophys. Res., vol. 98 (B4), 1993, pp. 6829–6844).