This study investigates the strong influence of a splitter plate on two- and three-dimensional wake transitions of a circular cylinder. Direct numerical simulations and Floquet analyses are conducted over a parameter space including Reynolds numbers (Re) of 10–480 and non-dimensional plate lengths (L/D) of 0–6. With the increase in L/D, the critical Re for the onset of vortex shedding (Recr2D) increases monotonically. The delayed onset of vortex shedding with elongation of the body is physically explained. The critical Re for the onset of three-dimensionality (Recr3D) and the three-dimensional wake instability modes and structures are also significantly altered by the splitter plate. Compared with an isolated cylinder, the Recr3D for L/D = 1 is significantly reduced via a long wavelength mode, whereas the Recr3D for L/D = 2–6 is significantly increased via other modes. For each L/D, with increasing Re over the wake transition process, the spanwise wavelength of the wake structure gradually decreases, and the wake structure becomes increasingly chaotic. The strong influence of the splitter plate on the formation of the primary vortices and three-dimensional wake structures alter the hydrodynamic characteristics strongly. In particular, optimal lift reduction is achieved at L/D ∼ 1. In addition, the existence/absence of a hysteresis effect at the onset of three-dimensionality is identified by three methods. Among which, the method involving the Landau equation may be contaminated by initial transients induced by stable Floquet modes and may thus lead to a false conclusion on the existence/absence of hysteresis.

We investigate the motion of weakly negatively buoyant spheres settling in surface gravity waves using laboratory experiments. The trajectories of the settling spheres are tracked over most of the water depth with simultaneous measurements of the background fluid flow. These experiments are conducted in the regime relevant for environmental and geophysical applications where both particle inertia and fluid inertia are important. Using these data, we show that the sphere motion is well described by the kinematic sum of the undisturbed fluid velocity and the particle terminal settling velocity as long as the fluid inertia is not too large. We show how this result can be understood in the context of an ad hoc Maxey–Riley–Gatignol-type equation where the drag on the particle is given by the Schiller–Naumann drag correlation. We also evaluate whether inertial particles experience enhanced settling in waves, finding that measurement uncertainties in the particle terminal settling velocity and the presence of Eulerian-mean flows do not allow the small percentage increase in the settling velocity to be measured. When the fluid inertia becomes large enough, we observe path instabilities caused by particle wake effects in both quiescent and wavy conditions. However, the particle velocity fluctuations associated with the path instabilities are unaffected by the background flow. The minimal influence of the wavy flow on the particle path instabilities is thought to be due to the large-scale separation between the waves and the particle.

We present results of three-dimensional direct numerical simulations of turbulent Rayleigh–Bénard convection of dilute polymeric solutions for Rayleigh number ($Ra$) ranging from $10^6$ to $ 10^{10}$, and Prandtl number $Pr=4.3$. The viscoelastic flow is simulated by solving the incompressible Navier–Stokes equations under the Boussinesq approximation coupled with the finitely extensible nonlinear elastic Peterlin constitutive model. The Weissenberg number ($Wi$) is either $Wi=5$ or $Wi=10$, with the maximum chain extensibility parameter $L=50$, corresponding to moderate fluid elasticity. Our results demonstrate that both heat transport and momentum transport are reduced by the presence of polymer additives in the studied parameter range. Remarkably, the specific parameters used in the current numerical study give similar heat transfer reduction values as observed in experiments. We demonstrate that polymers have different effects in different regions of the flow. The presence of polymers stabilises the boundary layer, which is found to be the primary cause of the overall heat transfer reduction. In the bulk region, the presence of polymers slows down the flow by increasing the effective viscosity, enhances the coherency of thermal plumes, and suppresses the small-scale turbulent fluctuations. For small $Ra$, the heat transfer reduction in the bulk region is associated with plume velocity reduction, while for larger $Ra$, it is caused by the competing effects of suppressed turbulent fluctuations and enhanced plume coherency.

Wall-resolved large-eddy simulation (LES) of a non-equilibrium turbulent boundary layer (TBL) is performed. The simulations are based on the experiments of Volino (2020a J. Fluid Mech. 897, A2), who reported profile measurements at several streamwise stations in a spatially developing zero pressure gradient TBL evolving through a region of favourable pressure gradient (FPG), a zero pressure gradient recovery and subsequently an adverse pressure gradient (APG) region. The pressure gradient quantified by the acceleration parameter $K$ was held constant in each of these three regions. Here, $K = -(\nu /\rho U_e^{3}) {\textrm d}P_e/{\textrm d}x$, where $\nu$ is the kinematic viscosity, $\rho$ is density, $U_e$ is the free stream velocity and ${\textrm d}P_e/{\textrm d}x$ is the streamwise pressure gradient at the edge (denoted by the subscript ‘$e$’) of the TBL. The simulation set-up is carefully designed to mimic the experimental conditions while keeping the computational cost tractable. The computational grid appropriately resolves the increasingly thinning and thickening of the TBL in the FPG and APG regions, respectively. The results are thoroughly compared with the available experimental data at several stations in the domain, showing good agreement. The results show that the computational set-up accurately reproduces the experimental conditions and the results demonstrate the accuracy of LES in predicting the complex flow field of the non-equilibrium TBL. The scaling laws and models proposed in the literature are evaluated and the response of the TBL to non-equilibrium conditions is discussed.

For Stokes waves in finite depth within the neighbourhood of the Benjamin–Feir stability transition, there are two families of periodic waves, one modulationally unstable and the other stable. In this paper we show that these two families can be joined by a heteroclinic connection, which manifests in the fluid as a travelling front. By shifting the analysis to the setting of Whitham modulation theory, this front is in wavenumber and frequency space. An implication of this jump is that a permanent frequency downshift of the Stokes wave can occur in the absence of viscous effects. This argument, which is built on a sequence of asymptotic expansions of the phase dynamics, is confirmed via energetic arguments, with additional corroboration obtained by numerical simulations of a reduced model based on the Benney–Roskes equation.

The passive flight of a thin wing or plate is an archetypal problem in flow–structure interactions at intermediate Reynolds numbers. This seemingly simple aerodynamic system displays an impressive variety of steady and unsteady motions that are familiar from fluttering leaves, tumbling seeds and gliding paper planes. Here, we explore the space of flight behaviours using a nonlinear dynamical model rooted in a quasisteady description of the fluid forces. Efficient characterisation is achieved by identification of the key dimensionless parameters, assessment of the steady equilibrium states and linear analysis of their stability. The structure and organisation of the stable and unstable flight equilibria proves to be complex, and seemingly related factors such as mass and buoyancy-corrected weight play distinct roles in determining the eventual flight patterns. The nonlinear model successfully reproduces previously documented unsteady states such as fluttering and tumbling while also predicting new types of motions, and the linear analysis accurately accounts for the stability of steady states such as gliding and diving. While the conditions for dynamic stability seem to lack tidy formulae that apply universally, we identify relations that hold in certain regimes and which offer mechanistic interpretations. The generality of the model and the richness of its solution space suggest implications for small-scale aerodynamics and related applications in biological and robotic flight.

This paper explores the construction of quadratic Lyapunov functions for establishing the conditional stability of shear flows described by truncated ordinary differential equations, addressing the limitations of traditional methods like the Reynolds–Orr equation and linear stability analysis. The Reynolds–Orr equation, while effective for predicting unconditional stability thresholds in shear flows due to the non-contribution of nonlinear terms, often underestimates critical Reynolds numbers. Linear stability analysis, conversely, can yield impractically high limits due to subcritical transitions. Quadratic Lyapunov functions offer a promising alternative, capable of proving conditional stability, albeit with challenges in their construction. Typically, sum-of-squares programs are employed for this purpose, but these can result in sizable optimisation problems as system complexity increases. This study introduces a novel approach using linear transformations described by matrices to define quadratic Lyapunov functions, validated through nonlinear optimisation techniques. This method proves particularly advantageous for large systems by leveraging analytical gradients in the optimisation process. Two construction methods are proposed: one based on general optimisation of transformation matrix coefficients, and another focusing solely on the system’s linear aspects for more efficient Lyapunov function construction. These approaches are tested on low-order models of subcritical transition and a two-dimensional Poiseuille flow model with degrees of freedom nearing 1000, demonstrating their effectiveness and efficiency compared with sum-of-squares programs.

We explore a reduced-order model (ROM) of plane Couette flow with a view to performing near-wall turbulence control. The ROM is derived through Galerkin projections of the incompressible Navier–Stokes system onto a basis of controllability modes. Such ROMs were found to reproduce key aspects of turbulence dynamics in Couette flow with only a few hundred degrees of freedom, and here we use them to devise a control strategy. We consider a ROM with an extra forcing term whose structure is given by a combination of eigenfunctions of a linear viscous diffusion equation, optimised in order to minimise the total fluctuation energy. The optimisation is performed at Reynolds numbers $Re=1000, 2000, 3000$, and produces a novel control mechanism wherein the optimal forcing leads the flow to laminarisation in all cases. The forcing acts by reducing the shear in a large portion of the channel, hindering the main energy input mechanism. The forced flow possesses a new laminar solution which is linearly stable at $Re=1000$ and unstable at higher $Re$, but whose transient growth of streaky structures is substantially lower than that of laminar Couette flow, leading the flow to full laminarisation when the forcing is removed. Forcings optimised in the ROM are subsequently applied in direct numerical simulations (DNS). The same control mechanisms are observed in the DNS, where laminarisation is also achieved. We show that the ROMs provide an effective framework to design turbulence control strategies, despite the high degree of truncation, which opens up interesting possibilities for turbulence control.

Turbulent flows over rough beds with macroroughness elements of low relative submergence are characteristic of natural river systems. These flows exhibit highly three-dimensional structures, including large-scale coherent patterns, complex nonlinear interactions and significant drag induced by immobile boulders. In this study, large-eddy simulations are conducted of the flow through an array of boulders on a rough bed, based on experiments by Papanicolaou et al. (2012) Acta Geophys. 60 (6), 1502–1546. The analysis includes the instantaneous flow dynamics, the parameterisation of hydrodynamic roughness on the averaged velocity profile and the application of the double-averaged methodology. These upscaling approaches reveal the combined influence of wake turbulence and secondary currents (SCs), and provide insights into momentum and energy conservation mechanisms, which are critical for transport processes in fluvial environments. Results indicate that the boulder array reduces total fluid stress at the rough bed surface to $0.5 \rho u_*^2$, which can have important implications for sediment transport. Form-induced stresses, primarily originating in the boulder wakes, reach up to 37 % of total fluid stress, with peak values comparable to turbulent stresses at mid-boulder elevation. Form-induced kinetic energy (DKE) is shown to have the same magnitude as the turbulent kinetic energy (TKE), highlighting energy transfers from mean flow drag to DKE, then to TKE, before final dissipation. This study underscores the critical role of macroroughness in stress distribution, and the importance of the joint action of SCs and wake turbulence in driving form-induced stresses, which partially counterbalance drag dissipation.

We investigate the dynamics of circular self-propelled particles in channel flow, modelled as squirmers using a two-dimensional lattice Boltzmann method. The simulations explore a wide range of parameters, including channel Reynolds numbers ($\textit{Re}_c$), squirmer Reynolds numbers ($\textit{Re}_s$) and squirmer-type factors ($\beta$). For a single squirmer, four motion regimes are identified: oscillatory motion confined to one side of the channel, oscillatory crossing of the channel centreline, stabilisation at a lateral equilibrium position with the squirmer tilted and stable upstream swimming near the channel centreline. For two squirmers, interactions produce not only these four corresponding regimes but also three additional ones: continuous collisions with repeated position exchanges, progressive separation and drifting apart and, most notably, the formation of a stable wedge-like conformation (regime D). A key finding is the emergence of regime D, which predominantly occurs for weak pullers ($\beta = 1$) and at moderate to high $\textit{Re}_c$ values. Hydrodynamic interactions align the squirmers with streamline bifurcations near the channel centreline, enabling stability despite transient oscillations. Additionally, the channel blockage ratio critically affects the range of $\textit{Re}_s$ values over which this regime occurs, highlighting the influence of geometric confinement. This study extends the understanding of squirmer dynamics, revealing how hydrodynamic interactions drive collective behaviours. The findings also offer insights into the design of self-propelled particles for biomedical applications and contribute to the theoretical framework for active matter systems. Future work will investigate three-dimensional effects and the stability conditions for spherical squirmers forming stable wedge-like conformations, further generalising these results.

Based on data from pore-resolved direct numerical simulation of turbulent flow over mono-disperse random sphere packs, we evaluate the budgets of the double-averaged turbulent kinetic energy (TKE) and the wake kinetic energy (WKE). While TKE results from temporal velocity fluctuations, WKE describes the kinetic energy in spatial variations of the time-averaged flow field. We analyse eight cases which represent sampling points within a parameter space spanned by friction Reynolds numbers $Re_\tau \in [150, 500]$ and permeability Reynolds numbers $Re_K \in [0.4, 2.8]$. A systematic exploration of the parameter space is possible by varying the ratio between flow depth and sphere diameter $h/D \in \{ 3, 5, 10 \}$. With roughness Reynolds numbers of $k_s^+ \in [20,200]$, the simulated cases lie within the transitionally or fully rough regime. Revisiting the budget equations, we identify a WKE production mechanism via viscous interaction of the flow field with solid surfaces. The scaling behaviour of different processes over $Re_K$ and $Re_\tau$ suggests that this previously unexplored mechanism has a non-negligible contribution to the WKE production. With increasing $Re_K$, progressively more WKE is transferred into TKE by wake production. A near-interface peak in the TKE production, however, primarily results from shear production and scales with interface-related scales. Conversely, further above the sediment bed, the TKE budget terms of cases with comparable $Re_\tau$ show similarity under outer-scaling. Most transport processes relocate energy in the near-interface region, whereas pressure diffusion propagates TKE and WKE into deeper regions of the sphere pack.

Sub-convective wall pressure fluctuations play a critical role in vibroacoustic and noise analyses of vehicle structures as they serve as the primary forcing function. However, measuring these fluctuations is challenging due to their weak pressure magnitudes, typically $10^{-3}{-}10^{-5}$ of convective fluctuations. This study introduces a non-intrusive measurement technique using an array of multi-pore Helmholtz resonator sensors to capture sub-convective fluctuations with high resolution. The array features large-area, spanwise-oriented sensors arranged linearly for optimal sampling. Results provide a continuous streamwise wavenumber–frequency spectrum, resolving sub-convective fluctuations with sufficient range and accuracy. Convergence analysis indicates that long sampling durations, $\mathcal{O}(10^6 \delta ^*/U_\infty )$, $\delta^*$ is the displacement thickness of the boundary layer. $U_\infty$ is the freestream velocity are necessary to capture true sub-convective levels. Comparisons with four existing wall pressure models, which account for sensor area averaging, reveal discrepancies in predicted levels, convection speed relations and convective ridge characteristics. Notably, the measured data align most closely with the Chase (1980, J. Sound Vib., vol.70, pp. 29–67) model at convective peak levels and in the sub-convective domain. However, the observed roll-off at wavenumbers exceeding the convective wavenumber decays more slowly than predicted, giving the convective ridge an asymmetric profile about the convective line. These findings underscore the need for improved wall pressure models that incorporate frequency-dependent convective speed relations, ridge asymmetry, and more accurate sub-convective levels. Further validation using a microphone array from Farabee & Geib (1991) confirms the accuracy of our measurements, which indicate sub-convective pressure levels lower than reported previously.

The transient dynamics of a wake vortex, modelled as a strong swirling $q$-vortex, is investigated with a focus on optimal transient growth driven by continuous eigenmodes associated with continuous spectra. The pivotal contribution of viscous critical-layer eigenmodes (Lee and Marcus, J. Fluid Mech. vol. 967, 2023, p. A2) amongst the entire eigenmode families to optimal perturbations is numerically confirmed, using a spectral collocation method for a radially unbounded domain that ensures correct analyticity and far-field behaviour. The consistency of the numerical method across different sensitivity tests supports the reliability of the results and provides flexibility for tuning. Both axisymmetric and helical perturbations with axial wavenumbers of order unity or less are examined through linearised theory and nonlinear simulations, yielding results that align with existing literature on energy growth curves and optimal perturbation structures. The initiation process of transient growth is also explored, highlighting its practical relevance. Inspired by ice crystals in contrails, the backward influence of inertial particles on the vortex flow, particularly through particle drag, is emphasised. In the pursuit of optimal transient growth, particles are initially distributed at the periphery of the vortex core to disturb the flow. Two-way coupled vortex–particle simulations reveal clear evidence of optimal transient growth during ongoing vortex–particle interactions, reinforcing the robustness and significance of transient growth in the original nonlinear vortex system over finite time periods.

Contactless manipulation of small objects is essential for biomedical and chemical applications, such as cell analysis, assisted fertilisation and precision chemistry. Established methods, including optical, acoustic and magnetic tweezers, are now complemented by flow control techniques that use flow-induced motion to enable precise and versatile manipulation. However, trapping multiple particles in fluid remains a challenge. This study introduces a novel control algorithm capable of steering multiple particles in flow. The system uses rotating disks to generate flow fields that transport particles to precise locations. Disk rotations are governed by a feedback control policy based on the optimising a discrete loss framework, which combines fluid dynamics equations with path objectives into a single loss function. Our experiments, conducted in both simulations and with the physical device, demonstrate the capability of the approach to transport two beads simultaneously to predefined locations, advancing robust contactless particle manipulation for biomedical applications.

Supersonic free jets are extensively employed across a range of applications, especially in high-tech industries such as semiconductor processing and aerospace propulsion. Due to the difficulties involved in flow measurement, previous research on supersonic free jets has primarily focused on investigating near-field shockwave structures, with quantitative experimental analysis of the far-field zone being relatively scarce. However, physical understanding of the far-field flow, particularly post-shockwave energy dissipation, holds significant importance for the application and utilisation of these jets in vacuum environments. Therefore, this study aims to provide a robust experimental foundation for a rarefied supersonic free jet through the analysis of the flow field in both the near- and far-field zones. Nanometre-sized tracer particles and molecules were utilised to measure the rarefied supersonic jet flow field using particle image velocimetry and acetone molecular tagging velocimetry, respectively. The experiments revealed that in rarefied conditions, the supersonic jet exhibits a one-barrel shockwave structure in the near field, and after passing the Mach disk, a long annular viscous layer develops downstream. Experimental data on the jet velocity profile and width demonstrated a transition to a laminar flow regime in the far-field zone. This transition aligns with the theoretically inferred flow regimes based on the complex Reynolds number. The velocity profile and potential core length of the laminar flow regime could be modelled using a bi-modal distribution, which represents the summation of symmetric Gaussian distributions.

The rupture of a liquid film, where a thin liquid layer between two other fluids breaks and forms holes, commonly occurs in both natural phenomena and industrial applications. The post-rupture dynamics, from initial hole formation to the complete collapse of the film, are crucial because they govern droplet formation, which plays a significant role in many applications such as disease transmission, aerosol formation, spray drying nanodrugs, oil spill remediation, inkjet printing and spray coating. While single-hole rupture has been extensively studied, the dynamics of multiple-hole ruptures, especially the interactions between neighbouring holes, are less well understood. Here, this study reveals that when two holes ‘meet’ on a curved film, the film evolves into a spinning twisted ribbon before breaking into droplets, distinctly different from what occurs on flat films. We explain the formation and evolution of the spinning twisted ribbon, including its geometry, orbits, corrugations and ligaments, and compare the experimental observations with models. We compare and contrast this phenomena with its counterpart on planar films. While our experiments are based on the multiple-hole ruptures in corona splash, the underlying principles are likely applicable to other systems. This study sheds light on understanding and controlling droplet formation in multiple-hole rupture, improving public health, climate science and various industrial applications.

We extend the perceived velocity gradient defined by a group of particles that was previously used to investigate the Lagrangian statistics of fluid turbulence to the study of inertial particle dynamics. Using data from direct numerical simulations, we observe the correlation between the strong compression in the particle phase and the instantaneous local fluid compression. Furthermore, the Lagrangian nature of the particle velocity gradient defined in this way allows an investigation of its evolution along particle trajectories, including the process after the caustic event, or the blow-up of the particle velocity gradient. Observations reveal that, for particles with Stokes number in the range $St \lesssim 1$, inertial particles experience the maximum compression by local fluid before the caustic event. Interestingly, data analyses show that, while the post-caustic process is mainly the relaxation of the particle motion and the particle relaxation time is the relevant time scale for the dynamics, the pre-caustic dynamics is controlled by the fluid–particle interaction and the proper time scale is determined by both the Kolmogorov time and the particle relaxation time.

The jet from a model-scale, internally mixed nozzle produced a loud howling when operated at jet Mach numbers between 0.80 and 1.00. Discrete tones dominated the noise radiated to the far field and powerful oscillations were present in the jet. To explain these observations, this paper leverages a blend of experimental acoustic and flow measurements and modal analyses thereof via the spectral proper orthogonal decomposition, computational fluid dynamics simulations and local, linear stability analyses of vortex-sheet models for the flow inside the nozzle. This blend of experiments, computations and theory makes clear the cause of the howling, what sets its characteristic frequency and how it may be suppressed. The flow around a small-radius, convex bend just upstream of the final-nozzle exit led to a pocket of locally supersonic flow that was terminated by a shock. The shock was strong enough to separate the boundary layer, but neither the attached nor separated states were stable. A periodic, shock-induced separation of the boundary layer resulted, and this shock-wave/boundary-layer interaction coupled with a natural acoustic mode of the nozzle’s interior in a feedback phenomenon of sorts. Acoustic tones and large flow oscillations were produced at the associated natural frequency of the nozzle’s interior.

Ventilated cavities in the wake of a two-dimensional bluff body are studied experimentally via time-resolved X-ray densitometry. With a systematic variation of flow velocity and gas injection rate, expressed as Froude number ($\textit{Fr}$) and ventilation coefficient ($C_{qs}$), four cavities with different closure types are identified. A regime map governed by $\textit{Fr}$ and $C_{qs}$ is constructed to estimate flow conditions associated with each cavity closure type. Each closure exhibits a different gas ejection mechanism, which in turn dictates the cavity geometry and the pressure in the cavity. Three-dimensional cavity closure is seen to exist for the supercavities at low $\textit{Fr}$. However, closure is nominally two-dimensional for supercavities at higher $\textit{Fr}$. At low $C_{qs}$, cavity closure is seen to be wake-dominated, while supercavities are seen to have interfacial perturbation near the closure at higher $C_{qs}$, irrespective of $\textit{Fr}$. With the measured gas fraction, a gas balance analysis is performed to quantify the gas ejection rate at the transitional cavity closure during its formation. For a range of $\textit{Fr}$, the transitional cavity closure is seen to be characterised by re-entrant flow, whose intensity depends on the flow inertia, dictating the gas ejection rates. Two different ventilation strategies were employed to systematically investigate the formation and maintenance gas fluxes. The interaction of wake and gas injection is suspected to dominate the cavity formation process and not the maintenance, resulting in ventilation hysteresis. Consequently, the ventilation gas flux required to maintain the supercavity is significantly less than the gas flux required to form the supercavity.