In the gap region of tandem cylinders, within the reattachment regime, bi-stability is seen to be cellular. Direct numerical simulations at a Reynolds number of 500, and a gap ratio of 3.0, show that shedding of gap vortices occurs in spanwise cells, with lengths between 0.3 and 2.7 cylinder diameters. These unstable vortex cells tend to be asymmetric with respect to the gap centreline, so that vortices are repeatedly shed from just one gap shear layer for several periods. The unstable cells appear within the basic spanwise cell structure dictated by the three-dimensional instability, and their cell lengths do not exceed those of the basic cells. Unstable cells intermittently become spanwise coherent, and this leads to a significant increase in drag amplitude. The mode change in the gap is associated with low-frequency variation of the reattachment and separation points on the downstream cylinder, causing low-frequency modulation of the vortex formation length.

The instability and vortex shedding in the bottom boundary layer caused by internal solitary waves of depression propagating along a shallow pycnocline of a fluid are computed by finite-volume code in two dimensions. The calculated transition to instability agrees very well with laboratory experiments (Carr et al., Phys. Fluids, vol. 20, issue 6, 2008, 06603) but disagrees with existing computations that give a very conservative instability threshold. The instability boundary expressed by the amplitude depends on the depth $d$ of the pycnocline divided by the water depth $H$, and decays by a factor of 2.2 when $d/H$ is 0.21, and by a factor of 1.6 when $d/H$ is 0.16, and the stratification Reynolds number increases by a factor of 32. The instability occurs at moderate amplitude at large scale. The calculated oscillatory bed shear stress is strong in the wave phase and increases with the scale. Its non-dimensional magnitude at stratification Reynolds number 650 000 is comparable to the turbulent stress that can be extracted from field measurements of internal solitary waves of similar nonlinearity, moving along a pycnocline of similar relative depth.

Resonant sloshing in circular cylinders was studied by Faltinsen et al. (J. Fluid Mech., vol. 804, 2016, pp. 608–645), whose theory was used to describe steady-state resonant waves due to a time-harmonic container's elliptic orbits. In the limit of longitudinal container motions, a symmetry breaking of the planar wave solution occurs, with clockwise and anti-clockwise swirling equally likely. In addition to this primary harmonic dynamics, previous experiments have unveiled that diverse super-harmonic dynamics are observable far from primary resonances. Among these, the so-called double-crest (DC) dynamics, first observed by Reclari et al. (Phys. Fluids, vol. 26, no. 5, 2014, 052104) for circular sloshing, is particularly relevant, as its manifestation is the most favoured by the spatial structure of the external driving. Following Bongarzone et al. (J. Fluid Mech., vol. 943, 2022, A28), in this work we develop a weakly nonlinear analysis to describe the system response to super-harmonic longitudinal forcing. The resulting system of amplitude equations predicts that a planar wave symmetry breaking via stable swirling may also occur under super-harmonic excitation. This finding is confirmed by our experimental observations, which identify three possible super-harmonic regimes, i.e. (i) stable planar DC waves, (ii) irregular motion and (iii) stable swirling DC waves, whose corresponding stability boundaries in the forcing frequency-amplitude plane quantitatively match the present theoretical estimates.

Bathymetric changes have been experimentally shown to affect the occurrence of rogue waves. We recently derived a non-homogeneous correction to the spectral analysis, allowing us to describe the evolution of the rogue wave probability over a shoal. Here, we extend this work to the evolution of the excess kurtosis of the surface elevation, that plays a central role in estimating rare event probabilities. Furthermore, we provide an upper bound to the excess kurtosis. In intermediate and deep water regimes, a shoal does not affect wave steepness nor bandwidth significantly, so that the vertical asymmetry between crests and troughs, the excess kurtosis and the exceedance probability of wave height stay rather constant. In contrast, in shallower water, a sharp increase in wave steepness increases the vertical asymmetry, resulting in a growth of both the tail of the exceedance probability and the excess kurtosis.

It is shown that thermal waves applied on the bounding surface of a horizontal slot generate a pumping effect. Reynolds stress developed by the change in the flow field due to the thermal inertia associated with the penetration of the wave into the fluid interior drives the net fluid movement in the horizontal direction. The induced flow rate increases with the wave speed, but excessive wave speed reduces it as convection becomes limited only in the near-wall area. The excessive increase of the wavelength and its excessive decrease reduce the flow rate. An increase in the wave amplitude increases the flow rate with a saturation limiting its growth. Judicious selection of the wave speed and wavelength of a thermal wave provides a means for optimizing the wave-induced pumping. The magnitude of the pumping increases with a reduction of Prandtl number.

We analyse a class of stochastic advection problems by conditionally averaging the passive tracer equation with respect to a given flow state. In doing so, we obtain expressions for the turbulent diffusivity as a function of the flow statistics spectrum. When flow statistics are given by a continuous-time Markov process with a finite state space, calculations are amenable to analytic treatment. When the flow statistics are more complex, we show how to approximate turbulent fluxes as hierarchies of finite state space continuous-time Markov processes. The ensemble average turbulent flux is expressed as a linear operator that acts on the ensemble average of the tracer. We recover the classical estimate of turbulent flux as a diffusivity tensor, the components of which are the integrated autocorrelation of the velocity field in the limit that the operator becomes local in space and time.

Wall-bounded turbulent shear flows are known to exhibit universal small-scale dynamics that are modulated by large-scale flow structures. Strong pressure gradients complicate this characterization, however. They can cause significant variation of the mean flow in the streamwise direction. For such situations, we perform asymptotic analysis of the Navier–Stokes equations to inform a model for the effect of mean flow growth on near-wall turbulence in a small domain localized to the boundary. The asymptotics are valid whenever the viscous length scale is small relative to the length scale over which the mean flow varies. To ensure the correct momentum environment, a dynamic procedure is introduced that accounts for the additional sources of mean momentum flux through the upper domain boundary arising from the asymptotic terms. Comparisons of the model's low-order, single-point statistics with those from direct numerical simulation and well-resolved large eddy simulation of adverse-pressure-gradient turbulent boundary layers indicate the asymptotic model successfully accounts for the effect of boundary layer growth on the small-scale near-wall turbulence.

The presence of thermocapillary (Marangoni) convection in microgravity may help to enhance the heat transfer rate of phase change materials (PCMs) in space applications. We present a three-dimensional numerical investigation of the nonlinear dynamics of a melting PCM placed in a cylindrical container filled with n-octadecane and surrounded by passive air. The heat exchange between the PCM and ambient air is characterized in terms of the Biot number, when the air temperature has a linear profile. The effect of thermocapillary convection on heat transfer and the topology of the melting front is studied by varying the applied temperature difference between the circular supports and the heat transfer through the interface. The evolution of Marangoni convection during the PCM melting leads to the appearance of hydrothermal instabilities. A new mathematical approach for the nonlinear analysis of emerging hydrothermal waves (HTWs) is suggested. Being applied for the first time to the examination of PCMs, this procedure allows us to explore the nature of the coupling between HTWs and heat gain/loss through the interface, and how it changes over time. We observe a variety of dynamics, including standing and travelling waves, and determine their dominant and secondary azimuthal wavenumbers. Coexistence of multiple travelling waves with different wavenumbers, rotating in the same or opposite directions, is among the most fascinating observations.