Preston's method of measuring skin friction in the turbulent boundary layer makes use of a circular Pitot tube resting on the wall. On the assumption of a velocity distribution in the wall region common to boundary layer and pipe flows the calibration curve for the Pitot tube can be obtained in fully developed pipe flow. Earlier experiments suggested that Preston's original calibration was in error, and a revised calibration curve has been obtained and is presented here.

From experiments in strong favourable and adverse pressure gradients, limits are assigned to the pressure-gradient conditions within which the calibration can be used with prescribed accuracy. It is shown that in sufficiently strong favourable gradients the ‘inner-law’ velocity distribution breaks down completely, and it is suggested that this breakdown is associated with reversion to laminar flow.

As an incidental result, values have been obtained for the constants occurring in the logarithmic expression for the inner-law velocity distribution.

Ensemble-averaged statistics at constant phase of the turbulent near-wake flow (Reynolds number ≈ 21400 around a square cylinder have been obtained from two-component laser-Doppler measurements. Phase was defined with reference to a signal taken from a pressure sensor located at the midpoint of a cylinder sidewall. The distinction is drawn between the near wake where the shed vortices are ‘mature’ and distinct and a base region where the vortices grow to maturity and are then shed. Differences in length and velocity scales and vortex celerities between the flow around a square cylinder and the more frequently studied flow around a circular cylinder are discussed. Scaling arguments based on the circulation discharged into the near wake are proposed to explain the differences. The relationship between flow topology and turbulence is also considered with vorticity saddles and streamline saddles being distinguished. While general agreement with previous studies of flow around a circular cylinder is found with regard to essential flow features in the near wake, some previously overlooked details are highlighted, e.g. the possibility of high Reynolds shear stresses in regions of peak vorticity, or asymmetries near the streamline saddle. The base region is examined in more detail than in previous studies, and vorticity saddles, zero-vorticity points, and streamline saddles are observed to differ in importance at different stages of the shedding process.

We consider horizontal static liquid layers on planar solid boundaries and analyse their instabilities. The layers are either evaporating, when the plates are heated, or condensing, when the plates are cooled. Vapour recoil, thermocapillary, and rupture instabilities are discussed, along with the effects of mass loss (or gain) and non-equilibrium thermodynamic effects. Particular attention is paid to the development of dryout. We derive long-wave evolution equations for the interface shapes that govern the two-dimensional nonlinear stability of the layers subject to the above coupled mechanisms. These equations are analysed and their predictions discussed. Previous theoretical and experimental results are reviewed and compared with the present results. Finally, we discuss limitations of the modelling and extend our derivation to the case of three-dimensional disturbances.

Data on turbulent mixing and other turbulent-flow phenomena suggest that a (mixing) transition, originally documented to occur in shear layers, also occurs in jets, as well as in other flows and may be regarded as a universal phenomenon of turbulence. The resulting fully-developed turbulent flow requires an outer-scale Reynolds number of Re = Uδ/v [gsim ] 1–2 × 104, or a Taylor Reynolds number of ReT = u′ λT/v [gsim ] 100–140, to be sustained. A proposal based on the relative magnitude of dimensional spatial scales is offered to explain this behaviour.

Consideration of the structure of wind-generated waves when the duration and fetch of the wind are large suggests that the smaller-scale components of the wave field may be in a condition of statistical equilibrium determined by the requirements for attachment of the crests of the waves. A dimensional analysis, based upon the idea of an equilibrium range in the wave spectrum, shows that for large values of the frequency ω, the spectrum Φ(ω) is of the form $\Phi (\omega) \sim \alpha g^2\omega^{-5}$ where α is an absolute constant. The instantaneous spatial spectrum Ψ (k) is proportional to k−4 for large wave numbers k, which is consistent with the observed occurrence of sharp crests in a well-developed sea, and the loss of energy from the wave system to turbulence and heat is proportional to $\rho _w u^3_*$, where ρw is the water density and u the friction velocity of the wind at the surface. This prediction of the form of Φ(ω) for large ω with α = 7·4×10−3, agrees well with measurements made by Burling (1955).

Stokes flow through a random, moderately dense bed of spheres is treated by a generalization of Brinkman's (1947) method, which is applicable to both stationary beds and suspensions. For stationary beds, Darcy's law with a permeability result similar to Brinkman's is derived. For suspensions an effective viscosity μ/(1–2·60ψ) is found, where ψ is the volume fraction of spheres. Also, an expression for the settling velocity is derived.

When two drops of radius R touch, surface tension drives an initially singular motion which joins them into a bigger drop with smaller surface area. This motion is always viscously dominated at early times. We focus on the early-time behaviour of the radius rm of the small bridge between the two drops. The flow is driven by a highly curved meniscus of length 2πrm and width Δ[Lt ]rm around the bridge, from which we conclude that the leading-order problem is asymptotically equivalent to its two-dimensional counterpart. For the case of inviscid surroundings, an exact two-dimensional solution (Hopper 1990) shows that Δ∝r3m and rm∼(tγ/πη) ln [tγ(ηR)]; and thus the same is true in three dimensions. We also study the case of coalescence with an external viscous fluid analytically and, for the case of equal viscosities, in detail numerically. A significantly different structure is found in which the outer-fluid forms a toroidal bubble of radius Δ∝r3/2m at the meniscus and rm∼(tγ/4πη) ln [tγ/(ηR)]. This basic difference is due to the presence of the outer-fluid viscosity, however small. With lengths scaled by R a full description of the asymptotic flow for rm(t)[Lt ]1 involves matching of lengthscales of order r2m, r3/2m, rm, 1 and probably r7/4m.

The commonly used linear K-l and K-ε models of turbulence are shown to be incapable of accurately predicting turbulent flows where the normal Reynolds stresses play an important role. By means of an asymptotic expansion, nonlinear K-l and K-ε models are obtained which, unlike all such previous nonlinear models, satisfy both realizability and the necessary invariance requirements. Calculations are presented which demonstrate that this nonlinear model is able to predict the normal Reynolds stresses in turbulent channel flow much more accurately than the linear model. Furthermore, the nonlinear model is shown to be capable of predicting turbulent secondary flows in non-circular ducts - a phenomenon which the linear models are fundamentally unable to describe. An additional application of this model to the improved prediction of separated flows is discussed briefly along with other possible avenues of future research.

Instability of the alternate-bar type in straight channels has long been identified as the cause of fluvial meandering. The condition of inerodible sidewalls, however, does not allow a meandering channel to develop. Herein a stability analysis of a sinuous channel with erodible banks allows for delineation of a ‘bend’ instability that does not occur in straight channels, and differs from the alternate-bar instability.

In the case of alluvial meanders, the two mechanisms are shown to operate at similar characteristic wavelengths. This provides a rationale for the continuous evolution of alternate bars into true bends such that each bend contains one alternate bar.

The same bend instability applies to incised meanders. A mechanism for incised alternate bars which differs from that for the alluvial case appears to operate at different characteristic wavelengths than that of bend instability. Analysis of data suggests that meandering in supraglacial meltwater streams is primarily due to the alternate bar mechanism, whereas the meandering of rills incised in cohesive material and of caves is likely due to the bend mechanism.

The meander wavelength of incised reaches of meandering streams is often longer than that of adjacent alluvial reaches. An explanation is offered in terms of bend instability.

It is shown that the dilatational terms that need to be modelled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of the compressible velocity field, which generates these dilatational terms, is explored by asymptotic analysis of the compressible Navier-Stokes equations in the case of homogeneous turbulence. The lowest-order equations for the compressible field are solved and explicit expressions for some of the associated one-point moments are obtained. For low Mach numbers, the compressible mode has a fast timescale relative to the incompressible mode. Therefore, it is proposed that, in moderate Mach number homogeneous turbulence, the compressible component of the turbulence is in quasi-equilibrium with respect to the incompressible turbulence. A non-dimensional parameter which characterizes this equilibrium structure of the compressible mode is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

The return to isotropy of homogeneous turbulence without mean velocity gradients is attacked by considering changes to be slow relative to turbulence time scales. This single assumption permits the problem to be cast as one of finding the form of three invariant functions. Examination of limiting behaviour for large Reynolds number and small anisotropy, as well as small Reynolds number and arbitrary anisotropy, places restrictions on the form of the functions. Realizability conditions (requiring that energies be non-negative) reduce the problem to two functions subject to further restrictions. A convenient interpolation form is found for the functions, satisfying all the restrictions, and it is shown that predictions based on this are in excellent agreement with all available data.

The basic property of equations describing dispersive waves is the existence of solutions representing uniform wave trains. In this paper a general theory is given for non-uniform wave trains whose amplitude, wave-number, etc., vary slowly in space and time, the length and time scales of the variation in amplitude, wave-number, etc., being large compared to the wavelength and period. Dispersive equations may be derived from a variational principle with appropriate Lagrangian, and the whole theory is developed in terms of the Lagrangian. Boussinesq's equations for long water waves are used as a typical example in presenting the theory.

The flow around a vertical circular pile exposed to a steady current is studied numerically and experimentally. The numerical model is a three-dimensional model. The model validation was achieved against new experimental data (which include two-component laser-Doppler anemometry (LDA) flow measurements and the hot-film bed shear stress measurements, and reported in the present paper) and the data of others, and a $k$–$\omega $ turbulence model was used for closure. The model does not have a free-surface facility and therefore is applicable only to cases where the Froude number is small ($\hbox{\it Fr}\,{<}\,O(0.2)$). The flow model was used to study the horseshoe vortex and lee-wake vortex flow processes around the pile. The influence on the horseshoe vortex of three parameters, namely the boundary-layer thickness, the Reynolds number and the bed roughness, was investigated. In the latter investigation, the steady solution of the model was chosen. A study of the influence of the unsteady solution on the previously mentioned flow processes was also carried out. The ranges of the parameters covered in the numerical simulations are: The boundary-layer-thickness-to-pile-diameter ratio is varied from 2$\,{\times}\, 10^{-2}$ to 10$^{2}$, the pile Reynolds number from 10$^{2}$ to $2\,{\times}\, 10^{6},$ and the pile diameter-to-roughness ratio from 2 to about 10$^{3}.$ The amplification of the bed shear stress around the pile (including the areas under the horseshoe vortex and the lee-wake region) was obtained for various values of the previously mentioned parameters. The steady-state flow model was coupled with a morphologic model to calculate scour around a vertical circular pile exposed to a steady current in the case of non-cohesive sediment. The morphologic model includes (i) a two-dimensional bed load sediment-transport description, and (ii) a description of surface-layer sand slides for bed slopes exceeding the angle of repose. The results show that the present numerical simulation captures all the main features of the scour process. The equilibrium scour depth obtained from the simulation agrees well with the experiments for the upstream scour hole. Some discrepancy (up to 30%) was observed, however, for the downstream scour hole. The calculations show that the amplification of the bed shear stress around the pile in the equilibrium state of the scour process is reduced considerably with respect to that experienced at the initial stage where the bed is plane.

Considerable discussion over the past few years has been devoted to the question of whether the logarithmic region in wall turbulence is indeed universal. Here, we analyse recent experimental data in the Reynolds number range of nominally $2\times 1{0}^{4} \lt {\mathit{Re}}_{\tau } \lt 6\times 1{0}^{5} $ for boundary layers, pipe flow and the atmospheric surface layer, and show that, within experimental uncertainty, the data support the existence of a universal logarithmic region. The results support the theory of Townsend (The Structure of Turbulent Shear Flow, Vol. 2, 1976) where, in the interior part of the inertial region, both the mean velocities and streamwise turbulence intensities follow logarithmic functions of distance from the wall.

The results of experiments involving instability, known as fingering, in a circular Hele Shaw cell with inward and outward flow are presented. The width of fingers in this situation is examined, and an approximate equation for the growth of fingers is proposed. The equation rα = cos (nθ) is shown to fit the shape of long fingers.

The problem of the response of a porous elastic bed to water waves is treated analytically on the basis of the three-dimensional consolidation theory of Biot (1941). Exact solutions for the pore-water pressure and the displacements of the porous medium are obtained in closed form for the case of waves propagating over the poro-elastic bed. The theoretical results indicate that the bed response to waves is strongly dependent on the permeability k and the stiffness ratio G/K’, where G is the shear modulus of the porous medium and K’ is the apparent bulk modulus of elasticity of the pore fluid. The earlier solutions for pore-water pressure by various authors are given as the limiting cases of the present solution. For the limits G/K′ → 0 or k→ ∞, the present solution for pressure approaches the solution of the Laplace equation by Putnam (1949). For the limit G/K′→ ∞, the present solution approaches the solution of the heat conduction equation by Nakamura et al. (1973) and Moshagen & Tørum (1975).

The theoretical results are compared with wave tank experimental data on pore-water pressure in coarse and fine sand beds which contain small amounts of air. Good agreement between theory and experiment is obtained.

The development of three-dimensional motions in a plane mixing layer was investigated experimentally. It is shown that superimposed on the primary, spanwise vortex structure there is a secondary, steamwise vortex structure. Three aspects of this secondary structure were studied. First, the spanwise vortex instability that generates the secondary structure was characterized by measurements of the critical Reynolds number and the spanwise wavelength at several flow conditions. While the critical Reynolds number was found to depend on the velocity ratio, density ratio and initial shear-layer-profile shape, the mean normalized wavelength is independent of these parameters. Secondly, flow visualization in water was used to obtain cross-sectional views of the secondary structure associated with the streamwise counter-rotating vortices. A model is proposed in which those vortices are part of a single vortex line winding back and forth between the high-speed side of a primary vortex and the low-speed side of the following one. Finally, the effect of the secondary structure on the spanwise concentration field was measured in a helium–nitrogen mixing layer. The spatial organization of the secondary structure produces a well-defined spanwise entrainment pattern in which fluid from each stream is preferentially entrained at different spanwise locations. These measurements show that the spanwise scale of the secondary structure increases with downstream distance.

In this paper, the change in energy dissipation due to a small hump on a body in a uniform steady flow is calculated. The result is used in conjunction with the variational methods of optimal control to obtain the optimality conditions for four minimum-drag problems of fluid mechanics. These conditions imply that the unit-area profile of smallest drag has a front end shaped like a wedge of angle 90°.

A model-based description of the scaling and radial location of turbulent fluctuations in turbulent pipe flow is presented and used to illuminate the scaling behaviour of the very large scale motions. The model is derived by treating the nonlinearity in the perturbation equation (involving the Reynolds stress) as an unknown forcing, yielding a linear relationship between the velocity field response and this nonlinearity. We do not assume small perturbations. We examine propagating helical velocity response modes that are harmonic in the wall-parallel directions and in time, permitting comparison of our results to experimental data. The steady component of the velocity field that varies only in the wall-normal direction is identified as the turbulent mean profile. A singular value decomposition of the resolvent identifies the forcing shape that will lead to the largest velocity response at a given wavenumber–frequency combination. The hypothesis that these forcing shapes lead to response modes that will be dominant in turbulent pipe flow is tested by using physical arguments to constrain the range of wavenumbers and frequencies to those actually observed in experiments. An investigation of the most amplified velocity response at a given wavenumber–frequency combination reveals critical-layer-like behaviour reminiscent of the neutrally stable solutions of the Orr–Sommerfeld equation in linearly unstable flow. Two distinct regions in the flow where the influence of viscosity becomes important can be identified, namely wall layers that scale with R+1/2 and critical layers where the propagation velocity is equal to the local mean velocity, one of which scales with R+2/3 in pipe flow. This framework appears to be consistent with several scaling results in wall turbulence and reveals a mechanism by which the effects of viscosity can extend well beyond the immediate vicinity of the wall. The model reproduces inner scaling of the small scales near the wall and an approach to outer scaling in the flow interior. We use our analysis to make a first prediction that the appropriate scaling velocity for the very large scale motions is the centreline velocity, and show that this is in agreement with experimental results. Lastly, we interpret the wall modes as the motion required to meet the wall boundary condition, identifying the interaction between the critical and wall modes as a potential origin for an interaction between the large and small scales that has been observed in recent literature as an amplitude modulation of the near-wall turbulence by the very large scales.

In this paper an attempt is made to formulate a model for the mechanism of wall turbulence that links recent flow-visualization observations with the various quantitative measurements and scaling laws established from anemometry studies. Various mechanisms are proposed, all of which use the concept of the horse-shoe, hairpin or ‘A’ vortex. It is shown that these models give a connection between the mean-velocity distribution, the broad-band turbulence-intensity distributions and the turbulence spectra. Temperature distributions above a heated surface are also considered. Although this aspect of the work is not yet complete, the analysis for this shows promise.