The collision of a gravitationally driven horizontal current with a barrier following release from a confining lock is investigated using a shallow water model of the motion, together with a sophisticated boundary condition capturing the local interaction. The boundary condition permits several overtopping modes: supercritical, subcritical and blocked flow. The model is analysed both mathematically and numerically to reveal aspects of the unsteady motion and to compute the proportion of the fluid trapped upstream of the barrier. Several problems are treated. Firstly, the idealised problem of a uniform incident current is analysed to classify the unsteady dynamical regimes. Then, the extreme regimes of a very close or distant barrier are tackled, showing the progression of the interaction through the overtopping modes. Next, the trapped volume of fluid at late times is investigated numerically, demonstrating regimes in which the volume is determined purely by volumetric considerations, and others where transient inertial effects are significant. For a Boussinesq gravity current, even when the volume of the confined region behind the barrier is equal to the fluid volume, $30\,\%$ of the fluid escapes the domain, and a confined volume three times larger is required for the overtopped volume to be negligible. For a subaerial dam-break flow, the proportion that escapes is in excess of $60\,\%$ when the confined volume equals the fluid volume, and a barrier as tall as the initial release is required for negligible overtopping. Finally, we compare our predictions with experiments, showing a good agreement across a range of parameters.

The sedimentation of two spherical solid objects in a viscous fluid has been extensively investigated and well understood. However, a pair of flat disks (in three dimensions) settling in the fluid shows more complex hydrodynamic behaviour. The present work aims to improve the understanding of this phenomenon by performing direct numerical simulation and physical experiments. The present results show that the sedimentation processes are significantly influenced by disk shape, characterized by a dimensionless moment of inertia I*, and Reynolds number Re of the leading disk. For the flatter disks with smaller I*, steady falling with enduring contact transits to periodic swinging with intermittent contacts as Re increases. The disks with larger I* tend to fall in a drafting-kissing-tumbling mode at low Re and to remain separated at high Re. Based on I* and Re, a phase diagram is created to classify the two-disk falling into ten distinctive patterns. The planar motion or three-dimensional motion of the disks is determined primarily by Re. Turbulent disturbance flows at a high Re contribute to the chaotic three-dimensional rotation of the disks. The chance for the two disks to contact is increased when I* and Re are reduced.

This paper investigates heat transport in penetrative convection with a marginally stable temporal-horizontal-averaged field or background field. Assuming that the background field is steady and is stabilised by the nonlinear perturbation terms, we obtain an eigenvalue problem with an unknown background temperature $\tau$ by truncating the nonlinear terms. Using a piecewise profile for $\tau$, we derived an analytical scaling law for heat transport in penetrative convection as $Ra\rightarrow \infty$: $Nu=(1/8)(1-T_M)^{5/3}Ra^{1/3}$ ($Nu$ is the Nusselt number; $Ra$ is the Rayleigh number and $T_M$ corresponds to the temperature at which the density is maximal). A conditional lower bound on $Nu$, under the marginal stability assumption, is then derived from a variational problem. All the solutions to the full system should deliver a higher heat flux than the lower bound if they satisfy the marginal stability assumption. However, data from the present direct numerical simulations and previous optimal steady solutions by Ding & Wu (J. Fluid Mech., vol. 920, 2021, A48) exhibit smaller $Nu$ than the lower bound at large $Ra$, indicating that these averaged fields are over-stabilised by the nonlinear terms. To incorporate a more physically plausible constraint to bound heat transport, an alternative approach, i.e. the quasilinear approach is invoked which delivers the highest heat transport and agrees well with Veronis's assumption, i.e. $Nu\sim Ra^{1/3}$ (Astrophys. J., vol. 137, 1963, p. 641). Interestingly, the background temperature $\tau$ yielded by the quasilinear approach can be non-unique when instability is subcritical.

Measurements of fluctuating wall pressure in a high-Reynolds-number flow over a body of revolution are described. With a strong axial pressure gradient and moderate lateral curvature, this non-equilibrium flow is relevant to marine applications as well as short-haul urban transportation. The wall-pressure spectrum and its scaling are discussed, along with its relation to the space–time structure. As the flow decelerates downstream, the root-mean-square level of the pressure drops together with the wall shear stress ($\tau _w$) and is consistently approximately 7$\tau _w$. While the associated dimensional spectra see a broadband reduction of over 15 dB per Hz, they appear to attain a single functional form, collapsing to within 2 dB when normalized with the wall-wake scaling where $\tau _w$ is the pressure scale and $U_e/\delta$ is the frequency scale. Here, $\delta$ is the boundary layer thickness and $U_e$ is the local free-stream velocity. The general success of the wall-wake scaling, including in the viscous $f^{-5}$ region, suggests that the large-scale motions in the outer layer play a predominant role in the near-wall turbulence and wall pressure. On investigating further, we find that the instantaneous wall-pressure fluctuations are characterized by a quasi-periodic feature that appears to convect downstream at speeds consistent with the outer peak in the turbulence stresses. The conditional structure of this feature, estimated through peak detection in the time series, resembles that of a roller, supporting the embedded shear layer hypothesis (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642; Balantrapu et al., J. Fluid Mech., vol. 929, 2021, A9). Therefore, the outer-region shear-layer-type motions may be important when devising strategies for flow control, drag and noise reduction for decelerating boundary layers.

We investigate the role of the laminar/turbulent interface in the interscale energy transfer in a boundary layer undergoing bypass transition with the aid of the Kármán–Howarth–Monin–Hill (KHMH) equation. A local binary indicator function is used to detect the interface and employed subsequently to define two-point intermittencies. These are used to decompose the standard-averaged interscale and interspace energy fluxes into conditionally averaged components. We find that the inverse cascade in the streamwise direction reported in an earlier work arises due to events across the downstream or upstream interfaces of a turbulent spot. However, the three-dimensional energy flux maps reveal significant differences between these two regions: in the downstream interface, inverse cascade is stronger and dominant over a larger range of streamwise and spanwise separations. We explain this finding by considering a propagating spot of simplified shape as it crosses a fixed streamwise location. We derive also the conditionally averaged KHMH equation, thus generalising similar equations for single-point statistics to two-point statistics. We compare the three-dimensional maps of the conditionally averaged production and total energy flux within turbulent spots against the maps of standard-averaged quantities within the fully turbulent region. The results indicate remarkable dynamical similarities between turbulent spots and the fully turbulent region for two-point statistics. This has been known only for single-point quantities, and we demonstrate here that the similarity extends to two-point quantities as well.

It has been reported that, in submarine landslides down a slope, only a small part of the landslide momentum was transferred to the tsunami while most was lost in the further propagation over the flat bed. The high momentum of landslides off a slope meant a long run-out and a great impact on underwater infrastructure. Nevertheless, little attention has been paid to the variation of momentum of deformable submarine landslides off a slope, i.e. the loss of momentum when the slides flow away from the slope over a flat bottom. In this paper, the translational momentum of deformable submarine granular landslides running down a non-erodible inclined bed is investigated with a two-phase smoothed particle hydrodynamics model. After flowing down the slope, the transport rate and flux of the landslide translational momentum along the propagation over the flat bottom are examined. The effects of physical variables of the slide, particularly the grain size, the initial compaction and the front intrusion angle on the variation of the translational momentum, are explored. Accordingly, scaling relations of the spatial-temporal maximum transport rate and flux of the landslide translational momentum, as well as those of the slide final run-out and the generated leading wave height, are proposed. These scaling relations, although based on numerical data of small-scale granular landslides, demonstrate a preliminary attempt to develop practical expressions in the framework of momentum for estimating the run-out of real submarine landslides and the impact on underwater infrastructure.

Effects of surface tension reduction on wind-wave growth are investigated using direct numerical simulations of air–water two-phase turbulent flows. The incompressible Navier–Stokes equations for air and water sides are solved using an arbitrary Lagrangian–Eulerian method with boundary-fitted moving grids. The wave growth of finite-amplitude and non-breaking gravity–capillary waves, whose wavelength is less than 0.07 m, is simulated for two cases of different surface tensions under a low-wind-speed condition of several metres per second. The results show that significant wave height for the smaller surface tension case increases faster than that for the larger surface tension case. Energy fluxes for gravity and capillary wave scales reveal that when the surface tension is reduced, the energy transfer from the significant gravity waves to capillary waves decreases, and the significant waves accumulate more energy supplied by wind. This results in faster wave growth for the smaller surface tension case. The effect on the scalar transfer across the air–water interface is also investigated. The results show that the scalar transfer coefficient on the water side decreases due to the surface tension reduction. The decrease is caused by suppression of turbulence in the water side. In order to support the conjecture, the surface tension effect is compared with laboratory experiments in a small wind-wave tank.