The ability of a large-eddy simulation to represent the large-scale motions in the interior of a turbulent flow is well established. However, concerns remain for the behaviour close to rigid surfaces where, with the exception of low-Reynolds-number flows, the large-eddy description must be matched to some description of the flow in which all except the larger-scale ‘inactive’ motions are averaged. The performance of large-eddy simulations in this near-surface region is investigated and it is pointed out that in previous simulations the mean velocity profile in the matching region has not had a logarithmic form. A number of new simulations are conducted with the Smagorinsky (1963) subgrid model. These also show departures from the logarithmic profile and suggest that it may not be possible to eliminate the error by adjustments of the subgrid lengthscale. An obvious defect of the Smagorinsky model is its failure to represent stochastic subgrid stress variations. It is shown that inclusion of these variations leads to a marked improvement in the near-wall flow simulation. The constant of proportionality between the magnitude of the fluctuations in stress and the Smagorinsky stresses has been empirically determined to give an accurate logarithmic flow profile. This value provides an energy backscatter rate slightly larger than the dissipation rate and equal to idealized theoretical predictions (Chasnov 1991).

We report on lattice-Boltzmann simulations of slow fluid flow past mono- and bidisperse random arrays of spheres. We have measured the drag force on the spheres for a range of diameter ratios, mass fractions and packing fractions; in total, we studied 58 different parameter sets. Our simulation data for the permeability agrees well with previous simulation results and the experimental findings. On the basis of our data for the individual drag force, we have formulated new drag force relations for both monodisperse and polydisperse systems, based on the Carman–Kozeny equations; the average deviation of our binary simulation data with the new relation is less than 5%. We expect that these new relations will significantly improve the numerical modelling of gas–solid systems with polydisperse particles, in particular with respect to mixing and segregation phenomena. For binary systems with large diameter ratios (1:4), the individual drag force on a particle, as calculated from our relations, can differ by up to a factor of five compared with predictions presently favoured in chemical engineering.

A visualization method is used to obtain the main features of the hydrodynamic field for flow past a circular cylinder moving at a uniform speed in a direction perpendicular to its generating lines in a tank filled with a viscous liquid. Photographs are presented to show the particular fineness of the experimental technique. More especially, the closed wake and the velocity distribution behind the obstacle are investigated; the changes in the geometrical parameters describing the eddies with Reynolds number (5 < Re < 40) and with the ratio λ between the diameters of the cylinder and tank are given. A comparison with existing numerical and experimental results is presented and some remarks are made about the calculation techniques proposed up to the present. The limits of the Reynolds-number range for which the twin vortices exist and adhere stably to the cylinder are determined.

A macroscopic equation of mass conservation is obtained by ensemble-averaging the basic conservation laws in a porous medium. In the long-time limit this ‘macro-transport’ equation takes the form of a macroscopic Fick's law with a constant effective diffusivity tensor. An asymptotic analysis in low volume fraction of the effective diffusivity in a bed of fixed spheres is carried out for all values of the Péclet number ℙ = Ua/Df, where U is the average velocity through the bed. a is the particle radius and Df is the molecular diffusivity of the solute in the fluid. Several physical mechanisms causing dispersion are revealed by this analysis. The stochastic velocity fluctuations induced in the fluid by the randomly positioned bed particles give rise to a convectively driven contribution to dispersion. At high Péclet numbers, this convective dispersion mechanism is purely mechanical, and the resulting effective diffusivities are independent of molecular diffusion and grow linearly with ℙ. The region of zero velocity in and near the bed particles gives rise to non-mechanical dispersion mechanisms that dominate the longitudinal diffusivity at very high Péclet numbers. One such mechanism involves the retention of the diffusing species in permeable particles, from which it can escape only by molecular diffusion, leading to a diffusion coefficient that grows as ℙ2. Even if the bed particles are impermeable, non-mechanical contributions that grow as ℙ ln ℙ and ℙ2 at high ℙ arise from a diffusive boundary layer near the solid surfaces and from regions of closed streamlines respectively. The results for the longitudinal and transverse effective diffusivities as functions of the Péclet number are summarized in tabular form in §6. Because the same physical mechanisms promote dispersion in dilute and dense fixed beds, the predicted Péclet-number dependences of the effective diffusivities are applicable to all porous media. The theoretical predictions are compared with experiments in densely packed beds of impermeable particles, and the agreement is shown to be remarkably good.

A theoretical method is given for the determination of the shape of a fluid drop in steady and unsteady flows by making an expansion in terms of the drop deformation. Effects of fluid viscosity and interfacial tension are taken into account. Examples given include the determination of the shape of a drop in shear and in hyperbolic flow when each is started impulsively from rest.

The primary objective of this investigation is to determine experimentally the sources of jet mixing noise. In the present study, four different approaches are used. It is reasonable to assume that the characteristics of the noise sources are imprinted on their radiation fields. Under this assumption, it becomes possible to analyse the characteristics of the far-field sound and then infer back to the characteristics of the sources. The first approach is to make use of the spectral and directional information measured by a single microphone in the far field. A detailed analysis of a large collection of far-field noise data has been carried out. The purpose is to identify special characteristics that can be linked directly to those of the sources. The second approach is to measure the coherence of the sound field using two microphones. The autocorrelations and cross-correlations of these measurements offer not only valuable information on the spatial structure of the noise field in the radial and polar angle directions, but also on the sources inside the jet. The third approach involves measuring the correlation between turbulence fluctuations inside a jet and the radiated noise in the far field. This is the most direct and unambiguous way of identifying the sources of jet noise. In the fourth approach, a mirror microphone is used to measure the noise source distribution along the lengths of high-speed jets. Features and trends observed in noise source strength distributions are expected to shed light on the source mechanisms. It will be shown that all four types of data indicate clearly the existence of two distinct noise sources in jets. One source of noise is the fine-scale turbulence and the other source is the large turbulence structures of the jet flow. Some of the salient features of the sound field associated with the two noise sources are reported in this paper.

The onset of periodic behaviour in two-dimensional laminar flow past bodies of various shapes is examined by means of finite-element simulations. The transition from steady to periodic flow is marked by a Hopf bifurcation, which we locate by solving an appropriate extended set of steady-state equations. The bodies considered are a circular cylinder, triangular prisms of various shapes, and flat plates and elliptical cylinders aligned over a range of angles to the direction of flow. Our results for the circular cylinder are in good agreement with experimental observations and with the results of time-dependent calculations.

The non-equilibrium behaviour of concentrated colloidal dispersions is studied using Stokesian Dynamics, a molecular-dynamics-like simulation technique for analysing suspensions of particles immersed in a Newtonian fluid. The simulations are of a monodisperse suspension of Brownian hard spheres in simple shear flow as a function of the Péclet number, Pe, which measures the relative importance of hydrodynamic and Brownian forces, over a range of volume fraction 0.316 [les ] ϕ [les ] 0.49. For Pe < 10, Brownian motion dominates the behaviour, the suspension remains well-dispersed, and the viscosity shear thins. The first normal stress difference is positive and the second negative. At higher Pe, hydrodynamics dominate resulting in an increase in the long-time self-diffusivity and the viscosity. The first normal stress difference changes sign when hydrodynamics dominate. Simulation results are shown to agree well with both theory and experiment.

This paper summarizes experimental studies of incompressible elliptic jets of different aspect ratios and initial conditions, and effects of excitations at selected frequencies and amplitudes. Elliptic jets are quite different from the extensively studied plane and circular jets - owing mainly to the fact that the azimuthal curvature variation of a vortical structure causes its non-uniform self-induction and hence complex three-dimensional deformation. Such deformation, combined with properly selected excitation can substantially alter entrainment and other turbulence phenomena, thus suggesting preference for the elliptic shape in many jet applications. The dominance of coherent structures in the jet far field is evident from the finding that switching over of the cross-section shape continues at least up to 100 equivalent diameters De. The locations and the number of switchovers are strongly dependent on the initial condition, on the aspect ratio, and, when excited, on the Strouhal number and the excitation level. We studied jets with constant exit momentum thickness θe, all around the perimeter, thus separating the effects of azimuthal variations of θe, (typical of elliptic jets) and of the shear-layer curvature. Also investigated are the instability characteristics, and enhanced entrainment caused by bifurcation as well as pairing of vortical structures. We discuss shear-layer and jet- column domains, and find the latter to be characterized by two modes : the preferred mode and the stable pairing mode - similar to those found in circular jets -both modes scaling on the newly-defined lengthscale De. The paper documents some time- average measurements and their comparison with those in circular and plane jets.

The paper describes studies of the turbulence of the liquid in a bubbly, grid-generated turbulent flow field. Laser-Doppler and hot-film anemometry are used for the experimental investigation. It is found that the turbulent kinetic energy increases strongly with the void fraction α. Roughly speaking, there exist two distinct regimes: the first one corresponds to low value of α, where hydrodynamic interactions between bubbles are negligible, and the second one to higher values, for which, owing to their mutual interactions, the bubbles transfer a greater amount of kinetic energy to the liquid. The Reynolds stress tensor shows that the quasi-isotropy is not altered. At low enough values of α, the difference between the turbulent kinetic energy in the liquid phase and the energy associated with the grid-generated turbulence proves to be approximately equal to the intensity of the pseudo-turbulence, defined as the fluctuating energy that would be induced by the motion of the bubbles under non-turbulent conditions. The one-dimensional spectra exhibit a large range of high frequencies associated with the wakes of the bubbles and the classical $-\frac{5}{3}$ power law is progressively replaced by a $-\frac{8}{3}$ dependence.

A formal algebraic analogy is drawn between meteorological parameters, such as the Richardson number, and the parameters describing the effect of rotation or streamline curvature on a turbulent flow. The analogy between the phenomena is a good first approximation. Semi-quantitative use of the analogy to apply meteorological data to curved shear layers shows that the effects of curvature on the apparent mixing length are appreciable if the shear-layer thickness exceeds roughly 1/300 of the radius of curvature; larger effects may occur in compressible flow. Application of the Monin-Oboukhov formula considerably improves the agreement between prediction and experiment in boundary layers on curved surfaces.

An apparatus is described for determining particle passages through any selected point on a tube cross-section. The method depends on the simultaneous blocking out of two mutually perpendicular light beams by a particle passing through their common region. The coincidence is registered and counted electronically. At higher particle concentrations coincidences are also registered arising from a pair-wise occupation of the light beams by two particles. An analysis is presented showing that these pair coincidences can be allowed for exactly in terms of experimentally measurable quantities.

The reliability and reproducibility of the method is discussed and illustrated by examples from sphere suspensions in Poiseuille flow.

In most work on the theory of stability of laminar flow, infinitesimal disturbances only have been considered, so that only the initial growth of the disturbance has been determined. It is the object of the present paper to extend the theory to larger amplitudes and to study the mechanics of disturbance growth with the inherent non-linearity of the hydrodynamical system taken into account.

The Reynolds stress (where averages are taken with respect to some suitable space coordinate) is the fundamental consequence of the non-linearity, and its effects can be anticipated as follows. Initially a disturbance grows exponentially with time according to the linear theory, but eventually it reaches such a size that the transport of momentum by the finite fluctuations is appreciable and the associated mean stress (the Reynolds stress) then has an appreciable effect on the mean flow. This distortion of the mean flow modifies the rate of transfer of energy from the mean flow to the disturbance and, since this energy transfer is the cause of the growth of the disturbance, there is a modification of the rate of growth of the latter.

It is suggested that, in many cases, an equilibrium state may be possible in which the rate of transfer of energy from the (distorted) mean flow to the disturbance balances precisely the rate of viscous dissipation of the energy of disturbance. A theory based on certain assumptions about the energy flow is given to describe both the growth of the disturbance and the final equilibrium state, and application is made to the cases of Poiseuille flow between parallel planes and flow between rotating cylinders. The distorted mean flow in the equilibrium state can be calculated and from this, in the latter case, the torque required to maintain the cylinders in motion. Good agreement is obtained with G. I. Taylor's measurements of the torque for the case when the inner cylinder rotates and the outer cylinder is at rest.

A non-inertial (zero Taylor number) viscoelastic instability is discovered for Taylor–Couette flow of dilute polymer solutions. A linear stability analysis of the inertialess flow of an Oldroyd-B fluid (using both approximate Galerkin analysis and numerical solution of the relevant small-gap eigenvalue problem) show the growth of an overstable (oscillating) mode when the Deborah number exceeds f(S) ε−½, where ε is the ratio of the gap to the inner cylinder radius, and f(S) is a function of the ratio of solvent to polymer contributions to the solution viscosity. Experiments with a solution of 1000 p.p.m. high-molecular-weight polyisobutylene in a viscous solvent show an onset of secondary toroidal cells when the Deborah number De reaches 20, for ε of 0.14, and a Taylor number of 10−6, in excellent agreement with the theoretical value of 21. The critical De was observed to increase as ε decreases, in agreement with the theory. At long times after onset of the instability, the cells become small in wavelength compared to those that occur in the inertial instability, again in agreement with our linear analysis. For this fluid, a similar instability occurs in cone-and-plate flow, as reported earlier. The driving force for these instabilities is the interaction between a velocity fluctuation and the first normal stress difference in the base state. Instabilities of the kind that we report here are likely to occur in many rotational shearing flows of viscoelastic fluids.

Direct numerical solutions of the three-dimensional time-dependent Navier-Stokes equations are presented for the evolution of three-dimensional finite-amplitude disturbances of plane Poiseuille and plane Couette flows. Spectral methods using Fourier series and Chebyshev polynomial series are used. It is found that plane Poiseuille flow can sustain neutrally stable two-dimensional finite-amplitude disturbances at Reynolds numbers larger than about 2800. No neutrally stable two-dimensional finite-amplitude disturbances of plane Couette flow were found.

Three-dimensional disturbances are shown to have a strongly destabilizing effect. It is shown that finite-amplitude disturbances can drive transition to turbulence in both plane Poiseuille flow and plane Couette flow at Reynolds numbers of order 1000. Details of the resulting flow fields are presented. It is also shown that plane Poiseuille flow cannot sustain turbulence at Reynolds numbers below about 500.

We study experimentally the behaviour of a drop deposited on a conical fibre. It is shown that for wetting liquids, such a drop spontaneously moves towards the region of lower curvature. The driving force is measured and shown to be a gradient of Laplace pressure, which allows us to characterize the dynamics of these self-propelling drops. We conclude by discussing the efficiency of this device for drying a solid initially coated with a liquid film.

A conceptual framework for analysing the energetics of density-stratified Boussinesq fluid flows is discussed. The concept of gravitational available potential energy is used to formulate an energy budget in which the evolution of the background potential energy, i.e. the minimum potential energy attainable through adiabatic motions, can be explicitly examined. For closed systems, the background potential energy can change only due to diabatic processes. The rate of change of background potential energy is proportional to the molecular diffusivity. Changes in the background potential energy provide a direct measure of the potential energy changes due to irreversible diapycnal mixing. For open systems, background potential energy can also change due to boundary fluxes, which can be explicitly measured. The analysis is particularly appropriate for evaluation of diabatic mixing rates in numerical simulations of turbulent flows. The energetics of a shear driven mixing layer is used to illustrate the analysis.

Experiments have been performed on the tip vortex trailing from a rectangular NACA 0012 half-wing. Preliminary studies showed the vortex to be insensitive to the introduction of a probe and subject only to small wandering motions. Meaningful velocity measurements could therefore be made using hot-wire probes.

Detailed analysis of the effects of wandering was performed to properly reveal the flow structure in the core region and to give confidence in measurements made outside the core. A theory has been developed to correct mean-velocity profiles for the effects of wandering and to provide complete quantitative estimates of its amplitude and contributions to Reynolds stress fields. Spectral decomposition was found to be the most effective method of separating these contributions from velocity fluctuations due to turbulence.

Outside the core the flow structure is dominated by the remainder of the wing wake which winds into an ever-increasing spiral. There is no large region of axisymmetric turbulence surrounding the core and little sign of turbulence generated by the rotational motion of the vortex. Turbulence stress levels vary along the wake spiral in response to the varying rates of strain imposed by the vortex. Despite this complexity, the shape of the wake spiral and its turbulent structure reach an approximately self-similar form.

On moving from the spiral wake to the core the overall level of velocity fluctuations greatly increases, but none of this increase is directly produced by turbulence. Velocity spectra measured at the vortex centre scale in a manner that implies that the core is laminar and that velocity fluctuations here are a consequence of inactive motion produced as the core is buffeted by turbulence in the surrounding spiral wake. Mean-velocity profiles through the core show evidence of a two-layered structure that dies away with distance downstream.

In this paper we present a rather personal view of the important developments in double-diffusive convection, a subject whose evolution has been the result of a close interaction between theoreticians, laboratory experimenters and sea-going oceano-graphers. More recently, applications in astrophysics, engineering and geology have become apparent. In the final section we attempt to draw some general conclusions and suggest that further progress will again depend on a close collaboration between fluid dynamicists and other scientists.

The flow around a circular cylinder placed at various heights above a plane boundary has been investigated experimentally. The cylinder spanned the test section of a wind tunnel and was aligned with its axis parallel to a long plate and normal to the free stream. It was placed 36 diameters downstream of the leading edge of the plate and its height above the plate was varied from zero, the cylinder lying on the surface, to 3·5 cylinder diameters. The thickness of the turbulent boundary layer on the plate at the cylinder position, but with it removed from the tunnel, was equal to 0·8 of the cylinder diameter. Distributions of mean pressure around the cylinder and along the plate were measured at a Reynolds number, based on cylinder diameter, of 4·5 × 104. Spectral analysis of hot-wire signals demonstrated that regular vortex shedding was suppressed for all gaps less than about 0·3 cylinder diameters. For gaps greater than 0·3 the Strouhal number was found to be remarkably constant and the only influence of the plate on vortex shedding was to make it a more highly tuned process as the gap was reduced. Flow-visualization experiments in a smoke tunnel revealed the wake structure at various gap-to-diameter ratios.