Direct numerical simulations of the flow past an impulsively started cylinder at high Reynolds numbers (25k–100k) reveal an intriguing portrait of unsteady separation. Vorticity generation and vortex shedding entails a cascade of separation events on the cylinder surface that are reminiscent of Kelvin–Helmholtz instabilities. Primary vortices roll up along the cylinder surface as a result of instabilities of the initially attached vortex sheets, followed by vortex eruptions, creation of secondary vorticity and formation of dipole structures that are subsequently ejected from the surface of the cylinder. We analyse the vortical structures and their relationship to the forces experienced by the cylinder. This striking cascade of vortex instabilities may serve as reference for reduced-order models of flow separation and as guide for flow control of separated flows at high Reynolds numbers.

A small spheroid settling in a quiescent fluid experiences an inertial torque that aligns it so that it settles with its broad side first. Here we show that an active particle experiences such a torque too, as it settles in a fluid at rest. For a spherical squirmer, the torque is $\boldsymbol {T}^\prime = -{\frac {9}{8}} m_f (\boldsymbol {v}_s^{(0)} \wedge \boldsymbol {v}_g^{(0)})$ where $\boldsymbol {v}_s^{(0)}$ is the swimming velocity, $\boldsymbol {v}_g^{(0)}$ is the settling velocity in the Stokes approximation and $m_f$ is the equivalent fluid mass. This torque aligns the swimming direction against gravity: swimming up is stable, swimming down is unstable.

Instabilities of fluid flows often generate turbulence. Using extensive direct numerical simulations, we study two-dimensional turbulence driven by a wavenumber-localised instability superposed on stochastic forcing, in contrast to previous studies of state-independent forcing. As the contribution of the instability forcing, measured by a parameter $\gamma$, increases, the system undergoes two transitions. For $\gamma$ below a first threshold, a regular large-scale vortex condensate forms. Above this threshold, shielded vortices (SVs) emerge within the condensate. At a second, larger value of $\gamma$, the condensate breaks down, and a gas of weakly interacting vortices with broken symmetry spontaneously emerges, characterised by preponderance of vortices of one sign only and suppressed inverse energy cascade. The latter transition is shown to depend on the damping mechanism. The number density of SVs in the broken symmetry state slowly increases via a random nucleation process. Bistability is observed between the condensate and mixed SV-condensate states. Our findings provide new evidence for a strong dependence of two-dimensional turbulence phenomenology on the forcing.

The well-known self-propulsion, or swimming, of a deformable body in Stokes flow (i.e. at low Reynolds number) can be understood and modelled from the variation in the configuration-dependent hydrodynamic resistance tensor throughout the period of deformation. Remarkably, at the other extreme of high Reynolds number, a deformable body may also self-propel without doing any net work on the fluid in potential flow. As a body deforms, the mass of fluid displaced – the so-called added mass – depends on the instantaneous body configuration, and a net displacement is possible over a period of deformation. This potential swimming takes a form identical to that for Stokes swimmers with the configuration-dependent added mass replacing the hydrodynamic resistance tensor. Analytical insight into the swimming of a deformable body is obtained through an expansion of the nonlinear spatial dependence of the hydrodynamic interactions and connections between previous studies of swimming in Stokes flow to those in potential flow are made.

A new two-equation model for gravity-driven liquid film flow based on the long-wave expansion has been derived. The novelty of the model consists in using a base velocity profile combining parabolic (Ruyer-Quil & Manneville, Eur. Phys. J. B, vol. 15, issue 2, 2000, pp. 357–369) and ellipse (Usha et al., Phys. Fluids, vol. 32, issue 1, 2020, 013603) profile functions in the wall-normal coordinate. The dependence on a free parameter $A$ related to the eccentricity of an ellipse serves as an adjustable parameter. The resulting models are consistent at $O(\varepsilon )$ for inertia terms and at $O(\varepsilon ^2)$ for viscous diffusion effects, and predict accurately the primary instability. Appropriate tuning of the adjustable parameter helps to recover accurate predictions for the asymptotic wave celerity of nonlinear solitary waves. Further, the model is shown to capture the closed separation vortices that can form underneath the troughs of precursory capillary ripples.

When a fluid-filled cube rotating rapidly about an axis passing through two opposite vertices is subjected to harmonic modulations of its rotation rate (librations) at a modulation frequency that is $2/\sqrt {3}$ times the mean rotation frequency, all walls of the cube have critical reflection slopes. As such, all inertial wave beams emitted from edges and vertices of the cube in response to the librations are trapped in thin oscillatory boundary layers for forcing amplitudes (Rossby numbers) below a critical value which depends on the Ekman number (ratio of rotation to viscous time scales). How the resulting oscillatory boundary layer flow, referred to as a boundary-confined wave, depends on Ekman and Rossby numbers is examined in detail over several decades. Of particular interest is how the mean flow grows with increasing forcing amplitude, leading to instability resulting from nonlinear interactions between the mean flow and waves in the oscillatory boundary layers, injecting intense small-scale structures throughout the cube.

Turbulent motions enhance the diffusion of large-scale flows and temperature gradients. Such diffusion is often parameterized by coefficients of turbulent viscosity ($\nu _{t}$) and turbulent thermal diffusivity ($\chi _{t}$) that are analogous to their microscopic counterparts. We compute the turbulent diffusion coefficients by imposing sinusoidal large-scale velocity and temperature gradients on a turbulent flow and measuring the response of the system. We also confirm our results using experiments where the imposed gradients are allowed to decay. To achieve this, we use weakly compressible three-dimensional hydrodynamic simulations of isotropically forced homogeneous turbulence. We find that the turbulent viscosity and thermal diffusion, as well as their ratio the turbulent Prandtl number, $\textit {Pr}_{t} = \nu _{t}/\chi _{t}$, approach asymptotic values at sufficiently high Reynolds and Péclet numbers. We also do not find a significant dependence of $\textit {Pr}_{t}$ on the microscopic Prandtl number $\textit {Pr} = \nu /\chi$. These findings are in stark contrast to results from the $k{-}\epsilon$ model, which suggests that $\textit {Pr}_{t}$ increases monotonically with decreasing $\textit {Pr}$. The current results are relevant for the ongoing debate on, for example, the nature of the turbulent flows in the very-low-$\textit {Pr}$ regimes of stellar convection zones.

We consider nonlinear dynamical systems driven by stochastic forcing. It has been largely evidenced in the literature that the linear response of non-normal systems (e.g. fluid flows) may exhibit a large variance amplification, even in a linearly stable regime. This linear response, however, is relevant only in the limit of vanishing forcing intensity. As the intensity increases, the frequency distribution and the variance of the response may be strongly modified due to nonlinear effects. To go beyond this limitation, we propose a theoretical approach to derive an amplitude equation governing the weakly nonlinear evolution of these systems. This approach does not rely on the presence of an eigenvalue close to the neutral axis, applying instead to any sufficiently non-normal operator, and the Fourier components of the response are allowed to be arbitrarily different from any eigenmode. The methodology is outlined for a generic nonlinear dynamical system, and the application case highlights a common non-normal mechanism in hydrodynamics: convective non-normal amplification in the flow past a backward-facing step.

An isolated charge-neutral droplet in a uniform electric field experiences no net force. However, a droplet pair can move in response to field-induced dipolar and hydrodynamic interactions. If the droplets are identical, the centre of mass of the pair remains fixed. Here, we show that if the droplets have different properties, the pair experiences a net motion due to non-reciprocal electrohydrodynamic interactions. We analyse the three-dimensional droplet trajectories using asymptotic theory, assuming spherical droplets and large separations, and numerical simulations based on a boundary integral method. The dynamics can be quite intricate depending on the initial orientation of the droplets line-of-centres relative to the applied field direction. Drops tend to migrate towards a configuration with line-of-centres either parallel or perpendicular to the applied field direction, while either coming into contact or indefinitely separating. We elucidate the conditions under which these different interaction scenarios take place. Intriguingly, we find that in some cases droplets can form a pair (tandem) that translates either parallel or perpendicular to the applied field direction.

We present an experimental study of Rayleigh–Bénard convection using liquid metal alloy gallium-indium-tin as the working fluid with a Prandtl number of $Pr=0.029$. The flow state and the heat transport were measured in a Rayleigh number range of $1.2\times 10^{4} \le Ra \le 1.3\times 10^{7}$. The temperature fluctuation at the cell centre is used as a proxy for the flow state. It is found that, as $Ra$ increases from the lower end of the parameter range, the flow evolves from a convection state to an oscillation state, a chaotic state and finally a turbulent state for $Ra>10^5$. The study suggests that the large-scale circulation in the turbulent state is a residual of the cell structure near the onset of convection, which is in contrast with the case of $Pr\sim 1$, where the cell structure is transiently replaced by high order flow modes before the emergence of the large-scale circulation in the turbulent state. The evolution of the flow state is also reflected by the heat transport characterised by the Nusselt number $Nu$ and the probability density function (p.d.f.) of the temperature fluctuation at the cell centre. It is found that the effective local heat transport scaling exponent $\gamma$, i.e. $Nu\sim Ra^{\gamma }$, changes continuously from $\gamma =0.49$ at $Ra\sim 10^4$ to $\gamma =0.25$ for $Ra>10^6$. Meanwhile, the p.d.f. at the cell centre gradually evolves from a Gaussian-like shape before the transition to turbulence to an exponential-like shape in the turbulent state. For $Ra>10^6$, the flow shows self-similar behaviour, which is revealed by the universal shape of the p.d.f. of the temperature fluctuation at the cell centre and a $Nu=0.19Ra^{0.25}$ scaling for the heat transport.

The shape of a sessile drop on a horizontal substrate depends upon the Bond number $Bo$ and the contact angle $\alpha$. Inspired by puddle approximations at large $Bo$ (Quéré, Rep. Prog. Phys., vol. 68, 2005, p. 2495), we address here the limit of small contact angles at fixed drop volume and arbitrary $Bo$. It readily leads to a pancake shape approximation, where the drop height and radius scale as $\alpha$ and $\alpha ^{-1/2}$, respectively, with capillary forces being appreciable only near the edge. The pancake approximation breaks down for $Bo=\textrm {ord}(\alpha ^{2/3})$. In that distinguished limit, capillary and gravitational forces are comparable throughout, and the drop height and radius scale as $\alpha ^{2/3}$ and $\alpha ^{-1/3}$, respectively. For $Bo\ll \alpha ^{2/3}$ these scalings remain, with the drop shape turning into a spherical cap. The asymptotic results are compared with a numerical solution of the exact problem.

It has been reported that a fully localized turbulent band in channel flow becomes sustained when the Reynolds number is above a threshold. Here we show evidence that turbulent bands are of a transient nature instead. When the band length is controlled to be fixed, the lifetime of turbulent bands appears to be stochastic and exponentially distributed, a sign of a memoryless transient nature. Besides increasing with the Reynolds number, the mean lifetime also strongly increases with the band length. Given that the band length always changes over time in real channel flow, this size dependence may translate into a time dependence, which needs to be taken into account when clarifying the relationship between channel flow transition and the directed percolation universality class.

We investigate the modulation of turbulence caused by the presence of finite-size dispersed particles. Bluff (isotropic) spheres versus slender (anisotropic) fibres are considered to understand the influence of the shape of the objects on altering the carrier flow. While at a fixed mass fraction – but different Stokes number – both objects provide a similar bulk effect characterized by a large-scale energy depletion, a scale-by-scale analysis of the energy transfer reveals that the alteration of the whole spectrum is intrinsically different. For bluff objects, the classical energy cascade shrinks in its extension but is unaltered in the energy content and its typical features, while for slender ones we find an alternative energy flux which is essentially mediated by the fluid–solid coupling.

Pattern-forming with externally imposed symmetry is ubiquitous in nature but little studied. We present experimental studies of pattern formation and selection by spatial periodic forcing in rapidly rotating convection. When periodic topographic structures are constructed on the heated boundary, they modulate the local temperature and velocity fields. Symmetric convection patterns in the form of regular vortex lattices are observed near the onset of convection, when the periodicity of the external forcing is set close to the intrinsic vortex spacing. We show that the new patterns arise as a dynamical process of imperfect bifurcation which is well described by a Ginzburg–Landau-like model. We explore the phase diagram of buoyancy strength and periodicity of external forcing to find the optimal experimental settings for which the vortex patterns best match that of the external forcing.

We consider direct statistical simulation (DSS) of a paradigm system of convection interacting with mean flows. In the Busse annulus model, zonal jets are generated through the interaction of convectively driven turbulence and rotation; non-trivial dynamics including the emergence of multiple jets and bursting ‘predator–prey’ type dynamics can be found. We formulate the DSS by expanding around the mean flow in terms of equal-time cumulants and arrive at a closed set of equations of motion for the cumulants. Here, we present results using an expansion terminated at the second cumulant (CE2); it is fundamentally a quasilinear theory. We focus on particular cases including bursting and bistable multiple jets and demonstrate that CE2 can reproduce the results of direct numerical simulation if particular attention is given to symmetry considerations.

Statistically homogeneous flows obey exact kinematic relations. The Betchov homogeneity constraints (Betchov, J. Fluid Mech., vol. 1, 1956, pp. 497–504) for the average principal invariants of the velocity gradient are among the most well-known and extensively employed homogeneity relations. These homogeneity relations have far-reaching implications for the coupled dynamics of strain and vorticity, as well as for the turbulent energy cascade. Whether the Betchov homogeneity constraints are the only possible ones or whether additional homogeneity relations exist has not been proven yet. Here we show that the Betchov homogeneity constraints are the only homogeneity constraints for incompressible and statistically isotropic velocity gradient fields. Our analysis also applies to compressible/perceived velocity gradients, and it allows the derivation of homogeneity relations involving the velocity gradient and other dynamically relevant quantities, such as the pressure Hessian and viscous stresses.

Symmetry relations are derived for the added mass and damping of structures where the shape is unchanged by rotation about the vertical axis through an angle $\theta = 2 {\rm \pi}/N$ with the integer $N\ge 3$. For this type of structure, the added mass and damping for horizontal translation are the same for all directions, as in the case of axisymmetric structures. The same symmetry applies to rotations about horizontal axes. The principal application is to offshore structures and other bodies floating on the free surface or submerged, but the same symmetry relations apply more generally to unsteady body motions in an ideal fluid.

A nonlinear ensemble-variational data assimilation is performed in order to estimate the unknown flow field over a slender cone at Mach 6, from isolated wall-pressure measurements. The cost functional accounts for discrepancies in wall-pressure spectra and total intensity between the experiment and the prediction using direct numerical simulations, as well as our relative confidence in the measurements and the estimated state. We demonstrate the robustness of the predicted flow by direct propagation of posterior statistics. The approach provides a unique first look at the flow beyond the sensor data, and rigorously accounts for the role of nonlinearity, unlike previous efforts that adopted ad hoc inflow syntheses. Away from the wall, two- and three-dimensional assimilated states both show rope-like structures, qualitatively similar to independent schlieren visualizations. Despite this resemblance, and even though the planar second modes are the most unstable upstream, three-dimensional waves must be included in the assimilation in order to accurately reproduce the wall-pressure measurements recorded in the AFRL Ludwieg Tube facility. The results highlight the importance of three-dimensionality of the field and of the base-state distortion on the instability waves in this experiment, and motivate future measurements that probe the three-dimensional nature of the flow field.

We derive from first principles analytic relations for the second- and third-order moments of $\boldsymbol{\mathsf{m}}$, the spatial gradient of fluid velocity $\boldsymbol{u}$, $\boldsymbol{\mathsf{m}} = \nabla \boldsymbol{u}$, in compressible turbulence, which generalize known relations in incompressible flows. These relations, although derived for homogeneous flows, hold approximately for a mixing layer. We also discuss how to apply these relations to determine all the second- and third-order moments of the velocity gradient experimentally for isotropic compressible turbulence.

We report on the phenomenon of detonation synchronization and demonstrate the existence of the Arnold tongues and devil's staircases in the problem of gaseous detonation in a periodically inhomogeneous reactive medium. Universal properties of these dynamical structures – the Farey tree, fractal dimension and period-doubling bifurcations – are revealed.