A numerical model is proposed for the flow in a deep turbulent boundary layer over water waves. The momentum equations are closed by the use of an isotropic eddy viscosity and the turbulent energy equation. For small amplitudes the results are similar to those of Townsend's (1972) linear model, but nonlinear effects become important as the ratio of wave height to wavelength increases. With uniform surface roughness zo, the predicted fractional rate of energy input per radian advance in phase, ζ, decreases slightly with increasing amplitude and is of the same order of magnitude as in Miles’ (1957, 1959) and Townsend's linear theories. If zo is allowed to vary with position along the wave, however, the fractional rate of energy input can be significantly increased for small amplitude waves. If the variation in zo is half the mean value and the maximum wave slope zak is 0.01, we find ζ ≈ 60 (ρair/ρwater) (uo/c)2, where uo is the friction velocity and c the wave phase speed. Comparison is also made with recent laboratory and field data.

The reported experimental coefficients describing the spreading rate and virtual origin of a plane turbulent jet exhibit substantial scatter. The hypothesis is made that the basic reason for these variations is that the growth of the jet is not linear on a large scale. Existing experimental results are used to support this hypothesis.

The bursting process in turbulent boundary layers provides new insight on turbulence phenomena, mechanics of sedimentation, and genesis of bedforms in natural geophysical flows. Recent visualization experiments suggest that the turbulent boundary layer can be divided into an inner zone, whose essential characteristics scale with inner (wall) variables, and an outer zone, whose properties scale with the fluid-dynamic variables of the entire flow. The inner zone is distinguished by (i) a viscous sublayer displaying spanwise alternations of high-and low-speed streaks and (ii) episodic disruption by lift-ups of low-speed streaks. Oscillatory growth and breakup stages of the Stanford model of bursting characterize the turbulent structure of the outer zone. The burst cycle exists in turbulent boundary layers of all natural flows except perhaps (i) open-channel flows in the upper part of the upper flow regime and (ii) wind-generated surface waves.

Fluid motions described as kolks and boils in incompressible open-channel flows correspond to the oscillatory growth stage and the late oscillatory growth and breakup stages, respectively, of the Stanford model of bursting. Supporting evidence includes (i) close similarity of gross fluid motions, (ii) equivalent scaling of boils and bursts, and (iii) intensification of boils and bursts in adverse pressure gradients and over rough beds. McQuivey's (1973) turbulence measurements show that the Eulerian integral time scale TE scales with the same outer variables as boil periodicity and burst periodicity. It is hypothesized that TE equals the mean duration of bursts at a point in the flow.

Bedforms governed by the turbulent structure of the inner zone (microforms) cannot form if the sublayer is disrupted by bed roughness. The conditions for the existence of two common microforms and their spacings scale with the inner variables. Grain roughness increases the vertical intensity of the turbulence (by enhancing lift-ups) within the inner zone, thereby explaining textural differences between the coarse ripple and fine ripple bed stages of Moss (1972).

Mesoforms respond to the fluid-dynamical regime in the outer zone and scale with the outer variables. The mean spacing of dunelike large-scale ripples in equilibrium open-channel flows is proportional to the boundary-layer thickness and equals the length scale formed by the product of the free-stream velocity and the boil period.

Strong upward flow in a burst provides the vertical anisotropy in the turbulence which is needed to suspend sediment. Bursting promotes the entrain-ment of more and coarser sediment than tractive forces alone can accomplish.

Approximations for shallow-water ship waves are sought which are valid near the critical speed U = (gh)1/2. For sufficiently thin bodies (struts or ships) the governing equation is dispersive. Simple analytic solutions are given which are valid for all F2 ≤ O(1). As the thickness increases, nonlinearity also enters. A soliton solution is discussed which applies to a sharp-nosed half-body at slightly supercritical speed.

An experiment has been carried out in which air bubbles were caused to grow isothermally at a vertical wall between two closely spaced horizontal plates. The experiment gives an approximate representation in two dimensions for the growth of a vapour bubble at a wall under some conditions of subcooled boiling in zero gravity. Although the effect of gravity was virtually eliminated in the experiment, it was found that a bubble could still detach itself from the wall, apparently owing to the effects of surface tension and inertia.

Also, bubbles were seen to detach from a wall despite the presence of a slight gravity force directed to oppose such detachment.

Computer simulations of fluid-element trajectories in mirror-symmetric and maximally helical turbulence are used to evaluate Moffatt's (1974) formulae for the magnetic diffusivity η(t) and the coefficient κ(t) of the alpha-effect. The passive-scalar diffusivity κ(t) and the mean response functions of scalar and magnetic field wave-vector modes are also computed. The velocity field is normal, stationary, homogeneous and isotropic with spectrum $E(k) = \frac{3}{2}v^2_0\delta(k - k_0)$ and time correlation exp [−1/2ω2/0(t-t)2]. The cases ω0 = O (frozen turbulence), ω = vo k0 and ω0 = 2v0 k0 are followed to t = 4/v0 k0. In the ω0 > 0 cases with maximal helicity, κ(t) and a(t) approach steady-state values of order vo/k0 and v0, respectively. They behave anomalously for ω0 = 0. In the mirror-symmetric: cases, q(t) and κ(t) differ very little from each other. At all the ω0 values, is bigger in the helical than in the mirror-symmetric case. The difference is marked for ω0 = 0. The simulation results imply that κ(t) becomes negative in non-normal mirror-symmetric turbulence with strong helicity fluctuations that persist over several correlation lengths and times. The computations of response functions indicate that asymptotic expressions for these functions, valid for k [Lt ]; k0, retain good accuracy for k ∼ k0. The mean-square magnetic field is found to grow exponentially, and its kurtosis also grows apidly with t, indicating rapid development of a highly intermittent distribution of magnetic field.

A higher-order lifting-line theory for a wing in unsteady motion is discussed. Apart from the addition of higher-order terms, it also differs from the theory derived by James (1975) in its emphasis on ‘physical’ interpretations. This emphasis has made it possible to shed some new light on Prandtl's classical lifting-line theory, as well as on Weissinger's 3/4-chord theory.

A new approximate theory is proposed for treating the flow past smoothly contoured two-dimensional bluff bodies in the intermediate Reynolds number range O(1) < Re < 0(102), where the displacement effect of the thick viscous layer near the surface of the body is large and a steady laminar wake is present. The theory is based on a new pressure hypothesis which enables one to take account of the displacement interaction and centrifugal effects in thick viscous layers using conventional first-order boundary-layer equations. The basic question asked is how the wall pressure gradient in ordinary boundary -layer theory must be modified if the pressure gradient along the displacement surface using the Prandtl pressure hypothesis is to be equal to the pressure gradient along this surface using a higher-order approximation to the Navier-Stokes equation in which centrifugal forces are considered. The inclusion of the normal pressure field with displacement interaction is shown to be equivalent to stretching the streamwise body co-ordinate in first-order boundary-layer theory such that the streamwise pressure gradient as a function of distance along the original and displacement body surfaces are equal.

While the new theory is of a non-rigorous nature, it yields results for the location of separation and detailed surface pressure and vorticity distribution which are in remarkably good agreement with the large body of available numerical Navier-Stokes solutions. A novel feature of the new boundary-value problem is the development of a simple but accurate approximate method for determining the inviscid flow past an arbitrary two-dimensional displacement body with its wake.

A study has been made of baroclinic instability in a differentially heated, rotating fluid annulus whose channel width varies azimuthally. Both laboratory experiments and an a.nalytica1 model employing a linear normal-mode analysis have been used. The experiments show three types of flow. For slow rotation the flow is 'symmetric’, whereas at high rotation speeds baroclinic waves occur at all azimuths. At intermediate rotation speeds it is possible to have a mixed flow which is ‘symmetric’ in the narrow part but has baroclinic waves in the wide part of the annulus. This result suggested the analytical investigation of the stability of a barocIinic flow whose meridional scale varies downstream. It was found that this model also permits three possible types of flow: everywhere stable, everywhere unstable, and also a mixed flow which is locally unstable where the meridional scale is largest but locally stable where the scale is smallest.

The selective withdrawal of a layered fluid from a reservoir has interesting properties which are not present if the stratification is continuous. The problem investigated is that of finding the ratios of the discharges from each flowing layer when the total discharge from all layers is regulated at the outlet. It is shown that the ratios are determined by the requirement that the flow be smooth at critical points in the flow, each critical point being a point at which the long-wave velocity on one of the interfaces is zero. If there are too many critical points, the flow is over-determined and becomes unsteady. It follows therefore that if the external conditions change slowly, such as in the slow draining of a reservoir, the flow alternates between intervals in which the flow ratios also change slowly and intervals in which the flow ratios oscillate unsteadily. An experiment is described in which the theoretical conclusions were tested for two-layer flow. Good agreement was obtained between theory and experiment.

The stability of laminar and turbulent channel flow is examined for cases where Coriolis forces are introduced by steady rotation about an axis perpendicular to the plane of mean flow. Linearized equations of motion are derived for small disturbances of the Taylor type. Conditions for marginal stability in laminar Couette and Poiseuille flow correspond, in part, to the analogous solutions of buoyancy-driven convection instabilities in heated fluid layers, and to those of Taylor instabilities in the flow between rotating cylinders. In plane Poiseuille flow with rotation, the critical disturbance mode occurs at a Reynolds number of Rec = 88.53 and rotation number Ro = 0.5. At higher Reynolds numbers, unstable conditions canexist over the range of rotation numbers given by 0 < Ro < 3, provided the undisturbed flow remains laminar. A two-layer model is devised to investigate the onset of longitudinal instabilities in turbulent flow. The linear disturbance equations are solved essentially in their laminar form, whereby the velocity gradient of laminar flow is replaced by a numerically computed profile for the gradient of the turbulent mean velocity. The turbulent stress levels in the stable and unstable flow regions are represented by integrated averages of the eddy viscosity. Onset of instability for Reynolds numbers between 6000 and 35 000 is predicted to occur at Ro = 0.022, a value in remarkable agreement with the experimentally observed appearance of roll instabilities in rotating turbulent channel flow.

To understand the turbulent generation of large-scale magnetic fields and to advance beyond purely kinematic approaches to the dynamo effect like that introduced by Steenbeck, Krause & Radler (1966)’ a new nonlinear theory is developed for three-dimensional, homogeneous, isotropic, incompressible MHD turbulence with helicity, i.e. not statistically invariant under plane reflexions. For this, techniques introduced for ordinary turbulence in recent years by Kraichnan (1971 a)’ Orszag (1970, 1976) and others are generalized to MHD; in particular we make use of the eddy-damped quasi-normal Markovian approximation. The resulting closed equations for the evolution of the kinetic and magnetic energy and helicity spectra are studied both theoretically and numerically in situations with high Reynolds number and unit magnetic Prandtl number.

Interactions between widely separated scales are much more important than for non-magnetic turbulence. Large-scale magnetic energy brings to equipartition small-scale kinetic and magnetic excitation (energy or helicity) by the ‘Alfvén effect’; the small-scale ‘residual’ helicity, which is the difference between a purely kinetic and a purely magnetic helical term, induces growth of large-scale magnetic energy and helicity by the ‘helicity effect’. In the absence of helicity an inertial range occurs with a cascade of energy to small scales; to lowest order it is a −3/2 power law with equipartition of kinetic and magnetic energy spectra as in Kraichnan (1965) but there are −2 corrections (and possibly higher ones) leading to a slight excess of magnetic energy. When kinetic energy is continuously injected, an initial seed of magnetic field will grow to approximate equipartition, at least in the small scales. If in addition kinetic helicity is injected, an inverse cascade of magnetic helicity is obtained leading to the appearance of magnetic energy and helicity in ever-increasing scales (in fact, limited by the size of the system). This inverse cascade, predicted by Frisch et al. (1975), results from a competition between the helicity and Alféh effects and yields an inertial range with approximately — 1 and — 2 power laws for magnetic energy and helicity. When kinetic helicity is injected at the scale linj and the rate $\tilde{\epsilon}^V$ (per unit mass), the time of build-up of magnetic energy with scale L [Gt ] linj is $t \approx L(|\tilde{\epsilon}^V|l^2_{\rm inj})^{-1/3}.$

In Poiseuille flow in a circular tube passive contaminant initially spread uniformly over the cross-section would be pulled out in a paraboloidal snout in the absence of any diffusive mechanism, and there would be a discontinuity in $\overline{C}$, the mean concentration over the cross-section, associated with the contaminant at the front of the snout. In reality molecular diffusion smooths out this snout in two ways: direct longitudinal diffusion and the interaction between lateral diffusion and advection. The effect of these two mechanisms is discussed, and determined for small values of κt/a2, where t is the time since injection, κ is the molecular diffusivity and a is the tube radius. For such values, important in many applications, the tube walls play no part in the smoothing process. It is shown that for $\kappa t/a^2 < 0.25(\overline{u}a/\kappa)^{-\frac{2}{3}}$, where $\overline{u}$ is the discharge velocity, the effect of longitudinal diffusion dominates over that of the interaction, which is, in turn, dominant for $\kappa t/a^2 > 2.5(\overline{u}a/\kappa)^{-\frac{2}{3}}$, when $\overline{C}$ is close to the form described by Lighthill (1966).

A simplified theoretical model of a collapsing mixed region in a stably stratified fluid, as studied experimentally by Wu (1969), is presented. Unlike previous work the model describes the mixed region itself in the ‘principal stage’ of collapse. An inertia-buoyancy balance assumption is invoked, resulting in good quantitative agreement with Wu's results.

The interaction between short internal gravity waves and a larger-scale mean flow in the ocean is analysed in the Wkbj approximation. The wave field determines the radiation-stress term in the momentum equation of the mean flow and a similar term in the buoyancy equation. The mean flow affects the propagation characteristics of the wave field. This cross-coupling is treated as a small perturbation. When relaxation effects within the wave field are considered, the mean flow induces a modulation of the wave field which is a linear functional of the spatial gradients of the mean current velocity. The effect that this modulation itself has on the mean flow can be reduced to the addition of diffusion terms to the equations for the mass and momentum balance of the mean flow. However, there is no vertical diffusion of mass and other passive properties. The diffusion coefficients depend on the frequency spectrum and the relaxation time of the internal-wave field and can be evaluated analytically. The vertical viscosity coefficient is found to be vv [ape ] 4 x 103cm2/s and exceeds values typically used in models of the general circulation by at least two orders of magnitude.

A turbulent boundary layer developing on a smooth heated uniform-temperature plate in a zero pressure gradient was set up. The origins of the layers were matched to remove the effect of an In heated starting length. Similarity proposals were tested. The mean flow field followed the usual law of the wall and defect law for both temperature and velocity. Broad-band measurements of stream wise velocity and temperature fluctuations were made, and wall similarity and Townsend's self-preserving flow similarity were found to be applicable, at least after a sufficient flow development.

Some initial attempts to arrive at a comparison between heat and momentum transport are presented. The results include conditionally sampled measurements of instantaneous heat and momentum fluxes and correlations between these two quantities. The fluxes were divided into quadrants. Conditional probabilities and weighted joint probability density functions were measured to determine whether there was a similarity in behaviour of these two fluxes. The concept of ‘hole size’ developed for momentum flux was extended to heat flux and events corresponding to bursts and sweeps in the momentum flux were found to be accompanied by corresponding events in the heat flux.

A surface-wave/internal-wave mode coupled model is constructed to describe the energy transfer from a linear surface wave field on the ocean to a linear internal wave field. Expressed in terms of action-angle variables the dynamic equations have a particularly useful form and are solved both numerically and in some analytic approximations. The growth time for internal waves generated by the resonant interaction of surface waves is calculated for an equilibrium spectrum of surface waves and for both the Garrett-Munk and two-layer models of the undersea environment. We find energy transfer rates as a function of undersea parameters which are much faster than those based on the constant Brunt-ViiisSila model used by Kenyon (1968) and which are consistent with the experiments of Joyce (1974). The modulation of the surface-wave spectrum by internal waves is also calculated, yielding a ‘mottled’ appearance of the ocean surface similar to that observed in photographs taken from an ERTS1 satellite (Ape1 et al. 1975b).

Conditions leading to the onset of air-flow separation over a mobile air-water interface are discussed. It is argued that, in a frame of reference in which the interfacial boundary assumes a steady shape, the occurrence of separation requires a stagnation point on the interface. In the case of air blowing over water waves, this corresponds to the onset of wave breaking. These arguments are strongly supported by flow visualization and pressure measurements carried out in a laboratory wind-wave flume. Furthermore, the pressure measurements show a greatly enhanced interfacial shear stress for a breaking wave compared with that over an unbroken wave of the same wavelength. The implications of these findings for wind-wave generation are discussed.

What transient motions occur as a fluid responds to gravitational forces in a rotating channel, and what equilibrium does the fluid adjust to? This problem is studied to illustrate how boundaries affect the process of adjustment to a geostrophic equilibrium. In particular, linear solutions are found for an infinitely long channel of constant width and constant depth when there is an initial discontinuity in the level of the free surface. The results are summarized in the figures, and can be described in terms of Poinearé waves and Kelvin waves. When the channel is wide compared with the Rossby radius, the final state involves a current of that width which follows the left-hand boundary (for Northern-Hemisphere rotation) to the position of the initial discontinuity, then crosses the channel and continues downstream along the right-hand wall.

Measurements are presented of both mean and fluctuating velocity components in a turbulent boundary layer subjected to a nearly homogeneous external turbulent shear flow. The Reynolds shear stress in the external shear flow is small compared with the wall shear stress. Its transverse mean velocity gradient λ (≃ 6 s−l) is also small compared with typical gradients based on outer variables (say Uw/δ, where Uwis the value of the linear velocity profile extrapolated to the wall and δ is the boundary-layer thickness), but is of the same order as Ut/δ (Ur is the friction velocity). The influence of both positive and negative transverse velocity gradients on the turbulent wall layer is investigated over a streamwise region where the normal Reynolds stresses in the external flow are approximately equal and constant in the streamwise direction. In this region, the integral length scale of the external flow is of the same order of magnitude as that of the wall layer. Measurements in the boundary layer are also given for an un-sheared external turbulent flow (λ = 0) with a turbulence level Tu of 1.5%, approximately the same as that for h = ± 6 s−1. (Tu, is defined as the ratio of the r.m.s. longitudinal velocity fluctuation to Uw.) The measurements are in good agreement with those available in the literature for a similar free-stream turbulence level and show that the external turbulence level and length scale exert a large influence on the turbulence structure in the boundary layer. The additional effect of the external shear on the mean velocity and turbulent energy budget distributions in the inner region of the boundary layer is found to be small. In the outer region, the ‘wake’ component of the mean velocity defect is lowered by the presence of free-stream turbulence and one extra effect due to the external shear is an increase in the Reynolds shear stress when h is positive and a decrease when h is negative. Another interesting effect due to the shear is the appearance near the edge of the layer of a small but distinct region where the local mean velocity is constant and the Reynolds shear stress is negligible.