The motion of finite-length cylindrical rods moving near a planar rigid surface is a scenario common across many engineering and natural settings. We study the low-Reynolds-number flow around finite rods that are allowed to rotate or translate in directions perpendicular or parallel to the plane. We develop a three-dimensional lubrication theory to characterize the pressure and hydrodynamic resistances of the cylinders through a special consideration of the cylinder's end effects. In addition, we use three-dimensional numerical simulations to solve these Stokes flows for cylinders of varying lengths and with varying gap sizes between the cylinder and plane, and the numerical results support the developed analytical descriptions. We also use visualizations of the flow to provide qualitative insights and rationalize the effect of the ends on the dynamics of the cylinders. The numerical simulations and theoretical predictions show good agreement in the long (isolated ends) and short (disk-like) limits.

Gravity induced film flow over a rigid smoothly corrugated substrate heated uniformly from below, is explored. This is achieved by reducing the governing equations of motion and energy conservation to a manageable form within the mathematical framework of the well-known long-wave approximation; leading to an asymptotic model of reduced dimensionality. A key feature of the approach and to solving the problem of interest, is proof that the leading approximation of the temperature field inside the film must be nonlinear to accurately resolve the thermodynamics beyond the trivial case of ‘a flat film flowing down a planar uniformly heated incline.’ Superior predictions are obtained compared with earlier work and reinforced via a series of corresponding solutions to the full governing equations using a purpose written finite element analogue, enabling comparisons to be made between free-surface disturbance and temperature predictions, as well as the streamline pattern and temperature contours inside the film. In particular, the free-surface temperature is captured extremely well at moderate Prandtl numbers. The stability characteristics of the problem are examined using Floquet theory, with the interaction between the substrate topography and thermo-capillary instability modes investigated as a set of neutral stability curves. Although there are no relevant experimental data currently available for the heated film problem, recent existing predictions and experimental data concerning the behaviour of corresponding isothermal flow cases are taken as a reference point from which to explore the effect of both heating and cooling.

A one-dimensional mechanism of deflagration to detonation transition is identified and investigated by an asymptotic analysis in the double limit of large activation energy and small Mach number of the laminar flame velocity. The unsteady analysis concerns the self-accelerating tip of an elongated flame in a smooth walled tube. The flame on the tip, considered as plane and orthogonal to the tube axis, is pushed from behind by the longitudinal flow resulting from the cumulative effect of the radial flows of burned gas issued from the lateral flame of the finger-like front (called backflow in the following). The analysis of the one-dimensional dynamics is performed by coupling the flame structure with the downstream-running compression waves propagating in the external flows. A critical elongation is identified from which the slightest increase in elongation leads to a pressure runaway producing the flame blow-off. The dynamics of the inner structure of the laminar flame on the tip which is accelerated by the self-induced backflow is characterized by a finite-time singularity of the reacting flow in the form of a dynamical saddle-node bifurcation.

When a lightning bolt darts across the sky, the thunderclap that reaches our ears a few seconds later is an example of a fluid dynamical shock: a wave across which flow properties such as pressure and density change almost discontinuously. In compressible fluids these shocks are associated with high-energy supersonic flows and so require specialist equipment to realise in steady state. But in granular media, shocks occur much more readily and at flow speeds easily obtainable in the laboratory. In the featured article, Khan et al. (J. Fluid Mech., vol. 884, 2020, R4) exploit this to explore a remarkable range of steady and oscillatory shocks and shock interactions, which demonstrate many of the unique rheological complexities of granular flow.

Geophysical turbulent flows, characterized by rapid rotation, quantified by small Rossby number, and stable stratification, often self-organize into a collection of coherent vortices, referred to as a vortex gas. The lowest-order asymptotic expansion in Rossby number is quasigeostrophy, which has purely horizontal velocities and cyclone–anticyclone antisymmetry. Ageostrophic effects are important components of many geophysical flows and, as such, these phenomena are not well modelled by quasigeostrophy. The next-order correction in Rossby number, which includes ageostrophic effects, is the so-called balanced dynamics. Balanced dynamics includes ageostrophic vertical velocity and breaks the geostrophic cyclone–anticyclone antisymmetry. Point-vortex solutions are well known in two-dimensional and quasigeostrophic dynamics and are useful for studying the vortex-gas regime of geophysical turbulence. Here, we find point-vortex solutions in fully three-dimensional continuously stratified QG$^{+1}$ dynamics, a particular formulation of balanced dynamics. Simulations of QG$^{+1}$ point vortices show several interesting features not captured by quasigeostrophic point vortices including significant vertical transport on long time scales. The ageostrophic component of QG$^{+1}$ point vortex point-vortex dynamics renders them useful in modelling flows where quasigeostrophy filters out important physical processes.

We analyse a class of stochastic advection problems by conditionally averaging the passive tracer equation with respect to a given flow state. In doing so, we obtain expressions for the turbulent diffusivity as a function of the flow statistics spectrum. When flow statistics are given by a continuous-time Markov process with a finite state space, calculations are amenable to analytic treatment. When the flow statistics are more complex, we show how to approximate turbulent fluxes as hierarchies of finite state space continuous-time Markov processes. The ensemble average turbulent flux is expressed as a linear operator that acts on the ensemble average of the tracer. We recover the classical estimate of turbulent flux as a diffusivity tensor, the components of which are the integrated autocorrelation of the velocity field in the limit that the operator becomes local in space and time.

Insights gained from modal analysis are invoked for predictive large-eddy simulation (LES) wall modelling. Specifically, we augment the law of the wall (LoW) by an additional mode based on a one-dimensional proper orthogonal decomposition (POD) applied to a two-dimensional turbulent channel. The constructed wall model contains two modes, i.e. the LoW-based mode and the POD-based mode, and the model matches with the LES at two, instead of one, off-wall locations. To show that the proposed model captures non-equilibrium effects, we perform a priori and a posteriori tests in the context of both equilibrium and non-equilibrium flows. The a priori tests show that the proposed wall model captures extreme wall-shear stress events better than the equilibrium wall model. The model also captures non-equilibrium effects due to adverse pressure gradients. The a posteriori tests show that the wall model captures the rapid decrease and the initial decrease of the streamwise wall-shear stress in channels subjected to suddenly imposed adverse and transverse pressure gradients, respectively, both of which are missed by currently available wall models. These results show promise in applying modal analysis for turbulence wall modelling. In particular, the results show that employing multiple modes helps in the modelling of non-equilibrium flows.

The stratified inner-shelf response to surf-zone-generated transient rip currents (TRC) is examined using idealized simulations with uniform initial thermal stratification $\mathrm {d}{T_0}/\mathrm {d}{z}$, with initial temperature T0 and vertical coordinate z, or initial squared buoyancy frequency $N_0^2$, varying from unstratified to highly stratified ($0.75^\circ \text {C}\ \text {m}^{-1}$). The TRC-induced depth-integrated cross-shore eddy kinetic energy flux is independent of $\mathrm {d}{T_0}/\mathrm {d}{z}$ and decays to near zero within $4L_{SZ}$ (surf-zone width $L_{SZ}=100$ m). Cross-shore inhomogeneous TRC mixing causes shoreward broadening isotherms, driving a near-field inner-shelf overturning circulation and a far-field geostrophic along-shore velocity that strengthen with $\mathrm {d}{T_0}/\mathrm {d}{z}$. TRC mixing mostly (90 %) increases background potential energy (BPE) and also available (APE, 10 %) potential energy, driving inner-shelf mean circulation. The specific BPE zero-crossing depth $d_{s}$ collapses the near-field $3$-layer cross-shore velocity, using an intrusive gravity current scaling $(d_{s} N_0)$. The approximately steady exchange flow exports low-stratification fluid ($N/N_0\approx 1/2$) at a depth $z/d_{s}\approx -1$, re-supplying the TRC region with stratified fluid from above/below, similar to localized mixing laboratory experiments. Offshore of $\approx 5L_{SZ}$, a self-similar far-field intrusion with characteristic isotherm slope $d_{s}/L_{R}$ (Rossby deformation radius $L_{R}\sim d_{s} f/N_0$) is in approximate geostrophic balance with the non-dimensional along-shore velocity. Inner-shelf near-field and far-field horizontal length scales vary as $x/d_{s}$ and $x/L_{R}$, respectively. The length scale $d_{s}$ is related to the work performed by TRC mixing using an idealized well-mixed wedge geometry. Idealized analytical scalings are qualitatively consistent with modelled BPE and APE distributions. Thus, the self-similar stratified inner-shelf response to TRC-driven mixing depends on key dimensional parameters $N_0$ and $d_{s}$.

An additional distant wall is known to highly alter the jetting scenarios of wall-proximal bubbles. Here, we combine high-speed photography and axisymmetric volume of fluid (VoF) simulations to quantitatively describe its role in enhancing the micro-jet dynamics within the directed jet regime (Zeng et al., J. Fluid Mech., vol. 896, 2020, A28). Upon a favourable agreement on the bubble and micro-jet dynamics, both experimental and simulation results indicate that the micro-jet velocity increases dramatically as $\eta$ decreases, where $\eta =H/R_{max}$ is the distance between two walls $H$ normalized by the maximum bubble radius $R_{max}$. The mechanism is related to the collapsing flow, which is constrained by the distant wall into a reverse stagnation-point flow that builds up pressure near the bubble's top surface and accelerates it into micro-jets. We further derive an equation expressing the micro-jet velocity $U_{jet}=87.94\gamma ^{0.5}(1+(1/3)(\eta -\lambda ^{1.2})^{-2})$, where ${\gamma =d/R_{max}}$ is the stand-off distance to the proximal wall with $d$ the distance between the initial bubble centre and the wall, $\lambda =R_{y,m}/R_{max}$ with $R_{y,m}$ the distance between the top surface and the proximal wall at the bubble's maximum expansion. Viscosity has a minimal impact on the jet velocity for small $\gamma$, where the pressure buildup is predominantly influenced by geometry.

The highly nonlinear evolution of the single-mode stratified compressible Rayleigh–Taylor instability (RTI) is investigated via direct numerical simulation over a range of Atwood numbers ($A_T=0.1$–$0.9$) and Mach numbers ($Ma=0.1$–$0.7$) for characterising the isothermal background stratification. After the potential stage, it is found that the bubble is accelerated to a velocity which is well above the saturation value predicted in the potential flow model. Unlike the bubble re-acceleration behaviour in quasi-incompressible RTI with uniform background density, the characteristics in the stratified compressible RTI are driven by not only vorticity accumulation inside the bubble but also flow compressibility resulting from the stratification. Specifically, in the case of strong stratification and high $A_T$, the flow compressibility dominates the bubble re-acceleration characters. To model the effect of flow compressibility, we propose a novel model to reliably describe the bubble re-acceleration behaviours in the stratified compressible RTI, via introducing the dilatation into the classical model that takes into account only vorticity accumulation.

When inertial particles are dispersed in a turbulent flow at sufficiently high concentrations, the continuous and dispersed phases are two-way coupled. Here, we show via laboratory measurements how, as the suspended particles modify the turbulence, their behaviour is also profoundly changed. In particular, we investigate the spatial distribution and motion of sub-Kolmogorov particles falling in homogeneous air turbulence. We focus on the regime considered in Hassaini & Coletti (J. Fluid Mech., vol. 949, 2022, A30), where the turbulent kinetic energy and dissipation rate were found to increase as the particle volume fraction increases from $10^{-6}$ to $5\times 10^{-5}$. This leads to strong intensification of the clustering, encompassing a larger fraction of the particles and over a wider range of scales. The settling rate is approximately doubled over the considered range of concentrations, with particles in large clusters falling even faster. The settling enhancement is due in comparable measure to the predominantly downward fluid velocity at the particle location (attributed to the collective drag effect) and to the larger slip velocity between the particles and the fluid. With increasing loading, the particles become less able to respond to the fluid fluctuations, and the random uncorrelated component of their motion grows. Taken together, the results indicate that the concentrated particles possess an effectively higher Stokes number, which is a consequence of the amplified dissipation induced by two-way coupling. The larger relative velocities and accelerations due to the increased fall speed may have far-reaching consequences for the inter-particle collision probability.

A slowly varying modes solution of Wentzel–Kramers–Brillouin type is derived for the problem of sound propagation in a slowly varying two-dimensional duct with homentropic inviscid sheared mean flow and acoustically lined walls of slowly varying impedance. The modal shape function and axial wavenumber are described by the Pridmore-Brown eigenvalue equation. The slowly varying modal amplitude is determined in the usual way by an equation resulting from a solvability condition. For a general mean flow, this equation can be solved in the form of an incomplete adiabatic invariant. Due to conservation of specific mean vorticity along streamlines, two simplifications prove possible for a linearly sheared mean flow: (i) an analytically exact approximation for the mean flow, and (ii) a complete adiabatic invariant for the acoustics. For this last configuration some example cases are evaluated numerically, where the Pridmore-Brown eigenvalue problem is solved by a Galerkin projection combined with an efficient nonlinear iteration.

Self-sustained, low-frequency, coherent flow unsteadiness over rigid, stationary aerofoils in the transonic regime is referred to as transonic buffet. This study examines the role of shock waves in sustaining this transonic phenomenon and its relation to low-frequency oscillations (LFO) that occur in flow over aerofoils in the incompressible regime (Zaman et al., J. Fluid Mech., vol. 202, 1989, pp. 403–442). This is investigated by performing large-eddy simulations of the flow over a NACA0012 profile for a wide range of flow conditions under free-transition conditions. At low Reynolds numbers, zero incidence angle and sufficiently high free-stream Mach numbers, $M$, transonic buffet occurs with shock waves present in the flow. However, when $M$ alone is lowered, self-sustained, periodic oscillations at a low frequency are observed even though shock waves are absent and the entire flow field remains subsonic at all times. At higher incidence angles, the oscillations are sustained at progressively lower $M$ and are present even at $M=0.3$, where compressibility effects are low. A spectral proper orthogonal decomposition (SPOD) shows that the spatial structure of these oscillations is consistent for all cases. The SPOD modes are topologically similar, suggesting a connection between transonic buffet and LFO in the incompressible regime. Comparisons with other studies examining transonic buffet on various aerofoils, under forced-transition and fully turbulent conditions support this hypothesis. Future studies using tools of global linear stability analysis, especially at high free-stream Reynolds numbers are required to examine whether the underlying mechanisms of transonic buffet and incompressible LFO are the same.