Convection occurs ubiquitously on and in rotating geophysical and astrophysical bodies. Prior spherical shell studies have shown that the convection dynamics in polar regions can differ significantly from the lower latitude, equatorial dynamics. Yet most spherical shell convective scaling laws use globally-averaged quantities that erase latitudinal differences in the physics. Here we quantify those latitudinal differences by analysing spherical shell simulations in terms of their regionalized convective heat-transfer properties. This is done by measuring local Nusselt numbers in two specific, latitudinally separate, portions of the shell, the polar and the equatorial regions, $Nu_p$ and $Nu_e$, respectively. In rotating spherical shells, convection first sets in outside the tangent cylinder such that equatorial heat transfer dominates at small and moderate supercriticalities. We show that the buoyancy forcing, parameterized by the Rayleigh number $Ra$, must exceed the critical equatorial forcing by a factor of ${\approx }20$ to trigger polar convection within the tangent cylinder. Once triggered, $Nu_p$ increases with $Ra$ much faster than does $Nu_e$. The equatorial and polar heat fluxes then tend to become comparable at sufficiently high $Ra$. Comparisons between the polar convection data and Cartesian numerical simulations reveal quantitative agreement between the two geometries in terms of heat transfer and averaged bulk temperature gradient. This agreement indicates that rotating spherical shell convection dynamics is accessible both through spherical simulations and via reduced investigatory pathways, be they theoretical, numerical or experimental.

In Rayleigh–Bénard convection, it has been found that the amount of heat passing through the fluid has a power-law dependence on the imposed temperature difference. Modifying this dependence, either enhancing or reducing the heat transfer capability of fluids, is important in many scientific and practical applications. Here, we present a simple means to control the vertical heat transfer in Rayleigh–Bénard convection by injecting heat through one lateral side of the fluid domain and extracting the same amount of heat from the opposite side. This horizontal heat flux regulates the large-scale circulation, and increases the heat transfer rate in the vertical direction. Our numerical and theoretical studies demonstrate how a classical Rayleigh–Bénard convection responds to such a perturbation when the system is near or well above the onset of convection.

We derive a limit on energy savings in controlled channel flow. For flow in a channel driven by pressure, shear or any combination of the two, and controlled via wall transpiration or spanwise wall motion, the uncontrolled laminar state requires the least net energy (accounting for the energetic cost of control). Thus, the optimal control solution is to laminarize the flow. Additionally, we raise the possibility of beating this limit. By simultaneously applying wall transpiration and spanwise wall motion, we show that it may be possible to attain sustained sub-laminar energy expenditure in a controlled flow. We provide a necessary design criterion for net energy savings.

The role of the Lagrangian mean flow, or drift, in modulating the geometry, kinematics and dynamics of rotational and irrotational deep-water surface gravity waves is examined. A general theory for permanent progressive waves on an arbitrary vertically sheared steady Lagrangian mean flow is derived in the Lagrangian reference frame and mapped to the Eulerian frame. A Lagrangian viewpoint offers tremendous flexibility due to the particle labelling freedom and allows us to reveal how key physical wave behaviour arises from a kinematic constraint on the vorticity of the fluid, inter alia the nonlinear correction to the phase speed of irrotational finite amplitude waves, the free surface geometry and velocity in the Eulerian frame, and the connection between the Lagrangian drift and the Benjamin–Feir instability. To complement and illustrate our theory, a small laboratory experiment demonstrates how a specially tailored sheared mean flow can almost completely attenuate the Benjamin–Feir instability, in qualitative agreement with the theory. The application of these results to problems in remote sensing and ocean wave modelling is discussed. We provide an answer to a long-standing question: remote sensing techniques based on observing current-induced shifts in the wave dispersion will measure the Lagrangian, not the Eulerian, mean current.

Active elements in active nematics can impose forces on immersed bodies and move them accordingly. We numerically investigate the vibrational motion of a cantilever beam placed in active nematics. The continuous energy transfer from vortices to the beam results in beam oscillation, whose direction and amplitude depend on the vortex strength, size and position. Referring to the kinetic-energy spectrum, we indicate that both the large- and small-scale vortices are the primary mechanism for the energy transfer between the fluid and beam, leading to the beam oscillatory motion, with the contribution from the large-scale vortices being higher. We investigate the effect of fluid properties such as activity, viscosity and elastic constant on the oscillation frequency. We show that the intensification of the activity increases peak frequency, and there is a linear correlation between the peak frequency and activity. We further demonstrate the reciprocal relationship between viscosity and peak frequency. Subsequently, we relate the increase and decrease in the peak frequency to the energy injection/dissipation by activity/viscosity. Moreover, we reveal the negligibly small dependency of beam peak frequency on the elastic constant and discuss free energy's role in accounting for this behaviour. The findings clearly demonstrate that active fluids can impose an oscillatory motion on flexible bodies, which might be used as a novel method for measuring the critical properties of active nematics.

Shock-containing supersonic jets undergoing resonance processes are challenging from both a measurement and simulation perspective. These jets are host to a broad range of complex fluid phenomena: intense acoustic waves, turbulence, wavepackets and strong shock waves. Strong shocks present a challenge to both the experimental and numerical researcher. In the paper of Léon et al. (J. Fluid Mech., vol. 947, 2022, A36), a novel optical technique based on multi-axis digital holographic interferometry is applied to the study of a highly underexpanded screeching jet, producing density measurements of unprecedented clarity and resolution. Where prior studies have been restricted to extrapolating the three-dimensional field from two-dimensional slices or projections, in this work the authors directly measure the three-dimensional helical structure of the wavepacket associated with jet screech.

We derive and investigate numerically a reduced model for wave–vortex interactions involving non-dispersive waves, which we study in a two-dimensional shallow water system with an eye towards applications in atmosphere–ocean fluid dynamics. The model consists of a coupled set of nonlinear partial differential equations for the Lagrangian-mean velocity and the wave-related pseudomomentum vector field defined in generalized Lagrangian-mean theory. It allows for two-way interactions between the waves and the balanced flow that is controlled by the distribution of Lagrangian-mean potential vorticity, and for strong solutions it features a desirable exact energy conservation law for the sum of wave energy and mean flow energy. Our model goes beyond standard ray tracing as we can derive weak solutions that contain discontinuities in the pseudomomentum field, using the theory of weakly hyperbolic systems. This allows caustics to form without predicting infinite wave amplitudes, as would be the case in the standard ray-tracing theory. Suitable wave forcing and dissipation terms are added to the model and a numerical scheme for the model is implemented as a coupled set of pseudo-spectral and finite-volume integrators. Idealized examples of interactions between wavepackets and simple vortex structures are presented to illustrate the model dynamics. The unforced and non-dissipative simulations suggest a heuristic rule of ‘greedy’ waves, i.e. in the long run the wave field always extracts energy from the mean flow.

Swirling waves frequently occur in a three-dimensional tank under a nearly resonant excitation with oblique excitation angle. The oblique excitation produces two excitation components and the secondary component triggers rotational waves. The detailed mechanism of the switch in circular direction is clearly described in this study. The evolution of the hydrodynamic forces (Fz) on the tank walls presents a beating phenomenon and the switch direction always occurs at the peak and trough of the envelope of the Fz history. The external excitation moment changes the circulation intensity, and finally the swirling wave becomes a single-direction wave, but only for a short time. The profile of the single-direction wave was determined by the dominant sloshing mode and the instantaneous tank motion was found to be a key factor to determine the consequent swirling flow circulating direction.

We present the results of a set of experiments designed to measure the dispersion of non-spherical particles in a wave–current flow. We released negatively buoyant discs, rods and unit-aspect-ratio cylinders into a flow both with and without waves and analysed their respective landing positions to quantify how much they had dispersed while in the flow. We found that the presence of waves significantly increased the dispersion of the particles, and that the magnitude of this increase depends on particle shape and volume. In particular, thinner rods and thinner discs have greater relative dispersion than thicker rods and thicker discs, respectively, and smaller particles have greater relative dispersion than larger particles. Although the particles travelled farther in the presence of waves, the increase in dispersion was much greater than could be explained solely by increased transport distance. These results indicate that models of microplastic transport must account for waves as well as particle characteristics.

The trajectory of surface gravity waves in the potential flow regime is affected by the gravitational acceleration, water density and sea bed depth. Although the gravitational acceleration and water density are approximately constant, the effect of water depth on surface gravity waves exponentially decreases as the water depth increases. In shallow water, cloaking an object from surface waves by varying the sea bed topography is possible, however, as the water depth increases, cloaking becomes a challenge because there is no physical parameter to be engineered and subsequently affects the wave propagation. In order to create an omnidirectional cylindrical cloaking device for finite-depth/deep-water waves, we propose an elastic composite plate that floats on the surface around a to-be-cloaked cylinder. The composite plate is made of axisymmetric, homogeneous and isotropic annular thin rings which provide adjustable degrees of freedom to engineer and affect the wave propagation. We first develop a pseudo-spectral method to efficiently determine the wave solution for a floating composite plate. Next, we optimise the physical parameters of the plate (i.e. flexural rigidity and mass of every ring) using an evolutionary algorithm to minimise the energy of scattered waves from the object and therefore cloak the inner cylinder from incident waves. We show that the optimised cloak reduces the energy of scattered waves as high as 99 % for the target wave number. We quantify the effectiveness of our cloak with different parameters of the plate and show that varying the flexural rigidity is essential to control wave propagation and the cloaking structure needs to be at least made of four rings with a radius of at least three times of the cloaked region. We quantify the wave drift force exerted on the structures and show that the optimised plate reduces the exerted force by 99.9 %. The proposed cloak, due to its structural simplicity and effectiveness in reducing the wave drift force, may have potential applications in cloaking offshore structures from water waves.

Motivated by modelling rotating turbulence in planetary fluid layers, we investigate precession-driven flows in ellipsoids subject to stress-free boundary conditions (SF-BC). The SF-BC could indeed unlock numerical constraints associated with the no-slip boundary conditions (NS-BC), but are also relevant for some astrophysical applications. Although SF-BC have been employed in the pioneering work of Lorenzani & Tilgner (J. Fluid Mech., vol. 492, 2003, pp. 363–379), they have scarcely been used due to the discovery of some specific mathematical issues associated with angular momentum conservation. We revisit the problem using asymptotic analysis in the low-viscosity regime, which is validated with numerical simulations. First, we extend the reduced model of uniform-vorticity flows in ellipsoids to account for SF-BC. We show that the long-term evolution of angular momentum is affected by viscosity in triaxial geometries, but also in axisymmetric ellipsoids when the mean rotation axis of the fluid is not the symmetry axis. In a regime relevant to planets, we analytically obtain the primary forced flow in triaxial geometries, which exhibits a second inviscid resonance. Then, we investigate the bulk instabilities existing in precessing ellipsoids. We show that using SF-BC would be useful to explore the non-viscous instabilities (e.g. Kerswell, Geophys. Astrophys. Fluid Dyn., vol. 72, 1993, pp. 107–144), which are presumably relevant for planetary applications but are often hampered in experiments or simulations with NS-BC.

The supersonic jet over a convex wall is numerically investigated using the delayed detached-eddy simulation method based on the two-equation shear-stress transport model. The current study focuses on instabilities, turbulent statistics and the influence of compressibility effects. A widely applicable data-driven modal decomposition approach, called dynamic mode decomposition is used to gain further insight into the dynamical behaviours of the flow. The results demonstrate that streamwise vortices caused by the centrifugal force play significant roles in shear layer instabilities. The spanwise modulation of the streamwise vortices induces inflection points in the flow, resulting in secondary shear layer instability. This instability, which is sustained by the side-to-side sway of the streamwise vortices to obtain energy from the mean flow, dominates the rapid growth of the shear layer and turbulent stresses in the growth region. In the self-similar region, there is not only self-similarity of velocity profiles, but also self-similarity of normalized turbulent stresses. The compressibility effect significantly inhibits the growth of the shear layer and the formation of large-scale streamwise vortices. The investigation of turbulent stresses in the self-similar region with increasing convective Mach number indicates that the compressibility effect enhances turbulence anisotropy.

An elastically mounted circular cylinder, immersed in a cross-current and free to move along a rectilinear path, is subjected to vortex-induced vibrations (VIV). These vibrations develop through a mechanism referred to as lock-in, where body motion and vortex shedding synchronize at a frequency that may deviate both from the oscillator natural frequency and from the vortex shedding frequency past a fixed cylinder. The present numerical study aims at extending the analysis to curved trajectories, by considering that the cylinder is free to translate along a circular path. The Reynolds number based on the body diameter ($D$) and current velocity ($U$) is set to $100$. A wide range of path radii, from $0.05D$ to $10D$, and values of the reduced velocity (inverse of the oscillator natural frequency non-dimensionalized by $D$ and $U$) up to $30$ are examined, for the concave and convex configurations, i.e. the circular path centre located upstream or downstream of the cylinder. Path curvature results in a major alteration of the flow–body system behaviour compared with rectilinear VIV, with substantially different evolutions in the concave and convex configurations. In addition to the typical lock-in mechanism, two subharmonic forms of synchronization, at half and one third of vortex formation frequency, are uncovered in the convex configuration. They coexist with a desynchronized regime where the body and the flow oscillate at incommensurable frequencies. The four interaction regimes exhibit contrasted trends in terms of structural response, spatiotemporal organization of the wake and associated forces. They particularly differ by their symmetry properties, which are closely linked to the possible reconfiguration of the oscillator due to mean fluid forcing.

This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction ${{\varPhi _b}}$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, ${{\varPhi _b}}= 0.2$, the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses.

Screech tones in supersonic jets are underpinned by resonance between downstream-travelling Kelvin–Helmholtz waves and upstream-travelling acoustic waves. Specifically, recent works suggest that the relevant acoustic waves are guided within the jet and are described by a discrete mode of the linearised Euler equations. However, the reflection mechanism that converts downstream-travelling waves into upstream-travelling waves, and vice versa, has not been thoroughly addressed, leading to missing physics within most resonance models. In this work, we investigate the reflection and transmission of waves generated by the interaction between a Kelvin–Helmholtz wave and a normal shock in an under-expanded jet using a mode-matching approach. Both vortex-sheet and finite-thickness shear-layer models are explored, quantifying the impact of the shear layer in the reflection process. This approach could enable more quantitative predictions of resonance phenomena in jets and other fluid systems.

Data-driven dimensionality reduction methods such as proper orthogonal decomposition and dynamic mode decomposition have proven to be useful for exploring complex phenomena within fluid dynamics and beyond. A well-known challenge for these techniques is posed by the continuous symmetries, e.g. translations and rotations, of the system under consideration, as drifts in the data dominate the modal expansions without providing an insight into the dynamics of the problem. In the present study, we address this issue for fluid flows in rectangular channels by formulating a continuous symmetry reduction method that eliminates the translations in the streamwise and spanwise directions simultaneously. We demonstrate our method by computing the symmetry-reduced dynamic mode decomposition (SRDMD) of sliding windows of data obtained from the transitional plane-Couette and turbulent plane-Poiseuille flow simulations. In the former setting, SRDMD captures the dynamics in the vicinity of the invariant solutions with translation symmetries, i.e. travelling waves and relative periodic orbits, whereas in the latter, our calculations reveal episodes of turbulent time evolution that can be approximated by a low-dimensional linear expansion.

Understanding the statistics of bedload particle motions is of great importance. To model the hop events which are defined as trajectories of particles moving successively from the start to the end of their motions, recently, Wu et al. (Water Resour. Res., vol. 56, 2020, p. e2019WR025116) have successfully performed individual-based simulations according to the Fokker–Planck equation for particle velocities. However, analytical solutions are still not available due to (i) difficulties in treating the velocity-dependent diffusivity, and (ii) a knowledge gap in incorporating the termination of particle motions for the equation. To tackle the above-mentioned challenges, we first specify a Robin boundary condition representing the deposition of particles. Second, for analytical solutions of hop statistics, a variable transformation is devised to deal with the velocity-dependent diffusivity. The original bedload transport problem is thus found to be governed by the classic equation for the solute transport in tube flows with a constant diffusivity after the transformation. Finally, through solving the spatial and temporal moments of the governing equation, we investigate the influence of the deposition rate on three key characteristics of particle hops. Importantly, we have related the deposition rate to the mean travel times and hop distances, enabling a direct determination of this physical parameter based on measured particle motion statistics. The analytical solutions are validated by experimental observations with different bedload particle diameters and transport conditions. Based on the limited experimental datasets, the deposition frequency is shown to decrease as the shear stress increases when the flow rate is not small.

Propagation of fluid-filled fractures by fluid buoyancy is important in a variety of settings, from magmatic dykes and veins to water-filled crevasses in glaciers. Industrial hydro-fracturing utilises fluid-driven fractures to increase the permeability of rock formations, but few studies have quantified the effect of buoyancy on fracture pathways in this context. Analytical approximations for the buoyant ascent rate facilitate observation-based inference of buoyant effects in natural and engineered systems. Such analysis exists for two-dimensional fractures, but real fractures are three-dimensional (3-D). Here we present novel analysis to predict the buoyant ascent speed of 3-D fractures containing a fixed-volume batch of fluid. We provide two estimates of the ascent rate: an upper limit applicable at early time, and an asymptotic estimate (proportional to $t^{-2/3}$) describing how the speed decays at late time. We infer and verify these predictions by comparison with numerical experiments across a range of scales and analogue experiments on liquid oil in solid gelatine. We find the ascent speed is a function of the fluid volume, density, viscosity and the elastic parameters of the host medium. Our approximate solutions predict the ascent rate of fluid-driven fractures across a broad parameter space, including cases of water injection in shale and magmatic dykes. Our results demonstrate that in the absence of barriers or fluid loss, both dykes and industrial hydro-fractures can ascend by buoyancy over a kilometre within a day. We infer that barriers and fluid loss must cause the arrest of ascending fractures in industrial settings.

An initially perturbed interface between two fluids of different densities is usually unstable when driven by an acceleration or a shock wave; it is known as a Rayleigh–Taylor instability or a Richtmyer–Meshkov instability. One of the most significant issues in these instabilities is the spatiotemporal development of fingers generated at the interface, which plays an important role in both scientific research (e.g. supernova explosion) and engineering applications (e.g. inertial confinement fusion). Accurate theoretical solution of these interfacial fingers remains as an unsolved and challenging problem since Taylor's seminal work more than seven decades ago. This paper reports a unified theory established for such phenomena by combining the classical potential-flow theory and a dual-source model to address the long-standing difficulty highlighted by the initial-value sensitivity and strong nonlinearity. It is the first time for a theory to accurately predict the long-time developments in both growth rate and shape curvature of interfacial fingers at all density ratios in two and three dimensions. Moreover, the new theory clearly reveals the nonlinear coupling mechanism for interfacial evolution, and especially explains the origin of overshot in the growth rate curve.

Conduits generated by the buoyant dynamics between two miscible Stokes fluids with high viscosity contrast, a type of core–annular flow, exhibit a rich nonlinear wave dynamics. However, little is known about the fundamental wave dispersion properties of the medium. In the present work, a pump is used to inject a time-periodic flow that results in the excitation of propagating small- and large-amplitude periodic travelling waves along the conduit interface. This wavemaker problem is used as a means to measure the linear and nonlinear dispersion relations and corresponding periodic travelling wave profiles. Measurements are favourably compared with predictions from a fully nonlinear, long-wave model (the conduit equation) and the analytically computed linear dispersion relation for two-Stokes flow. A critical frequency is observed, marking the threshold between propagating and non-propagating (spatially decaying) waves. Measurements of wave profiles and the wavenumber–frequency dispersion relation quantitatively agree with wave solutions of the conduit equation. An upshift from the conduit equation's predicted critical frequency is observed and is explained by incorporating a weak recirculating flow into the full two-Stokes flow model. When the boundary condition corresponds to the temporal profile of a nonlinear periodic travelling wave solution of the conduit equation, weakly nonlinear and strongly nonlinear, cnoidal-type waves are observed that quantitatively agree with the conduit nonlinear dispersion relation and wave profiles. This wavemaker problem is an important precursor to the experimental investigation of more general boundary value problems in viscous fluid conduit nonlinear wave dynamics.