In this work, an attempt is made to gain a better understanding of the breakage of low-viscosity drops in turbulent flows by determining the dynamics of deformation of an inviscid drop in response to a pressure variation acting on the drop surface. Known scaling relationships between wavenumbers and frequencies, and between pressure fluctuations and velocity fluctuations in the inertial subrange are used in characterizing the pressure fluctuation. The existence of a maximum stable drop diameter dmax follows once scaling laws of turbulent flow are used to correlate the magnitude of the disruptive forces with the duration for which they act.

Two undetermined dimensionless quantities, both of order unity, appear in the equations of continuity, motion, and the boundary conditions in terms of pressure fluctuations applied on the surface. One is a constant of proportionality relating root-mean-square values of pressure and velocity differences between two points separated by a distance l. The other is a Weber number based on turbulent stresses acting on the drop and the resisting stresses in the drop due to interfacial tension. The former is set equal to 1, and the latter is determined by studying the interaction of a drop of diameter equal to dmax with a pressure fluctuation of length scale equal to the drop diameter. The model is then used to study the breakage of drops of diameter greater than dmax and those with densities different from that of the suspending fluid.

It is found that, at least during breakage of a drop of diameter greater than dmax by interaction with a fluctuation of equal length scale, a satellite drop is always formed between two larger drops. When very large drops are broken by smaller-length-scale fluctuations, highly deformed shapes are produced suggesting the possibility of further fragmentation due to instabilities. The model predicts that as the dispersed-phase density increases, dmax decreases.

We report experimental and numerical observations on the way initially symmetric and time-periodic fluid oscillations in baffled channels develop in complexity. Experiments are carried out in a spatially periodic baffled channel with a sinusoidal oscillatory flow. At modest Reynolds number the observed vortex structure is symmetric and time periodic. At higher values the flow progressively becomes three-dimensional, asymmetric and aperiodic. A two-dimensional simulation of incompressible Newtonian flow is able to follow the flow pattern at modest oscillatory Reynolds number. At higher values we report the development of both asymmetry and a period-doubling cascade leading to a chaotic flow regime. A bifurcation diagram is constructed that can describe the progressive increase in complexity of the flow.

The multi-layer exchange equations for gravitationally driven flows between two basins with stable Boussinesq type of stratification in discrete layers, specified far upstream on either side of a connecting strait, result in a hydraulic control condition that must be satisfied at the narrowest part of the contraction, the control point. If one stagnant layer is present at the control point, the control condition that applies to all layers collectively may be separated into two such conditions that apply independently to two groups of layers going in opposite directions separated by the stagnant layer. Such bidirectional flow regimes exist if the structure of the prespecified density profiles permits each of the opposing groups to vertically reduce their thickness by the ratio 2/3 relative to their upstream thicknesses, leaving space for the stagnant layer to protrude through the contraction. Under these restrictions, the bidirectional flow is controlled by the fastest propagating wave mode and the stationary solution then relies on the superposition of two previously known unidirectional self-similar flow regimes that are completely decoupled. Techniques for their numerical computation are presented. The transition into loosely coupled and fully coupled flow is discussed. The decoupling principle also applies when several non-adjacent stagnant layers are simultaneously present at control in which case multiple groups of decoupled layers flow in alternating directions.

This paper examines whether G. K. Batchelor's (1988) theory of the propagation of planar concentration disturbances and the occurrence of instabilities in uniform fluidized beds can be applied to bubbly flows. According to this theory the propagation of long weakly nonlinear gas volume concentration waves is governed by the Burgers equation. Experiments on the propagation of weak concentration shock waves and small, but finite, amplitude periodic waves are presented; good agreement is found with classic solutions of Burgers’ equation. For example, the phenomenon of amplitude saturation, familiar from nonlinear acoustics, is established here for concentration waves. Batchelor's instability conditions are given for bubbly flows, and his model for the bulk modulus of elasticity of the dispersed phase is used to obtain estimates of the critical volume concentration at which a uniform bubbly flow becomes unstable to planar disturbances. Observations of the onset of instabilities of bubbly flow in a pipe are described, and compared with what would be expected from Batchelor's theory.

A weakly nonlinear analysis is presented of the small oscillations of nearly inviscid liquid bridges subjected to almost resonant axial vibrations of the disks. An amplitude equation is derived for the evolution of the complex amplitude of the oscillations that exhibits hysteresis and period doublings. We also analyse the steady streaming in the bulk forced by the oscillatory boundary layers near the disks; the boundary layer near the free surface produces no forcing terms. In particular some experimentally observed patterns are explained, and some new, non-observed ones are predicted. We conclude that the structure of this steady flow is not the appropriate one to counterbalance steady thermocapillary convection, but our results indicate how to get the desired counterbalancing effect.

We examine the linear stability of elliptical columns of uniform potential vorticity subject to two-dimensional (horizontal) straining within a rapidly rotating, stratified (quasi-geostrophic) fluid. We find that horizontal straining can promote the exponential growth of three-dimensional disturbances when the vortex height-to-width aspect ratio exceeds, qualitatively, three times the ratio of the Coriolis parameter to the buoyancy frequency. This instability is not related to the usual baroclinic instability which operates on shallow vortex columns whose potential vorticity changes sign with height. The nonlinear development of these instabilities is investigated numerically using a high-resolution contour surgery algorithm. Simulations are conducted for both a Boussinesq (ocean-like) fluid and a compressible (atmospheric-like) fluid having exponentially decreasing density with height. The simulations reveal a generic nonlinear development that results in a semi-ellipsoidal baroclinic vortex dome at the lower surface and, in the case of a Boussinesq fluid, another such dome at the upper surface.

The related problem of two interacting vortex columns is also examined. A generic three-dimensional instability and nonlinear development occurs no matter how great the distance between the vortex columns, provided that they are sufficiently tall.

Our results may bear upon the observed structure of many atmospheric and oceanic vortices, whose height-to-width aspect ratios are consistent with our findings. Remarkably, even strongly ageostrophic vortices, such as tropical cyclones, fit the pattern. Our results furthermore re-open questions about the long-time nature of freely decaying quasi-geostrophic turbulence, for which recent simulations indicate a progressive two-dimensionalization by vortex alignment, while earlier simulations have indicated long-lived baroclinic vortices, not unlike what we find here.

The upper bound on momentum transport in the turbulent regime of plane Couette flow is considered. Busse (1970) obtained a bound from a variational formulation based on total energy conservation and the mean momentum equation. Two-dimensional asymptotic solutions of the resulting Euler-Lagrange equations for the system were obtained in the large-Reynolds-number limit. Here we make a toroidal poloidal decomposition of the flow and impose an additional power integral constraint, which cannot be satisfied by two-dimensional flows. Nevertheless, we show that the additional constraint can be met by only small modifications to Busse's solution, which leaves his momentum transport bound unaltered at lowest order. On the one hand, the result suggests that the addition of further integral constraints will not significantly improve bound estimates. On the other, our optimal solution, which possesses a weak spanwise roll in the outermost of Busse's nested boundary layers, appears to explain the three-dimensional structures observed in experiments. Only in the outermost boundary layer and in the main stream is the solution three-dimensional. Motion in the thinner layers remains two-dimensional characterized by streamwise rolls.

The development of a two-dimensional viscous incompressible flow generated from an infinitesimally thin flat plate, impulsively started or uniformly accelerated normal to the free stream is studied computationally. An adaptive numerical scheme, based on vortex methods, is used to integrate the vorticity–velocity formulation of the Navier–Stokes equations. The results of the computations complement relevant experimental works while providing us with quantities such as the vorticity field and the unsteady forces experienced by the body. For the uniformly accelerated plate the present simulations capture the development of a number of centers of vorticity along the primary separating shear layer. This phenomenon has been observed in experimental works but has not been predicted by inviscid models. The present simulations suggest that this Kelvin–Helmholtz-type instability is driven by the interaction of primary and secondary vorticity near the tips of the plate and depends on the acceleration of the plate.

The phenomenon of pattern formation of free-surface waves of a fluid confined in an annulus the inner wall of which vibrates radially, is investigated both theoretically and experimentally. Although the waves are excited by harmonic axisymmetric deformations of the inner shell, depending on the vibration frequency both axisymmetric and non-symmetric wave patterns may arise.

Experimental observations have revealed that waves are excited in two different resonance regimes. The first type corresponds to forced resonance, in which axisymmetric patterns are realized with eigenfrequencies equal to the frequency of excitation. The second kind is parametric resonance, in which case the waves are ‘transverse’, with their crests and troughs aligned perpendicular to the vibrating wall. These so-called cross-waves have frequencies equal to half of that of the wavemaker.

Both kinds of resonance were investigated theoretically using Lamé's method of superposition. It was shown experimentally that the pure forced resonant standing waves are not realized when the amplitude of excitation is beyond the threshold of parametric resonance for non-symmetric waves. The experimental observations agree very well with the theoretical results.

We examine results of direct numerical simulations (DNS) of homogeneous turbulence in the presence of stable stratification. We focus on the effects of stratification on eddy diffusion, and the distribution of pairs of particles released in the flow. DNS results are presented over a range of stratification, and at Reynolds numbers compatible with aliased free spectral results for a resolution of 128 mesh points. We compare results for particle dispersion to simple analytic theories such as that proposed by Csanady (1964) and Pearson et al. (1983) by adapting the basic Langevin model to decaying turbulence at low Reynolds numbers. Stable stratification is found to arrest both single particle displacements and pair separation in the direction of stratification, but it leaves these quantities nearly unaltered in the transverse direction. With respect to the dynamics of stratified flows, we find that regions of strong viscous dissipation are intermittently spaced, and are associated with large horizontal vorticity, consistent with recent experimental results by Fincham et al. (1994).

The axisymmetric flows arising in a rotating annulus with a superimposed forced flow are investigated with a pseudo-spectral numerical method. The flow enters the annulus at the inner radius with a radial velocity, then develops into a geostrophic flow azimuthally directed and flanked by two Ekman (nonlinear) boundary layers, and finally exits the outer radius, with a radially directed velocity. In this study the rotation rate of the cavity is fixed and very high. When the forced flow is weak, the flow is steady. On increasing the mass flow rate, the flow evolves to a chaotic temporal behaviour through several bifurcations, which perturbs the basic spatial configuration of the flow. The first bifurcation drives the steady state into an oscillatory regime, associated with a break of symmetry with respect to the midheight of the annulus. The entry flow travels radially through the cavity as in the steady flow, but it wavers and then is alternately sucked towards each Ekman layer. The frequency of this oscillation is close to the rotation rate frequency of the cavity, which is characteristic of inertial waves in rotating flows. A second transition to a quasi-periodic regime is characterized by the appearance of a second frequency. Further increases in the flow rate lead to a period-five state, via a locking of both frequencies, and then to a chaotic motion. This second frequency is of the order of the inverse of the Ekman spin-up characteristic time, suggesting that this instability is originated by the relaxation of the perturbations in the flow field. These perturbations of the unsteady flow field are corotating vortices along the rigid boundary walls. They are excited by the entry flow and their strength diminishes with increasing radius due to the low value of the Reynolds number. The parameters characterizing the unstable flows are also consistent with this explanation. The conclusion is that in this configuration, the origin of the described dynamical behaviour is not the instability of the Ekman boundary layers, as could be expected, but the instability of the entry flow. The reason is the importance of the nonlinear inertial terms in cavities with small radius of curvature.

We present theoretical results on spontaneous stratification of sedimenting suspensions and colloids caused by a lateral temperature gradient. Fluid motion is treated in the Stokes approximation, and motion of suspended particles is described by Burgers equation with convection. The internal structure and interaction of shocks at convection roll boundaries is studied numerically using a reduced one-dimensional model based on a Galerkin approach. Qualitative comparison is made to experimental data.

New experiments have been carried out in a large laboratory channel to explore the structure of turbulent motion in the water layer beneath surface gravity waves. These experiments involve pure wind waves as well as wind-ruffled mechanically generated waves. A submersible two-component LDV system has been used to obtain the three components of the instantaneous velocity field along the vertical direction at a single fetch of 26 m. The displacement of the free surface has been determined simultaneously at the same downstream location by means of wave gauges. For both types of waves, suitable separation techniques have been used to split the total fluctuating motion into an orbital contribution (i.e. a motion induced by the displacement of the surface) and a turbulent contribution. Based on these experimental results, the present paper focuses on the structure of the water turbulence. The most prominent feature revealed by the two sets of experiments is the enhancement of both the turbulent kinetic energy and its dissipation rate with respect to values found near solid walls. Spectral analysis provides clear indications that wave–turbulence interactions greatly affect energy transfers over a significant frequency range by imposing a constant timescale related to the wave-induced strain. For mechanical waves we discuss several turbulent statistics and their modulation with respect to the wave phase, showing that the turbulence we observed was deeply affected at both large and small scales by the wave motion. An analysis of the phase variability of the bursting suggests that there is a direct interaction between the waves and the underlying turbulence, mainly at the wave crests. Turbulence budgets show that production essentially takes place in the wavy region of the flow, i.e. above the wave troughs. These results are finally used to address the nature of the basic mechanisms governing wave–turbulence interactions.

It is now well-known that the wake transition regime for a circular cylinder involves two modes of small-scale three-dimensional instability (modes A and B), depending on the regime of Reynolds number (Re), although almost no understanding of the physical origins of these instabilities, or indeed their effects on near-wake formation, have hitherto been made clear. We address these questions in this paper. In particular, it is found that the two different modes A and B scale on different physical features of the flow. Mode A has a larger spanwise wavelength of around 3–4 diameters, and scales on the larger physical structure in the flow, namely the primary vortex core. The wavelength for mode A is shown to be the result of an ‘elliptic instability’ in the nearwake vortex cores. The subsequent nonlinear growth of vortex loops is due to a feedback from one vortex to the next, involving spanwise-periodic deformation of core vorticity, which is then subject to streamwise stretching in the braid regios. This mode gives an out-of-phase streamwise vortex pattern.

In contrast, mode-B instability has a distinctly smaller wavelength (1 diameter) which scales on the smaller physical structure in the flow, the braid shear layer. It is a manifestation of an instability in a region of hyperbolic flow. It is quite distinct from other shear flows, in that it depends on the reverse flow of the bluff-body wake; the presence of a fully formed streamwise vortex system, brought upstream from a previous half-cycle, in proximity to the newly evolving braid shear layer, leads to an in-phase stream-wise vortex array, in strong analogy with the ‘Mode 1’ of Meiburg & Lasheras (1988) for a forced unseparated wake. In mode B, we also observe amalgamation of streamwise vortices from a previous braid with like-sign vortices in the subsequent braid.

It is deduced that the large scatter in previous measurements concerning mode A is due to the presence of vortex dislocations. Dislocations are triggered at the sites of some vortex loops of mode A, and represent a natural breakdown of the periodicity of mode A instability. By minimizing or avoiding the dislocations which occur from end contamination or which occur during wake transition, we find an excellent agreement of both critical Re and spanwise wavelength of mode A with the recent secondary stability analysis of Barkley & Henderson (1996).

Wake transition is further characterized by velocity and pressure measurements. It is consistent that, when mode-A instability and large-scale dislocations appear, one finds a reduction of base suction, a reduction of (two-dimensional) Reynolds stress level, a growth in size of the formation region, and a corresponding drop in Strouhal frequency. Finally, the present work leads us to a new clarification of the possible flow states through transition. Right through this regime of Re, there exist two distinct and continuous Strouhal frequency curves: the upper one corresponds with purley small- scale instabilities (e.g. denoted as mode A), while the lower curve corresponds with a combination of small-scale plus dislocation structures (e.g. mode A*). However, some of the flow states are transient or ‘unstable’, and the natural transitioning wake appears to follow the scenario: (2D→A*→B).