We study experimentally the dynamics of a water droplet on a tilted and vertically oscillating rigid fibre. As we vary the frequency and amplitude of the oscillations the droplet transitions between different modes: harmonic pumping, subharmonic pumping, a combination of rocking and pumping modes, and a combination of pumping and swinging modes. We characterize these responses and report how they affect the sliding speed of the droplet along the fibre. The swinging mode is explained by a minimal model making an analogy between the droplet and a forced elastic pendulum.

We study numerically the characteristics of fluid–fluid displacement in simple mixed-wet porous micromodels using a dynamic pore network model. The porous micromodel consists of distinct water-wet and oil-wet regions, whose fractions are varied systematically to yield a variety of displacement patterns over a wide range of capillary numbers. We find that the impact of mixed-wettability is most prominent at low capillary numbers, and it depends on the complex interplay between wettability fraction and the intrinsic contact angle of the water-wet regions. For example, the fractal dimension of the displacement pattern is a monotonically increasing function of wettability fraction in flow cells with strongly water-wet clusters, but it becomes non-monotonic with respect to wettability fraction in flow cells with weakly water-wet clusters. Additionally, mixed-wettability also manifests itself in the injection pressure signature, which exhibits fluctuations especially at low wettability fraction. Specifically, preferential filling of water-wet regions leads to reduced effective permeability and higher injection pressure, even at vanishingly small capillary numbers. Finally, we demonstrate that scaling analyses based on a weighted average description of the overall wetting state of the mixed-wet system can effectively capture the variations in observed displacement pattern morphology.

We study the onset of the three-dimensional mode A instability in the near wake behind a circular cylinder. We show that long-wavelength perturbations organise in a time-shifting pattern such that the in-plane velocity in each streamwise slice corresponds to the base flow solution at shifted times. This observation introduces an additional unifying characteristic for certain mode A type instabilities. We then analyse the mechanisms which control the growth or decay of these perturbations and highlight the crucial role played by the tilting mechanism which operates via non-local interactions in a manner similar to Biot–Savart induction. We characterise its domain of influence using a Green's function-based approach which allows us to rationalise the non-trivial dependence of the growth rate on the spanwise wavenumber. We connect this behaviour to the subtle balance between the local growth of the perturbations as they are swept along by the flow and the feedback on the perturbations that are generated during the next period of the time-periodic base flow. Finally, we discuss generalisations of our findings to other types of flows.

The traditional approximation neglects the cosine components of the Coriolis acceleration, and this approximation has been widely used in the study of geophysical phenomena. However, the justification of the traditional approximation is questionable under a few circumstances. In particular, dynamics with substantial vertical velocities or geophysical phenomena in the tropics have non-negligible cosine Coriolis terms. Such cases warrant investigations with the non-traditional setting, i.e. the full Coriolis acceleration. In this manuscript, we study the effect of the non-traditional setting on an isothermal, hydrostatic and compressible atmosphere assuming a meridionally homogeneous flow. Employing linear stability analysis, we show that, given appropriate boundary conditions, i.e. a bottom boundary condition that allows for a vertical energy flux and non-reflecting boundary at the top, the atmosphere at rest becomes prone to a novel unstable mode. The validity of assuming a meridionally homogeneous flow is investigated via scale analysis. Numerical experiments were conducted, and Rayleigh damping was used as a numerical approximation for the non-reflecting top boundary. Our three main results are as follows: (i) experiments involving the full Coriolis terms exhibit an exponentially growing instability, yet experiments subjected to the traditional approximation remain stable; (ii) the experimental instability growth rate is close to the theoretical value; (iii) a perturbed version of the unstable mode arises even under sub-optimal bottom boundary conditions. Finally, we conclude our study by discussing the limitations, implications and remaining open questions. Specifically, the influence on numerical deep-atmosphere models and possible physical interpretations of the unstable mode are discussed.

Wave interaction with graded metamaterials exhibits the phenomenon of rainbow reflection, in which broadband wave signals slow down and separate into their frequency components before being reflected. This phenomenon has been qualitatively understood by describing the wave field in the metamaterial using the local Bloch wave approximation (LBWA), which locally represents the wave field as a superposition of propagating wave solutions in the cognate infinite periodic media (so-called Bloch waves). We evaluate the performance of the LBWA quantitatively in the context of two-dimensional linear water-wave scattering by graded arrays of surface-piercing vertical barriers. To do this, we implement the LBWA numerically so that the Bloch waves in one region of the graded array are coupled to Bloch waves in adjacent regions. This coupling is computed by solving the scattering of Bloch waves across the interface between two semi-infinite arrays of vertical barriers, where the barriers in each semi-infinite array can have different submergence depths. Our results suggest that the LBWA accurately predicts the free surface amplitude across a wide range of frequencies, except those just above the cutoff frequencies associated with each of the vertical barriers in the array. This highlights the importance of decaying Bloch modes above the cutoff in rainbow reflection.

Compound surfaces, consisting of periodic arrays of solid patches and free surfaces, exhibit hydrodynamic slipperiness which is quantified by their slip length. The limit of small solid fractions, where the slip length diverges, is of fundamental interest. This paper addresses longitudinal liquid flows over a periodic array of grooves which are partially invaded by the liquid. Assuming that the slats separating the grooves are infinitely thin, the solid fraction $\epsilon$ is set by the invasion depth. Inspired by the singular small-solid-fraction limit for non-invaded grooves (Schnitzer, Phys. Rev. Fluids, vol. 1, 2016, 052101R), we consider the idealised geometry of $90^{\circ }$ protrusion angles, where an integral force balance in the limit $\epsilon \to 0$ implies a slip length that scales as $\epsilon ^{-1}$. The problem exhibits a nested structure, where the liquid domain is conceptually decomposed into four distinct regions: a unit-cell region on the scale of the period, where the wetted portion of the slat appears as a point singularity; two regions on the scale of the wetted slat, where the flow essentially varies in one dimension; and a transition region about the tip of the slat. Analysing these regions using matched asymptotic expansions and conformal mappings yields the ratio of slip length to semi-period as $2\epsilon ^{-1} - (10/{\rm \pi} )\ln 2 + \cdots$.

Axisymmetric steady solutions of Taylor–Couette flow at high Taylor numbers are studied numerically and theoretically. As the axial period of the solution shortens from approximately one gap length, the Nusselt number goes through two peaks before returning to laminar flow. In this process, the asymptotic nature of the solution changes in four stages, as revealed by the asymptotic analysis. When the aspect ratio of the roll cell is approximately unity, the solution has the Nusselt number proportional to the quarter power of the Taylor number, and captures quantitatively the characteristics of the classical turbulence regime. By shortening the axial period the Nusselt number can even reach the experimental value around the onset of the ultimate turbulence regime. However, at higher Taylor numbers, the theoretical predictions eventually underestimate the experimental values. An important consequence of the asymptotic analyses is that the mean angular momentum should become uniform in the core region unless the axial wavelength is too short. The theoretical scaling laws deduced for the steady solutions can be carried over to Rayleigh–Bénard convection.

We demonstrate that the relative motion of horizontal parallel plates can be generated using patterned heating. This movement is driven by nonlinear thermal streaming associated with a pitchfork bifurcation. The propulsive effect is strongest when all the heating energy is concentrated in a single Fourier mode of the spatial heating pattern; it increases with a decrease in the Prandtl number and increases with the addition of a uniform heating component.

The formation of swirling vortex rings and their early time evolution, resulting from the controlled discharge of an incompressible, Newtonian fluid into a stationary equivalent fluid bulk, is explored for weak to moderate swirl number $S \in [0, 1]$. Two practically realisable inlet conditions are investigated with swirl simultaneously superposed onto a linear momentum discharge; the corresponding circulation based Reynolds number is 7500. The results obtained reveal that for $S > 1/2$, the addition of swirl promotes the breakdown of the leading primary vortex ring structure, giving rise to the striking feature of significant negative azimuthal vorticity generation in the region surrounding the primary vortex ring core, whose strength scales with ${S}^2$. Through a nonlinear interaction with the vortex breakdown, the radius of the primary toroidal vortex core is rapidly increased; consequently, the self-induced propagation velocity of the leading ring decreases with $S$ and vortex stretching along the circular primary vortex core increases counteracting viscous diffusion effects. The latter governs the evolution of the peak vorticity intensity and the swirl velocity magnitude in the primary ring core, the circulation growth rate of the primary ring, as well as the vorticity intensity of the trailing jet and hence its stability. This combination of effects leads to an increased dimensionless kinetic energy for the primary ring with increasing $S$ and results in an almost linearly decreasing circulation based formation number, $F$.

The log law of the wall, joining the inner, near-wall mean velocity profile (MVP) in wall-bounded turbulent flows to the outer region, has been a permanent fixture of turbulence research for over hundred years, but there is still no general agreement on the value of the prefactor, the inverse of the Kármán ‘constant’ $\kappa$, or on its universality. The choice diagnostic tool to locate logarithmic parts of the MVP is to look for regions where the indicator function $\varXi$ (equal to the wall-normal coordinate $y^+$ times the mean velocity derivative $\mathrm {d} U^+/\mathrm {d} y^+$) is constant. In pressure-driven flows, however, such as channel and pipe flows, $\varXi$ is significantly affected by a term proportional to the wall-normal coordinate, of order $O({Re}_{\tau }^{-1})$ in the inner expansion, but moving up across the overlap to the leading $O(1)$ in the outer expansion. Here we show that, due to this linear overlap term, ${Re}_{\tau }$ values well beyond $10^5$ are required to produce one decade of near constant $\varXi$ in channels and pipes. The problem is resolved by considering the common part of the inner asymptotic expansion carried to $O({Re}_{\tau }^{-1})$, and the leading order of the outer expansion. This common part contains a superposition of the log law and a linear term $S_0 \,y^+{Re}_{\tau }^{-1}$, and corresponds to the linear part of $\varXi$, which, in channel and pipe, is concealed up to $y^+ \approx 500\unicode{x2013}1000$ by terms of the inner expansion. A new and robust method is devised to simultaneously determine $\kappa$ and $S_0$ in pressure-driven flows at currently accessible ${Re}_{\tau }$ values, yielding $\kappa$ values which are consistent with the $\kappa$ values deduced from the Reynolds number dependence of centreline velocities. A comparison with the zero-pressure-gradient turbulent boundary layer, further clarifies the issues and improves our understanding.

We examine the equilibrium configurations of a nematic liquid crystal with an immersed body in two dimensions. A complex variables formulation provides a means for finding analytical solutions in the case of strong anchoring. Local tractions, forces and torques on the body are discussed in a general setting. For weak (finite) anchoring strengths, an effective boundary technique is proposed which is used to determine asymptotic solutions. The energy-minimizing locations of topological defects on the body surface are also discussed. A number of examples are provided, including circular and triangular bodies, and a Janus particle with hybrid anchoring conditions. Analogies to classical results in fluid dynamics are identified, including d'Alembert's paradox, Stokes’ paradox and the Kutta condition for circulation selection.

In this paper, we study the linear stability of a two-dimensional shear flow of an air layer overriding a water layer of finite depth. The air layer is considered to be of an infinite extent with an exponential velocity profile. Three different background conditions are considered in the finite-depth water layer: a quiescent background, a linear velocity profile and a quadratic velocity profile. It is known that the cases of the quiescent water layer and the linear velocity profile allow for analytical treatment. We further provide an analytical solution for the case of the quadratic velocity field: we specifically consider a flow-reversal profile, although the result could be generalized to other quadratic profiles as well. The role of water layer depth on the growth rate of the Miles and rippling instabilities is studied in each of the three cases. Using asymptotic analysis, with the air–water density ratio being a small parameter, we obtain an analytical expression for the growth rate of the Miles mode and discuss the condition for the existence of a long-wave cutoff for these profiles. We provide analytical expressions for the stability boundary in the parameter space of inverse squared Froude number and wavenumber. In scenarios where a long-wave cutoff does not exist, we have carried out a long-wave asymptotic study to obtain the growth rate behaviour in that regime.

The present paper numerically investigates the fluid dynamics associated with the flow over an elastically mounted circular cylinder at a Reynolds number of 500. A slit normal to the flow direction is placed at the cylinder's centre. The study covers the large parametric set of calculations for the combined mass-damping ratio $m^{*}\zeta = 0.2$, including the slit offset for six different offset angles ($-30^\circ \leq \alpha \leq +30^\circ$) measured from the vertical axis at the centre point of the cylinder and the effect of slit widths ($0.1 \leq s/D \leq 0.3$) on the aerodynamic loading, vibration response and associated flow characteristics. Furthermore, a wide range of reduced velocities ($3\leq U_r \leq 7$) are examined for the complete closure of studying the effect on the vortex-induced vibration response. The results demonstrate that adding the normal slit increases the periodic suction-blowing phenomena, strengthening the vortex shedding. The results suggest that the normal slit can suppress or increase vortex-induced vibration depending on the slit-offset angle. Placing it toward the front stagnation point results in the increased oscillation amplitude (with a wider wake behind the cylinder) while shifting it toward the rear stagnation point diminishes the cylinder's oscillations. The paper reveals that this dual behaviour of the normal slit (based on its offset placement) is closely related to the phase difference between the lift force and the oscillation amplitude. Various vortex-shedding patterns associated with the different slit-offset angles are duly reported in the paper. Furthermore, a magnet is attached to the slit cylinder, and its effect on the energy harvesting capability via a coil-magnet arrangement is explored in detail for a wide range of reduced velocities.

An exact analytical solution is obtained for the dynamical system derived in Part 1 of this series (Moffatt & Kimura, J. Fluid Mech, vol. 861, 2019a, pp. 930–967), which describes the approach of two initially circular vortices of finite but small cross-section symmetrically located on inclined planes. This exact solution, applicable in the inviscid limit, allows determination of the amplification $\mathcal {A}_{\omega }$ of the axial vorticity within the finite time $T$ during which the basic assumptions of the model continue to apply. It is first shown that, for arbitrarily prescribed $\mathcal {A}_{\omega }$, it is possible to specify smooth initial conditions of finite energy such that, in the inviscid limit, this amplification is achieved within the time $T$. When viscosity is included, an estimate is provided for the minimum vortex Reynolds number that is sufficient for the same result to hold. The predictions are broadly compatible with results from direct numerical simulations at moderate Reynolds numbers. Moreover, it is shown that one may come arbitrarily close to a finite-time singularity of the Navier–Stokes equation by appropriate choice of an initial, smooth, finite-energy velocity field; however, this approach to a singularity is ultimately thwarted through breach of the assumptions on which the dynamical system is based. Thus we make no claim here concerning realisation of a Navier–Stokes singularity. Moreover, we find that the conditions required to attain a large amplification $\mathcal {A}_{\omega }\gg 1$ during the time $T$ are far beyond those that can be realised in either experiment or direct numerical simulation.

Shock buffet on wings is a phenomenon caused by strong shock-wave/boundary-layer interaction resulting first in self-sustained flow unsteadiness and eventually in a detrimental structural response called buffeting. While it is an important aspect of wing design and aircraft certification, particularly for modern transonic air transport, not all of the underlying multidisciplinary physics is thoroughly understood. Building upon a single-discipline shock-buffet stability study, this work now investigates the impact of an elastic structure in these extreme flow conditions. Specifically, a triglobal stability analysis of a fluid–structure coupled system is presented, utilising the implicitly restarted Arnoldi method with a sparse iterative Krylov solver and novel preconditioner. Asymmetry resulting from a static aeroelastic simulation based on a finite-element model of the underlying geometry in a wind tunnel modifies the global modes of the earlier fluid-only symmetric full-span analysis. A flutter stability analysis at wind-tunnel flow conditions below shock-buffet onset finds no instability in the structural degrees-of-freedom, whereas in shock-buffet flow with globally unstable fluid modes additional marginally unstable structural (and fluid) modes emerge. The developed stability tool for coupled analysis is instrumental in identifying those physically relevant and strongly coupled modes where a standard pk-type (p being eigenvalue and k reduced frequency) flutter analysis fails. With the complementary computation of adjoint eigenmodes, the core of the instability is pinpointed to a relatively small wing area which may help to effect the control and delay of this detrimental transonic unsteadiness. We contribute to the question on how the presence of the elastic wing structure impacts on the otherwise pure aerodynamic three-dimensional shock-buffet dynamics.

Power minimisation of fluid transport in branched fluidic networks has become of paramount importance for microfluidics, additive manufacturing and hierarchical functional materials. For fully developed laminar flow of Newtonian fluids, Murray's theory provides a solution for the channel and network dimensions that minimise power consumption. However, design and optimisation of networks that transport complex fluids is still challenging. Here, we generalise Murray's theory towards fluid rheologies, including non-Newtonian (power-law) and yield-stress fluids (Bingham, Herschel–Bulkley, Casson). A straightforward graphical approach is presented that provides the optimal radii in a branching network, and the angles between these branches. The wall shear stress is found to be uniform over the entire network, and the velocity profile is self-similar. Furthermore, the effect of non-optimal channel radii on the power consumption of the network is investigated. Finally, examples illustrate how this approach applies to a wide variety of systems.

The development of novel drug delivery systems, which are revolutionizing modern medicine, is benefiting from studies on microorganisms’ swimming. In this paper we consider a model microorganism (a squirmer) enclosed in a viscous droplet to investigate the effects of medium heterogeneity or geometry on the propulsion speed of the caged squirmer. We first consider the squirmer and droplet to be spherical (no shape effects) and derive exact solutions for the equations governing the problem. For a squirmer with purely tangential surface velocity, the squirmer is always able to move inside the droplet (even when the latter ceases to move as a result of large fluid resistance of the heterogeneous medium). Adding radial modes to the surface velocity, we establish a new condition for the existence of a co-swimming speed (where squirmer and droplet move at the same speed). Next, to probe the effects of geometry on propulsion, we consider the squirmer and droplet to be in Newtonian fluids. For a squirmer with purely tangential surface velocity, numerical simulations reveal a strong dependence of the squirmer's speed on shapes, the size of the droplet and the viscosity contrast. We found that the squirmer speed is largest when the droplet size and squirmer's eccentricity are small, and the viscosity contrast is large. For co-swimming, our results reveal a complex, non-trivial interplay between the various factors that combine to yield the squirmer's propulsion speed. Taken together, our study provides several considerations for the efficient design of future drug delivery systems.

Boundary-layer flow over a realistic porous wall might contain both the effects of wall-permeability and wall-roughness. These two effects are typically examined in the context of a rough-wall flow, i.e. by defining a ‘roughness’ length or equivalent to capture the effect of the surface on momentum deficit/drag. In this work, we examine the hypothesis of Esteban et al. (Phys. Rev. Fluids, vol. 7, no. 9, 2022, 094603), that a turbulent boundary layer over a porous wall could be modelled as a superposition of the roughness effects on the permeability effects by using independently obtained information on permeability and roughness. We carry out wind tunnel experiments at high Reynolds number ($14\ 400 \leq Re_{\tau } \leq 33\ 100$) on various combinations of porous walls where different roughnesses are overlaid over a given permeable wall. Measurements are also conducted on the permeable wall as well as the rough walls independently to obtain the corresponding length scales. Analysis of mean flow data across all these measurements suggests that an empirical formulation can be obtained where the momentum deficit ($\Delta U^+$) is modelled as a combination of independently obtained roughness and permeability length scales. This formulation assumes the presence of outer-layer similarity across these different surfaces, which is shown to be valid at high Reynolds numbers. Finally, this decoupling approach is equivalent to the area-weighted power-mean of the respective permeability and roughness length scales, consistent with the approach recently suggested by Hutchins et al. (Ocean Engng, vol. 271, 2023, 113454) to capture the effects of heterogeneous rough surfaces.

The propagation of a finite blister of fluid beneath an elastic sheet is controlled by the local dynamics at the peeling periphery of the current. Previous works have described constant volume elastic blisters by considering the peeling due to curvature at these edges. Here, we show that an along-slope component of gravity fundamentally changes the dynamics by removing the role of curvature at the trailing, upslope edge. The local dynamics of this trailing edge is instead controlled by shear stress in the sheet, as in the elastic Landau–Levich problem, and thereby allows for a receding edge, in contrast to propagation by peeling for which only an advancing contact line is possible. Using an asymptotic analysis, we show that this receding edge condition allows for a new, nearly translating regime in which the body of the blister moves at an approximately constant speed, leaving behind a thin layer of fluid. This prediction is verified by detailed numerical modelling of the two-dimensional downslope spreading. We conclude by discussing the applicability of these results in the rapid spreading of subglacial meltwater.

A spectral method using associated Legendre functions with algebraic mapping is developed for a linear stability analysis of wake vortices. These functions serve as Galerkin basis functions, capturing correct analyticity and boundary conditions for vortices in an unbounded domain. The incompressible Euler or Navier–Stokes equations linearised on a swirling flow are transformed into a standard matrix eigenvalue problem of toroidal and poloidal streamfunctions, solving perturbation velocity eigenmodes with their complex growth rate as eigenvalues. This reduces the problem size for computation and distributes collocation points adjustably clustered around the vortex core. Based on this method, strong swirling $q$ vortices with linear perturbation wavenumbers of order unity are examined. Without viscosity, neutrally stable eigenmodes associated with the continuous eigenvalue spectrum having critical-layer singularities are successfully resolved. The inviscid critical-layer eigenmodes numerically tend to appear in pairs, implying their singular degeneracy. With viscosity, the spectra pertaining to physical regularisation of critical layers stretch out toward an area, referring to potential eigenmodes with wavepackets found by Mao & Sherwin (J. Fluid Mech., vol. 681, 2011, pp. 1–23). However, the potential eigenmodes exhibit no spatial similarity to the inviscid critical-layer eigenmodes, doubting that they truly represent the viscous remnants of the inviscid critical-layer eigenmodes. Instead, two distinct continuous curves in the numerical spectra are identified for the first time, named the viscous critical-layer spectrum, where the similarity is noticeable. Moreover, the viscous critical-layer eigenmodes are resolved in conformity with the $Re^{-1/3}$ scaling law. The onset of the two curves is believed to be caused by viscosity breaking the singular degeneracy.