The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘$-1$’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient ($A_1$) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than $4 \times 10^{7}$, the inferred $A_1$ converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘$-1$’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to $A_1$.

This study presents a novel investigation into the vortex dynamics of flow around a near-wall rectangular cylinder based on direct numerical simulation at $Re=1000$, marking the first in-depth exploration of these phenomena. By varying aspect ratios ($L/D = 5$, $10$, $15$) and gap ratios ($G/D = 0.1$, $0.3$, $0.9$), the study reveals the vortex dynamics influenced by the near-wall effect, considering the incoming laminar boundary layer flow. Both $L/D$ and $G/D$ significantly influence vortex dynamics, leading to behaviours not observed in previous bluff body flows. As $G/D$ increases, the streamwise scale of the upper leading edge (ULE) recirculation grows, delaying flow reattachment. At smaller $G/D$, lower leading edge (LLE) recirculation is suppressed, with upper Kelvin–Helmholtz vortices merging to form the ULE vortex, followed by instability, differing from conventional flow dynamics. Larger $G/D$ promotes the formation of an LLE shear layer. An intriguing finding at $L/D = 5$ and $G/D = 0.1$ is the backward flow of fluid from the downstream region to the upper side of the cylinder. At $G/D = 0.3$, double-trailing-edge vortices emerge for larger $L/D$, with two distinct flow behaviours associated with two interactions between gap flow and wall recirculation. These interactions lead to different multiple flow separations. For $G/D = 0.9$, the secondary vortex (SV) from the plate wall induces the formation of a tertiary vortex from the lower side of the cylinder. Double-SVs are observed at $L/D = 5$. Frequency locking is observed in most cases, but is suppressed at $L/D = 10$ and $G/D = 0.9$, where competing shedding modes lead to two distinct evolutions of the SV.

The transition to chaos in the subcritical regime of counter-rotating Taylor–Couette flow is investigated using a minimal periodic domain capable of sustaining coherent structures. Following a Feigenbaum cascade, the dynamics is found to be remarkably well approximated by a simple discrete map that admits rigorous proof of its chaotic nature. The chaotic set that arises for the map features densely distributed periodic points that are in one-to-one correspondence with unstable periodic orbits (UPOs) of the Navier–Stokes system. This supports the increasingly accepted view that UPOs may serve as the backbone of turbulence and, indeed, we demonstrate that it is possible to reconstruct every statistical property of chaotic fluid flow from UPOs.

Undulatory swimming is among the most common swimming forms found in nature across various length scales. In this study, we analyse the inertial effects of both the fluid and the swimmer on the transient motion of undulatory swimming using Taylor’s waving sheet model. We derive the transient velocity of the sheet for combined longitudinal and transverse waves in the Laplace domain, identifying three contributions to the velocity: the ‘slip’ velocity, fluid convection and a hydrodynamic force contribution. By numerically inverting the Laplace transform, we obtain the time history of the velocity for swimmers with varying swimming parameters and initial configurations. The acceleration performance of two types of swimmers is analysed by considering three dimensionless parameters: the acceleration rate $1/T$, sheet mass $M$, and Reynolds number $Re$, representing the effects of unsteady, convective and swimmer inertia, respectively. Under a relatively strong inertia effect, the start-up time scales as $\sim TM^2\,Re$ and $\sim TM^2$ for longitudinal and transverse waving sheets, respectively. Under weak inertia effects, the start-up time approximately reaches a constant for longitudinal waves, while it scales as $\sim T$ for transverse waves. Additionally, the transverse waving may induce a velocity overshoot, and enhances the burst swimming performance.

In rotating convection, analysis of heat transfer reveals a distinct shift in behaviour as the system transitions from a steep scaling regime near the onset of convection to a shallower scaling at higher Rayleigh numbers ($Ra$), irrespective of whether the top and bottom plates have stress-free, no-slip or no boundaries (homogeneous convection). However, while most research on this transition focuses on no-slip boundary conditions, geophysical and astrophysical flows commonly involve stress-free and homogeneous convection models as well. This study delves into the transition from the rapidly rotating regime to the non-rotating one with both stress-free and homogeneous models, leveraging direct numerical simulations (DNS) and existing literature data. Our findings unveil that for stress-free boundary conditions, the transitional Rayleigh number ($Ra_T$) exhibits a relationship $Ra_T\sim Ek^{-12/7}$, whereas for homogeneous rotating convection, $Ra_T\sim Ek^{-2}\, Pr$, where $Ek$ denotes the Ekman number, and $Pr$ denotes the Prandtl number. Both of these relationships align with the data obtained through DNS.

Head-on collisions between elliptic vortex rings (EVRs) and walls were studied experimentally using planar laser-induced fluorescence visualisations and time-resolved particle image velocimetry. Aspect ratios of $AR=2$ and 4 EVRs at a Reynolds number of $Re=4000$ were used. Collision locations were based on four key axis-switching stages of freely translating EVRs, which would shed light upon how axis-switching behaviour and aspect ratio variations affect the collision outcomes. Results show that non-uniform circumferential induced velocities in both colliding EVRs produce different behaviours along major and minor planes, where vortex-stretching/compression and hence circumferential flows play key roles in the vortex dynamics. Non-uniform formations of secondary/tertiary EVRs also lead to varied entanglements around the primary EVRs. As such, secondary vortex rings form vortex loops that may congregate along the collision axis, depending on the exact collision location. Vortex-core trajectories show the net primary/secondary vortex-core movements result from a balance between EVR diameter expansion due to collision and EVR segment motions associated with the axis-switching stage at the point of collision. Confinement effects are also observed to dominate over aspect ratio effects when the collision occurs closest to the orifice. While increasing the aspect ratio leads to different vortex-stretching/compression behaviour and more varied vortex-core trajectories due to the greater non-uniform induced velocities, they could still be understood by the preceding interpretations. Finally, three-dimensional vortex flows are reconstructed based on the experimental results to further explain the flow mechanisms.

A reactive control strategy is implemented to attenuate the streaks formed on a wing boundary layer due to free-stream turbulence (FST). Numerical simulations are performed on a section of a NACA0008 profile, considering its leading edge, while forced by FST with turbulence intensities of 0.5 % and 2.5 %. The controller is composed of localised sensors and actuators, with the control law consisting of a linear quadratic Gaussian regulator designed on a reduced-order model based only on the impulse responses of the system. Three configurations are evaluated by considering three different numbers of sensors/actuators along the spanwise direction. It is found that all configurations are effective in damping the streaks inside the boundary layer, whose effect is sustained downstream of the objective function location. However, distinct behaviours are observed when comparing the capability of the controllers with delay transition, where the best performance is attained for the case with larger number of sensors/actuators. This is attributed to the effectiveness of the controller in damping the streaks that will later break down, which in this case are associated with relatively short spanwise wavelength. This observation is confirmed by analysing the stability of the flow before the appearance of turbulent spots. Our results suggest that for an effective transition delay, efforts should not only be put into control of streaks with average spanwise wavelength, but also in the short spanwise wavelength associated with breakdown.

We investigate the deformation, dynamics and rheology of a single and a suspension of elastic capsules in inertial shear flow using high-fidelity particle-resolved simulations. For a single capsule in the shear flow, we elucidate the interplay of flow inertia and viscosity ratio, revealing the mechanism behind the stretching of capsule surface during tank-treading motion and the sign changes in normal stress differences with increasing inertia. When examining capsule suspensions, we thoroughly discuss the impact of volume fraction on average deformation, diffusion and rheology. Notably, we observe the formation of bridge structures due to hydrodynamic interactions, which enhance the inhomogeneity of the microstructure and alter the surface stress distribution within the suspension. We identify a critical Reynolds number range that marks the transition of capsule diffusion from non-inertial to inertial regimes. Furthermore, we reveal close connections between the behaviour of individual capsules and dense suspensions, particularly regarding capsule deformation and dynamics. Additionally, we propose multiple new empirical correlations for predicting the deformation factor of a single capsule and the relative viscosity of the suspension. These findings provide valuable insights into the complex behaviour of elastic capsules in inertial flows, informing the design of more accurate and efficient inertial microfluidic systems.

Sedimenting flows occur in a range of society-critical systems, such as circulating fluidised bed reactors and pyroclastic density currents (PDCs), the most hazardous volcanic process. In these systems, mass loading is sufficiently high ($\gg \mathcal {O}(1)$) and momentum coupling between the phases gives rise to mesoscale behaviour, such as formation of coherent structures capable of generating and sustaining turbulence in the carrier phase and directly impacting large-scale quantities of interest, such as settling time. While contemporary work has explored the physical processes underpinning these multiphase phenomena for monodispersed particles, polydispersed behaviour has been largely understudied. Since all real-world flows are polydisperse, understanding the role of polydispersity in gas–solid systems is critical for informing closures that are accurate and robust. This work characterises the sedimentation behaviour of two polydispersed gas–solid flows, with properties of the particles sampled from historical PDC ejecta. Highly resolved data at two volume fractions (1 % and 10 %) are collected using an EulerLagrange framework and is compared with monodisperse configurations of particles with diameters equivalent to the arithmetic mean of the polydisperse configurations. From these data, we find that polydispersity has an important impact on cluster formation and structure and that this is most pronounced for dilute flows. At higher volume fraction, the effect of polydispersity is reduced. We also propose a new metric for predicting the degree of clustering, termed ‘surface loading’, and a model for the coefficient of drag that accurately captures the settling velocity observed in the high-fidelity data.

Experiments are conducted over smooth and rough walls to explore the influence of pressure-gradient histories on skin friction and mean flow of turbulent boundary layers. Different pressure-gradient histories are imposed on the boundary layer through an aerofoil mounted in the free stream. Hot-wire measurements are taken at different free-stream velocities downstream of the aerofoil where the flow has locally recovered to zero pressure gradient but retains the history effects. Direct skin friction measurements are also made using oil film interferometry for smooth walls and a floating-element drag balance for rough walls. The friction Reynolds number, $Re_\tau$, varies between $3000$ and $27\,000$, depending both on the surface conditions and the free-stream velocity ensuring sufficient scale separation. Results align with previous findings, showing that adverse pressure gradients just upstream of the measurement location increase wake strength and reduce the local skin friction while favourable pressure gradients suppress the wake and increase skin friction. The roughness length scale, $y_0$, remains constant across different pressure-gradient histories for rough wall boundary layers. Inspired by previous works, a new correlation is proposed to infer skin friction based on the mean flow. The difference in skin friction by matching the turbulence profiles and flow structure between an arbitrary pressure-gradient history and zero pressure-gradient condition can be predicted using only the local wake strength parameter ($\Pi$), and the variations in wake strength for different histories are related to a weighted integral of the pressure-gradient history normalised by local quantities. This allows us to develop a general correlation that can be used to infer skin friction for turbulent boundary layers experiencing arbitrary pressure-gradient histories.

The quasi-steady shock refraction at a diffusive air–SF$_6$ interface (fast–slow type) is investigated numerically and theoretically. A new refraction pattern where both shock and expansion waves are simultaneously present in the reflected waves (named RRR-E) is first observed at the diffusive interface. The new refraction pattern is a regular pattern that is not expected to occur in classical shock refraction at a sharp fast–slow interface. Through the shock polar method, continuous refraction processes occur within the diffusion layer to satisfy the kinematic relationship between the reflected wave and the transmitted shock, which results in the RRR-E formation. Subsequently, the conditions for the RRR-E occurrence are obtained theoretically and verified numerically. In the phase diagram of the refraction patterns, the presence of RRR-E results in the transition boundaries of different refraction patterns at the sharp fast–slow interface no longer being valid. Specifically, the appearance of RRR-E delays the Mach reflection refraction (MRR) process, which is of great significance for the design of scramjet engines.

Geophysical and astrophysical fluid flows are typically driven by buoyancy and strongly constrained at large scales by planetary rotation. Rapidly rotating Rayleigh–Bénard convection (RRRBC) provides a paradigm for experiments and direct numerical simulations (DNS) of such flows, but the accessible parameter space remains restricted to moderately fast rotation rates (Ekman numbers ${ {Ek}} \gtrsim 10^{-8}$), while realistic ${Ek}$ for geo- and astrophysical applications are orders of magnitude smaller. On the other hand, previously derived reduced equations of motion describing the leading-order behaviour in the limit of very rapid rotation ($ {Ek}\to 0$) cannot capture finite rotation effects, and the physically most relevant part of parameter space with small but finite ${Ek}$ has remained elusive. Here, we employ the rescaled rapidly rotating incompressible Navier–Stokes equations (RRRiNSE) – a reformulation of the Navier–Stokes–Boussinesq equations informed by the scalings valid for ${Ek}\to 0$, recently introduced by Julien et al. (2024) – to provide full DNS of RRRBC at unprecedented rotation strengths down to $ {Ek}=10^{-15}$ and below, revealing the disappearance of cyclone–anticyclone asymmetry at previously unattainable Ekman numbers (${Ek}\approx 10^{-9}$). We also identify an overshoot in the heat transport as ${Ek}$ is varied at fixed $\widetilde { {Ra}} \equiv {Ra}{Ek}^{4/3}$, where $Ra$ is the Rayleigh number, associated with dissipation due to ageostrophic motions in the boundary layers. The simulations validate theoretical predictions based on thermal boundary layer theory for RRRBC and show that the solutions of RRRiNSE agree with the reduced equations at very small ${Ek}$. These results represent a first foray into the vast, largely unexplored parameter space of very rapidly rotating convection rendered accessible by RRRiNSE.

Previous studies on the scaling of pressure fluctuations in wall-bounded turbulent flows have typically employed the same frameworks as those used for mean flow, with inner scaling based on frictional velocity and viscous length scales, and outer scaling relying on boundary layer thickness or displacement thickness. These traditional scales primarily reflect the characteristics of the mean streamwise velocity profile and momentum balance. In this work, we propose novel scaling frameworks for pressure fluctuations in turbulent channel and pipe flows, derived from the Poisson equation for pressure fluctuations. Applying the scaling patch approach, we analyse the rapid and slow terms in the Poisson equation, and introduce new scaling for pressure fluctuation variance in both the inner and outer regions. These new scales are designed to better capture the influence of Reynolds stresses by incorporating their peak values. Additionally, we establish a strong correlation between the root mean square (r.m.s.) of pressure fluctuations and the Reynolds shear stress, resulting in an empirical equation that accurately predicts their ratio. This equation provides a practical method for estimating the r.m.s. of pressure fluctuations in the flow, which remains challenging to measure in experimental investigations.

In this paper, we discuss the transport of sediment and the formation of bedforms in turbulent river flows, under flow conditions typical of flooding events. Through the implementation of an immersed boundary method, a wall model and a morphological model, we were able to simulate complex and mobile geometries under high Reynolds numbers at an affordable computational cost. In particular, we examined the evolution of bedforms on a loose sediment bed under turbulent flow conditions, using input parameters obtained from laboratory measurements. Over time, the bedforms become more three-dimensional and irregular in shape, leading to changes in the shear layer, crest angle and separation patterns. The bedforms continue to evolve until a quasi-steady equilibrium is reached. Our simulations highlight the crucial role played by the small-scale bedforms, which significantly affect the flow dynamics: an increase in the total drag is observed, related to the form drag generated by the local recirculation and the increased size of the large-scale recirculation bubble. Furthermore, a stronger turbulent activity ensues from the shear layers forming on the crests of the small-scale bedforms. Finally, a wider shedding angle of the shear layer is caused by the irregular crest line.

Motivated by microfluidic applications, we investigate drag reduction in laminar pressure-driven flows in channels with streamwise-periodic superhydrophobic surfaces (SHSs) contaminated with soluble surfactant. We develop a model in the long-wave and weak-diffusion limit, where the streamwise SHS period is large compared with the channel height and the Péclet number is large. Using asymptotic and numerical techniques, we determine the influence of surfactant on drag reduction in terms of the relative strength of advection, diffusion, Marangoni effects and bulk–surface exchange. In scenarios with strong exchange, the drag reduction exhibits a complex dependence on the thickness of the bulk-concentration boundary layer and surfactant strength. Strong Marangoni effects immobilise the interface through a linear surfactant distribution, whereas weak Marangoni effects yield a quasi-stagnant cap. The quasi-stagnant cap has an intricate structure with an upstream slip region, followed by intermediate inner regions and a quasi-stagnant region that is mediated by weak bulk diffusion. The quasi-stagnant region differs from the immobile region of a classical stagnant cap, observed for instance in surfactant-laden air bubbles in water, by displaying weak slip. As exchange weakens, the bulk and interface decouple: the surfactant distribution is linear when the surfactant is strong, whilst it forms a classical stagnant cap when the surfactant is weak. The asymptotic solutions offer closed-form predictions of drag reduction across much of the parameter space, providing practical utility and enhancing understanding of surfactant dynamics in flows over SHSs.

Carbon storage in saline aquifers is a prominent geological method for reducing CO2 emissions. However, salt precipitation within these aquifers can significantly impede CO2 injection efficiency. This study examines the mechanisms of salt precipitation during CO2 injection into fractured matrices using pore-scale numerical simulations informed by microfluidic experiments. The analysis of varying initial salt concentrations and injection rates revealed three distinct precipitation patterns, namely displacement, breakthrough and sealing, which were systematically mapped onto regime diagrams. These patterns arise from the interplay between dewetting and precipitation rates. An increase in reservoir porosity caused a shift in the precipitation pattern from sealing to displacement. By incorporating pore structure geometry parameters, the regime diagrams were adapted to account for varying reservoir porosities. In hydrophobic reservoirs, the precipitation pattern tended to favour displacement, as salt accumulation occurred more in larger pores than in pore throats, thereby reducing the risk of clogging. The numerical results demonstrated that increasing the gas injection rate or reducing the initial salt concentration significantly enhanced CO2 injection performance. Furthermore, identifying reservoirs with high hydrophobicity or large porosity is essential for optimising CO2 injection processes.