An experimental investigation of the flow generated by a homogeneous population of bubbles rising in water is reported for three different bubble diameters (d = 1.6, 2.1 and 2.5 mm) and moderate gas volume fractions (0.005 ≤ α ≤ 0.1). The Reynolds numbers, Re = V0d/ν, based on the rise velocity V0 of a single bubble range between 500 and 800. Velocity statistics of both the bubbles and the liquid phase are determined within the homogeneous bubble swarm by means of optical probes and laser Doppler anemometry. Also, the decaying agitation that takes place in the liquid just after the passage of the bubble swarm is investigated from high-speed particle image velocimetry measurements. Concerning the bubbles, the average velocity is found to evolve as V0α−0.1 whereas the velocity fluctuations are observed to be almost independent of α. Concerning the liquid fluctuations, the probability density functions adopt a self-similar behaviour when the gas volume fraction is varied, the characteristic velocity scaling as V0α0.4. The spectra of horizontal and vertical liquid velocity fluctuations are obtained with a resolution of 0.6 mm. The integral length scale Λ is found to be proportional to V02/g or equivalently to d/Cd0, where g is the gravity acceleration and Cd0 the drag coefficient of a single rising bubble. Normalized by using Λ, the spectra are independent on both the bubble diameter and the volume fraction. At large scales, the spectral energy density evolves as the power −3 of the wavenumber. This range starts approximately from Λ and is followed for scales smaller than Λ/4 by a classic −5/3 power law. Although the Kolmogorov microscale is smaller than the measurement resolution, the dissipation rate is however obtained from the decay of the kinetic energy after the passage of the bubbles. It is found to scale as α0.9V03/Λ. The major characteristics of the agitation are thus expressed as functions of the characteristics of a single rising bubble. Altogether, these results provide a rather complete description of the bubble-induced turbulence.

In this paper, the results of experiments on unsteady disturbances in the boundary-layer flow over a disk rotating in otherwise still air are presented. The flow was perturbed impulsively at a point corresponding to a Reynolds number R below the value at which transition from laminar to turbulent flow is observed. Among the frequencies excited are convectively unstable modes, which form a three-dimensional wave packet that initially convects away from the source. The wave packet consists of two families of travelling convectively unstable waves that propagate together as one packet. These two families are predicted by linear-stability theory: branch-2 modes dominate close to the source but, as the packet moves outwards into regions with higher Reynolds numbers, branch-1 modes grow preferentially and this behaviour was found in the experiment. However, the radial propagation of the trailing edge of the wave packet was observed to tend towards zero as it approaches the critical Reynolds number (about 510) for the onset of radial absolute instability. The wave packet remains convectively unstable in the circumferential direction up to this critical Reynolds number, but it is suggested that the accumulation of energy at a well-defined radius, due to the flow becoming radially absolutely unstable, causes the onset of laminar–turbulent transition. The onset of transition has been consistently observed by previous authors at an average value of 513, with only a small scatter around this value. Here, transition is also observed at about this average value, with and without artificial excitation of the boundary layer. This lack of sensitivity to the exact form of the disturbance environment is characteristic of an absolutely unstable flow, because absolute growth of disturbances can start from either noise or artificial sources to reach the same final state, which is determined by nonlinear effects.

Material pitting from cavitation bubble collapse is investigated numerically including two-way fluid–structure interaction (FSI). A hybrid numerical approach which links an incompressible boundary element method (BEM) solver and a compressible finite difference flow solver is applied to capture non-spherical bubble dynamics efficiently and accurately. The flow codes solve the fluid dynamics while intimately coupling the solution with a finite element structure code to enable simulation of the full FSI. During bubble collapse high impulsive pressures result from the impact of the bubble re-entrant jet on the material surface and from the collapse of the remaining bubble ring. A pit forms on the material surface when the impulsive pressure is large enough to result in high equivalent stresses exceeding the material yield stress. The results depend on bubble dynamics parameters such as the size of the bubble at its maximum volume, the bubble standoff distance from the material wall, and the pressure driving the bubble collapse. The effects of these parameters on the re-entrant jet, the following bubble ring collapse pressure, and the generated material pit characteristics are investigated.

A novel data-driven modal decomposition of fluid flow is proposed, comprising key features of proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD). The first mode is the normalized real or imaginary part of the DMD mode that minimizes the time-averaged residual. The $N$th mode is defined recursively in an analogous manner based on the residual of an expansion using the first $N-1$ modes. The resulting recursive DMD (RDMD) modes are orthogonal by construction, retain pure frequency content and aim at low residual. Recursive DMD is applied to transient cylinder wake data and is benchmarked against POD and optimized DMD (Chen et al., J. Nonlinear Sci., vol. 22, 2012, pp. 887–915) for the same snapshot sequence. Unlike POD modes, RDMD structures are shown to have purer frequency content while retaining a residual of comparable order to POD. In contrast to DMD, with exponentially growing or decaying oscillatory amplitudes, RDMD clearly identifies initial, maximum and final fluctuation levels. Intriguingly, RDMD outperforms both POD and DMD in the limit-cycle resolution from the same snapshots. Robustness of these observations is demonstrated for other parameters of the cylinder wake and for a more complex wake behind three rotating cylinders. Recursive DMD is proposed as an attractive alternative to POD and DMD for empirical Galerkin models, in particular for nonlinear transient dynamics.

Wind, blowing over a water surface, induces a thin layer of high vorticity in which the wind stress is supported by molecular viscosity; the magnitude of the surface drift, the velocity difference across the layer, being of the order of 3% of the wind speed. When long waves move across the surface, there is a nonlinear augmentation of the surface drift near the long-wave crests, so that short waves, superimposed on the longer ones, experience an augmented drift in these regions. This is shown to reduce the maximum amplitude that the short waves can attain when they are at the point of incipient breaking.

Theoretical estimates of the reduction are compared with measurements in wind-wave tanks by the authors and by Mitsuyasu (1966) in which long mechanically generated waves are superimposed on short wind-generated waves. The reductions measured in the energy density of the short waves by increasing the slope of the longer ones at constant wind speed are generally consistent with the predictions of the theory in a variety of cases.

The failure of local isotropy to describe the experimentally obtained derivative moments in turbulent shear flows has previously been well-documented, but is briefly reviewed. The same data are then used to evaluate the hypothesis that the turbulence is locally axisymmetric. Locally axisymmetric turbulence is defined herein as turbulence which is locally invariant to rotations about a preferred axis.

The derivative moment relations are derived from the general form of the two-point velocity correlation tensor near the origin for axisymmetric turbulence. These are used to derive relations for the rate of dissipation of kinetic energy, the mean-square vorticity, and the components of each. Almost all of the experimental derivative moment data are shown to be consistent with these equations, and thus with local axisymmetry.

The present paper contains an analysis of the model of a porous material proposed in part 1, and carries out calculations which allow comparison between theory and the experiments described therein. The relevant boundary conditions to be applied at an interface between a fluid and such a material are considered.

Direct numerical simulation is used to examine the interaction of turbulence with a wall in the absence of mean shear. The influence of a solid wall on turbulence is analysed by first considering two ‘simpler’ types of boundaries. The first boundary is an idealized permeable wall. This boundary isolates and elucidates the viscous effects created by the wall. The second boundary is an idealized free surface. This boundary complements the first by allowing one to isolate and investigate the kinematic effects that occur near boundaries. The knowledge gained from these two simpler flows is then used to understand how turbulence is influenced by solid walls where both viscous and kinematic effects occur in combination.

Examination of the instantaneous flow fields confirms the presence of previously hypothesized structures (splats), and reveals an additional class of structures (antisplats). Statistical analysis of the Reynolds stresses and Reynolds stress transport equations indicates the relative importance of dissipation, intercomponent energy transfer, and energy transport. It is found that it is not the structures themselves, but the imbalance between structures which leads to intercomponent energy transfer. Remarkably, this imbalance (and hence near-wall intercomponent energy transfer) is controlled by viscous processes such as dissipation and diffusion. The analysis presented herein is a departure from past notions of how boundaries influence turbulence. The efficacy of these qualitative physical concepts is demonstrated in Part 2 where improved near-wall turbulence models are derived based on these ideas.

In this paper we describe the design and operating characteristics of a computer-controlled four-roll mill for investigations of particle and drop dynamics in two-dimensional linear flows. The control system is based upon the use of: a video camera to visualize the instantaneous position of the drop or particle; a PDP 11/23 computer, with a pipeline processor acting as an interface between the camera and computer, to calculate the position and implement a control strategy, and d.c. stepping motors to convert an electronic signal to angular velocities of the four rollers. The control objective is to keep the centre of mass of the drop/particle at the centre of the region between the rollers where there is a stagnation point in the undisturbed flow, while maintaining the shear-rate and the ratio of vorticity to strain rate in the flow at fixed values. The resulting system is suitable for studies of: the rotational motions of single solid particles; the deformation and burst of single droplets; or the hydrodynamic interactions of two particles or drops, one of which is held with its centre-of-mass fixed at the stagnation point of the undisturbed flow. In all cases, the flow can be varied from pure rotation to pure strain, and the shear rate can be either steady or changing as a prescribed function of time.

We present closures for the drag and virtual mass force terms appearing in a two-fluid model for flow of a mixture consisting of uniformly sized gas bubbles dispersed in a liquid. These closures were deduced through computational experiments performed using an implicit formulation of the lattice Boltzmann method with a BGK collision model. Unlike the explicit schemes described in the literature, this implicit implementation requires iterative calculations, which, however, are local in nature. While the computational cost per time step is modestly increased, the implicit scheme dramatically expands the parameter space in multiphase flow calculations which can be simulated economically. The closure relations obtained in our study are limited to a regular array of uniformly sized bubbles and were obtained by simulating the rise behaviour of a single bubble in a periodic box. The effect of volume fraction on the rise characteristics was probed by changing the size of the box relative to that of the bubble. While spherical bubbles exhibited the expected hindered rise behaviour, highly distorted bubbles tended to rise cooperatively. The closure for the drag force, obtained in our study through computational experiments, captured both hindered and cooperative rise. A simple model for the virtual mass coefficient, applicable to both spherical and distorted bubbles, was also obtained by fitting simulation results. The virtual mass coefficient for isolated bubbles could be correlated with the aspect ratio of the bubbles.

In Part 1 an extension of the attached eddy hypothesis was developed and applied to equilibrium pressure gradient turbulent boundary layers. In this paper the formulation is applied to data measured by the authors from non-equilibrium layers and agreement with the extended theory is encouraging. Also power spectra of the Reynolds stresses as developed from the extended theory compare favourably with experiment. The experimental data include a check of cone-angle effects by using a flying hot wire.

Investigations of the Richtmyer–Meshkov instability carried out in shock tubes have traditionally used membranes to separate the two gases. The use of membranes, in addition to introducing other experimental difficulties, impedes the use of advanced visualization techniques such as planar laser-induced fluorescence (PLIF). Jones & Jacobs (1997) recently developed a new technique by which a perturbed, membrane-free gas–gas interface can be created in a shock tube. The gases enter the shock tube from opposite ends and exit through two small slots on opposite sides of the test section, forming a stagnation point flow at the interface location. A gentle rocking motion of the shock tube then provides the initial perturbation in the form of a standing wave. The original investigation using this technique utilized dense fog seeding for visualization, which allowed large-scale effects to be observed, but was incapable of resolving smaller-scale features. PLIF visualization is used in the present study to investigate the instability generated by two incident shock strengths (Ms = 1.11 and 1.21), yielding very clear digital images of the flow. Early-time growth rate measurements obtained from these experiments are found to be in excellent agreement with incompressible linear stability theory (appropriately adjusted for a diffuse interface). Very good agreement is also found between the late-time amplitude measurements and the nonlinear models of Zhang & Sohn (1997) and Sadot et al. (1998). Comparison of images from the Ms = 1.11 and 1.21 sequences reveals a significant increase in the amount of turbulent mixing in the higher-Mach-number experiments, suggesting that a mixing transition has occurred.

The derivation of a new vectorial bedload formulation for the transport of coarse sediment by fluid flow is presented in the first part of the paper. This relation has been developed for slopes up to the angle of repose both in the streamwise and transverse directions. The pressure distribution is assumed to be hydrostatic. The bed shear stress for the onset of particle motion and mean particle velocity are obtained from the mean force balance on a particle. A new generalized Bagnold hypothesis is introduced to calculate the sediment content of the bedload layer. The new formulation possesses two innovative features. It is fully nonlinear and vectorial in nature, in addition, it behaves smoothly up to the angle of repose.

A mathematical model of the time evolution of straight river channels is presented in the second half of the paper. This study focuses on the evolution process due to bank erosion in the presence of bedload only. The bed and bank material is taken to be coarse, non-cohesive and uniform in size. The sediment continuity and the fluid momentum conservation equations describe the time evolution of the bed topography and flow field. These equations are coupled through the fluid shear stress acting on the bed. This bed shear stress distribution is predicted with the aid of a simple algebraic turbulent closure model. As regards the computation of the sediment flux, the new fully nonlinear vectorial formulation is found to perform well and renders the evolution model fully mechanistic.

The formation of an erosional front in the time development of straight river channels has been so far obscured in physical experiments. Herein, with the help of the new bedload formulation, the existence and migration speed of the front of erosion are inferred from the analysis of the sediment continuity equation.

The model successfully describes the time relaxation of an initially trapezoidal channel toward an equilibrium cross-sectional shape, as evidenced by comparison with experimental data. This equilibrium is characterized by a constant width, vanishing sediment transport in the transverse direction, and a small but non-vanishing streamwise transport rate of bed sediment.

The two-dimensional motion of a fluid confined between two long horizontal planes, heated and salted from below, is examined. By a combination of perturbation analysis and direct numerical solution of the governing equations, the possible forms of large-amplitude motion are traced out as a function of the four non-dimensional parameters which specify the problem: the thermal Rayleigh number RT, the saline Rayleigh number ES, the Prandtl number σ and the ratio of the diffusivities τ. A branch of time-dependent asymptotic solutions is found which bifurcates from the linear oscillatory instability point. In general, for fixed σ, τ and RS, as RT increases three further abrupt transitions in the form of motion are found to take place independent of the initial conditions. At the first transition, a rather simple oscillatory motion changes into a more complicated one with different structure, at the second, the motion becomes aperiodic and, at the third, the only asymptotic solutions are time independent. Disordered motions are thus suppressed by increasing RT. The time-independent solutions exist on a branch which, it is conjectured, bifurcates from the time-independent linear instability point. They can occur for values of RT less than that at which the third transition point occurs. Hence for some parameter ranges two different solutions exist and a hysteresis effect occurs if solutions obtained by increasing RT and then decreasing RT are followed. The minimum value of RT for which time-independent motion can occur is calculated for fourteen different values of σ, τ and RS. This minimum value is generally much less than the critical value of time-independent linear theory and for the larger values of σ and RS and the smaller values of τ, is less than the critical value of time-dependent linear theory.

An explanation is proposed for the dependence on Reynolds number and other parameters of the number of waves which appear on vortex rings formed by pushing fluid out of a tube. It is shown that the number of waves can be sensitive to the vorticity distribution in the core of the ring. The process of ring formation is discussed and it is concluded that peaked vorticity distributions, limited by viscosity, will occur. Quantitative estimates of the number of waves are made. Agreement with observation is satisfactory.

The evolution of long, finite amplitude Rossby waves in a horizontally sheared zonal current is studied. The wave evolution is described by the Korteweg–de Vries equation or the modified Korteweg-de Vries equation depending on the atmospheric stratification. In either case, the cross-stream modal structure of these waves is given by the long-wave limit of the neutral eigensolutions of the barotropic stability equation. Both non-singular and singular eigensolutions are considered and the appropriate analysis is developed to yield a uniformly valid description of the motion in the critical-layer region where the wave speed matches the flow velocity. The analysis demonstrates that coherent, propagating, eddy structures can exist in stable shear flows and that these eddies have peculiar interaction properties quite distinct from the traditional views of turbulent motion.

The global linear stability of incompressible, two-dimensional shear flows is investigated under the assumptions that far-field pressure feedback between distant points in the flow field is negligible and that the basic flow is only weakly non-parallel, i.e. that its streamwise development is slow on the scale of a typical instability wavelength. This implies the general study of the temporal evolution of global modes, which are time-harmonic solutions of the linear disturbance equations, subject to homogeneous boundary conditions in all space directions. Flow domains of both doubly infinite and semi-infinite streamwise extent are considered and complete solutions are obtained within the framework of asymptotically matched WKBJ approximations. In both cases the global eigenfrequency is given, to leading order in the WKBJ parameter, by the absolute frequency ω0(Xt) at the dominant turning point Xt of the WKBJ approximation, while its quantization is provided by the connection of solutions across Xt. Within the context of the present analysis, global modes can therefore only become time-amplified or self-excited if the basic flow contains a region of absolute instability.

Improved models for the movement of fluid by cilia are presented. A theory which models the cilia of an organism by an array of flexible long slender bodies distributed over and attached at one end to a plane surface is developed. The slender bodies are constrained to move in similar patterns to the cilia of the microorganisms Opalina, Paramecium and Pleurobrachia.

The velocity field is represented by a distribution of force singularities (Stokes flow) along the centre-line of each slender body. Contributions to the velocity field from all the cilia distributed over the plane are summed, to give a streaming effect which in turn implies propulsion of the organism. From this we have been able to model the mean velocity field through the cilia sublayer for the three organisms. We find that, in a frame of reference situated in the organism, the velocity near the surface of the organism is very small – up to one half the length of the cilium – but it increases rapidly to near the velocity of propulsion from then on. This is because of the beating pattern of the cilia; they beat in a near rigid-body rotation during the effective (‘power’) stroke, but during the recovery stroke move close to the wall. Backflow (‘reflux’) is found to occur in the organisms exhibiting antiplectic metachronism (i.e. Paramecium and Pleurobrachia). The occurrence of gradient reversal, but not backflow, has recently been confirmed experimentally (Sleigh & Aiello 1971).

Other important physical values that are obtained from this analysis are the force, bending moment about the base of a cilium and the rate of working. It is found, for antiplectic metachronism, that the force exerted by a cilium in the direction of propulsion is large and positive during the effective stroke whereas it is small and negative during the recovery stroke. However, the duration of the recovery stroke is longer than the effective stroke so the force exerted over one cycle of a ciliary beat is very small. The bending moment follows a similar pattern to the component of force in the direction of propulsion, being larger in the effective stroke for antiplectic metachronism. In symplectic metachronism (i.e. Opalina) the force and bending moment are largest in magnitude when the bending wave is propagated along the cilium. The rate of working indicates that more energy is consumed in the effective stroke for Paramecium and Pleurobrachia than in the recovery stroke, whereas in Opalina it is found to be large during the propagation of the bending wave.

We investigate the detailed nature of the ‘mixing transition’ through which turbulence may develop in both homogeneous and stratified free shear layers. Our focus is upon the fundamental role in transition, and in particular the associated ‘mixing’ (i.e. small-scale motions which lead to an irreversible increase in the total potential energy of the flow) that is played by streamwise vortex streaks, which develop once the primary and typically two-dimensional Kelvin–Helmholtz (KH) billow saturates at finite amplitude.

Saturated KH billows are susceptible to a family of three-dimensional secondary instabilities. In homogeneous fluid, secondary stability analyses predict that the stream-wise vortex streaks originate through a ‘hyperbolic’ instability that is localized in the vorticity braids that develop between billow cores. In sufficiently strongly stratified fluid, the secondary instability mechanism is fundamentally different, and is associated with convective destabilization of the statically unstable sublayers that are created as the KH billows roll up.

We test the validity of these theoretical predictions by performing a sequence of three-dimensional direct numerical simulations of shear layer evolution, with the flow Reynolds number (defined on the basis of shear layer half-depth and half the velocity difference) Re = 750, the Prandtl number of the fluid Pr = 1, and the minimum gradient Richardson number Ri(0) varying between 0 and 0.1. These simulations quantitatively verify the predictions of our stability analysis, both as to the spanwise wavelength and the spatial localization of the streamwise vortex streaks. We track the nonlinear amplification of these secondary coherent structures, and investigate the nature of the process which actually triggers mixing. Both in stratified and unstratified shear layers, the subsequent nonlinear amplification of the initially localized streamwise vortex streaks is driven by the vertical shear in the evolving mean flow. The two-dimensional flow associated with the primary KH billow plays an essentially catalytic role. Vortex stretching causes the streamwise vortices to extend beyond their initially localized regions, and leads eventually to a streamwise-aligned collision between the streamwise vortices that are initially associated with adjacent cores.

It is through this collision of neighbouring streamwise vortex streaks that a final and violent finite-amplitude subcritical transition occurs in both stratified and unstratified shear layers, which drives the mixing process. In a stratified flow with appropriate initial characteristics, the irreversible small-scale mixing of the density which is triggered by this transition leads to the development of a third layer within the flow of relatively well-mixed fluid that is of an intermediate density, bounded by narrow regions of strong density gradient.

The orientations of fibres in a semi-dilute, index-of-refraction-matched suspension in a Newtonian fluid were observed in a cylindrical Couette device. Even at the highest concentration (nL3 = 45), the particles rotated around the vorticity axis, spending most of their time nearly aligned in the flow direction as they would do in a Jeffery orbit. The measured orbit-constant distributions were quite different from the dilute orbit-constant distributions measured by Anczurowski & Mason (1967b) and were described well by an anisotropic, weak rotary diffusion. The measured ϕ-distributions were found to be similar to Jeffery's solution. Here, ϕ is the meridian angle in the flow-gradient plane. The shear viscosities measured by Bibbo (1987) compared well with the values predicted by Shaqfeh & Fredrickson's theory (1990) using moments of the orientation distribution measured here.