The influence of outer large-scale motions (LSMs) on near-wall structures in compressible turbulent channel flows is investigated. To separate the compressibility effects, velocity fluctuations are decomposed into solenoidal and dilatational components using the Helmholtz decomposition method. Solenoidal velocity fluctuations manifest as near-wall streaks and outer large-scale structures. The spanwise drifting of near-wall solenoidal streaks is found to be driven by the outer LSMs, while LSMs have a trivial influence on the spanwise density of solenoidal streaks, consistent with the outer LSM impacts found in incompressible flows (Zhou et al., J. Fluid Mech., vol. 940, 2022, p. A23). Dilatational motions are characterized by the near-wall small-scale travelling-wave packets and the large-scale parts in the outer region. The streamwise advection velocity of the near-wall structures remains at $16 \sim 18u_{\tau }$, hardly influenced by Mach numbers, Reynolds numbers and wall temperatures. The spanwise drifting of near-wall dilatational structures, quantified by the particle image velocimetry method, follows a mechanism distinct from solenoidal streaks. This drifting velocity is notably larger than those of the solenoidal streaks, and the influence of outer LSMs is not the primary trigger for this drifting.

We present an experimental study of reactive control of turbulent jets, in which we target axisymmetric coherent structures, known to play a key role in the generation of sound. We first consider a forced jet, in which coherent structures are amplified above background levels, facilitating their detection, estimation and control. We then consider the more challenging case of an unforced jet. The linear control targets coherent structures in the region just downstream of the nozzle exit plane, where linear models are known to be appropriate for description of the lowest-order azimuthal modes of the turbulence. The control law is constructed in frequency space, based on empirically determined transfer functions. And the Wiener–Hopf formalism is used to enforce causality and to provide an optimal controller, as opposed to the sub-optimal control laws provided by simpler wave-cancellation methods. Significant improvements are demonstrated in the control of both forced and unforced jets. In the former case, order-of-magnitude reductions are achieved; and in the latter, turbulence levels are reduced by up to 60 %. The results open new perspectives for the control of turbulent flow at high Reynolds number.

The stimulation of instability and transport in the bottom boundary layer by internal solitary waves has been documented for over twenty years. However, the challenge of shallow slopes and a disparity of scales between the large-scale wave and the small-scale boundary layer has proven challenging for simulations. We present laboratory scale simulations that resolve the three-dimensionalisation in the boundary layer during the entire shoaling process. We find that the late stage, in which the incoming wave fissions into boluses, provides the most consistent source of three-dimensionalisation. In the early stage of shoaling, three-dimensionalisation occurs not so much due to separation bubble instability, but due to the interaction of vortices shed from the separation bubble with the overlying pycnocline. This interaction overturns the pycnocline, and creates bursts in kinetic energy and viscous dissipation, suggesting that the shed vortices induce turbulent motion and sediment resuspension in the water column above and behind the separation bubble.

The dynamics of a single highly elastic fibre settling under gravity in a very viscous fluid is studied numerically. We employ the bead model and multipole expansion of the Stokes equations, corrected for lubrication that is implemented in the precise Hydromultipole numerical codes. Four attracting regular dynamical modes of highly elastic fibres are found: two stationary shapes (one translating and the other rotating and translating), and two periodic oscillations around such shapes. The phase diagram of these modes is presented. It illustrates that the existence of each mode depends not only on the elasto-gravitation number but also on the fibre aspect ratio. Characteristic time scales, fibre deformation patterns and motion in the different modes are determined.

Steep wave breaking on a vertical cylinder (a typical foundation supporting offshore wind turbines) will induce slam loads. Many questions on the important violent wave loading and the associated secondary load cycle remain unanswered. We use laboratory experiments with unidirectional waves to investigate the fluid loading on vertical cylinders. We use a novel three-phase decomposition approach that allows us to separate different types of nonlinearity. Our findings reveal the existence of an additional quasi-impulsive loading component that is associated with the secondary load cycle and occurs in the backwards direction against that of the incoming waves. This quasi-impulsive force occurs at the end of the secondary load cycle and close to the passage of the downward zero-crossing point of the undisturbed wave. Wavelet analysis showed that the impulsive force exhibits superficially similar behaviour to a typical wave-slamming event but in the reverse direction. To monitor the scattered wave field and extract run-up on the cylinder, we installed a four-camera synchronised video system and found a strong temporal correlation between the arrival time of the Type-II scattered wave onto the cylinder and the occurrence of this quasi-impulsive force. The temporal characteristics of this quasi-impulsive force can be approximated by the Goda wave impact model, taking the collision of the Type-II scattered waves at the rear stagnation point as the impact source.

The yield stress and shear thinning properties of mucus are identified as critical for ciliary coordination and mucus transport in human airways. We use here numerical simulations to explore the hydrodynamic coupling of cilia and mucus with these two properties using the Herschel–Bulkley model, in a lattice Boltzmann solver for the fluid flow. Three mucus flow regimes, i.e. a poorly organized regime, a swirly regime, and a fully unidirectional regime, are observed and analysed by parametric studies. We systematically investigate the effects of ciliary density, interaction length, Bingham number and flow index on the mucus flow regime formation. The underlying mechanism of the regime formation is analysed in detail by examining the variation of two physical quantities (polarization and integral length) and the evolution of the flow velocity, viscosity and shear-rate fields. Mucus viscosity is found to be the dominant parameter influencing the regime formation when enhancing the yield stress and shear thinning properties. The present model is able to reproduce the solid body rotation observed in experiments (Loiseau et al., Nat. Phys., vol. 16, 2020, pp. 1158–1164). A more precise prediction can be achieved by incorporating non-Newtonian properties into the modelling of mucus as proposed by Gsell et al. (Sci. Rep., vol. 10, 2020, 8405).

In this experimental study, we explore the dynamics of the thermal boundary layer in liquid metal Rayleigh–Bénard convection, covering the parameter ranges of $0.026 \leq$ Prandtl numbers $(Pr) \leq 0.033$ and Rayleigh numbers ($Ra$) up to $2.9\times 10^9$. Our research focuses on characterising the thermal boundary layer near the top plate of a cylindrical convection cell with an aspect ratio of 0.5, distinguishing between two distinct regions: the shear-dominated region around the centre of the top plate and a location near the side wall where the boundary layer is expected to be affected by the impact or ejection of thermal plumes. The dependencies of the boundary layer thickness on $Ra$ at these positions reveal deviating scaling exponents with the difference diminishing as $Ra$ increases. We find stronger fluctuations in the boundary layer and increasing deviation from the Prandtl–Blasius–Pohlhausen profile with increasing $Ra$, as well as in the measurements outside the centre region. Our data illustrate the complex interplay between flow dynamics and thermal transport in low-$Pr$ convection.

We consider an internally heated fluid between parallel plates with fixed thermal fluxes. For a large class of heat sources that vary in the direction of gravity, we prove that $\smash { \smash {{\langle {\delta T} \rangle _h}} } \geq \sigma R^{-1/3} - \mu$, where $\smash { \smash {{\langle {\delta T} \rangle _h}} }$ is the average temperature difference between the bottom and top plates, $R$ is a ‘flux’ Rayleigh number and the constants $\sigma,\mu >0$ depend on the geometric properties of the internal heating. This result implies that mean downward conduction (for which $\smash { \smash {{\langle {\delta T} \rangle _h}} }< 0$) is impossible for a range of Rayleigh numbers smaller than a critical value $R_0:=(\sigma /\mu )^{3}$. The bound demonstrates that $R_0$ depends on the heating distribution and can be made arbitrarily large by concentrating the heating near the bottom plate. However, for any given fixed heating profile of the class we consider, the corresponding value of $R_0$ is always finite. This points to a fundamental difference between internally heated convection and its limiting case of Rayleigh–Bénard convection with fixed-flux boundary conditions, for which $\smash {{\langle {\delta T} \rangle _h}}$ is known to be positive for all $R$.

Two-dimensional free-surface flow over localised topography is examined, with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a dip, respectively. Steady hydraulic-fall solutions to the full incompressible, irrotational Euler equations are computed, and their linear and nonlinear stability is analysed by computing eigenspectra of the pertinent linearised operator and by solving an initial value problem. The computations are carried out numerically using a specially developed computational framework based on the finite-element method. The Hamiltonian structure of the problem is demonstrated, and stability is determined by computing eigenspectra of the pertinent linearised operator. It is found that a hydraulic-fall flow over a bump is spectrally stable. The corresponding flow over a dip is found to be linearly unstable. In the latter case, time-dependent simulations show that ultimately, the flow settles into a time-periodic motion that corresponds to an invariant solution in an appropriately defined phase space. Physically, the solution consists of a localised large-amplitude wave that pulsates above the dip while simultaneously emitting nonlinear cnoidal waves in the upstream direction and multi-harmonic linear waves in the downstream direction.

Consider the motion of a thin layer of electrically conducting fluid, between two closely spaced parallel plates, in a classical Hele-Shaw geometry. Furthermore, let the system be immersed in a uniform external magnetic field (normal to the plates) and let electrical current be driven between conducting probes immersed in the fluid layer. In the present paper, we analyse the ensuing fluid flow at low Hartmann numbers. Physically, the system is particularly interesting because it allows for circulation in the flow, which is not possible in the standard pressure-driven Hele-Shaw cell. We first elucidate the mechanism of flow generation both physically and mathematically. After formulating the problem using complex variables, we present mathematical solutions for a class of canonical multiply connected geometries in terms of the prime function framework developed by Crowdy (Solving Problems in Multiply Connected Domains, SIAM, 2020). We then demonstrate how recently developed fast numerical methods may be applied to accurately determine the flow field in arbitrary geometries.