In their transition from a laminar state to turbulence, some unstable flows pass through a set of well-defined stages involving different and distinct processes. This is so, in particular, for Kelvin–Helmholtz instability, although details of its transition still retain many mysterious aspects. Billows develop in the primary stage of this stratified shear flow instability, separated by thin braids in which the shear is relatively high. Fluid is statically unstable within the billows and consequently potentially prone to convective instability. Numerical studies by Mashayek & Peltier (J. Fluid Mech., this issue, vol 708, 2012a,b, pp. 5–44 and 45–70) have discovered several new types of secondary instability in the braids and billow cores that may hasten the eventual transition to turbulence. The instabilities are illustrated by the authors in colour figures, remarkable for their beauty and (recalling William Blake’s ‘The Tyger’) ‘fearful symmetry’. But are they helpful in establishing the subsequent turbulence in the natural environment?

We study the competition between various secondary instabilities that co-exist in a preturbulent stratified parallel flow subject to Kelvin–Helmholtz instability. In particular, we investigate whether a secondary braid instability might emerge prior to the overturning of the statically unstable regions that develop in the cores of the primary Kelvin–Helmholtz billows. We identify two groups of instabilities on the braid. One group is a shear instability which extracts its energy from the background shear and is suppressed by the straining contribution of the background flow. The other group, which seems to have no precedent in the literature, includes phase-locked modes which grow at the stagnation point on the braid and are almost entirely driven by the straining contributions of the background flow. For the latter group, the braid shear has a negative influence on the growth rate. Our analysis demonstrates that the probability of finite-amplitude growth of both braid instabilities is enhanced with increasing Reynolds number and Richardson number. We also show that the possibility of emergence of braid instabilities decreases with the Prandtl number for the shear modes and increases for the stagnation point instabilities. Through detailed non-separable linear stability analysis, we show that both braid instabilities are fundamentally three dimensional with the shear modes being of small wavenumbers and the stagnation point modes dominating at large wavenumber.

The linear stability analyses described in Mashayek & Peltier (J. Fluid Mech., vol. 708, 2012, 5–44, hereafter MP1) are extended herein in an investigation of the influence of stratification on the evolution of secondary instabilities to which an evolving Kelvin–Helmholtz (KH) wave is susceptible in an initially unstable parallel stratified shear layer. We show that over a wide range of background stratification levels, the braid shear instability has a higher probability of emerging at early stages of the flow evolution while the secondary convective instability (SCI), which occurs in the eyelids of the individual Kelvin ‘cats eyes’, will remain a relevant and dominant instability at high Reynolds numbers. The evolution of both modes is greatly influenced by the background stratification. Various other three-dimensional secondary instabilities are found to exist over a wide range of stratification levels. In particular, the stagnation point instability (SPI), which was discussed in detail in MP1, may be of great potential importance providing alternate routes for transition of an initially two-dimensional KH wave into fully developed turbulence. The energetics of the secondary instabilities revealed by our simulations are analysed in detail and the preturbulent mixing properties are studied.

We extend the classical energy criterion for stability, the Lagrange–Dirichlet theorem, to rotating non-smooth complex fluids. The stability test so developed is very general and may be applied to most rotating non-smooth systems where the spectral method is inapplicable. In the process, we rigourously define an appropriate coordinate system in which to investigate stability – this happens to be the well-known Tisserand mean axis of the body – as well as systematically distinguish perturbations that introduce angular momentum and/or jumps in the stress state from those that do not. With a view to future application to planetary objects, we specialize the stability test to freely rotating self-gravitating ellipsoids. This is then employed to investigate the stability to homogeneous perturbations of rotating inviscid fluid ellipsoids. We recover results consistent with earlier predictions, and, in the process, also reconcile some contradictory conclusions about the stability of Maclaurin spheroids. Finally, we consider the equilibrium and stability of freely rotating self-gravitating Bingham fluid ellipsoids. We find that the equilibrium shapes of most such ellipsoids are secularly stable to homogeneous perturbations that preserve angular momentum, but not otherwise. We also touch upon the effect of shear thinning on stability.

We investigate and compare the boundary conditions that are to be applied to free-surface problems involving inlet and outlets of Newtonian fluid, typically found in coating processes. The flux of fluid is a priori known at an inlet, but unknown at an outlet, where it is governed by the local behaviour near the film-forming meniscus. In the limit of vanishing capillary number it is well known that the flux scales with , but this classical result is non-uniform as the contact angle approaches . By examining this limit we find a solution that is uniformly valid for all contact angles. Furthermore, by considering the far-field behaviour of the free surface we show that there exists a critical capillary number above which the problem at an inlet becomes over-determined. The implications of this result for the modelling of coating flows are discussed.

We analyse numerically the linear stability of a liquid-metal flow in a rectangular duct with perfectly electrically conducting walls subject to a uniform transverse magnetic field. A non-standard three-dimensional vector stream-function/vorticity formulation is used with a Chebyshev collocation method to solve the eigenvalue problem for small-amplitude perturbations. A relatively weak magnetic field is found to render the flow linearly unstable as two weak jets appear close to the centre of the duct at the Hartmann number In a sufficiently strong magnetic field, the instability following the jets becomes confined in the layers of characteristic thickness located at the walls parallel to the magnetic field. In this case the instability is determined by which results in both the critical Reynolds number and wavenumber scaling as Instability modes can have one of the four different symmetry combinations along and across the magnetic field. The most unstable is a pair of modes with an even distribution of vorticity along the magnetic field. These two modes represent strongly non-uniform vortices aligned with the magnetic field, which rotate either in the same or opposite senses across the magnetic field. The former enhance while the latter weaken one another provided that the magnetic field is not too strong or the walls parallel to the field are not too far apart. In a strong magnetic field, when the vortices at the opposite walls are well separated by the core flow, the critical Reynolds number and wavenumber for both of these instability modes are the same: and The other pair of modes, which differs from the previous one by an odd distribution of vorticity along the magnetic field, is more stable with an approximately four times higher critical Reynolds number.

In this work, experiments using a pendulum apparatus were conducted for two particles engaged in oblique, wetted collisions over a range of impact angles, impact velocities, coating thicknesses, liquid viscosities, particle materials, and particle radii. From previous studies on normal or head-on collisions, the two particles bounce apart if the Stokes number (a ratio of particle inertia to viscous forces) exceeds a critical value, whereas they stick together if the Stokes number is below this critical value. However, for oblique collisions, an additional outcome is observed at moderate Stokes numbers and impact angles, in which the spheres initially stick together, rotate as a doublet, and then separate due to centrifugal forces. We refer to this outcome as ‘stick–rotate–separate’. For subcritical Stokes numbers exhibiting this new outcome, the experimental results for the apparent coefficient of normal restitution and angle of rotation from impact to separation show only weak dependence on the fluid viscosity and thickness and the dry restitution coefficient, whereas they both decrease with increasing particle radius. These results are in contrast with those for supercritical Stokes numbers in which the spheres bounce upon impact. An accompanying theory based on lubrication forces, the glass transition of the liquid layer, and solid deformation and rebound agrees well with experimental results and gives insight into the observed trends.

Streamwise rolls and accompanying streamwise streaks are ubiquitous in wall-bounded shear flows, both in natural settings, such as the atmospheric boundary layer, as well as in controlled settings, such as laboratory experiments and numerical simulations. The streamwise roll and streak structure has been associated with both transition from the laminar to the turbulent state and with maintenance of the turbulent state. This close association of the streamwise roll and streak structure with the transition to and maintenance of turbulence in wall-bounded shear flow has engendered intense theoretical interest in the dynamics of this structure. In this work, stochastic structural stability theory (SSST) is applied to the problem of understanding the dynamics of the streamwise roll and streak structure. The method of analysis used in SSST comprises a stochastic turbulence model (STM) for the dynamics of perturbations from the streamwise-averaged flow coupled to the associated streamwise-averaged flow dynamics. The result is an autonomous, deterministic, nonlinear dynamical system for evolving a second-order statistical mean approximation of the turbulent state. SSST analysis reveals a robust interaction between streamwise roll and streak structures and turbulent perturbations in which the perturbations are systematically organized through their interaction with the streak to produce Reynolds stresses that coherently force the associated streamwise roll structure. If a critical value of perturbation turbulence intensity is exceeded, this feedback results in modal instability of the combined streamwise roll/streak and associated turbulence complex in the SSST system. In this instability, the perturbations producing the destabilizing Reynolds stresses are predicted by the STM to take the form of oblique structures, which is consistent with observations. In the SSST system this instability exists together with the transient growth process. These processes cooperate in determining the structure of growing streamwise roll and streak. For this reason, comparison of SSST predictions with experiments requires accounting for both the amplitude and structure of initial perturbations as well as the influence of the SSST instability. Over a range of supercritical turbulence intensities in Couette flow, this instability equilibrates to form finite amplitude time-independent streamwise roll and streak structures. At sufficiently high levels of forcing of the perturbation field, equilibration of the streamwise roll and streak structure does not occur and the flow transitions to a time-dependent state. This time-dependent state is self-sustaining in the sense that it persists when the forcing is removed. Moreover, this self-sustaining state rapidly evolves toward a minimal representation of wall-bounded shear flow turbulence in which the dynamics is limited to interaction of the streamwise-averaged flow with a perturbation structure at one streamwise wavenumber. In this minimal realization of the self-sustaining process, the time-dependent streamwise roll and streak structure is maintained by perturbation Reynolds stresses, just as is the case of the time-independent streamwise roll and streak equilibria. However, the perturbation field is maintained not by exogenously forced turbulence, but rather by an endogenous and essentially non-modal parametric growth process that is inherent to time-dependent dynamical systems.

Kinetic theory provides a rigorous foundation for calculating the dynamics of gas flow at arbitrary degrees of rarefaction, with solutions of the Boltzmann equation requiring numerical methods in many cases of practical interest. Importantly, the near-continuum regime can be examined analytically using asymptotic techniques. These asymptotic analyses often assume steady flow, for which analytical slip models have been derived. Recently, developments in nanoscale fabrication have stimulated research into the study of oscillatory non-equilibrium flows, drawing into question the applicability of the steady flow assumption. In this article, we present a formal asymptotic analysis of the unsteady linearized Boltzmann–BGK equation, generalizing existing theory to the oscillatory (time-varying) case. We consider the near-continuum limit where the mean free path and oscillation frequency are small. The complete set of hydrodynamic equations and associated boundary conditions are derived for arbitrary Stokes number and to second order in the Knudsen number. The first-order steady boundary conditions for the velocity and temperature are found to be unaffected by oscillatory flow. In contrast, the second-order boundary conditions are modified relative to the steady case, except for the velocity component tangential to the solid wall. Application of this general asymptotic theory is explored for the oscillatory thermal creep problem, for which unsteady effects manifest themselves at leading order.

We present a theoretical and numerical study of the effective mean force exerted on an oceanic mean flow due to the presence of small-amplitude internal waves that are forced by the oscillatory flow of a barotropic tide over undulating topography and are also subject to dissipation. This extends the classic lee-wave drag problem of atmospheric wave–mean interaction theory to a more complicated oceanographic setting, because now the steady lee waves are replaced by oscillatory internal tides and, most importantly, because now the three-dimensional oceanic mean flow is defined by time averaging over the fast tidal cycles rather than by the zonal averaging familiar from atmospheric theory. Although the details of our computation are quite different, we recover the main action-at-a-distance result from the atmospheric setting, namely that the effective mean force that is felt by the mean flow is located in regions of wave dissipation, and not necessarily near the topographic wave source. Specifically, we derive an explicit expression for the effective mean force at leading order using a perturbation series in small wave amplitude within the framework of generalized Lagrangian-mean theory, discuss in detail the range of situations in which a strong, secularly growing mean-flow response can be expected, and then compute the effective mean force numerically in a number of idealized examples with simple topographies.

The incompressible impulsive time scale for inviscid liquid sloshing in open rigid containers suddenly put into motion is defined as the intermediate time scale in between the acoustic time scale and the gravitational time scale. Surge and sway boundary-value problems for incompressible impulsive sloshing in some realistic container shapes are solved analytically to the leading order in a small-time expansion. A solution is provided for two types of horizontal cylinders: a triangular cylindrical wedge and a half-filled circular cylinder. The surface velocity and the hydrodynamic force with its corresponding virtual fluid mass are calculated. The cases of constant impulsive velocity and constant impulsive acceleration are linked by transformation equations. Flows with waterline singularities are discussed, being leading-order outer flows in terms of matched asymptotic expansions.

A numerical model is developed to analyse the interaction of artificial cilia with the surrounding fluid in a three-dimensional setting in the limit of vanishing fluid inertia forces. The cilia are modelled using finite shell elements and the fluid is modelled using a boundary element approach. The coupling between both models is performed by imposing no-slip boundary conditions on the surface of the cilia. The performance of the model is verified using various reference problems available in the literature. The model is used to simulate the fluid flow due to magnetically actuated artificial cilia. The results show that narrow and closely spaced cilia create the largest flow, that metachronal waves along the width of the cilia create a significant flow in the direction of the cilia width and that the recovery stroke in the case of the out-of-plane actuation of the cilia strongly depends on the cilia width.

A linear spatial stability analysis is performed on the velocity profiles measured in the wake of an actively flexible robotic elliptical fin to find the frequency of maximum spatial growth, that is, the hydrodynamic resonant frequency of the time-averaged jet. It is found that: (i) optima in propulsive efficiency occur when the driving frequency of a flapping fin matches the resonant frequency of the jet profile; (ii) there can be multiple wake resonant frequencies and modes corresponding to multiple peaks in efficiency; and (iii) some wake structures transition from one pattern to another when the wake instability mode transitions. A theoretical framework, termed wake resonance theory, is developed and utilized to explain the mechanics and energetics of unsteady self-propulsion. Experimental data are used to validate the theory. The analysis, although one-dimensional, captures the performance exhibited by a three-dimensional propulsor, showing the robustness and broad applicability of the technique.

Force and particle image velocimetry measurements were conducted on a NACA 0012 aerofoil undergoing small-amplitude high-frequency plunging oscillation at low Reynolds numbers and angles of attack in the range 0–. For angles of attack less than or equal to the stall angle, at high Strouhal numbers, significant bifurcations are observed in the time-averaged lift coefficient resulting in two lift-coefficient branches. The upper branch is associated with an upwards deflected jet, and the lower branch is associated with a downwards deflected jet. These branches are stable and highly repeatable, and are achieved by increasing or decreasing the frequency in the experiments. Increasing frequency refers to starting from stationary and increasing the frequency very slowly (while waiting for the flow to reach an asymptotic state after each change in frequency); decreasing frequency refers to impulsively starting at the maximum frequency and decreasing the frequency very slowly. For the latter case, angle of attack, starting position and initial acceleration rate are also parameters in determining which branch is selected. The bifurcation behaviour is closely related to the properties of the trailing-edge vortices. The bifurcation was therefore not observed for very small plunge amplitudes or frequencies due to insufficient trailing-edge vortex strength, nor at larger angles of attack due to greater asymmetry in the strength of the trailing-edge vortices, which creates a preference for a downward deflected jet. Vortex strength and asymmetry parameters are derived from the circulation measurements. It is shown that the most appropriate strength parameter in determining the onset of deflected jets is the circulation normalized by the plunge velocity.

Closed and open hollow wakes are considered as analytic models for the two-dimensional inviscid steady flow past a plate normal to the stream. It is shown that only open configurations which satisfy the Kutta condition exist. The main argument is based on considering a plate located on the edge of a step with varying height. It is shown that solutions for open wakes exist for backward-, null and forward-facing steps, while closed wakes only exist for backward-facing steps. The occurrence of secondary separation has been modelled by adding a hollow region attached to the downstream corner. Peculiar accuracy issues of the problem are pointed out which may explain other contradictory results from the literature. It is shown how the Kirchhoff wake is a limiting solution for certain values of the governing parameters.

A finite-core heton is a baroclinic -plane modon of a special type: it is composed of two patches of uniform quasi-geostrophic potential vorticity (PV) residing in different layers of a two-layer rotating fluid. This paper focuses on numerical construction of steadily translating, doubly symmetric, finite-core hetons and testing their stability. Such a heton, which possesses symmetry about the translation axis and the transverse axis, is a stationary solution to the equations of PV conservation in each of the layers when considered in a comoving frame of reference. When constructing the heton solutions and examining their bifurcations, we identify a heton by a pair of independent non-dimensional parameters, the half-length (in the translation direction) of a PV patch and the distance of the front point of the upper patch from the translation axis. The advantage of this method over other tried approaches is that it allows one to obtain solutions of new, previously unknown types. The results of testing the heton stability are presented on the plane made by a mean radius of a PV patch and the (horizontal) separation between the centroids of the patches. Two kinds of stability are tested separately, the stability to arbitrary perturbations that do not preserve the symmetry of the initial state and the stability to so-called symmetric perturbations that do not violate the initial symmetry. The hetons comparable in size with the Rossby radius, and smaller, are always stable in both senses. However, when some critical size is exceeded, the heton stability becomes dependent on the separation, and the larger the heton, the higher the separation required for stability. The separation guaranteeing the stability to symmetric perturbations is smaller than that required for the stability to arbitrary perturbations. Interrelations between instabilities and bifurcations are briefly discussed.

A two-phase Newtonian surface fluid is modelled as a surface Cahn–Hilliard–Navier–Stokes equation using a stream function formulation. This allows one to circumvent the subtleties in describing vectorial second-order partial differential equations on curved surfaces and allows for an efficient numerical treatment using parametric finite elements. The approach is validated for various test cases, including a vortex-trapping surface demonstrating the strong interplay of the surface morphology and the flow. Finally the approach is applied to a Rayleigh–Taylor instability and coarsening scenarios on various surfaces.

This article deals with open-loop control of open-cavity flows with harmonic forcings. Two-dimensional laminar open-cavity flows usually undergo a supercritical Hopf bifurcation at some critical Reynolds number: a global mode becomes unstable and its amplitude converges towards a limit cycle. Such behaviour may be accurately captured by a Stuart–Landau equation, which governs the amplitude of the global mode. In the present article, we study the effect on such a flow of a forcing characterized by its frequency , its amplitude and its spatial structure . The system reacts like a forced Van der Pol oscillator. In the general case, such a forcing modifies the linear dynamics of the global mode. It is then possible to predict preferred forcing frequencies , at which the global mode may be stabilized with the smallest possible forcing amplitude . In the case of a forcing frequency close to the frequency of the global mode, a locking phenomenon may be observed if the forcing amplitude is sufficiently high: the frequency of the flow on the limit cycle may be modified with a very small forcing amplitude . In each case, we compute all harmonics of the flow field and all coefficients that enter the amplitude equations. In particular, it is possible to find preferred forcing structures that achieve strongest impact on the flow field. In the general case, these are the optimal forcings, which are defined as the forcings that trigger the strongest energy response. In the case of a forcing frequency close to the frequency of the global mode, a forcing structure equal to the adjoint global mode ensures the lowest forcing amplitude . All predictions given by the amplitude equations are checked against direct numerical simulations conducted at a supercritical Reynolds number. We show that a global mode may effectively be stabilized by a high-frequency harmonic forcing, which achieves suppression of the perturbation frequencies that are lower than the forcing frequency, and that a near-resonant forcing achieves locking of the flow onto the forcing frequency, as predicted by the amplitude equations.

When a drop impacts at low velocity onto a pool surface, a hemispheric air layer cushions and can delay direct contact. Herein we use ultra-high-speed video to study the rupture of this layer, to explain the resulting variety of observed distribution of bubbles. The size and distribution of micro-bubbles is determined by the number and location of the primary punctures. Isolated holes lead to the formation of bubble necklaces when the edges of two growing holes meet, whereas bubble nets are produced by regular shedding of micro-bubbles from a sawtooth edge instability. For the most viscous liquids the air film contracts more rapidly than the capillary–viscous velocity through repeated spontaneous ruptures of the edge. From the speed of hole opening and the total volume of micro-bubbles we conclude that the air sheet ruptures when its thickness approaches .

The dynamics of solitary gravity–capillary water waves propagating on the surface of a three-dimensional fluid domain is studied numerically. In order to accurately compute complex time-dependent solutions, we simplify the full potential flow problem by using surface variables and taking a particular cubic truncation possessing a Hamiltonian with desirable properties. This approximation agrees remarkably well with the full equations for the bifurcation curves, wave profiles and the dynamics of solitary waves for a two-dimensional fluid domain, and with higher-order truncations in three dimensions. Fully localized solitary waves are then computed in the three-dimensional problem and the stability and interaction of both line and localized solitary waves are investigated via numerical time integration of the equations. There are many solitary wave branches, indexed by their finite energy as their amplitude tends to zero. The dynamics of the solitary waves is complex, involving nonlinear focusing of wavepackets, quasi-elastic collisions, and the generation of propagating, spatially localized, time-periodic structures akin to breathers.