In confined systems, the entrapment of a gas volume with an equivalent spherical diameter greater than the dimension of the channel can form extended bubbles that obstruct fluid circuits and compromise performance. Notably, in sealed vertical tubes, buoyant long bubbles cannot rise if the inner tube radius is below a critical value near the capillary length. This critical threshold for steady ascent is determined by geometric constraints related to the matching of the upper cap shape with the lubricating film surrounding the elongated part of the bubble. Developing strategies to overcome this threshold and release stuck bubbles is essential for applications involving narrow liquid channels. Effective strategies involve modifying the matching conditions with an external force field to facilitate bubble ascent. However, it is unclear how changes in acceleration conditions affect the motion onset of buoyancy-driven long bubbles. This study investigates the mobility of elongated bubbles in sealed tubes with an inner radius near the critical value inhibiting bubble motion in a vertical setting. Two strategies are explored to tune bubble motion, leveraging variations in axial and transversal accelerations: tube rotation around its axis and tube inclination relative to gravity. By revising the geometrical constraints of the simple vertical setting, the study predicts new thresholds based on rotational speed and tilt angle, respectively, providing forecasts for the bubble rising velocity under modified apparent gravity. Experimental measurements of motion threshold and rising velocity compare well with theoretical developments, thus suggesting practical approaches to control and tune bubble motion in confined environments.

We examine the mechanisms responsible for the onset of the three-dimensional mode B instability in the wake behind a circular cylinder. We show that it is possible to explicitly account for the stabilising effect of spanwise viscous diffusion and then demonstrate that the remaining mechanisms involved in this short-wavelength instability are preserved in the limit of zero wavelength. Using the resulting simplified equations, we show that perturbations in different fluid particles interact only through the in-plane viscous diffusion which turns out to have a destabilising effect. We also show that in the presence of viscous diffusion, the closed trajectories which had been conjectured to play a crucial role in the onset of the mode B instability are not actually a prerequisite for the growth of mode B type perturbations. We combine these observations to identify the three essential ingredients for the development of the mode B instability: (i) the amplification of perturbations in the braid regions due to the stretching mechanism; and the spreading of perturbations through (ii) viscous diffusion, and (iii) cross-flow advection which transports fluid between the two braid regions on either side of the cylinder. Finally, we develop a simple criterion that allows the prediction of the regions where three-dimensional short-wavelength perturbations are amplified by the stretching mechanism. The approach used in our study is general and has the potential to give insights into the onset of three-dimensionality via short-wavelength instabilities in other flows.

Aqueous foams coarsen with time due to gas diffusion through the liquid between the bubbles. The mean bubble size grows, and small bubbles vanish. However, coarsening is little understood for foams with an intermediate liquid content, particularly in the presence of surfactant-induced attractive forces between the bubbles, measured by the interface contact angle where thin films meet the bulk liquid. Rigorous bubble growth laws have yet to be developed, and the evolution of bulk foam properties is unclear. We present a quasistatic numerical model for coarsening in two-dimensional wet foams, focusing on growth laws and related bubble properties. The deformation of bubble interfaces is modelled using a finite-element approach, and the gas flow through both films and Plateau borders is approximated. We give results for disordered two-dimensional wet foams with $256$ to $1024$ bubbles, at liquid fractions from $2\,\%$ to $25\,\%$, beyond the zero-contact-angle unjamming transition, and with contact angles up to $10^\circ$. Simple analytical models for the bubble pressures, film lengths and coarsening growth rates are developed to aid interpretation. If the contact angle is non-zero, we find that a prediction of the coarsening rate approaches a non-zero value as the liquid fraction is increased. We also find that an individual bubble's effective number of neighbours determines whether it grows or shrinks to a good approximation.

This study explores heat and turbulent modulation in three-dimensional multiphase Rayleigh–Bénard convection using direct numerical simulations. Two immiscible fluids with identical reference density undergo systematic variations in dispersed-phase volume fractions, $0.0 \leq \varPhi \leq 0.5$, and ratios of dynamic viscosity, $\lambda _{\mu }$, and thermal diffusivity, $\lambda _{\alpha }$, within the range $[0.1\unicode{x2013}10]$. The Rayleigh, Prandtl, Weber and Froude numbers are held constant at $10^8$, $4$, $6000$ and $1$, respectively. Initially, when both fluids share the same properties, a 10 % Nusselt number increase is observed at the highest volume fractions. In this case, despite a reduction in turbulent kinetic energy, droplets enhance energy transfer to smaller scales, smaller than those of single-phase flow, promoting local mixing. By varying viscosity ratios, while maintaining a constant Rayleigh number based on the average mixture properties, the global heat transfer rises by approximately 25 % at $\varPhi =0.2$ and $\lambda _{\mu }=10$. This is attributed to increased small-scale mixing and turbulence in the less viscous carrier phase. In addition, a dispersed phase with higher thermal diffusivity results in a 50 % reduction in the Nusselt number compared with the single-phase counterpart, owing to faster heat conduction and reduced droplet presence near walls. The study also addresses droplet-size distributions, confirming two distinct ranges dominated by coalescence and breakup with different scaling laws.

We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity $\tilde {w}_b$ compared with the rise velocity in quiescent flow $w_b$. The decrease in mean rise velocity of bubbles $\tilde {w}_b/w_b$ is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number $Fr=u'/\sqrt {gd}$, where $u'$ is the root-mean-square velocity fluctuations, $g$ is gravity and $d$ is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For $Fr > 0.5$, we confirm the scaling $\tilde {w}_b / w_b \propto 1 / Fr$, as proposed previously by Ruth et al. (J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability.

In convergent geometry, the effect of convergence and compression on the Rayleigh–Taylor instability (RTI) and Richtmyer–Meshkov instability (RMI) modifies the growth rate and behaviour of the instabilities. In order to better understand how compression/expansion caused by axial strain rates (i.e. strain rates normal to the interface) change the instability dynamics, axial strain rates are applied to RMI in planar geometry, isolated from the effects of convergence. Potential flow theory for the linear regime shows the growth rate of the instability is modified to include the background velocity difference of the instability's width. Resolved two-dimensional simulations of single-mode RMI showed the potential flow model is accurate whilst the amplitude is small compared with the wavelength. The application of strain rate to an RMI-induced mixing layer was investigated using three-dimensional implicit large eddy simulations (ILES) of the quarter-scale $\theta$-group case by Thornber et al. (Phys. Fluids, vol. 29, 2017, 105107). Whilst the background strain rate contributed to the mixing layer's growth, it was to a smaller extent than expected. The shear production of axial turbulent kinetic energy from the strain rate modified the rate of bulk entrainment, affecting the mixing layer's growth and mixedness, such that the strained simulations no longer attained the same self-similar state. The capability of the buoyancy-drag model by Youngs & Thornber (Physica D, vol. 410, 2020, 132517) to predict the mixing layer width was investigated, using a model calibrated to the unstrained case. New terms were introduced into the buoyancy-drag model, which correspond to the shear production of turbulent kinetic energy.

We investigate theoretically the steady incompressible viscoelastic flow in a rigid axisymmetric cylindrical pipe with varying cross-section. We use the Oldroyd-B viscoelastic constitutive equation to model the fluid viscoelasticity. First, we derive exact general formulae: for the total average pressure-drop as a function of the wall shear rate and the viscoelastic axial normal extra-stress; for the viscoelastic extra-stress tensor and the Trouton ratio as functions of the fluid velocity on the axis of symmetry; and for the viscoelastic extra-stress tensor along the wall in terms of the shear rate at the wall. Then we exploit the classic lubrication approximation, valid for small values of the square of the aspect ratio of the pipe, to simplify the original governing equations. The final equations are solved analytically using a regular perturbation scheme in terms of the Deborah number, De, up to eighth order in De. For a hyperbolically shaped pipe, we reveal that the reduced pressure-drop and the Trouton ratio can be recast in terms of a modified Deborah number, Dem, and the polymer viscosity ratio, η, only. Furthermore, we enhance the convergence and accuracy of the eighth-order solutions by deriving transformed analytical formulae using Padé diagonal approximants. The results show the decrease of the pressure drop and the enhancement of the Trouton ratio with increasing Dem and/or increasing η. Comparison of the transformed solutions with numerical simulations of the lubrication equations using pseudospectral methods shows excellent agreement between the results, even for high values of Dem and all values of η, revealing the robustness, validity and efficiency of the theoretical methods and techniques developed in this work. Last, it is shown that the exact solution for the Trouton ratio gives a well-defined and finite solution for any value of Dem and reveals the reason for the failure of the corresponding high-order perturbation series for Dem > 1/2.

Wave-assisted propulsion (WAP) systems directly convert wave energy into thrust using elastically mounted hydrofoils. The wave conditions as well as the design of the hydrofoil drive the fluid–structure interaction of the hydrofoil and, consequently, its performance. We employ simulations using a sharp-interface immersed boundary method to examine the effect of three key parameters on the flow physics, the fluid–structure interaction as well as thrust performance of these systems – the stiffness of the torsional spring, the location of the pitch axis and the Strouhal number. We demonstrate the utility of ‘maps’ of energy exchange between the flow and the hydrofoil system, as a way to understand and predict these characteristics. The force-partitioning method (FPM) is used to decompose the pressure forces into interpretable components and to quantify the mechanisms associated with thrust generation. Based on the results from FPM, a phenomenological model for the thrust generated by the WAP foil is presented. The parameters associated with this model are estimated based on data from over 450 distinct simulations. The predictions of the model are compared with the simulations and the use of this model for guiding WAP design is discussed.

We investigate the dynamics of heavy inertial particles in a flow field due to an isolated, non-axisymmetric vortex. For our study, we consider a canonical elliptical vortex – the Kirchhoff vortex and its strained variant, the Kida vortex. Contrary to the anticipated centrifugal dispersion of inertial particles, which is typical in open vortical flows, we observe the clustering of particles around co-rotating attractors near the Kirchhoff vortex due to its non-axisymmetric nature. We analyse the inertia-modified stability characteristics of the fixed points, highlighting how some of the fixed points migrate in physical space, collide and then annihilate with increasing particle inertia. The introduction of external straining, the Kida vortex being an example, introduces chaotic tracer transport. Using a Melnikov analysis, we show that particle inertia and external straining can compete, where chaotic transport can be suppressed beyond a critical value of particle inertia.

Quantifying transport by strongly stratified turbulence in low Prandtl number ($Pr$) fluids is critically important for the development of better models for the structure and evolution of stellar and planetary interiors. Motivated by recent numerical simulations showing strongly anisotropic flows suggestive of a scale-separated dynamics, we perform a multiscale asymptotic analysis of the governing equations. We find that, in all cases, the resulting slow–fast systems take a quasilinear form. Our analysis also reveals the existence of several distinct dynamical regimes depending on the emergent buoyancy Reynolds and Péclet numbers, $Re_b = \alpha ^2 Re$ and $Pe_b = Pr Re_b$, respectively, where $\alpha$ is the aspect ratio of the large-scale turbulent flow structures, and $Re$ is the outer-scale Reynolds number. Scaling relationships relating the aspect ratio, the characteristic vertical velocity and the strength of the stratification (measured by the Froude number $Fr$) naturally emerge from the analysis. When $Pe_b \ll \alpha$, the dynamics at all scales is dominated by buoyancy diffusion, and our results recover the scaling laws empirically obtained from direct numerical simulations by Cope et al. (J. Fluid Mech., vol. 903, 2020, A1). For $Pe_b \ge O(1)$, diffusion is negligible (or at least subdominant) at all scales and our results are consistent with those of Chini et al. (J. Fluid Mech., vol. 933, 2022) for strongly stratified geophysical turbulence at $Pr =O(1)$. Finally, we have identified a new regime for $\alpha \ll Pe_b \ll 1$, in which slow, large scales are diffusive while fast, small scales are not. We conclude by presenting a map of parameter space that clearly indicates the transitions between isotropic turbulence, non-diffusive stratified turbulence, diffusive stratified turbulence and viscously dominated flows, and by proposing parameterisations of the buoyancy flux, mixing efficiency and turbulent diffusion coefficient for each regime.

We investigate exact nonlinear waves on surfaces locally approximating the rotating sphere for two-dimensional inviscid incompressible flow. Our first system corresponds to a $\beta$-plane approximation at the equator, and the second to a $\gamma$ approximation, with the latter describing flow near the poles. We find exact wave solutions in the Lagrangian reference frame that cannot be written down in closed form in the Eulerian reference frame. The wave particle trajectories, contours of potential vorticity and Lagrangian mean velocity take relatively simple forms. The waves possess a non-trivial Lagrangian mean flow that depends on the amplitude of the waves and on a particle label that characterizes values of constant potential vorticity. The mean flow arises due to potential vorticity conservation on fluid particles. Solutions over the entire space are generated by assuming that the flow far from the origin is zonal and there is a region of uniform potential vorticity between this zonal flow and the waves. In the $\gamma$ approximation, a class of waves is found that, based on analogous solutions on the plane, we call Ptolemaic vortex waves. The mean flow of some of these waves, which we can describe in highly nonlinear scenarios due to the exact nature of the solutions, resembles polar jet streams. Several illustrative solutions are used as initial conditions in the fully spherical rotating Navier–Stokes equations, where integration is performed via the numerical scheme presented in Salmon & Pizzo (Atmosphere, vol. 14, issue 4, 2023, 747). The potential vorticity contours found from these numerical experiments vary between stable permanent progressive form and fully turbulent flows generated by wave breaking.

We present a thin-film viscoplastic fluid model of a mountain range, where a uniform fluid layer is deformed by a vertical backstop moving at constant speed. This represents a simplification of the geometry in subduction zones, where an overlying tectonic plate scrapes sediments off the underthrusting plate, thereby forming an accretionary wedge. By using a viscoplastic rheology, we aim to generalise Newtonian models that capture the effective viscous behaviour of rock at large length scales. The model system is characterised by the dimensionless Bingham number, which is the ratio of the yield stress to characteristic shear stress. At low and high Bingham numbers and at early and late times the system is found to be asymptotically self-similar, which we confirm by solving the governing equations numerically. In addition, we test the high-Bingham-number results experimentally using ultrasound gel as the working fluid. The experiments reproduce many features of the theoretically predicted behaviour. The size of the observed wedge grows in the manner predicted, and the fluid surface profile is found to collapse to a universal shape. The viscoplastic fluid wedge exhibits features of both viscous continuum and Coulomb models of accretionary wedges, suggesting that a viscoplastic rheology may provide quantitative insights into the dynamics of real accretionary wedges found in convergent tectonic settings.

We develop a model for the interaction of a fluid flowing above an otherwise static particle bed, with generally the particles being entrained or detrained into the fluid from the upper surface of the particle bed, and thereby forming a fully two phase fluidized cloud above the particle bed. The flow in this large-scale fluidized region is treated as a two-phase flow, whilst the key processes of entrainment and detrainment from the particle bed are treated by examining the local dynamical force balances on the particles in a thin transition layer at the interface between the fully fluidized region and the static particle bed. This detailed consideration leads to the formation of an additional macroscopic boundary condition at this interface, which closes the two-phase flow problem in the bulk fluidized region above. We then introduce an elementary model of the well-known helicopter brownout problem, and use the theory developed in the first part of the paper to fully analyse this model, both analytically and numerically.

Supersonic shear layers experience instabilities that generate significant adverse effects; in complex configurations, these instabilities have global impacts as they foster compounding complications with other independent flow features. We consider the flow near the exit of a dual-stream rectangular nozzle. The supersonic core and sonic bypass streams mix downstream of a splitter plate trailing edge (SPTE) just above an adjacent deck. We perform two-dimensional direct numerical simulations with laminar inflow conditions to parametrically explore the influence of active flow control, considering various actuation angles and locations. The goal is to alleviate the prominent tone associated with vortices shed at the SPTE; these vortices initiate an unsteady shock system that affects the entire flow field through a shock-induced separation and the downstream evolution of plume shear layers. Resolvent analysis is performed on the baseline flow. The identified optimal response guides the placement of steady-blowing actuators. Since the resolvent analysis fundamentally investigates the input–output dynamics of a system, it is also utilized to uncover actuation-induced changes in the forcing–response dynamics. Spectral analysis shows that the baseline flow fluctuating energy is concentrated in the shedding instability. Actuating at optimal angles based on location disperses this energy into various flow features; this affects the shedding itself, and the structure and unsteadiness of the shock system, and thus the response of the deck and nozzle wall boundary layers and the plume. The resolvent analysis indicates, and Navier–Stokes solutions confirm, that favourable control is obtained by either indirectly or directly mitigating the baseline instability.

A hydrodynamic theory of premixed flame propagation within closed vessels is developed assuming the flame is much thinner than all other fluid dynamic lengths. In this limit, the flame is confined to a surface separating the unburned mixture from burned combustion products, and propagates at a speed determined from the analysis of its internal structure. Unlike freely propagating flames that propagate under nearly isobaric conditions, combustion in a closed vessel results in continuous increases in pressure, burning rate and flame temperature, and a progressive decrease in flame thickness. The flame speed is shown to depend on the voluminal stretch rate, which measures the deformation of a volume element of the flame zone, and on the rate of pressure rise. Both effects are modulated by pressure-dependent Markstein numbers that depend on heat release and mixture properties while capturing the effects of temperature-dependent transport and stoichiometry. The model applies to flames of arbitrary shape propagating in general flows, laminar or turbulent, within vessels of general configurations. The main limitation of hydrodynamic flame theories is the assumption that variations inside the flame zone due to chemistry or turbulence, which could potentially alter its internal structure, are physically unresolved. Nonetheless, the theory, deduced from physical first principles, identifies the various mechanisms involved in the combustion process as demonstrated in detailed discussions of planar flames propagating in rectangular channels and spherically expanding flames in spherical vessels. It also enables the construction of instructive models to numerically simulate the evolution of multi-dimensional and corrugated flames under confinement.

In the generalized quasilinear approximation (GQLA) (Marston et al., Phys. Rev. Lett., vol. 116, 2016, 214501), a threshold wavenumber ($k_0$) in the direction of translational symmetry segregates the total into large- ($l$) and small-scale ($h$) fields. While the governing equation for the large-scale field is fully nonlinear, that for the small scales is linearized with respect to the large-scale field. In addition, some nonlinear triad interactions are omitted in the GQLA. Herein, the GQLA is applied to two-dimensional planar Rayleigh–Bénard convection (RBC). A scale separation between the large-scale convection rolls and small-scale turbulent fluctuations is typical in RBC. The present work explores the efficacy of GQLA in capturing the scale-by-scale energy transfer processes in RBC. The initial condition for the GQLA simulations was either the statistically stationary state obtained in direct numerical simulation (DNS) or random fluctuations superimposed on the linear conductive temperature profile and $\boldsymbol {u}=0$. The GQLA simulations can capture the convection rolls for $k_0$ larger than or equal to the dominant wavenumber for thermal driving of the flow ($k_0 \ge k_{\hat {Q}}$). Additionally, the GQLA emulates the fully nonlinear dynamics for $k_0$ larger than or equal to the first harmonic of the convection-roll wavenumber ($k_0 \ge 2k_{roll}$). In the intermediate regime with $2k_{roll} > k_0 \ge k_{\hat {Q}}$, the dynamics captured in GQLA simulations is different from the DNS. In DNS, two primary energy transfer processes dominate: (i) the energy transfer to/from the convection rolls and (ii) the scale-by-scale inverse kinetic energy and forward thermal energy cascades mediated by the convection rolls. The fully nonlinear dynamics is emulated by GQLA when these energy transfer processes are faithfully reproduced. Utilizing the framework of altering the triad interactions in GQLA, an additional intrusive calculation, including target triad interactions, is performed here to study their influence. This intrusive calculation shows that the convection rolls are not captured in GQLA for $k_0 < k_{\hat {Q}}$ because of the exclusion of the $h \to l \to l$ and $l \to h \to h$ triad interactions in GQLA. The inclusion of these triad interactions in the intrusive calculation yields the convection rolls, and the reproduced dynamics is similar to that of the intermediate GQLA regime with $2k_{roll} > k_0 \ge k_{\hat {Q}}$.

Shock refraction in a gas–liquid interface is ubiquitous in nature and engineering. This study investigates the shock refraction phenomena in air–water interfaces for various inclination angles. The interface inclination angles are achieved using a tiltable vertical shock tube. The time-resolved schlieren images are compared with numerical simulations performed using the BlastFoam solver in the OpenFOAM software. The stiffened gas equation of state is used to model water in the simulations. The shock polar analysis using modified shock relations for a stiffened gas is used to elucidate the refraction patterns. A regular refraction pattern with a reflected shock wave and a bound precursor refraction with a regular reflection are observed experimentally for the first time in an air–water system. Further, a new free precursor refraction pattern with a Mach reflection is observed. The transition criteria and the corresponding boundaries for each refraction pattern are demarcated in the ($M_S, \theta _w^c$)-plane. The refraction sequence and the range for various incident shock strength regimes are also identified.

Shear-induced self-diffusion is a fundamental mode of transport in granular flows. Yet its critical behaviour and dependence on the particle solid fraction are still unclear. Here, we rationalize these dependencies by performing two-dimensional pressure-imposed numerical simulations of dense non-Brownian frictional suspensions. Our results, combined with existing numerical data on inertial granular flows, show that the shear-induced diffusion coefficients of both systems can be captured by a single function of the distance to jamming. They further show that the grain diffusive behaviour is underpinned by a specific random walk process, having a constant elementary step length driven at a frequency that increases with the solid fraction. The proposed scaling laws pave the way for a better understanding of mixing processes in granular media.

We examine the fluid flow forced by precession of a rotating cylindrical container using numerical simulations and experimental flow measurements with ultrasonic Doppler velocimetry. The analysis is based on the decomposition of the flow field into contributions with distinct azimuthal symmetry or analytically known inertial modes and the corresponding calculation of their amplitudes. We show that the predominant fraction of the kinetic energy of the precession-driven fluid flow is contained only within a few large-scale modes. The most striking observation shown by simulations and experiments is the transition from a flow dominated by large-scale structures to a more turbulent behaviour with the small-scale fluctuations becoming increasingly important. At a fixed rotation frequency (parametrized by the Reynolds number, $Re$) this transition occurs when a critical precession ratio is exceeded and consists of a two-stage collapse of the directly driven flow going along with a massive modification of the azimuthal circulation (the zonal flow) and the appearance of an axisymmetric double-roll mode limited to a narrow range of precession ratios. A similar behaviour is found in experiments which make it possible to follow the transition up to Reynolds numbers of $Re\approx 2\times 10^6$. We find that the critical precession ratio decreases with rotation, initially showing a particular scaling ${\propto }Re^{-({1}/{5})}$ but developing an asymptotic behaviour for $Re\gtrsim 10^5$ which might be explained by the onset of turbulence in boundary layers.

We discuss the applicability of quasilinear-type approximations for a turbulent system with a large range of spatial and temporal scales. We consider a paradigm fluid system of rotating convection with vertical and horizontal temperature gradients. In particular, the interaction of rotation with the horizontal temperature gradient drives a ‘thermal wind’ shear flow whose strength is controlled by the horizontal temperature gradient. Varying this parameter therefore systematically alters the ordering of the shearing time scale, the convective time scale and the correlation time scale. We demonstrate that quasilinear-type approximations work well when the shearing time scale or the correlation time scale is sufficiently short. In all cases, the generalised quasilinear approximation systematically outperforms the quasilinear approximation. We discuss the consequences for statistical theories of turbulence interacting with mean gradients.