This work investigates the spatio-temporal evolution of coherent structures in the wake of a generic high-speed train, based on a three-dimensional database from large eddy simulation. Spectral proper orthogonal decomposition (SPOD) is used to extract energy spectra and energy ranked empirical modes for both symmetric and antisymmetric components of the fluctuating flow field. The spectrum of the symmetric component shows overall higher energy and more pronounced low-rank behaviour compared with the antisymmetric one. The most dominant symmetric mode features periodic vortex shedding in the near wake, and wave-like structures with constant streamwise wavenumber in the far wake. The mode bispectrum further reveals the dominant role of self-interaction of the symmetric component, leading to first harmonic and subharmonic triads of the fundamental frequency, with remarkable deformation of the mean field. Then, the stability of the three-dimensional wake flow is analysed based on two-dimensional local linear stability analysis combined with a non-parallelism approximation approach. Temporal stability analysis is first performed for both the near-wake and the far-wake regions, showing a more unstable condition in the near-wake region. The absolute frequency of the near-wake eigenmode is determined based on spatio-temporal analysis, then tracked along the streamwise direction to find out the global mode growth rate and frequency, which indicate a marginally stable global mode oscillating at a frequency very close to the most dominant SPOD mode. The global mode wavemaker is then located, and the structural sensitivity is calculated based on the direct and adjoint modes derived from a local spatial analysis, with the maximum value localized within the recirculation region close to the train tail. Finally, the global mode shape is computed by tracking the most spatially unstable eigenmode in the far wake, and the alignment with the SPOD mode is computed as a function of streamwise location. By combining data-driven and theoretical approaches, the mechanisms of coherent structures in complex wake flows are well identified and isolated.

This work presents models for the behaviour of both upstream- and downstream-travelling waves in screeching elliptical jets. Proper orthogonal decomposition is performed on experimental velocity data in both the major and minor axis planes, for an aspect ratio $AR=2$ converging elliptical jet operating at nozzle pressure ratios of $2.6$ and $3.4$. From this decomposition, the radial and axial structure of the guided-jet mode (GJM) and the Kelvin–Helmholtz instability are educed. Linear-stability analysis (LSA) is performed using both the experimentally obtained mean flow, and one obtained using Reynolds-averaged Navier–Stokes (RANS) at matched conditions. It is shown that the wavenumber predicted by LSA for both waves are within the range of experimentally observed wavenumbers. Furthermore, the model accurately predicts the structure of these waves at multiple axial locations, using either the experimental or RANS mean flow. Most critically, it is demonstrated that the GJM is only predicted to be neutrally stable at the screech frequency for a relatively limited streamwise domain, the size and location of which is dependent on the nozzle pressure ratio. A comparison with the amplitude envelope for the GJM extracted from the experimental measurements indicates that the maximum fluctuations associated with the GJM are collocated with this region of the flow that is predicted to support the GJM. While there have been extensive discussions about the frequency dependence of the GJM, this is the first demonstration that its existence is highly dependent on streamwise position within the flow.

Recognising that surfactants can impede the drag reduction resulting from superhydrophobic surfaces (SHS), we investigate the impact of spatio–temporal fluctuations in surfactant concentration on the drag-reduction properties of SHS. We model the unsteady transport of soluble surfactant in a channel flow bounded by two SHS. The flow is laminar, pressure driven and the SHS are periodic in the streamwise and spanwise directions. We assume that the channel length is much longer than the streamwise period, the streamwise period is much longer than the channel height and spanwise period, and bulk diffusion is sufficiently strong for cross-channel concentration gradients to be small. By combining long-wave and homogenisation theories, we derive an unsteady advection–diffusion equation for surfactant-flux transport over the length of the channel, which is coupled to a quasi-steady advection–diffusion equation for surfactant transport over individual plastrons. As diffusion over the length of the channel is typically small, the surfactant flux is governed by a nonlinear wave equation. In the fundamental case of the transport of a bolus of surfactant, we predict its propagation speed and describe its nonlinear evolution via interaction with the SHS. The propagation speed can fall below the average streamwise velocity as the surfactant adsorbs and rigidifies the plastrons. Smaller concentrations of surfactant are advected faster than larger ones, so that wave-steepening effects can lead to shock formation in the surfactant-flux distribution. Our asymptotic results reveal how unsteady surfactant transport can affect the spatio–temporal evolution of the slip velocity, drag reduction and effective slip length in SHS channels.

Reynolds-averaged models for solving the Navier–Stokes equations are implicitly based on Kolmogorov's theory for describing energy transfers between the different turbulent scales, which means that all the energy produced at large scales is transferred at a constant rate to the smallest turbulent scales where it is dissipated. As a result, these models use a single scale to describe the turbulence spectrum, which in cases of non-equilibrium turbulence does not provide an adequate description of the transfers actually observed. This is particularly the case for wall-bounded flows at high Reynolds numbers, such as turbulent channel flows. Taking up an approach developed by Schiestel (2007 Modeling and Simulation of Turbulent Flows, ISTE Ltd and John Wiley & Sons), which aims to define a Reynolds-averaged Navier–Stokes model transporting several scales of turbulence, a two-scale Reynolds stress model (RSM) was developed in order to take into account the interactions between the inner and outer regions of wall-bounded flows. The results obtained with the model are compared with the direct numerical simulations (DNS) of Lee & Moser (J. Fluid Mech., vol. 860, 2019, pp. 886–938) in a turbulent channel for several friction Reynolds numbers up to $Re_{\tau }=5200$, for which partial integrations in spectral space were carried out, highlighting distinct behaviours between small and large scales of turbulence. The model developed provides an accurate description of the contributions at small and large scales and thus reproduces the high-Reynolds-number effects observed in DNS data. In addition, comparisons with the DNS data served to validate a large part of the closure relations used for the various terms in the two-scale RSM.

Colloidal gels have strong industrial relevance as they can behave as liquids or solids. The latter allows them to support a buoyant weight against gravity. However, the system is intrinsically out of equilibrium, which means that the colloids must eventually settle out of the suspension. The process of settling has been captured theoretically, but the presence of a delay time during which the gel appears relatively unaffected by gravity has not. Here, we modify existing frameworks to capture this delay, by treating the gel as a continuum with viscous response that is based on the local bond density. We can solve our model numerically to obtain the evolution of the colloid density profile and recover qualitatively the accumulation of a dense layer on top of the settling gel, as is observed experimentally in depletion gels. This numerical study is complemented by a theoretical analysis that allows us to identify an emergent time and length scale that set the dynamics of the gel. Our model provides a solid foundation for future studies that incorporate hydrodynamic erosion and tackle industrially relevant geometries.

We consider the direct numerical simulation of the flow over a forward-facing step protruding in a turbulent boundary layer. Proper orthogonal decomposition (POD) is applied to the velocity field in different regions using Fourier modes in the spanwise direction. The upstream flow is characterized by a structure with a spanwise modulation of the order of the step height, the origin of which is consistent with a centrifugal instability. The structure is associated with ejections over the step of low-speed fluid from the upstream recirculation, and organized into streaks through the action of strong spanwise motions along the step face. The spanwise-averaged instantaneous momentum deficit created by the ejections is directly related to the maximal shear rate at the edge of the step, and is well correlated with the dynamics of the downstream recirculation. The most energetic patterns consist of three-dimensional motions with a large spanwise wavelength located in the shear layer developing at the edge of the step, as well as two-dimensional fluctuations downstream of the reattachment. A linear model based on the interaction of the mean flow with the dominant POD modes is able to recover the main frequencies of the fluctuations at these wavenumbers. Including the time variations of the ejections into the model yields temporal spectra that resemble qualitatively those computed from the simulation. The results suggest that the global dynamics of the flow are at least partly driven by linear mechanisms and depend on the characteristic structure identified in the upstream region close to the step.

Rotating convection is considered on the tilted $f$-plane where gravity and rotation are not aligned. For sufficiently large rotation rates, $\Omega$, the Taylor–Proudman effect results in the gyroscopic alignment of anisotropic columnar structures with the rotation axis giving rise to rapidly varying radial length scales that vanishes as $\Omega ^{-1/3}$ for $\Omega \rightarrow \infty$. Compounding this phenomenon is the existence of viscous (Ekman) layers adjacent to the impenetrable bounding surfaces that scale as $\Omega ^{-1/2}$. In this investigation, these constraints are relaxed upon utilising a non-orthogonal coordinate representation of the fluid equations where the upright coordinate aligns with rotation axis. This exposes the problem to asymptotic perturbation methods that permit: (i) relaxation of the constraints of gyroscopic alignment; (ii) the filtering of Ekman layers through the uncovering of parameterised velocity pumping boundary conditions; and (iii) the development of reduced quasi-geostrophic systems valid in the limit $\Omega \rightarrow \infty$. Linear stability investigations reveal excellent quantitative agreement between results from parameterised or unapproximated mechanical boundary conditions. For no-slip boundaries, it is demonstrated that the associated Ekman pumping alters convective onset through an enhanced destabilisation of large spatial scales. The range of unstable modes at a fixed thermal forcing is thus significantly extended with a direct dependence on $\Omega$. This holds true even for geophysical and astrophysical regimes characterised by extreme values of the non-dimensional Ekman number $E$. The nonlinear regime is explored via the global heat and momentum transport of single-mode solutions to the quasi-geostrophic systems which indicate $O(1)$ changes which do not scale with the size of $E$.

A cylindrical liquid thread readily destabilizes into a series of drops due to capillary instability, which is also responsible for undesirable bead-on-fibre structures observed when coating a thin fibre. In this experimental study, we show how a falling liquid thread can be stabilized by internally distorting the cross-sectional shape using two vertically hung fibres. Below a critical flow rate $Q_c$, the dual-fibre system deforms the falling thread into a smooth column with a non-circular cross-section, thereby suppressing instability. Above $Q_{{c}}$, the cylindrical thread is left undeformed by the fibres and destabilizes into beads connected by a stable, non-cylindrical film. An empirical stability threshold is identified showing that flow delays the onset of capillary instability when compared with a marginally stable quasi-static coating. When the flow is unstable $Q>Q_c$, the bead velocity $v$ obeys a simple scaling law that is well supported by our experiments over a large parameter range. This suppression technique can be extended to other slender geometries, such as a ribbon, which shows similar qualitative results but exhibits a different stability threshold due to spontaneous dewetting about its short edge.