A potential-flow based analytical model is developed for the temporal development of the spectra of ripples and dunes generated on an initially flattened bed by an open-channel flow. The analysis predicts the occurrence of two distinct spectral peaks at small times; the physical origin of each is explained. One of the peaks has spatial frequency equal to that given by the Airy relation for a small amplitude stationary free-surface wave. These results are confirmed by Jain's (1971) experimental data. Jain's data also are used to determine quantitativevalues of the lag distance andHayashi's (1970) inclination factor, α, which appear in the sediment-transport relation used in the analysis. A variance-cascade model is used to derive the ‘minus-three-power law’ that describes the spatial spectra of rippIes and dunes at higher wavenumbers. Finally, the relation between the Froude number and dominant wavelength implied by Jain & Kennedy's (1971) non-dimensionalized spectra is discussed and compared with data presented by Jain (1971) and Nordin (1971).

Forced nonlinear acoustic oscillations near the resonant frequency of closed and open tubes are studied experimentally. In particular, the motion in tubes terminated with different orifice plates is studied, and comparison is made with second- and third-order theories of the motion which contain an adjustable end-wall reflexion coefficient.

It is found that oscillations at resonance in an open tube exhibit remarkably large amplitudes despite the fact that in some cases shock waves are emitted from the open end. For oscillations at resonance in a closed tube, the effect of substituting an orifice plate for the solid end wall is to reduce the amplitude and thicken the compressive portion of the shock waves which occur under these conditions. In both the open-tube and closed-tube experiments the reflexion coefficients which are evaluated by fitting theory to experiment are found to increase with increasing amplitude, in agreement with the observations of previous investigators (Ingard & Ising 1967). In fact, for the open end the same linear dependence upon amplitude is observed, but the constant of proportionality is different. Qualitative differences are observed between the reflexion coefficients of a given orifice at the open-end and the closed-end resonant frequencies; at the open-end frequency the reflexion from the given orifice is less ideal than at the closed-end frequency. The implications of reflexion coefficients dependent on the wave forms are discussed.

The form taken by a finite amplitude internal seiche in a shallow two-layer fluid within a long narrow container is akin to a shock wave in a gas-filled tube; an undular bore or surge is present which reflects back and forth along the length of the container. The internal surge in Loch Ness is sometimes observed to retain its amplitude over several seiche periods, and a resonance with the wind appears possible. This idea is explored by developing a theoretical model and by making a laboratory experiment, but difficulties are encountered in estimating the size of the parameters which characterize the natural phenomenon.

The perturbation of a turbulent flow by an organized wavelike disturbance is examined using a dynamical, rather than phenomenological, approach. On the basis of the assumption that an infinitesimal perturbation results in a linear change in the statistics of the turbulence, and that the turbulence is either weak or that the turbulent perturbations are quasi-Gaussian, a method for predicting the perturbation turbulent Reynolds stresses is developed. The novel aspect of the analysis is that all averaging is delayed until the dynamical equations have been solved rather than attempting to find, apriori, equations for averaged quantities. When applied to long-wave perturbations the analysis indicates that the perturbation shear stress is of primary dynamical importance, and that this stress is determined by the principal component of mean shear through a relation which depends on the spectrum of the turbulent velocity component parallel to the gradient of the undisturbed mean velocity (the component perpendicular to the wall in a turbulent boundary layer). Theoretical arguments and observations are used to estimate the form of this spectrum in a constant stress shear layer. This results in a prediction of the constitutive law relating turbulent stress and the mean flow. The law is visco-elastic in nature, and is in agreement with the known constitutive relation for stress perturbations to a constant stress boundary layer; it resembles the eddy viscosity relation used successfully by others in describing perturbations in turbulent flows. The details of the constitutive law depend on how well the turbulence obeys Taylor's hypothesis that phase velocity equals mean flow velocity, and some insight into this question is given.

Measurements of the mean flow, intermittent structure and turbulent fluctuations were made in a cold-wall boundary layer at a stream Mach number of 9·4 and Reynolds number based on momentum thickness of 36 800. For these conditions, the r.m.s. sublayer thickness was 32 times smaller than that of the boundary layer proper, and the interfacial standard deviation of the latter was about three times proportionately smaller than has been found at low speeds. The mean flow data, which extended well into the sublayer, revealed a large increase in static pressure from the layer edge to the wall and a quadratic law relation between the total temperature and velocity. While the transformed velocity profile was in good agreement with the incompressible law of the wake, no indication of a linear variation of velocity in the sublayer was detected.

Hot-wire fluctuation data, interpreted with the use of appropriate assumptions concerning the nature of the sound field, indicated that the turbulence is dominated by high-frequency pressure fluctuations whose magnitude at the wall and beyond the layer edge agree with extrapolation of data acquired at supersonic speeds. The static temperature fluctuations agreed with expectations from adiabatic, supersonicdata apparently because they were suppressed by the cooled-wall condition. The fluctuations in the longitudinal velocity component were generally small and differed little from lower Mach number results. The high turbulence Reynolds numbers found generated an inertial-subrange spectral decay, while the longitudinal integral scales were found independent of turbulence mode and about one-fifth the boundary-layer thickness.

The equations governing the nonlinear development of a centred three-dimensional disturbance to plane parallel flow at slightly supercritical Reynolds numbers are obtained, In contrast to the corresponding equation for two-dimensional disturbances, two slowly varying functions are needed to describe the development: the amplitude function and a function related to the secular pressure gradient produced by the disturbance. These two functions satisfy a pair of coupled partial differential equations. The equations derived in Hocking, Stewartson & Stuart (1972) are shown to be incorrect, Some of the properties of the governing equations are discussed briefly.

The sharp interface which exists in unbounded turbulent shear flows was studied using a linear array of twenty hot-wire probes. The probe arrangement was such that the location of the interface could be monitored instantaneously. The particular flow examined was a two-dimensional turbulent wall jet. It was found that the interface is a highly contorted surface which exhibits a significant amount of folding. Quantitative methods for characterizing this behaviour are presented, together with pertinent measurements. In addition, measurements of the mean surface area of the interface, and space-time corrrelations of the width of the turbulence were obtained. The latter were used to find characteristic scales and convection velocities of the interface.

The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and non-axisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.

Axisymmetric solutions become unstable to non-axisymmetric perturbations at values of the Rayleigh number that depend on Prandtl number and shell thickness. However, even when stable, axisymmetric solutions are not a preferred solution to the full equations; steady non-axisymmetric solutions are obtained for the same parameter values. Initial conditions determine the characteristics of the finite-amplitude solutions, including, in the cases of non-axisymmetry, whether or not a steady state is achieved. Transitions in horizontal flow structure can occur, accompanied by a transition in functional dependence of heat flux on Rayleigh number. The dominant modes in the solutions are usually the modes most unstable to the onset of convection, but not always.

A theoretical analysis has been made of turbulent viscous interaction on iso-thermal surfaces at hypersonic speeds. The important parameters governing the effects of incidence and displacement have been obtained under both strong and weak interaction conditions for flat-plate flows. A more general expression relating boundary-layer growth to the external pressure field and effective body shape has been obtained. The method is applied to the wedge compression corner problem and the results compared with some experimental data.

In this paper we assume the existence of a nonlinear boundary layer centred on the critical point, and explore its effect on the development of unstable parallel shear flows. A velocity matching condition derived in a qualitative discussion suggests a growth of harmonics which differs from that predicted by previous theories; however, the prediction is in excellent agreement with experimental data. A hyperbolic-tangent velocity profile, subjected to perturbations with wavenumbers and frequencies close to marginal values, is then chosen as a mathematical model of the nonlinear development, both temporal and spatial instability growth being considered.

A singularity in the analysis which has been treated in previous theories by the introduction of viscosity is dealt with in the present work by the introduction of a growth boundary layer. The asymptotics are non-uniform and the time-dependent solution does not resemble the steady viscous solutions, even as the growth rate tends to zero. The theory suggests that the instability will develop as a series of temporally growing spiral vortices, a description differing from that of a cat's-eye pattern predicted by existing theories, but in accord with experimental and field observations.

Detailed measurements of pressure distributions, mean velocity profiles and Reynolds stresses were made in the thick axisymmetric turbulent boundary layer near the tail of a body of revolution. The results indicate a number of important differences between the behaviour of a thick and a thin boundary layer. The thick boundary layer is characterized by significant variations in static pressure across it and an abnormally low level of turbulence. The static-pressure variation is associated with a strong interaction between the boundary layer and the potential flow outside it, while the changes in the turbulence structure appear to be a consequence of the transverse surface curvature. In order to predict the behaviour of the flow in the tail region of a body of revolution it is not therefore possible to use conventional thin-boundary-layer calculation procedures.

The steady flow of an ideal gas past a conical body is investigated by the method of matched asymptotic expansions, with particular emphasis on the flow near the singular ray occurring in linearized theory. The first-order problem governing the flow in this region is formulated, leading to the equation of Kuo, and an approximate solution is obtained in the case of compressive flow behind the main front. This solution is compared with the results of previous investigations with a view to assessing the applicability of the Lighthill-Whitham theories.

Laplace transforms and regular double asymptotic expansions are used to solve the problem of ordinary chemical mass transport in a permeable tube, where there is small radial convection through the membrane wall and where the length-to-diameter ratio is large. The system is taken to be dilute and Newtonian and the solution is found to higher order in two small parameters. Results indicate that the exit concentration decreases markedly as the diameter, membrane permeability and tube length increase, and that changes in mass transport owing to variations in radial convection are much more significant than those due to the same order of magnitude changes in the resistance of the chemical solute to passage through the membrane (transmittance). In addition, the maximum effects of changes in the radial convection and transmittance are not at the membrane itself (r = 1), but rather roughly at radial values of 0·6 and 0, respectively.

The inviscid stability of swirling flows with mean velocity profiles similar to that obtained by Batchelor (1964) for a trailing vortex from an aircraft is studied with respect to infinitesimal non-axisymmetric disturbances. The flow is characterized by a swirl parameter q involving the ratio of the magnitude of the maximum swirl velocity to that of the maximum axial velocity. It is found that, as the swirl is continuously increased from zero, the disturbances die out quickly for a small value of q if n = 1 (n is the azimuthal wavenumber of the Fourier disturbance of type exp{i(αx + nϕ − αct)}); but for negative values of n, the amplification rate increases and then decreases, falling to negative values at q slightly greater than 1·5 for n = −1. The maximum amplification rate increases for increasingly negative n up to n = −6 (the highest mode investigated), and corresponds to q ≃ 0·85. The applicability of these results to attempts at destabilizing vortices is briefly discussed.

Mean velocity and mean turbulent field measurements are performed for the case of a three-dimensional turbulent boundary layer on an axially rotated cylinder. The cylinder model consists of two parts: a stationary section followed by a spinning afterbody. Techniques of hot-wire anemometry are employed, which yield complete mean velocity and turbulence measurements in skewed flows. The general behaviour of the three-dimensional boundary layer is first discussed: two asymptotic layers analogous to the two-dimensional wall and defect layers are observed; they are shown to evolve from the equations of mean motion. The hypothesis of scalar eddy viscosity is investigated in the light of these results. Using conventional length scale assumptions together with the Reynolds stress tensor equations, a prediction of curvature effects in the law of the wall region is developed; a result in the present case is a smaller slope of the semi-logarithmic portion of the law of the wall, No assumptions over and above those necessary for plane, two-dimensional flow are required for this analysis. The geometry of the model is such that a rapid change in mean rate of strain occurs along the streamlines. From the history of the components of the $\overline{u_iu_j}$ tensor, it is possible to draw some fundamental conclusions concerning the dynamics of the energy dissipation, diffusion and redistribution processes.

We consider the two-dimensional problem of a linearly stratified salt solution contained between two infinite vertical plates. The fluid and the plates are initially at the same temperature. At t = 0, one of the plates is given a step increase in temperature, while the other is maintained at the initial temperature. A time-dependent basic flow is thus generated. The stability of such a time-dependent flow is analysed using an initial value problem approach to the linear stability equations. The method consists of initially distributing small random disturbances of given vertical wavelength throughout the fluid. The disturbances may be in the vorticity, temperature or salinity. The linearized field equations are integrated numerically. The growth or decay of the kinetic energy of the perturbation delineates unstable and stable states. Results have been obtained for a wide range of gap widths. The critical wavelength and the critical Rayleigh number compare favourably with those obtained previously in both physical and numerical experiments.

The limitations of existing one-dimensional experiments on double-diffusive convection are discussed, and a variety of new two-dimensional phenomena are described. We have used the sugar-salt system and shadowgraph photography to make exploratory studies of motions which can arise in a fluid with two smooth, opposing, vertical concentration gradients, with and without horizontal gradients. Many different effects have been observed, the most important of which are the following, (a) In the ‘finger’ case, local disturbances can propagate rapidly as wave motions, which cause a simultaneous breakdown to convection over large horizontal distances. (b) Layers formed in the’ diffusive’ sense overturn locally to produce fingers, but propagate more slowly, as convective rather than wave motions, (c) A series of layers, separated by diffusive interfaces, can become unstable, in the sense that successive layers merge in time as their densities become equal, (d) The presence of horizontally separated sources of water of similar density but different T, S characteristics can lead to the development of strong vertical gradients and extensive quasi-horizontal layering.

Most of our results are qualitative, but it is hoped that they will stimulate further quantitive work on each of the new processes described. It is already clear that much more needs to be done before the mechanism of formation of layers observed in the ocean can be regarded as properly understood.

Finite-amplitude disturbances in plane Poiseuille flow are studied by a method involving Fourier expansion with numerical solution of the resulting partial differential equations in the coefficient functions. A number of solutions are developed which extend to relatively long times so that asymptotic stability or instability can be established with a degree of confidence. The amplitude for neutral stability is established for a fixed wavenumber for two values of the Reynolds number. Details of the neutral velocity fluctuation are presented. These and earlier results are expressed in terms of the asymptotic amplitude and compared with results obtained by prior workers. The results indicate that the expansion methods used by prior workers may be valid only for amplitudes considerably smaller than 0·01.

Davis has shown by means of a lengthy calculation that, for two-dimensional oscillations in a canal of width 2a, the mth eigenvalue has the form \[ N_m = {\textstyle\frac{1}{2}}m\pi - \frac{\lambda_1+\lambda_2}{4m\pi}+o\bigg(\frac{1}{m}\bigg), \] where λ1/a and λ2/a are the curvatures of the bounding cross-sectional curve C at its vertical intersections with the free surface. Here the same result is obtained more simply.