A global transient linear stability analysis of the three-dimensional time-dependent flow around an aerofoil undergoing small-amplitude pitching motion is performed using the optimally time-dependent (OTD) framework. The most salient linear instabilities associated with the instantaneous basic state are computed and tracked over time. The resulting OTD modes reflect the variations in the basic state and can be used as predictors of its spatial and temporal evolution, including the formation of a laminar separation bubble and its gradual spanwise modulation via primary global instability, leading to secondary instability and finally rapid breakdown to turbulence. The study confirms and expands upon earlier stability analyses of the same case based on the local properties of spanwise averaged velocity profiles in the bubble that predicted the onset of absolute instability soon followed by rapid breakdown of the separation bubble. The three-dimensional structure of the most unstable OTD mode is extracted, which compares well with both the locally absolutely unstable mode and the evolution of the basic state itself.

Microorganisms are ubiquitous in nature and technology. They inhabit diverse environments, ranging from small river tributaries and lakes, to oceans, as well as wastewater treatment plants and food manufacturing. In many of these environments, microorganisms coexist with settling particles. Here, we investigate the effects of microbial activity (swimming E. coli) on the settling dynamics of passive colloidal particles using particle tracking methods. Our results reveal the existence of two distinct regimes in the correlation length scale ($L_u$) and the effective diffusivity of the colloidal particles ($D_{eff}$), with increasing bacterial concentration ($\phi _b$). At low $\phi _b$, the parameters $L_u$ and $D_{eff}$ increase monotonically with increasing $\phi _b$. Beyond critical $\phi _b$, a second regime is found where both $D_{eff}$ and $L_u$ are independent of $\phi _b$. We demonstrate that the transition between these regimes is characterized by the emergence of bioconvection. We use experimentally measured particle-scale quantities $L_u$ and $D_{eff}$ to predict the critical bacterial concentration for the diffusion–bioconvection transition.

When a sphere rotates near a rigid boundary coated with a thin layer of viscous liquid, ‘tracks’ are generated both behind and over the sphere. This paper describes a theory for the simplest one-track configuration which can occur under particular experimental conditions. The theoretical predictions are in good agreement with experimental observations.

This study proposes a self-learning algorithm for closed-loop cylinder wake control targeting lower drag and lower lift fluctuations with the additional challenge of sparse sensor information, taking deep reinforcement learning (DRL) as the starting point. The DRL performance is significantly improved by lifting the sensor signals to dynamic features (DFs), which predict future flow states. The resulting DF-based DRL (DF-DRL) automatically learns a feedback control in the plant without a dynamic model. Results show that the drag coefficient of the DF-DRL model is 25 % less than the vanilla model based on direct sensor feedback. More importantly, using only one surface pressure sensor, DF-DRL can reduce the drag coefficient to a state-of-the-art performance of approximately 8 % at Reynolds number $(Re) = 100$ and significantly mitigates lift coefficient fluctuations. Hence, DF-DRL allows the deployment of sparse sensing of the flow without degrading the control performance. This method also exhibits strong robustness in flow control under more complex flow scenarios, reducing the drag coefficient by 32.2 % and 46.55 % at $Re =500$ and 1000, respectively. Additionally, the drag coefficient decreases by 28.6 % in a three-dimensional turbulent flow at $Re =10\,000$. Since surface pressure information is more straightforward to measure in realistic scenarios than flow velocity information, this study provides a valuable reference for experimentally designing the active flow control of a circular cylinder based on wall pressure signals, which is an essential step toward further developing intelligent control in a realistic multi-input multi-output system.

We study the stability of a zero-pressure gradient boundary layer subjected to free-stream disturbances by means of local stability analysis. The dataset under study corresponds to a direct numerical simulation (DNS) of a flat plate with a sharp leading edge in realistic wind tunnel conditions, with a turbulence level of 3.45 % at the leading edge. We present a method to track the convective evolution of the secondary instabilities of streaks by performing sequential stability calculations following the wave packet, connecting successive unstable eigenfunctions. A scattered nature, in time and space, of secondary instabilities is seen in the stability calculations. These instabilities can be detected before they reach finite amplitude in the DNS, preceding the nucleation of turbulent spots, and whose appearance is well correlated to the transition onset. This represents further evidence regarding the relevance of secondary instabilities of streaks in the bypass transition in realistic flow conditions. Consistent with the spatio-temporal nature of this problem, our approach allows us to integrate directly the local growth rates to obtain the spatial amplification ratio of the individual instabilities, where it is shown that instabilities reaching an $N$-factor in the range [2.5,4] can be directly correlated to more than 65 % of the nucleation events. Interestingly, it is found that high amplification is not only attained by modes with high growth rates, but also by instabilities with sustained low growth rates for a long time.

We present data from direct numerical simulation (DNS) of the fully turbulent flow through nominally two-dimensional channels containing longitudinal, surface-mounted, rectangular ribs whose widths ($W$) are either one third of or equal to the gap ($S-W$) between consecutive ribs across the domain, where $S$ is the span (centre-to-centre spacing) of the ribs. A range of the ratio of channel half-height ($H$) to span ($S$) is considered, covering $0.25\le H/S\le 2.5$. In each case, a fixed rib height ($h$) of $0.1H$ was used, but a number of cases with much smaller heights, $h/H=0.025$ or 0.05, were also studied. The secondary flows resulting from the presence of the ribs are examined, along with their sources in terms of the axial vorticity transport equation, which highlights the effects of spanwise inhomogeneity in the Reynolds stresses. We show that the strength of the secondary flows depends strongly on $H/S$ (and, correspondingly, on $W/S$) and that the major sources of axial vorticity arise near the top corners of the ribs, with convection of that vorticity dominating its spread. We also show that for smaller ribs, the secondary flow strengths are similar to those predicted by Zampino et al. (J. Fluid Mech., vol. 944, 2022, A4) using a linearised model of the Reynolds-averaged equations, which does not include the vorticity convection process; the behaviour of secondary flow topology and strength with varying $W/H$ is thus noticeably different.

In 1945–1949, Lars Onsager made an exact analysis of the high-Reynolds-number limit for individual turbulent flow realisations modelled by incompressible Navier–Stokes equations, motivated by experimental observations that dissipation of kinetic energy does not vanish. I review here developments spurred by his key idea that such flows are well described by distributional or ‘weak’ solutions of ideal Euler equations. 1/3 Hölder singularities of the velocity field were predicted by Onsager and since observed. His theory describes turbulent energy cascade without probabilistic assumptions and yields a local, deterministic version of the Kolmogorov $4/5$th law. The approach is closely related to renormalisation group methods in physics and envisages ‘conservation-law anomalies’, as discovered later in quantum field theory. There are also deep connections with large-eddy simulation modelling. More recently, dissipative Euler solutions of the type conjectured by Onsager have been constructed and his $1/3$ Hölder singularity proved to be the sharp threshold for anomalous dissipation. This progress has been achieved by an unexpected connection with work of John Nash on isometric embeddings of low regularity or ‘convex integration’ techniques. The dissipative Euler solutions yielded by this method are wildly non-unique for fixed initial data, suggesting ‘spontaneously stochastic’ behaviour of high-Reynolds-number solutions. I focus in particular on applications to wall-bounded turbulence, leading to novel concepts of spatial cascades of momentum, energy and vorticity to or from the wall as deterministic, space–time local phenomena. This theory thus makes testable predictions and offers new perspectives on large-eddy simulation in the presence of solid walls.

The initial flow past an impulsively started rotating and translating circular cylinder is asymptotically analysed using a Brinkman penalization method on the Navier–Stokes equation. In our previous study (J. Fluid Mech., vol. 929, 2021, A31), the asymptotic solution was obtained within the second approximation with respect to the small parameter, $\epsilon$, that is of the order of $1 / \lambda$. Here, $\lambda$ is the penalization parameter. In addition, the Reynolds number based on the cylinder radius and the translating velocity is assumed to be of the order of $\epsilon$. The previous study asymptotically analysed the initial flow past an impulsively started translating circular cylinder and investigated the influence of the penalization parameter $\lambda$ on the drag coefficient. It was concluded that the drag coefficient calculated from the integration of the penalization term exhibits a half-value of the results of Bar-Lev & Yang (J. Fluid Mech., vol. 72, 1975, pp. 625–647) as $\lambda \to \infty$. Furthermore, the derivative of vorticity in the normal direction was found to be discontinuous on the cylinder surface, which is caused by the tangential gradient of the pressure on the cylinder surface. The present study hence aims to investigate the variance on the drag coefficient against the result of Bar-Lev & Yang (1975). First, we investigate the problem of an impulsively started rotating circular cylinder. In this situation, the moment coefficient is independent of the pressure on the cylinder surface so that we can elucidate the role of the pressure to the hydrodynamic coefficients. Then, the problem of an impulsively started rotating and translating circular cylinder is investigated. In this situation, the pressure force induced by the unsteady flow far from the cylinder is found to play a key role on the drag force for the agreement with the result of Bar-Lev & Yang (1975), whereas the variance still exists on the lift force. To resolve the variance, an alternative formula to calculate the hydrodynamic force is derived, assuming that there is the pressure jump between the outside and inside of the cylinder surface. The pressure jump is obtained in this analysis asymptotically. Of particular interest is the fact that this pressure jump can cause the variance on the lift force calculated by the integration of the penalization term.

We investigate the global instability mechanism of the flow past a three-dimensional stepped cylinder. A comprehensive study is performed for different diameter ratios of the two joined cylinders ($r=D/d$) ranging from $r=1.1$ to $r=4$. Independently of $r$, the spectrum of the linearised Navier–Stokes operator reveals a pair of complex conjugate eigenvalues, with Strouhal number $St \approx 0.11$. The initial transition is triggered by a two-dimensional mechanism of the larger cylinder only, not affected by the presence of the junction and the smaller cylinder ($Re_{D,cr}\approx 47$). The structural sensitivity analysis is used to identify where the instability mechanism acts. The onset of transition is solely localised in the large cylinder wake (L cell), where the wavemaker has two symmetric lobes across the separation bubble. When the Reynolds number increases, a second and a third unstable pair of complex conjugate eigenvalues appears. They are localised in the small cylinder (S) wake and modulation (N) region. For any $r$, the appearance of unstable eigenmodes resembling the three cells S–N–L in the wake is observed. The nonlinear simulation results support this finding, in contrast with the previous classification of the laminar vortex shedding in direct (L–S) and indirect (L–N–S) modes interaction Lewis & Gharib (Phys. Fluids, vol. 4, 1992, pp. 104–117). This result indicates that each cell undergoes a supercritical Hopf bifurcation for any $r$. As $r$ approaches $1$, the modal linear stability results also show an unstable eigenmode in the wake of the small cylinder resembling a new modulation cell, named N2, similar to the N cell but mirrored with respect to the junction plane.

Spray formed by a myriad of secondary droplets generated by the impact of raindrops on a deep-water pool is studied with a laboratory rain facility. Experiments are performed with two rain rates and raindrops fall on the water surface at a nearly constant velocity. The secondary droplets at various heights above the pool's water surface are recorded with a cinematic digital in-line holographic technique that consists of a high-speed camera, a pulsed Nd:YLF laser and associated optics. The experimental results show that in the heat-map scatter plots of radius versus velocity near the water surface of the pool, the droplets are distributed into three regions, corresponding to distinct physical mechanisms of droplet generation. It is found that the diameter distribution of the droplets in the rain field changes with height above the pool's water surfaces. Both numerical simulation and experimental data reveal that the liquid water content, due to the presence of secondary droplets, in the atmospheric surface layer decreases exponentially with increasing height.

The present study investigates the modal stability of the steady incompressible flow inside a toroidal pipe for values of the curvature $\delta$ (ratio between pipe and torus radii) approaching zero, i.e. the limit of a straight pipe. The global neutral stability curve for $10^{-7} \leq \delta \leq ~10^{-2}$ is traced using a continuation algorithm. Two different families of unstable eigenmodes are identified. For curvatures below $1.5 \times 10^{-6}$, the critical Reynolds number ${{Re}}_{cr}$ is proportional to $\delta ^{-1/2}$. Hence, the critical Dean number is constant, ${{De}}_{cr} = 2\,{{Re}}_{cr}\,\sqrt {\delta } \approx 113$. This behaviour confirms that the Hagen–Poiseuille flow is stable to infinitesimal perturbations for any Reynolds number and suggests that a continuous transition from the curved to the straight pipe takes place as far as it regards the stability properties. For low values of the curvature, an approximate self-similar solution for the steady base flow can be obtained at a fixed Dean number. Exploiting the proposed semi-analytic scaling in the stability analysis provides satisfactory results.

Wave turbulence provides a conceptual framework for weakly nonlinear interactions in dispersive media. Dating from five decades ago, applications of wave turbulence theory to oceanic internal waves assigned a leading-order role to interactions characterized by a near equivalence between the group velocity of high-frequency internal waves with the phase velocity of near-inertial waves. This scale-separated interaction leads to a Fokker–Planck (generalized diffusion) equation. More recently, starting four decades ago, this scale-separated paradigm has been investigated using ray tracing methods. These ray methods characterize spectral transport of energy by counting the amplitude and net velocity of wave packets in phase space past a high-wavenumber gate prior to ‘breaking’. This explicitly advective characterization is based on an intuitive assignment and lacks theoretical underpinning. When one takes an estimate of the net spectral drift from the wave turbulence derivation and makes the corresponding assessment, one obtains a prediction of spectral transport that is an order of magnitude larger than either observations or reported ray tracing estimates. Motivated by this contradiction, we report two parallel derivations for transport equations describing the refraction of high-frequency internal waves in a sea of random inertial waves. The first uses standard wave turbulence techniques and the second is an ensemble-averaged packet transport equation characterized by the dispersion of wave packets about a mean drift in the spectral domain. The ensemble-averaged transport equation for ray tracing differs in that it contains the intuitively motivated advective term. We conclude that the aforementioned contradiction between theory, numerics and observations needs to be taken at face value and present a pathway for resolving this contradiction.

Surfactant-like impurities are omnipresent in multiphase emulsions and may substantially affect the motion of small droplets by altering their interfacial properties. Usually these surfactants are soluble in the bulk and undergo adsorption–desorption onto the interface which modifies their surface concentration and hence their overall influence on droplet motion. Yet, the impact of the bulk solubility and transport of surfactants on droplet dynamics, especially in the presence of bounding walls, remains poorly understood. As such, in this article, we assess the impact of bulk soluble surfactants on the settling of a spherical drop towards a plane wall. We consider coupled bulk and interfacial transport of surfactants, mediated by adsorption–desorption processes and construct a semi-analytical framework for arbitrary values of ‘bulk interaction parameter’, which dictates the strength of adsorption–desorption kinetics compared with bulk diffusion. Our results indicate that while mass exchange between the bulk and the interface can remobilize the drop, a finite bulk diffusion rate restricts this process and therefore slows down the drop. This also results in bulk concentration depletion near the south pole and accumulation near the north pole, the extent of which becomes strongly asymmetric with an enhanced intensity of depletion, as the drop approaches the wall. Presence of the wall and bulk solubility are found to aid each other towards remobilizing the drop by aptly modifying the interfacial concentration. Our results may provide fundamental insights into the kinetics of surfactant-laden drops, with potential applications in food and pharmaceutical industries, separation processes, etc.

A theory is presented for wave-driven propulsion of floating bodies driven into oscillation at the fluid interface. By coupling the equations of motion of the body to a quasipotential flow model of the fluid, we derive expressions for the drift speed and propulsive thrust of the body which in turn are shown to be consistent with global momentum conservation. We explore the efficacy of our model in describing the motion of SurferBot (Rhee et al., Bioinspir. Biomim., vol. 17, issue 5, 2022), demonstrating close agreement with the experimentally determined drift speed and oscillatory dynamics. The efficiency of wave-driven propulsion is then computed as a function of driving oscillation frequency and the forcing location, revealing optimal values for both of these parameters which await confirmation in experiments. A comparison with other modes of locomotion and applications of our model with competitive water sports is discussed in conclusion.

We study self-similar viscous fingering for the case of divergent flow within a wedge-shaped Hele-Shaw cell. Previous authors have conjectured the existence of a countably infinite number of selected solutions, each distinguished by a different value of the relative finger angle. Interestingly, the associated solution branches have been posited to merge and disappear in pairs as the surface tension decreases. For the first time, we demonstrate how the selection mechanism can be derived based on exponential asymptotics. Asymptotic predictions of the finger-to-wedge angle are additionally given for different sized wedges and surface-tension values. The merging of solution branches is explained; this feature is qualitatively different to the case of classic Saffman–Taylor viscous fingering in a parallel channel configuration. Moreover, because the asymptotic framework does not highly depend on specifics of the wedge geometry, the proposed theory for branch merging in our self-similar problem likely relates much more widely to tip-splitting instabilities in time-dependent flows in circular and other geometries, where the viscous fingers destabilise and divide in two.