Waves are formed on the surface of a sessile drop driven through substrate vibrations oriented at a slanting angle from the normal. A mathematical model is derived, which leads to an infinite system of coupled Mathieu equations governing the wave dynamics that are solved using Floquet theory. The spatial structure of the waves is described by the mode number pair $[\ell,m]$, where $\ell$ and $m$ are the polar and azimuthal mode numbers, respectively. Limiting cases corresponding to horizontal and vertical vibrations are discussed with predictions agreeing well with prior literature. We focus our results on three drop motions – (1) harmonic $[1,1]$ rocking mode, (2) harmonic $[2,0]$ pumping mode, and (3) subharmonic rocking $[1,1]$ mode – as they depend upon the slanting angle, static contact angle, and contact-line conditions, which we assume to be either pinned or freely moving with fixed contact angle. New theoretical predictions are tested through experiments over a range of parameters, showing good agreement.

In the absence of large-scale coherent structures, a widely used statistical theory of two-dimensional turbulence developed by Kraichnan, Leith, and Batchelor (KLB) predicts a power-law scaling for the energy, $E(k)\propto k^\alpha$ with an integral exponent $\alpha ={-3}$, in the inertial range associated with the direct cascade. A power-law scaling is also observed in the presence of coherent structures, but the scaling exponent becomes fractal and often differs substantially from the value predicted by the KLB theory. Here we present a dynamical theory that sheds new light on the relationship between the spatial and temporal structure of the large-scale flow and the scaling of small-scale structures representing filamentary vorticity. Specifically, we find hyperbolic regions of the large-scale flow to play a key role in the flux of enstrophy between scales. Small-scale vorticity in these regions can be described by dynamically self-similar solutions of the Euler equation, which explains the power-law scaling. Furthermore, we find that correlations between different hyperbolic regions are responsible for the emergence of fractal scaling exponents.

The first experimental results on pattern transitions in the co-rotation regime (i.e. the rotation ratio $\varOmega = \omega _o/\omega _i > 0$, where $\omega _i$ and $\omega _o$ are the angular speeds of the inner and outer cylinders, respectively) of the Taylor–Couette flow (TCF) are reported for a neutrally buoyant suspension of non-colloidal particles, up to a particle volume fraction of $\phi = 0.3$. While the stationary Taylor vortex flow (TVF) is the primary bifurcating state in dilute suspensions ($\phi \leq ~0.05$), the non-axisymmetric oscillatory states, such as the spiral vortex flow (SVF) and the ribbon (RIB), appear as primary bifurcations with increasing particle loading, with an overall de-stabilization of the primary bifurcating states (TVF/SVF/RIB) being found with increasing $\phi$ for all $\varOmega \geq ~0$. At small co-rotations ($\varOmega \sim 0$), the particles play the dual role of stabilization ($\phi < 0.1$) and destabilization ($\phi \geq ~0.1$) on the secondary/tertiary oscillatory states. The distinctive features of the ‘particle-induced’ spiral vortices are identified and contrasted with those of the ‘fluid-induced’ spirals that operate in the counter-rotation regime.

The asymptotic analysis of steady azimuthally invariant electromagnetically driven flows occurring in a shallow annular layer of electrolyte undertaken in Part 1 of this study (McCloughan & Suslov, J. Fluid Mech., vol. 980, 2024, A59) predicted the existence of a two-tori flow state that has not been detected previously. In Part 2 of the study we confirm its existence by numerical time integration of the governing equations. We observe a hysteresis, where the type of solution obtained for the same set of governing parameters depends on the choice of the initial conditions and the way the governing parameters change, which is fully consistent with the analytic results of Part 1. Subsequently, we perform a linear stability analysis of the newly obtained steady state and deduce that the experimentally observed anti-cyclonic free-surface vortices appear on its background as a result of a centrifugal (Rayleigh-type) instability of the interface separating two counter-rotating toroidal structures that form the newly found flow solution. The quantitative characteristics of such instability structures are determined. It is shown that such structures can only exist in sufficiently thin layers with the depth not exceeding a certain critical value.

Acoustic resonance is an important factor that contributes to aeroengine compressor failure. In this study, a plane cascade of compressor blades was designed to reproduce acoustic resonance via a low-speed wind tunnel test. A high-frequency hot-wire, microphone and strain gauge were used to synchronously measure the fluid, acoustic and structural parameters. We analysed the variation in the amplitude and frequency of the multi-field parameters with increasing mean flow velocity and explored the multi-field interaction mechanism that induces the acoustic resonance of the plane cascade. The plane cascade effectively reproduced the acoustic resonance phenomenon. The first-order acoustic-mode frequency of the plane cascade flow duct, second-order torsional vibration mode frequency of the blade and shedding mode frequency of the tip clearance leakage vortex were equal under acoustic resonance. The fluid, acoustic and structural fields showed a strong interaction effect, achieving the maximum blade vibration amplitude and causing fatigue cracks of torsional vibration at the blade root. The frequency lock-in region of the compressor plane cascade was divided into an ‘acoustic–structure’ interaction region, a ‘fluid–acoustic–structure’ interaction region and a first-order acoustic-mode dominant region with increasing mean flow velocity, which demonstrates an interesting phenomenon in which the fluid–acoustic–structure modes compete: acoustic mode > blade vibration mode > vortex shedding mode. The results demonstrate a unique approach to the study of acoustic resonance that provides insight into the acoustic resonance mechanism in a cascade of compressor blades.

During stroke reversals, insect wings interact with their own wake flow from the preceding half-stroke, resulting in an unsteady aerodynamic mechanism known as ‘wing–wake interaction’ or ‘wake capture’. To better elucidate this mechanism, we numerically solved the incompressible Navier–Stokes equations at Reynolds numbers $10^2$ and $10^3$. Simulations were conducted for wing planforms defined using the beta function distribution with varying aspect ratios ($AR=2\unicode{x2013}6$) and radial centroid locations ($\hat {r}_1=0.4\unicode{x2013}0.6$), whilst employing representative normal hovering kinematics. The wake development from the considered flapping wing planforms was investigated, and the wake capture contribution to aerodynamic force production was quantified by comparing the force generation between the fifth and first stroke cycles at multiple sections along the wingspan. Our results revealed that on the inboard wing region experiencing an attached leading-edge vortex (LEV) structure, wing–wake interaction is dominated by an unsteady downwash effect, resulting in a reduction in local force production. However, in regions closer to the wingtip experiencing detachment of the LEV, wing–wake interaction is dominated by an unsteady upwash effect, leading to an increase in local force production. Consequently, the global wake capture force production is controlled by the extent of LEV detachment, which primarily increases with the increase of wing aspect ratio. This suggests that for normal hovering flapping wings, the typical loss in translational force production due to wingtip stall is partially mitigated by wake capture effects.

We consider steady surface waves in an infinitely deep two-dimensional ideal fluid with potential flow, focusing on high-amplitude waves near the steepest wave with a 120$^{\circ }$ corner at the crest. The stability of these solutions with respect to coperiodic and subharmonic perturbations is studied, using new matrix-free numerical methods. We provide evidence for a plethora of conjectures on the nature of the instabilities as the steepest wave is approached, especially with regards to the self-similar recurrence of the stability spectrum near the origin of the spectral plane.

The background potential energy (BPE) is the only reservoir that double diffusive instabilities can tap their energy from when developing from an unforced motionless state with no available potential energy (APE). Recently, Middleton and Taylor linked the extraction of BPE into APE to the sign of the diapycnal component of the buoyancy flux, but their criterion can predict only diffusive convection instability, not salt finger instability. Here, we show that the problem can be corrected if the sign of the APE dissipation rate is used instead, making it emerge as the most fundamental criterion for double diffusive instabilities. A theory for the APE dissipation rate for a two-component fluid relative to its single-component counterpart is developed as a function of three parameters: the diffusivity ratio, the density ratio, and a spiciness parameter. The theory correctly predicts the occurrence of both salt finger and diffusive convection instabilities in the laminar unforced regime, while more generally predicting that the APE dissipation rate for a two-component fluid can be enhanced, suppressed, or even have the opposite sign compared to that for a single-component fluid, with important implications for the study of ocean mixing. Because negative APE dissipation can also occur in stably stratified single-component and doubly stable two-component stratified fluids, we speculate that only the thermodynamic theory of exergy can explain its physics; however, this necessitates accepting that APE dissipation is a conversion between APE and the internal energy component of BPE, in contrast to prevailing assumptions.

Recently, subcritical transition to turbulence in the quasi-two-dimensional (quasi-2-D) shear flow with strong linear friction (Camobreco et al., J. Fluid Mech., vol. 963, 2023, R2) has been demonstrated by the 2-D mechanism at $Re = 71\,211$, and the nonlinear Tollmien–Schlichting (TS) waves related to the edge state were approached independently of initial optimal disturbances. For 2-D plane Poiseuille flow, transition to the fully developed turbulence requires that the Reynolds number is several times larger than the critical Reynolds number $Re_c$ (Markeviciute & Kerswell, J. Fluid Mech., vol. 917, 2021, A57). In this paper, we observed the subcritical transitional flow in 2-D plane Poiseuille flow driven by the nonlinear TS waves by both linear and nonlinear optimal disturbances ($Re < Re_c$) with different quantitative edge states. The nonlinear optimal disturbances could trigger the sustained subcritical transitional flow for $Re \geqslant 2400$. The initial energy for nonlinear optimal disturbance is more efficient than the linear optimal disturbance in reaching the subcritical transitional flow for $2400 \leqslant Re \leqslant 5000$. Moreover, the initial energy of linear optimal disturbance is larger than the energy of its edge state. The nonlinear TS waves along the edge state are formed by the nonlinear optimal disturbances to trigger transitional flow, which agrees well with the main conclusions of Camobreco et al. (J. Fluid Mech., vol. 963, 2023, R2), while the required $Re$ of 2-D plane Poiseuille flow is much smaller.

Resistive tearing instabilities are common in fluids that are highly electrically conductive and carry strong currents. We determine the effect of stable stratification on the tearing instability under the Boussinesq approximation. Our results generalise previous work that considered only specific parameter regimes, and we show that the length scale of the fastest-growing mode depends non-monotonically on the stratification strength. We confirm our analytical results by solving the linearised equations numerically, and we discuss whether the instability could operate in the solar tachocline.

Recent years have seen the emergence of new technologies that exploit nanoscale evaporation, ranging from nanoporous membranes for distillation to evaporative cooling in electronics. Despite the increasing depth of fundamental knowledge, there is still a lack of simulation tools capable of capturing the underlying non-equilibrium liquid–vapour phase changes that are critical to these and other such technologies. This work presents a molecular kinetic theory model capable of describing the entire flow field, i.e. the liquid and vapour phases and their interface, while striking a balance between accuracy and computational efficiency. In particular, unlike previous kinetic models based on the isothermal assumption, the proposed model can capture the temperature variations that occur during the evaporation process, yet does not require the computational resources of more complicated mean-field kinetic approaches. We assess the present kinetic model in three test cases: liquid–vapour equilibrium, evaporation into near-vacuum condition, and evaporation into vapour. The results agree well with benchmark solutions, while reducing the simulation time by almost two orders of magnitude on average in the cases studied. The results therefore suggest that this work is a stepping stone towards the development of an accurate and efficient computational approach to optimising the next generation of nanotechnologies based on nanoscale evaporation.

Whether colliding drops will merge with or bounce off each other is critical to numerous processes, and the physics involved is notoriously complex. In particular, experiments show that both sufficiently slow and fast head-on drop collisions lead to merging, but that there is often an intermediate regime in which bouncing is observed; these transitions in behaviour were recently discovered to be surprisingly sensitive to the radius of the drops and the ambient gas pressure. We show here that these transitions between bouncing and merging are governed by nanoscale phenomena; namely, gas-kinetic and disjoining pressure effects. To capture these crucial effects, a novel, open-source computational model is developed for the simulation of colliding drops. The model uses a hybrid approach, based on solving the Navier–Stokes equations in the drop with a lubrication approach for the unconventional physics of the gas film. Our simulations show remarkably good agreement with experiments of head-on collisions and also provide new experimentally verifiable predictions.

A high-resolution numerical simulation of an air–water turbulent upward bubbly flow in a pipe is performed to investigate the turbulence characteristics and bubble interaction with the wall. We consider three bubble equivalent diameters and three total bubble volume fractions. The bulk and bubble Reynolds numbers are $Re_{bulk}= u_{bulk} D/\nu _w = 5300$ and $Re_{bub}= (\langle u_{bub}\rangle - u_{bulk}) d_{eq}/\nu _w = 533\unicode{x2013}1000$, respectively, where $u_{bulk}$ is the water bulk velocity, $\langle u_{bub}\rangle$ is the overall bubble mean velocity, $D$ is the pipe diameter and $\nu _w$ is the water kinematic viscosity. The mean water velocity near the wall significantly increases due to bubble interaction with the wall, and the root-mean-square water velocity fluctuations are proportional to $\bar {\psi }(r)^{0.4}$, where $\bar {\psi } (r)$ is the mean bubble volume fraction. For the cases considered, the bubble-induced turbulence suppresses the shear-induced turbulence and becomes the dominant flow characteristic at all radial locations including near the wall. Rising bubbles near the wall mostly bounce against the wall rather than slide along the wall or hang around the wall without collision. Low-speed streaks observed in the near-wall region in the absence of bubbles nearly disappear due to the bouncing bubbles. These bouncing bubbles generate counter-rotating vortices in their wake, and increase the skin friction by sweeping high-speed water towards the wall. We also suggest an algebraic Reynolds-averaged Navier–Stokes model considering the interaction between shear-induced and bubble-induced turbulence. This model provides accurate predictions for a wide range of liquid bulk Reynolds numbers.