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.

Periodic gravity-capillary waves on a fluid of finite depth with constant vorticity are studied theoretically and numerically. The classical Stokes expansion method is applied to obtain the wave profile and the interior flow up to the fourth order of approximation, which thereby extends the works of Barakat & Houston (1968) J. Geophys. Res. 73 (20), 6545–6554 and Hsu et al. (2016) Proc. R. Soc. Lond. A 472, 20160363. The classical perturbation scheme possesses singularities for certain wavenumbers, whose variations with depth are shown to be affected by the vorticity. This analysis also reveals that for any given value of the physical depth, there exists a threshold value of the vorticity above which there are no singularities in the theoretical solution. The validity of the third- and fourth-order solutions is examined by comparison with exact numerical results, which are obtained with a method based on conformal mapping and Fourier series expansions of the wave surface. The outcomes of this comparison are surprising as they report important differences in the internal flow structure, when compared with the third-order predictions, even though both approximations predict almost perfectly the phase velocity and the surface profiles. Usually, this occurs when the wavenumber is far enough from a critical value and the steepness is not too large. In these non-resonant cases, it is found that the fourth-order theory is more consistent with the exact numerical results. With negative vorticity the improvement is noticeable both beneath the crest and the trough, whereas with positive vorticity the fourth-order theory does a better job either beneath the crest or beneath the trough, depending of the type of the wave.

Two-dimensional simulations incorporating detailed chemistry are conducted for detonation initiation induced by dual hot spots in a hydrogen/oxygen/argon mixture. The objective is to examine the transient behaviour of detonation initiation as facilitated by dual hot spots, and to elucidate the underlying mechanisms. Effects of hot spot pressure and distance on the detonation initiation process are assessed; and five typical initiation modes are identified. It is found that increasing the hot spot pressure promotes detonation initiation, but the impact of the distance between dual hot spots on detonation initiation is non-monotonic. During the initiation process, the initial hot spot autoignites, and forms the cylindrical shock waves. Then, the triple-shock structure, which is caused by wave collisions and consists of the longitudinal detonation wave, transverse detonation wave and cylindrical shock wave, dominates the detonation initiation behaviour. A simplified theoretical model is proposed to predict the triple-point path, whose curvature quantitatively indicates the diffraction intensity of transient detonation waves. The longitudinal detonation wave significantly diffracts when the curvature of the triple-point path is large, resulting in the failed detonation initiation. Conversely, when the curvature is small, slight diffraction effects fail to prevent the transient detonation wave from developing. The propagation of the transverse detonation wave is affected not only by the diffraction effects but also by the mixture reactivity. When the curvature of the triple-point trajectory is large, a strong cylindrical shock wave is required to compress the mixture, enhancing its reactivity to ensure the transverse detonation wave can propagate without decoupling.

The Edgerton crown is an iconic manifestation of drop impact splashing, with its prominent cylindrical edge decorated with detaching droplets. Herein, we identify the formation of an intriguing double-crown, when a high-viscosity drop impacts on a shallow pool of a lower-viscosity immiscible liquid. High-speed imaging shows that after the initial fine horizontal ejecta sheet, the first inner crown emerges vertically from the film liquid. This is followed by the second crown which forms near the outer base of the first crown, as the tip of the horizontally spreading viscous drop approaches the outer free surface. Axisymmetric numerical simulations, using the volume-of-fluid method with adaptive grid refinement, show that the flow squeezed out between the viscous drop and the solid surface, generates two counter-rotating vortex rings, which travel radially outwards together and drive out the second crown through the free surface. The bottom vortex emerges from the separated boundary layer at the solid wall, while the top one detaches from the underside of the viscous drop. We map out the narrow parameter regime, where this ephemeral structure emerges, in terms of viscosity ratio, impact velocity and film thickness.

We investigate the drag reduction effects by two representative blowing/suction-based control methods having different drag reduction mechanisms, i.e. the opposition control and uniform blowing (UB), in a bump-installed turbulent channel flow through direct numerical simulations. We consider two different bulk Reynolds numbers ${\textit {Re}}_b = 5600$ and $12\,600$, and bump heights $h^+ \approx 20$ and $40$. In the opposition-controlled case, the friction drag reduction effect in the case with a bump is similar to that in the case without a bump, while the control effect on the pressure drag is hardly observed. The total drag reduction rate decreases for the higher bump height because the ratio of the pressure drag to the total drag increases as the bump height. In the UB case, UB at $0.1\,\%$ or $0.5\,\%$ of the bulk-mean velocity is imposed on the lower wall with a bump, while the same amount of uniform suction (US) is applied on the upper flat wall to keep the mass flow rate. Although the total friction drag increases due to a detrimental effect of US on the upper wall, the wall-normal motions due to the existence of a bump on the lower wall are suppressed by the UB, so that the pressure drag is decreased, unlike the opposition-controlled case. Due to the difference in the inherent drag reduction mechanisms, the flow separation in the region behind the bump is enhanced by the opposition control, while suppressed by UB.

The critical points of vorticity in a two-dimensional viscous flow are essential for identifying coherent structures in the vorticity field. Their bifurcations as time progresses can be associated with the creation, destruction or merging of vortices, and we analyse these processes using the equation of motion for these points. The equation decomposes the velocity of a critical point into advection with the fluid and a drift proportional to viscosity. Conditions for the drift to be small or vanish are derived, and the analysis is extended to cover bifurcations. We analyse the dynamics of vorticity extrema in numerical simulations of merging of two identical vortices at Reynolds numbers ranging from 5 to 1500 in the light of the theory. We show that different phases of the merging process can be identified on the basis of the balance between advection and drift of the critical points, and identify two types of merging, one for low and one for high values of the Reynolds number. In addition to local maxima of positive vorticity and minima of negative vorticity, which can be considered centres of vortices, minima of positive vorticity and maxima of negative vorticity can also exist. We find that such anti-vortices occur in the merging process at high Reynolds numbers, and discuss their dynamics.

The hydrodynamic interactions between a sedimenting microswimmer and a solid wall have ubiquitous biological and technological applications. A plethora of gravity-induced swimming dynamics near a planar no-slip wall provide a platform for designing artificial microswimmers that can generate directed propulsion through their translation–rotation coupling near a wall. In this work, we provide exact solutions for a squirmer (a model swimmer of spherical shape with a prescribed slip velocity) facing either towards or away from a planar wall perpendicular to gravity. These exact solutions are used to validate a numerical code based on the boundary integral method with an adaptive mesh for distances from the wall down to 0.1 % of the squirmer radius. This boundary integral code is then used to investigate the rich gravity-induced dynamics near a wall, mapping out the detailed bifurcation structures of the swimming dynamics in terms of orientation and distance to the wall. Simulation results show that a squirmer may traverse the wall, move to a fixed point at a given height with a fixed orientation in a monotonic way or in an oscillatory fashion, or oscillate in a limit cycle in the presence of wall repulsion.

In the present study, we investigate the modulation effects of particles on compressible turbulent boundary layers at a Mach number of 6, employing high-fidelity direct numerical simulations based on the Eulerian–Lagrangian point-particle approach. Our findings reveal that the mean and fluctuating velocities in particle-laden flows exhibit similarities to incompressible flows under compressibility transformations and semi-local viscous scaling. With increasing particle mass loading, the reduction in Reynolds shear stress and the increase in particle feedback force constitute competing effects, leading to a non-monotonic variation in skin friction, particularly in turbulence over cold walls. Furthermore, dilatational motions near the wall, manifested as travelling-wave structures, persist under the influence of particles. However, these structures are significantly weakened due to the suppression of solenoidal bursting events and the negative work exerted by the particle feedback force. These findings align with the insight of Yu et al. (J. Fluid. Mech., vol. 984, 2024, A44), who demonstrated that dilatational motions are generated by the vortices associated with intense bursting events, rather than acting as evolving perturbations beneath velocity streaks. The attenuation of travelling-wave structures at higher particle mass loadings also contributes to the reduction in the intensities of wall shear stress and heat flux fluctuations, as well as the probability of extreme events. These results highlight the potential of particle-laden flows to mitigate aerodynamic forces and thermal loads in high-speed vehicles.

We analyse the small-scale characteristics, such as enstrophy, total strain and normality/non-normality, in the three-dimensional, separated flow around a NACA 0018 wing using direct numerical simulations. The angle of attack is $10^\circ$ and the Reynolds number (based on the chord length) is $Re_c=5000$. The role of non-normality is investigated by performing Schur decomposition of the velocity gradient tensor. We also apply the Schur decomposition to derive new expressions for the production of enstrophy and total strain arising from the mean flow inhomogeneity. We focus on two sections of the flow, across the recirculating zone and along the transitioning shear layer, and compare our results with homogeneous isotropic turbulence (HIT). Within the recirculating region, the non-normality index is approximately 0 (and close to the HIT value), indicating almost equal normal and non-normal contributions. However, in the separating layer non-normal effects strongly dominate, especially in the region of kinetic energy growth. Only in the decay region do the values of the non-normality index gradually approximate HIT values. The production of enstrophy due to vortex stretching is dominated by the mixed (interaction) term, where normal strain stretches non-normal vorticity. The same component also dominates the strain self-amplification term. The contributions of different QR regions to the production terms are also examined. Production due to mean strain rate is triggered upstream compared with production due to fluctuating strain fields.

The settling velocity of frozen hydrometeors in the atmospheric surface layer depends on their inertial and drag properties, and on the intensity of ambient turbulence. Thin, solid and perforated circular disks have been investigated through high-speed imaging, under laboratory conditions, to reproduce the settling of snow plates and dendrites in quiescent and turbulent flows. Different perforations made it possible to test the parameterisation of the fall speed in quiescent air, based on the geometric description of the solidity of the disk cross-sectional area. Interestingly, different falling styles, ranging from stable horizontal to fluttering and tumbling, were observed to depend significantly on the perforation geometry, which resulted in the stabilisation of the particle rotation and in a modulation of the drag coefficient. Ambient turbulence is observed to primarily induce cross-flow drag on the disks settling in the nonlinear regime, causing a reduction of the settling velocity in all cases investigated. Turbulence also manifests a secondary effect on the disk rotational dynamics, in particular when tumbling and stable falling styles co-exist, based on the phase space defined by the Reynolds number $Re$ and the inertia ratio $I^*$. The interaction between ambient turbulence, particle anisotropy and permeability and the likelihood of tumbling is inferred to be the main reason for the observed settling velocity variability of snow dendrites in nature.