A model for galloping detonations is conceived as a sequence of very fast re-ignitions followed by long periods of evolution with quenched reactions. Numerical simulations of the one-dimensional Euler equations are conducted in this limit. While the mean speed and structure is found in reasonable agreement with Chapman–Jouguet theory, very strong pulsations of the lead shock appear, along with a train of rear-facing N-waves. These dynamics are analysed using characteristics. A closed-form solution for the lead shock dynamics is formulated, which is found in excellent agreement with numerics. The model relies on the presence of a single time scale of the process, the pulsation period, which controls the shock dynamics via the shock change equations and establishes a shock decay with a single time constant. These long periods of shock decay with known dynamics are punctuated by energy release events, with ‘kicks’ in the shocked speed controlled by the pressure increase and resulting lead shock amplification. Model predictions are found in excellent agreement with previous numerical results of pulsating detonations far from the stability limit.

Elastic turbulence can lead to increased flow resistance, mixing and heat transfer. Its control – either suppression or promotion – has significant potential, and there is a concerted ongoing effort by the community to improve our understanding. Here we explore the dynamics of uncertainty in elastic turbulence, inspired by an approach recently applied to inertial turbulence in Ge et al. (J. Fluid Mech., vol. 977, 2023, A17). We derive equations for the evolution of uncertainty measures, yielding insight on uncertainty growth mechanisms. Through numerical experiments, we identify four regimes of uncertainty evolution, characterised by (i) rapid transfer to large scales, with large-scale growth rates of $\tau ^{6}$ (where $\tau$ represents time), (ii) a dissipative reduction of uncertainty, (iii) exponential growth at all scales and (iv) saturation. These regimes are governed by the interplay between advective and polymeric contributions (which tend to increase uncertainty), viscous, relaxation and dissipation effects (which reduce uncertainty) and inertial contributions. In elastic turbulence, reducing Reynolds number increases uncertainty at short times, but does not significantly influence the growth of uncertainty at later times. At late times, the growth of uncertainty increases with Weissenberg number, with decreasing polymeric diffusivity and with the logarithm of the maximum length scale, as large flow features adjust the balance of advective and relaxation effects. These findings provide insight into the dynamics of elastic turbulence, offering a new approach for the analysis of viscoelastic flow instabilities.

In this study, we experimentally investigate the stress field around a gradually contaminated bubble as it moves straight ahead in a dilute surfactant solution with an intermediate Reynolds number ($20 \lt {{\textit{Re}}} \lt 220$) and high Péclet number. Additionally, we investigate the stress field around a falling sphere unaffected by surface contamination. A newly developed polarisation measurement technique, highly sensitive to the stress field in the vicinity of the bubble or the sphere, was employed in these experiments. We first validated this method by measuring the flow around a solid sphere sedimenting in a quiescent liquid at a terminal velocity. The measured stress field was compared with established numerical results for ${{\textit{Re}}} = 120$. A quantitative agreement with the numerical results validated this technique for our purpose. The results demonstrated the ability to determine the boundary layer. Subsequently we measured a bubble rising in a quiescent surfactant solution. The drag force on the bubble, calculated from its rise velocity, was set to transiently vary from that of a clean bubble to a solid sphere within the measurement area. With the intermediate drag force between clean bubble and solid sphere, the stress field in the vicinity of the bubble front was observed to be similar to that of a clean bubble, and the structure near the rear was similar to that of a solid sphere. Between the front and rear of the bubble, the phase retardation exhibited a discontinuity around the cap angle at which the boundary conditions transitioned from no slip to slip, indicating an abrupt change in the flow structure. A reconstruction of the axisymmetric stress field from the phase retardation and azimuth obtained from polarisation measurements experimentally revealed that stress spikes occur around the cap angle. The cap angle (stress jump position) shifted as the drag on the bubble increased owing to surfactant accumulation on its surface. Remarkably, the measured cap angle as a function of the normalised drag coefficient quantitatively agreed with the numerical results at intermediate ${{\textit{Re}}} = 100$ of Cuenot et al. (1997 J. Fluid Mech. 339, 25–53), exhibiting only a slight deviation from the curve predicted by the stagnant cap model at low ${\textit{Re}}$ (creeping flow) proposed by Sadhal & Johnson (1983 J. Fluid Mech. 126, 237–250).

This numerical investigation focuses on the mechanisms, flow topology and onset of Kelvin–Helmholtz instabilities (KHIs), that drive the leading-edge shear-layer destabilisation in the wake of wall-mounted long prisms. Large-eddy simulations are performed at ${\textit{Re}} = 2.5\times 10^3, 5\times 10^3$ and $1\times 10^4$ for prisms with a range of aspect ratio (AR, height-to-width) between $0.25$ and $1.5$, and depth ratios (DR, length-to-width) of $1{-}4$. Results show that shear-layer instabilities enhance flow irregularity and modulate spanwise vortex structures. The onset of KHI is strongly influenced by depth ratio, such that long prisms (${\textit{DR}}= 4$) experience earlier initiation compared with shorter ones (${\textit{DR}}= 1$). At higher Reynolds numbers, the onset of KHI shifts upstream towards the leading-edge, intensifying turbulence kinetic energy and increasing flow irregularity, especially for long prisms. The results further show that in this configuration, energy transfer from the secondary recirculation region contributes to the destabilisation of the leading-edge shear layer by reinforcing low-frequency modes. A feedback mechanism is identified wherein energetic flow structures propagate upstream through reverse boundary-layer flow, re-energising the leading-edge shear layer. Quantification using probability density functions reveals rare, intense upstream energy convection events, driven by this feedback mechanism. These facilitate the destabilisation process regardless of Reynolds number. This study provides a comprehensive understanding of the destabilisation mechanisms for leading-edge shear layers in the wake of wall-mounted long prisms.

We study the dynamics of salt fingers in the regime of slow salinity diffusion (small inverse Lewis number) and strong stratification (large density ratio), focusing on regimes relevant to Earth’s oceans. Using three-dimensional direct numerical simulations in periodic domains, we show that salt fingers exhibit rich, multiscale dynamics in this regime, with vertically elongated fingers that are twisted into helical shapes at large scales by mean flows and disrupted at small scales by isotropic eddies. We use a multiscale asymptotic analysis to motivate a reduced set of partial differential equations that filters internal gravity waves and removes inertia from all parts of the momentum equation except for the Reynolds stress that drives the helical mean flow. When simulated numerically, the reduced equations capture the same dynamics and fluxes as the full equations in the appropriate regime. The reduced equations enforce zero helicity in all fluctuations about the mean flow, implying that the symmetry-breaking helical flow is generated spontaneously by strictly non-helical fluctuations.

Undulatory slender objects have been a central theme in the hydrodynamics of swimming at low Reynolds number, where the slender body is usually assumed to be inextensible, although some microorganisms and artificial microrobots largely deform with compression and extension. Here, we theoretically study the coupling between the bending and compression/extension shape modes, using a geometrical formulation of kinematic microswimmer hydrodynamics to deal with the non-commutative effects between translation and rotation. By means of a coarse-grained minimal model and systematic perturbation expansions for small bending and compression/extension, we analytically derive the swimming velocities and report three main findings. First, we revisit the role of anisotropy in the drag ratio of the resistive force theory, and generally demonstrate that no motion is possible for uniform compression with isotropic drag. We then find that the bending–compression/extension coupling generates lateral and rotational motion, which enhances the swimmer’s manoeuvrability, as well as changes in progressive velocity at a higher order of expansion, while the coupling effects depend on the phase difference between the two modes. Finally, we demonstrate the importance of often-overlooked Lie bracket contributions in computing net locomotion from a deformation gait. Our study sheds light on compression as a forgotten degree of freedom in swimmer locomotion, with important implications for microswimmer hydrodynamics, including understanding of biological locomotion mechanisms and design of microrobots.

Cross-shelf transport in the inner continental shelf is governed by wind, wave and tidal interactions, but the role of Langmuir circulation (LC), induced by wave–current interaction and modulated by tides, has remained under-studied in this setting. We develop a Reynolds-averaged Navier–Stokes (RANS) model incorporating the Craik–Leibovich vortex force to resolve LC, coupled with a mass-conserving undertow and oscillating along-shelf tidal currents, and compare results against field data from the Martha’s Vineyard Coastal Observatory (MVCO). Under strong wave forcing (significant wave height $H_{\textit{sig}} = 2.12\,\mathrm{m}$ and significant wave period $T_w = 5.8\,\mathrm{s}$), LC persists throughout the tidal cycle, reducing vertical shear in the tidally averaged cross-shelf velocity profile compared with simulations excluding LC. During peak tidal velocity (reaching $25\,\mathrm{cm\,s^{-1}}$ with period of $ 12.42\,\mathrm{h}$), LC is temporarily suppressed but reforms rapidly as tidal energy declines, sustaining high vertical mixing. Conversely, under weak wave forcing ( $H_{\textit{sig}} = 0.837\,\mathrm{m}$, $T_w = 4.3\,\mathrm{s}$), tidal currents persistently suppress LC, resulting in a cross-shelf undertow profile with greater vertical shear compared with strong-wave conditions. Model–observation comparisons show that only simulations including both the Craik–Leibovich vortex force and tidal forcing reproduce the observed undertow structure at MVCO. These results demonstrate that accurate prediction of cross-shelf transport at tidal and subtidal time scales requires resolving both the generation and disruption of LC by tides.

Understanding the interplay between buoyancy and fluid motions within stably stratified shear layers is crucial for unravelling the contribution of flow structures to turbulent mixing. In this study, we examine statistically the local relationship between stratification and fluid deformation rate in wave and turbulent regimes, using experimental datasets obtained from a stratified inclined duct (SID) containing fluids of different densities that form an exchange flow. We introduce rotational and shear components of varying strength within the vorticity and a family of coherent gradient Richardson numbers ($Ri_C$), ratios related to the buoyancy frequency and the strength of either the rotational or shearing motion. Conditional statistical analysis reveals that both shear and stratification intensity affect the probability distribution of the $Ri_C$, with extreme events occurring more frequently in areas of weak stratification. In the wave regime, we identify the persistence of fast-spin vortices within the strongly stratified density interface. However, scouring of the density interface is primarily driven by shearing motions, with baroclinic torque making a notable contribution to enstrophy transport. In the turbulent regime, rigid-body rotations occur at significantly shorter time scales than that associated with the local buoyancy frequency, making them more disruptive to stratification than shear. Additionally, correlation analysis reveals that irrotational strain distorts stable stratification similarly to shearing motions, but is weaker than both shearing and rotational motions and less likely to have a time scale longer than that related to the buoyancy frequency. Moreover, we observed that the interplay between rotational and shearing motions intensifies as stratification increases. Finally, a comparison of length scales along the shear layers highlights the $Ri_C$ as a valuable measure of the relative sizes of different motions compared with the Ozmidov scale and shows that stratification can influence sub-Ozmidov scales through baroclinic torque. This study highlights the critical impact of the type, strength and location of fluid deformations on localised mixing, providing new insights into the role of rotational motions in shear-driven stratified flows.

Interactions of turbulent boundary layers with a compliant surface are investigated experimentally at Reτ = 3300–8900. Integrating tomographic particle tracking with Mach–Zehnder interferometry enables simultaneous mapping of the compliant wall deformation and the three-dimensional velocity and pressure fields. Our initial study (J. Fluid. Mech. vol. 980, R2) shows that the flow–deformation correlations decrease with increasing Reτ, despite an order of magnitude increase in deformation amplitude. To elucidate the mechanisms involved, the same velocity, pressure and kinetic energy fields are decomposed to ‘wave-coherent’ and ‘stochastic’ parts using a Hilbert projection method. The phase dependent coherent variables, especially the pressure, are highly correlated with the wave, but decrease with increasing Reτ. While the coherent energy is 6 %–10 % of the stochastic level, the pressure root mean square is comparable near the wall. The energy flux between the coherent and stochastic parts and the pressure diffusion reverse sign at the critical layer. To explain the Reτ dependence, the characteristic deformation wavelength (three times the thickness) is compared with the scales of the energy-containing eddies in the boundary layer represented by the k−1 range in the energy spectrum. When the deformation wavelength is matched with the kxEuu peak at the present lowest Reτ, the flow–deformation correlations and coherent pressure become strong, even for submicron deformations. In this case, the flow and wall motion become phase locked, suggesting resonant behaviours. As Reτ increases, the wall wavelengths and spectral range of attached eddies are no longer matched, resulting in reduced correlations and lower coherent energy and pressure, despite larger deformation.

This work combines Navier–Stokes–Korteweg dynamics and rare event techniques to investigate the transition pathways and times of vapour bubble nucleation in metastable liquids under homogeneous and heterogeneous conditions. The nucleation pathways deviate from classical theory, showing that bubble volume alone is an inadequate reaction coordinate. The nucleation mechanism is driven by long-wavelength fluctuations with densities slightly different from the metastable liquid. We propose a new strategy to evaluate the typical nucleation times by inferring the diffusion coefficients from hydrodynamics. The methodology is validated against state-of-the-art nucleation theories in homogeneous conditions, revealing non-trivial, significant effects of surface wettability on heterogeneous nucleation. Notably, homogeneous nucleation is detected at moderate hydrophilic wettabilities despite the presence of a wall, an effect not captured by classical theories but consistent with atomistic simulations. Hydrophobic surfaces, instead, anticipate the spinodal. The proposed approach is fairly general and, despite the paper discussing results for a prototypical fluid, it can be easily extended, also in complex geometries, to any real fluid provided the equation of state is available, paving the way to model complex nucleation problems in real systems.

The paper uses three-dimensional large eddy simulation (LES) to investigate the structure and propagation of dam break waves of non-Newtonian fluids described by a power-law rheology. Simulations are also conducted for the limiting case of a dam-break wave of Newtonian fluid (water). Turbulent dam-break waves are found to have a two-layer structure and to generate velocity streaks beneath the region in which the flow is strongly turbulent and lobes at the front. The bottom part of the wave resembles a boundary layer and contains a log-law sublayer, while the streamwise velocity is close to constant inside the top layer. The value of the von Kármán constant is found to reach the standard value (i.e. $\kappa$ ≈ 0.4) associated with turbulent boundary layers of Newtonian fluids only inside the strongly turbulent region near the front of Newtonian dam-break waves. Much higher values of the slope of the log law are predicted for non-Newtonian dam-break waves (i.e. $\kappa$ ≈ 0.28) and in the regions of weak turbulence of Newtonian waves. LES shows that a power-law relationship can well describe the temporal evolution of the front position during the acceleration and deceleration phases, and that increasing the shear-thinning behaviour of the fluid increases the speed of the front. The numerical experiments are then used to investigate the predictive abilities of shallow water equation (SWE) models. The paper also proposes a novel one-dimensional (1-D) SWE model which accounts for the bottom friction by employing a friction coefficient regression valid for power-law fluids in the turbulent regime. An analytical approximate solution is provided by splitting the current into an outer region, where the flow is considered inviscid and friction is neglected, and an inner turbulent flow region, close to the wave front. The SWE numerical and analytical solutions using a turbulent friction factor are found to be in better agreement with LES compared with the agreement shown by an SWE numerical model using a laminar friction coefficient. The paper shows that inclusion of turbulence effects in SWE models used to predict high-Reynolds-number Newtonian and non-Newtonian dam break flows results is more accurate predictions.

The dynamics of self-propelled colloidal particles is strongly influenced by their environment through hydrodynamic and, in many cases, chemical interactions. We develop a theoretical framework to describe the motion of confined active particles by combining the Lorentz reciprocal theorem with a Galerkin discretisation of surface fields, yielding an equation of motion that efficiently captures self-propulsion without requiring an explicit solution for the bulk fluid flow. Applying this framework, we identify and characterise the long-time behaviours of a Janus particle near rigid, permeable and fluid–fluid interfaces, revealing distinct motility regimes, including surface-bound skating, stable hovering and chemo-hydrodynamic reflection. Our results demonstrate how the solute permeability and the viscosity contrast of the surface influence a particle’s dynamics, providing valuable insights into experimentally relevant guidance mechanisms for autophoretic particles. The computational efficiency of our method makes it particularly well suited for systematic parameter sweeps, offering a powerful tool for mapping the phase space of confined active particles and informing high-fidelity numerical simulations.

We report our finding from direct numerical simulations that polygonal cell structures are formed by inertial particles in turbulent Rayleigh–Bénard convection in a large aspect ratio channel at Rayleigh numbers of $10^6, 10^7$ and $10^8$, and Prandtl number of 0.7. The settling of small particles modified the flow structures only through momentum interactions. From the simulations performed for various sizes and mass loadings of particles, we discovered that for small- and intermediate-sized particles, cell structures such as square, pentagonal or hexagonal cells were observed, whereas a roll structure was formed by large particles. As the mass loading increased, the sizes of the cells or rolls decreased for all particle sizes. The Nusselt number increased with the mass loading of intermediate and large particles, whereas it decreased with the mass loading of small particles compared with the value for particle-free convection. A detailed investigation of the effective feedback forces of the settling particles inside the hot and cold plumes near the walls revealed that the feedback forces break the up–down symmetry between the hot and cold plumes near the surfaces. This enhances the hot plume ascent while not affecting the cold plume, which leads to the preferred formation of cellular structures. The energy budget analysis provides a detailed interaction between particles and fluid, revealing that the net energy is transferred from the fluid to particles when the particles are small, while settling intermediate and large particles drag the fluid so strongly that energy is transferred from particles to fluid.

Thixotropic fluids with a non-monotonic flow curve display viscosity bifurcations at certain stresses. It has been proposed that these transitions can introduce interfaces (or shear bands) into thin films that can destabilize inertialess flows over inclined planes. This proposition is confirmed in the present paper by formulating a thin-film model, then using this model to construct sheet-like base flows and test their linear stability. It is also found that viscosity bifurcations, and the associated interfaces, are not necessary for instability, but that the time-dependent relaxation of the microstructure responsible for thixotropy within the bulk of the film can promote instability instead. Computations with the thin-film model demonstrate that instabilities saturate supercritically into steadily propagating nonlinear waves that travel faster than the mean flow.

Within the context of machine learning-based closure mappings for Reynolds-averaged Navier Stokes turbulence modelling, physical realisability is often enforced using ad hoc postprocessing of the predicted anisotropy tensor. In this study, we address the realisability issue via a new physics-based loss function that penalises non-realisable results during training, thereby embedding a preference for realisable predictions into the model. Additionally, we propose a new framework for data-driven turbulence modelling which retains the stability and conditioning of optimal eddy viscosity-based approaches while embedding equivariance. Several modifications to the tensor basis neural network to enhance training and testing stability are proposed. We demonstrate the conditioning, stability and generalisation of the new framework and model architecture on three flows: flow over a flat plate, flow over periodic hills and flow through a square duct. The realisability-informed loss function is demonstrated to significantly increase the number of realisable predictions made by the model when generalising to a new flow configuration. Altogether, the proposed framework enables the training of stable and equivariant anisotropy mappings, with more physically realisable predictions on new data. We make our code available for use and modification by others. Moreover, as part of this study, we explore the applicability of Kolmogorov–Arnold networks to turbulence modelling, assessing its potential to address nonlinear mappings in the anisotropy tensor predictions and demonstrating promising results for the flat plate case.

Riparian vegetation along riverbanks and seagrass along coastlines interact with water currents, significantly altering their flow. To characterise the turbulent fluid motion along the streamwise-edge of a region covered by submerged vegetation (canopy), we perform direct numerical simulations of a half-channel partially obstructed by flexible stems, vertically clamped to the bottom wall. An intense streamwise vortex forms along the canopy edge, drawing high-momentum fluid into the side of the canopy and ejecting low-momentum fluid from the canopy tip, in an upwelling close to the canopy edge. This mechanism has a profound impact on the mean flow and on the exchange of momentum between the fluid and the structure, which we thoroughly characterise. The signature of the canopy-edge vortex is also found in the dynamical response of the stems, assessed for two different values of their flexibility. Varying the flexibility of the stems, we observe how different turbulent structures over the canopy are affected, while the canopy-edge vortex does not exhibit major modifications. Our results provide a better understanding of the flow in fluvial and coastal environments, informing engineering solutions aimed at containing the water flow and protecting banks and coasts from erosion.

The Lamb–Oseen vortex is a model for practical vortical flows with a finite vortex core. Vortices with a Lamb–Oseen vortex velocity profile are stable according to the Rayleigh criterion in an infinite domain. Practical situations introduce boundary conditions over finite domains. Direct numerical simulations are performed on the evolution of perturbations to a viscous Lamb–Oseen vortex with uniform inlet axial velocity in a pipe of finite length. Linear stability boundaries are determined in the $(\textit{Re},\omega )$ plane. For a given swirl ratio $\omega$, the flow is found to become linearly unstable when the Reynolds number $\textit{Re}$ is above a critical value. The complete evolution history of the flow is followed until it reaches its final state. For small swirl ratios, the axisymmetric mode is linearly unstable and evolves to a final steady axisymmetric but non-columnar accelerated flow state after nonlinear saturation. For large swirl ratios, the spiral mode is linearly unstable. The spiral mode is found to force growth of an axisymmetric component due to nonlinear interaction. The flow evolves to a final unsteady spiral vortex breakdown state after it undergoes nonlinear saturation. The energy transfer between the mean flow and perturbations is studied by the Reynolds–Orr equation. The pressure work at the exit of the finite pipe is a major source of energy production. Finite-domain boundary conditions also modify the perturbation mode shapes, which can render the vortex core from absorbing energy to producing energy, and thus lead to instabilities. As the pipe length increases, the stability behaviour of the flow is found to approach that predicted by the classical Rayleigh criterion.

This paper provides direct experimental evidence for the coexistence of both a laminar separation bubble and a secondary vortex on the advancing side of a rotating sphere when subjected to the inverse Magnus effect. Detailed experiments were conducted in a wind tunnel on two spheres of varying surface roughness to investigate both ordinary and inverse Magnus effects. Experiments took place for $0.5\times 10^{5}\leqslant {\textit{Re}}\leqslant 3\times 10^{5}$ and rotation rates $0\leqslant \alpha \leqslant 0.45$, where the spheres were rotated via a shaft that was oriented perpendicularly to the free stream flow. Static pressure measurements were made on the non-shaft hemisphere using a spline of taps spanning from the equator to the pole. The ordinary Magnus effect was generally observed at the lowest ${\textit{Re}}$ tested, with a transition to the inverse Magnus effect occurring as ${\textit{Re}}$ increased. Time-averaged pressure coefficient distributions across the equatorial plane were obtained for the smooth and rough spheres. Cross-flow particle image velocimetry was used to visualise the downstream wake velocity field. A pair of counter-rotating wing-tip-like vortices were detected when the sphere experienced the ordinary Magnus effect, generated by flow leakage from the advancing to the retreating side. When the sphere experienced the inverse Magnus effect, the polarity of the counter-rotating vortex pair reversed. This is the first experimental observation of the vortex polarity reversal associated with the inverse Magnus effect in the wake of a rotating sphere. The results provide qualitative visualisation of the complex fluid dynamics and inform future applications of the Magnus effect.

Steady flow at low Reynolds (Re) number through a planar channel with converging or diverging width is investigated in this study. Along the primary direction of flow, the small dimension of the channel cross-section remains constant while the sidewalls bounding the larger dimension are oriented at a constant angle. Due in part to ease of manufacturing, parallel-plate geometries such as this have found widespread use in microfluidic devices for mixing, heat exchange, flow control and flow patterning at small length scales. Previous analytical solutions for flows of this nature have required the converging or diverging aspect of the channel to be gradual. In this work, we derive a matched asymptotic solution, validated against numerical modelling results, that is valid for any sidewall angle, without requiring the channel width to vary gradually. To accomplish this, a cylindrical coordinate system defined by the angle of convergence between the channel sidewalls is considered. From the mathematical form of the composite expansion, a delineation between two secondary flow components emerges naturally. The results of this work show how one of these two components, originating from viscous shear near the channel sidewalls, corresponds to convective mixing, whereas the other component impresses the sidewall geometry on streamlines in the outer flow.