Drop-on-demand printing applications involve a drop connected to a fluid reservoir between which volume can be exchanged, a situation that can be idealized as a sessile drop with prescribed volume flux across the drop/reservoir boundary. Here we compute the frequency spectrum for these pressure disturbances, as it depends upon the static contact-angle $\alpha$ (CA) and an empirical constant $\chi$ relating the reservoir pressure to volume exchanged, for either (i) pinned or (ii) free contact-lines. Mode shapes are characterized by the mode number pair $[k,\ell ]$ with property $k+\ell =\mathbb {Z}^{+}_{odd}$ that can be associated with the symmetry properties of the Rayleigh drop modes for the free sphere. We report instabilities to the axisymmetric $[1,0]$ and non-axisymmetric rocking $[2,1]$ modes that are related to centre-of-mass motions, and show how the spectral degeneracy of the Rayleigh drop modes breaks with the model parameters.

A simple model for the motion of shape-changing swimmers in Poiseuille flow was recently proposed and numerically explored by Omori et al. (J. Fluid Mech., vol. 930, 2022, A30). These explorations hinted that a small number of interacting mechanics can drive long-time behaviours in this model, cast in the context of the well-studied alga Chlamydomonas and its rheotactic behaviours in such flows. Here, we explore this model analytically via a multiple-scale asymptotic analysis, seeking to formally identify the causal factors that shape the behaviour of these swimmers in Poiseuille flow. By capturing the evolution of a Hamiltonian-like quantity, we reveal the origins of the long-term drift in a single swimmer-dependent constant, whose sign determines the eventual behaviour of the swimmer. This constant captures the nonlinear interaction between the oscillatory speed and effective hydrodynamic shape of deforming swimmers, driving drift either towards or away from rheotaxis.

Motivated by the variation of local shear produced by internal waves in the ocean, we use direct numerical simulations to investigate the effect of a time-dependent shear forcing on the evolution and mixing of turbulence produced by Kelvin–Helmholtz instability (KHI) at high Reynolds number. The forcing is implemented using a tilting coordinate system which causes the background shear to accelerate and decelerate periodically. We demonstrate that, with suitable timing between development of instability and the shear oscillation cycle, turbulence produced by KHI with a decelerating shear mixes in a distinctly different way from the flow with constant background shear, specifically with the energy for turbulent motions extracted from alternative sources. As a result, the total amount of mixing as measured by the change in background potential energy can in fact be significantly larger for flows in which the shear is decelerated, despite the fact that the total kinetic energy in the flow is significantly smaller. The mixing has characteristics more in common with convectively driven rather than shear-driven flows, supporting the argument for an underlying change in the mechanisms triggering the turbulence.

Analysis of the classic problem of shallow film flow on an inclined plane is revisited for a Brownian suspension. The particle phase, considered in a two-fluid model, is predicted to cause pronounced changes to the instability characteristics of the flow. These are due to an indirect effect of the non-Newtonian rheology, with the normal stresses causing migration and a viscosity stratification which strongly alters the base state from its Newtonian counterpart. Both the short- and long-wavelength inertial stability boundaries are altered, and, more strikingly, an instability at zero Reynolds number is found.

We investigate the counter-intuitive initial decrease and subsequent increase in the Nusselt number $Nu$ with increasing wall Reynolds number $Re_w$ in the sheared Rayleigh–Bénard (RB) system by studying the energy spectra of convective flux and turbulent kinetic energy for Rayleigh number $Ra = 10^{7}$, Prandtl number $Pr=1.0$ and inverse Richardson numbers $0 \leq 1/Ri \leq 10$. These energy spectra show two distinct high-energy regions corresponding to the large-scale superstructures in the bulk and small-scale structures in the boundary layer (BL) regions. A greater separation between these scales at the thermal BL height correlates to a higher $Nu$ and indicates that the BLs are more turbulent. The minimum $Nu$, which occurs at $1/Ri=1.0$, is accompanied by the smallest separation between the large- and small-scale structures at the thermal BL height. At $1/Ri=1.0$, we also observe the lowest value of turbulent kinetic energy normalized with the square of friction velocity within the thermal BL. Additionally, we find that the domain size has a limited effect on the heat and momentum transfer in the sheared RB system as long as the domain can accommodate the small-scale convective structures at the thermal BL height, signifying that capturing the large-scale superstructures is not essential to obtain converged values of $Nu$ and shear Reynolds number $Re_{\tau }$. When the domain is smaller than these small-scale convective structures, the overall heat and momentum transfer reduces drastically.

Stochastic models of varying complexity have been proposed to describe the dispersion of particles in turbulent flows, from simple Brownian motion to complex temporally and spatially correlated models. A method is needed to compare competing models, accounting for the difficulty in estimating the additional parameters that more complex models typically introduce. We employ a data-driven method, Bayesian model comparison, which assigns probabilities to competing models based on their ability to explain observed data. We focus on the comparison between the Brownian and Langevin dynamics for particles in two-dimensional isotropic turbulence, with data that consist of sequences of particle positions obtained from simulated Lagrangian trajectories. We show that, while on sufficiently large time scales the models are indistinguishable, there is a range of time scales on which the Langevin model outperforms the Brownian model. While our set-up is highly idealised, the methodology developed is applicable to more complex flows and models of particle dynamics.

We seek the best stroke sequences of a finite-size swimming predator chasing a non-motile point or finite-size prey at low Reynolds number. We use optimal control to seek the globally optimal solutions for the former and reinforcement learning (RL) for general situations. The predator is represented by a squirmer model that can translate forward and laterally, rotate and generate a stresslet flow. We identify the predator's best squirming sequences to achieve the time-optimal (TO) and efficiency-optimal (EO) predation. For a point prey, the TO squirmer executing translational motions favours a two-fold $L$-shaped trajectory that enables it to exploit the disturbance flow for accelerated predation; using a stresslet mode expedites significantly the EO predation, allowing the predator to catch the prey faster yet with lower energy consumption and higher predatory efficiency; the predator can harness its stresslet disturbance flow to suck the prey towards itself; compared to a translating predator, its compeer combining translation and rotation is less time-efficient, and the latter occasionally achieves the TO predation via retreating in order to advance. We also adopt RL to reproduce the globally optimal predatory strategy of chasing a point prey, qualitatively capturing the crucial two-fold attribute of a TO path. Using a numerically emulated RL environment, we explore the dependence of the optimal predatory path on the size of prey. Our results might provide useful information that help in the design of synthetic microswimmers such as in vivo medical microrobots capable of capturing and approaching objects in viscous flows.

A rapid predictive tool based on the linearised Reynolds-averaged Navier–Stokes equations is proposed in this work to investigate secondary currents generated by streamwise-independent surface topography modulations in turbulent channel flow. The tool is derived by coupling the Reynolds-averaged momentum equation to the Spalart–Allmaras transport equation for the turbulent eddy viscosity, using a nonlinear constitutive relation for the Reynolds stresses to capture correctly secondary motions. Linearised equations, describing the steady flow response to arbitrary surface modulations, are derived by assuming that surface modulations are shallow. Since the equations are linear, the superposition principle holds and the flow response induced by an arbitrary modulation can be obtained by combining appropriately the elementary responses obtained over sinusoidal modulations at multiple spanwise length scales. The tool permits a rapid exploration of large parameter spaces characterising structured surface topographies previously examined in the literature. Here, channels with sinusoidal walls and with longitudinal rectangular ridges are considered. For sinusoidal walls, a large response is observed at two spanwise wavelengths scaling in inner and outer units respectively, mirroring the amplification mechanisms in turbulent shear flows observed from transient growth analysis. For longitudinal rectangular ridges, the model suggests that the analysis of the response and the interpretation of the topology of secondary structures is facilitated when the ridge width and the gap between ridges are used instead of other combinations proposed in the literature.

This study examines the impact of small-scale irregular topographic features on the dynamics and evolution of large-scale barotropic flows in the ocean. A multiscale theory is developed, which makes it possible to represent large-scale effects of the bottom roughness without explicitly resolving small-scale variability. The analytical model reveals that the key mechanism of topographic control involves the generation of a small-scale eddy field associated with considerable Reynolds stresses. These eddy stresses are inversely proportional to the large-scale velocity and adversely affect mean circulation patterns. The multiscale model is applied to the problem of topography-induced spin-down of a large circularly symmetric vortex and is validated by corresponding topography-resolving simulations. The small-scale bathymetry chosen for this configuration conforms to the Goff–Jordan statistical spectrum. While the multiscale model formally assumes a substantial separation between the scales of interacting flow components, it is remarkably accurate even when scale separation is virtually non-existent.

Linear modal instabilities of flow over untapered wings with aspect ratios $AR=4$ and 8, based on the NACA 0015 profile, have been investigated numerically over a range of angles of attack, $\alpha$, and angles of sweep, $\varLambda$, at chord Reynolds numbers $100\le Re\le 400$. Laminar base flows have been generated using direct numerical simulation and selective frequency damping, as appropriate. Several families of unstable three-dimensional linear global (TriGlobal) eigenmodes have been identified and their dependence on geometric parameters has been examined in detail at $Re=400$. The leading global mode A is associated with the peak recirculation in the three-dimensional laminar separation bubble formed on the wing and becomes unstable when recirculation reaches $\textit {O}(10\,\%)$. On unswept wings, this mode peaks in the midspan region of the wake and moves towards the wing tip with increasing $\varLambda$, following the displacement of peak recirculation; its linear amplification leads to wake unsteadiness. Additional amplified modes exist at nearly the same and higher frequencies compared to mode A. The critical $Re$ has been identified and it is shown that amplification increases with increasing sweep, up to $\varLambda \approx 10^\circ$. At higher $\varLambda$, all global modes become less amplified and are ultimately stable at $\varLambda =30^\circ$. An increase in amplification of the leading mode with sweep was not observed over the $AR=4$ wing, where tip vortex effects were shown to dominate.

This paper derives approximate ‘sunset’ similarity solutions for receding plane strain and radially symmetric hydraulic fractures in permeable elastic media close to the point of closure. Local analysis is used to show that a receding hydraulic fracture has a linear aperture asymptote $\hat {w}\sim \hat {s}$ in the fracture tip, where $\hat {s}$ is the distance from the fracture front. Due to the regularity of the linear asymptote, it is possible to determine similarity solutions in the form of power series expansions, which, for integers $N\ge 2$ and values of the radius decay exponent $\gamma =1/N$, can be shown to terminate to yield polynomial solutions for the fracture aperture of degree $N$. Of this countable infinity of polynomial solutions, the final aperture profile as the fracture approaches closure is associated with the second-degree polynomial with $\gamma =1/2$ called the sunset solution. For the reverse time $t^{\prime }$ measured from closure, the sunset solution is characterized by $w\sim t^{\prime }$ and $R\sim t^{\prime 1/2}$. Of all the admissible polynomial similarity solutions, the sunset solution is shown to form an attractor, as $t^{\prime }\rightarrow 0$, for receding hydraulic fractures associated with a wide variety of points in parametric space. Using the sunset solution, it is possible to estimate the duration of recession, assuming that the fracture aperture and radius at the start of recession are given, and determine how it scales with a dimensionless shut-in parameter. As the fracture approaches closure, the term responsible for coupling the elastic force balance and fluid conservation becomes subdominant to the other terms in the lubrication equation, which reduces to a local kinematic relation between the decaying fracture aperture and the leak-off velocity. This fundamental decoupling of dynamics from kinematics results in the sunset solution being dependent on only a single material parameter – namely the leak-off coefficient. This isolation of the leak-off coefficient by the sunset solution opens the possibility to determine this parameter from laboratory or field measurements.

When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because of the small air–water density ratio, $\rho _{{a}}/\rho _{{w}}\equiv \epsilon \ll 1$, the wind and the ripples are weakly coupled, and the eigenvalue problem can be solved perturbatively. At the leading order, the eigenvalue is equal to the phase speed $c_0$ of surface waves. At order $\epsilon$, the eigenvalue has a finite imaginary part, which implies growth. Miles (J. Fluid Mech., vol. 3, 1957, pp. 185–204) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level $z=z_c$, where the wind speed is equal to $c_0$ and the waves extract energy from the wind. Here, we construct uniform asymptotic approximations of the leading-order eigenfunction for long waves, which we use to calculate the growth rate as a function of $k$. In the strong wind limit, we find that the fastest growing wave is such that the aerodynamic pressure is in phase with the wave slope. The results are confirmed numerically.

The effect of gravity during the water entry of two-dimensional and axisymmetric bodies is investigated analytically and numerically. An extension to the Wagner model of water impact is proposed in order to take into account the effect of gravity. For this purpose, the free-surface condition is modified. The pressure is computed using the modified Logvinovich model of Korobkin (Eur. J. Appl. Maths, vol. 6, 2004, pp. 821–838). The model has been implemented and validated through comparisons with fully nonlinear potential flow simulations of different two-dimensional and axisymmetric water entry problems. Our investigation shows that it is equally important to account for gravity when computing the pressure distribution and to account for gravity when computing the size of the wetted surface in order to obtain accurate force results with the Wagner model. Simulations of wedges and cones with different values of deadrise angle ($\beta$) entering water at constant speed ($V$) demonstrate the accuracy of the semi-analytical model and show that the effect of gravity in such water impacts is governed by the effective Froude number defined as $Fr^*=V/(\sqrt {gh}\sqrt {\tan \beta })$, with $g$ the acceleration due to gravity and $h$ the penetration depth. The accuracy of the semi-analytical model for decelerated water entries is also demonstrated by investigating the water entry of a wedge and a cone with a $15^\circ$ deadrise angle with deceleration until full stop. The semi-analytical model is able to accurately predict the effect of gravity during both two-dimensional and axisymmetric water entry problems with deceleration.

The unsteady wind loading on high-rise buildings has the potential to influence strongly their structural performance in terms of serviceability, habitability and occupant comfort. This paper investigates numerically the flow structures around a canonical high-rise building immersed in an atmospheric boundary layer, using wall-resolved large eddy simulations. The switching between two vortex shedding modes is explored, and the influence of the atmospheric boundary layer on suppressing symmetric vortex shedding is identified. It is shown that the antisymmetric vortex shedding mode is prevalent in the near wake behind the building, with strong coherence between the periodic fluctuations of the building side force and the antisymmetric vortex shedding mode demonstrated. Two feedback control strategies, exploiting this idea, are designed to alleviate the aerodynamic side-force fluctuations, using pressure sensing on just a single building wall. The sensor response to synthetic jet actuation along the two ‘leading edges’ of the building is characterised using system identification. Both the designed linear controller and the least mean square adaptive controller attenuate successfully the side-force fluctuations when implemented in simulations. The linear controller exhibits a better performance, and its effect on the flow field is to delay the formation of dominant vortices and increase the extent of the recirculation region. Feedback control that requires a smaller sensing area is then explored, with a comparable control effect achieved in the attenuation of the unsteady loading. This study could motivate future attempts to understand and control the unsteady loading of a high-rise building exposed to oncoming wind variations.

The Rayleigh collapse time is the time it would take to shrink an empty spherical bubble in an infinite liquid domain to zero size, which is a function of ambient pressure and initial bubble radius. If a solid boundary is located in the vicinity of the shrinking or collapsing bubble, then liquid flow is hindered, such that the collapse time is prolonged. This can be quantified with the Rayleigh prolongation factor $k$. Here, we provide $k$ for intermediate to smallest bubble to wall stand-off distances. It is measured with single laser-induced cavitation bubbles in water close to a solid boundary. Maximum bubble radii are determined from microscopic high-speed imaging at one million frames per second. Collapse times are measured acoustically via the acoustic transients emitted during bubble seeding and collapse. The experimental findings are compared, with good agreement, to numerical simulations based on a volume of fluid method. As a result, a polynomial fit of $k$ versus stand-off distance is given for the near-wall bubble collapse in water. Then the influence of the viscosity on $k$ is studied numerically in the near-wall regime.

A combination of time-resolved tomographic particle image velocimetry, refractive index matching technique and machine vision algorithms was used to measure the translational and rotational motion of freely moving, nearly neutrally buoyant spheres in a fully developed turbulent boundary layer (TBL). Located in the buffer and logarithmic layers, the hydrogel spheres ($\sim$70 inner wall units in diameter) were refractive index matched with the water and tagged by ‘spokes’. Besides translational motion, the spheres exhibited significant rotation. The spheres were surrounded by typical coherent structures observed in TBLs, among them hairpin packets and transverse and longitudinal vortices that induced ejections and sweeps. While the majority of instantaneous sphere Reynolds numbers did not exceed 100, and vortex shedding was not observed, the results showed that the spheres may affect the evolution of hairpin packets in TBLs due to their finite size. The instantaneous rotation-, wall- and shear-induced lift forces, as well as the drag forces, acting on the spheres were estimated using available correlations for the lift and drag coefficients. Results hinted at negative shear-induced lift due to flow separation at a smaller critical Reynolds number than incorporated in the correlations that do not include the effect of ambient turbulence. The results indicated further that the drag force aided by the rotation-induced lift force was instrumental in keeping one of the spheres aloft. For the wall-ward moving spheres, lift forces opposed sphere motion. As a result, the spheres approached the wall with velocities lower than their quiescent settling velocity.

We investigate experimentally the coalescence cascade process for a confined swarm of deformable bubbles immersed in a bidimensional vertical cell filled with water. For different gas volume fractions, air bubbles of size $D_0$ larger than the cell thickness are injected at the bottom of the cell. The bubbles swarms transformation is explored using high-speed visualizations. The time evolution of each bubble in the swarm is determined using a specifically developed algorithm, enabling bubble tracking and coalescence detection. We determine the evolution of the bubble size distribution downstream from the injection point, and show that the stages of the coalescence cascade are characterized by the diameter, $D_{V90}$, representative of the largest bubbles. The collision frequency of pairs of bubbles of sizes $D_k$ and $D_{k'}$, $h(D_k, D_{k'})$, and their coalescence efficiency, $\lambda$, are obtained from the experiments. The efficiency is nearly constant, independently of the bubble sizes and of the gas volume fraction. Concerning collision frequency, our results reveal the existence of two different coalescence regimes depending on the capability of the bubbles to deform. Models describing $h(D_k, D_{k'})$ for both regimes are provided. They take into account the specific response of the bubble pair, which depends on the reduced diameter $D_p = 2 D_k D_{k'} / (D_k + D_{k'})$, to the global swarm-induced agitation governed by $D_{V90}$ and the gas volume fraction. In the first regime, occurring for smaller $D_p$, bubbles are brought together by agitation and rapidly coalesce, while for sufficiently large $D_p$, both bubbles are able to deform and spend more time adapting mutually their shapes before coalescing.

The influence of the wing root junction on the laminar separation bubble forming on the suction surface of a semispan NACA 0018 wing cantilevered from the wind tunnel test section wall is studied using surface flow visualisations, particle image velocimetry and surface pressure measurements at a chord Reynolds number of 125 000 and an angle of attack of 6$^\circ$. The test section wall boundary layer upstream of the wing is turbulent, and the spanwise influence of the junction on the separation bubble extends well beyond the test section wall boundary layer thickness. Substantial three-dimensionality is seen in the separation bubble flowfield near the wing root, where earlier transition and a reduction in separation bubble thickness is observed. In contrast with the wing tip, earlier transition and a reduction in separation bubble length occurs near the wing root. Outside of the junction affected region, the separation bubble is similar to separation bubbles forming on two-dimensional geometries, and displays mild spanwise waviness. The transition process away from the end affected regions is characterised by the formation of spanwise roll-up vortices that are shed in a nearly two-dimensional manner across the span. The analysis of the results shows that, near the wing root, the increased level of perturbations leads to earlier vortex roll-up and spanwise flow contributes to more rapid vortex breakdown. The results in the wing root region are complemented by the analysis of data from Toppings and Yarusevych (J. Fluid Mech., vol. 929, 2021, A39) in the wing tip region to provide a more holistic outlook on the laminar separation bubble topology and dynamics on the entire finite wing.

This article deals with the symmetry breaking of azimuthal thermoacoustic modes in annular combustors. Using a nominally symmetric annular combustor, we present experimental evidence of a predicted spontaneous reflectional symmetry breaking, and also an unexpected explicit rotational symmetry breaking in the neighbourhood of the Hopf bifurcation which separates linearly stable azimuthal thermoacoustic modes from self-oscillating modes. We derive and solve a multidimensional Fokker–Planck equation to unravel a unified picture of the phase space topology. We demonstrate that symmetric probability density functions of the thermoacoustic state vector are elusive, because the effect of asymmetries, even imperceptible ones, is magnified close to the bifurcation. This conclusion implies that the thermoacoustic oscillations of azimuthal modes in real combustors will systematically exhibit a statistically dominant orientation of the mode in the vicinity of the Hopf bifurcation.

Transonic buffet refers to the self-sustained periodic motion of shock waves observed in transonic flows over wings and can limit the flight envelope of aircraft. Based on the boundary layer characteristics at the shock foot, buffet has been classified as laminar or turbulent and the mechanisms underlying the two have been proposed to be different (Dandois et al., J. Fluid Mech., vol. 18, 2018, pp. 156–178). The effect of various flow parameters (freestream Mach and Reynolds numbers and sweep and incidence angles) on laminar transonic buffet on an infinite wing (Dassault Aviation's supercritical V2C aerofoil) is reported here by performing large-eddy simulations (LES) for a wide range of parameters. A spectral proper orthogonal decomposition identified the presence of a low-frequency mode associated with buffet and high-frequency wake modes related to vortex shedding. A flow reconstruction based only on the former shows periodic boundary-layer separation and reattachment accompanying shock wave motion. A modal reconstruction based only on the wake mode suggests that the separation bubble breathing phenomenon reported by Dandois et al. is due to this mode. Together, these results indicate that the physical mechanisms governing laminar and turbulent buffet are the same. Buffet was also simulated at zero incidence. Shock waves appear on both aerofoil surfaces and oscillate out of phase with each other indicating the occurrence of a Type I buffet (Giannelis et al., Aerosp. Sci. Technol., vol. 18, 2018, pp. 89–101) on a supercritical aerofoil. These results suggest that the mechanisms underlying different buffet types are the same.