The maximum spreading of drops impacting on smooth and rough surfaces is measured from low to high impact velocity for liquids with different surface tensions and viscosities. We demonstrate that dynamic wetting plays an important role in the spreading at low velocity, characterized by the dynamic contact angle at maximum spreading. In the energy balance, we account for the dynamic wettability by introducing the capillary energy at zero impact velocity, which relates to the spreading ratio at zero impact velocity. Correcting the measured spreading ratio by the spreading ratio at zero velocity, we find a correct scaling behaviour for low and high impact velocity and, by interpolation between the two, we find a universal scaling curve. The influence of the liquid as well as the nature and roughness of the surface are taken into account properly by rescaling with the spreading ratio at zero velocity, which, as demonstrated, is equivalent to accounting for the dynamic contact angle.

Low-Reynolds-number effects are observed in the inner region of a fully developed turbulent channel flow, using data obtained either from experiments or by direct numerical simulations. The Reynolds-number influence is observed on the turbulence intensities and to a lesser degree on the average production and dissipation of the turbulent energy. In the near-wall region, the data confirm Wei & Willmarth's (1989) conclusion that the Reynolds stresses do not scale on wall variables. One of the reasons proposed by these authors to account for this behaviour, namely the ‘geometry’ effect or direct interaction between inner regions on opposite walls, was investigated in some detail by introducing temperature at one of the walls, both in experiment and simulation. Although the extent of penetration of thermal excursions into the opposite side of the channel can be significant at low Reynolds numbers, the contribution these excursions make to the Reynolds shear stress and the spanwise vorticity in the opposite wall region is negligible. In the inner region, spectra and co-spectra of the velocity fluctuations u and v change rapidly with the Reynolds number, the variations being mainly confined to low wavenumbers in the u spectrum. These spectra, and the corresponding variances, are discussed in the context of the active/inactive motion concept and the possibility of increased vortex stretching at the wall. A comparison is made between the channel and the boundary layer at low Reynolds numbers.

Measurements of the surface-height profile and the vertical distributions of velocity and total head were made behind a two-dimensional fully submerged hydrofoil moving horizontally at constant speed and angle of attack. These measurements were used to resolve the drag on the foil into two parts: one associated with the turbulent breaking region that is sometimes present on the forward face of the first wave, and the other associated with the remaining non-breaking wavetrain. It was found that at ‘incipient breaking’ the first wave existed in either a breaking or a non-breaking state depending on the starting procedures. It was possible to induce steady breaking when the wave slope was 17° or higher. The wake survey measurements showed that the drag associated with breaking reached more than three times the maximum drag that could theoretically be obtained with non-breaking waves. The drag associated with breaking was found to be proportional to the downslope component of the weight of the breaking region.

A new theoretical formulation is presented for mass transport across the dynamic interface associated with a spherical bubble undergoing volume oscillations. As a consequence of the changing internal pressure of the bubble that accompanies volume oscillations, the concentration of the dissolved gas in the liquid at the interface undergoes large-amplitude oscillations. The convection-diffusion equations governing transport of dissolved gas in the liquid are written in Lagrangian coordinates to account for the moving domain. The Henry's law boundary condition is split into a constant and an oscillating part, yielding the smooth and the oscillatory problems respectively. The solution of the oscillatory problem is valid everywhere in the liquid but differs from zero only in a thin layer of the liquid in the neighbourhood of the bubble surface. The solution to the smooth problem is also valid everywhere in the liquid; it evolves via convection-enhanced diffusion on a slow timescale controlled by the Péclet number, assumed to be large. Both the oscillatory and smooth problems are treated by singular perturbation methods: the oscillatory problem by boundary-layer analysis, and the smooth problem by the method of multiple scales in time. Using this new formulation, expressions are developed for the concentration field outside a bubble undergoing arbitrary nonlinear periodic volume oscillations. In addition, the rate of growth or dissolution of the bubble is determined and compared with available experimental results. Finally, a new technique is described for computing periodically driven nonlinear bubble oscillations that depend on one or more physical parameters. This work extends a large body of previous work on rectified diffusion that has been restricted to the assumptions of infinitesimal bubble oscillations or of threshold conditions, or both. The new formulation represents the first self-consistent, analytical treatment of the depletion layer that accompanies nonlinear oscillating bubbles that grow via rectified diffusion.

When the dimensions of a convective system in a saturated porous medium are sufficiently great, diffusion effects can be neglected except in regions where the gradients of fluid properties are very large. A boundary-layer theory is developed for vertical plane flows in such regions. In special cases, the theory is equivalent to that for laminar incompressible flow in a two-dimensional half-jet, or in a plane jet or round jet, for which similarity solutions are well known.

A number of experiments have been performed using a Hele-Shaw cell immersed in water, with a source of potassium permanganate solution located between the plates. At very small values of the source strength, a flow analogous to that of a plane jet from a slit is obtained. The distance advanced by the jet front, or cap, is measured as a function of time, and the velocity is found to be nearly proportional to the velocity of the fluid on the axis of the steady jet behind the cap, as given by the similarity law of Schlichting and Bickley. At large values of the source strength, a two-dimensional ‘broad jet’ of homogeneous solution descending under gravity is produced; the shape of the flow region can be calculated with little error from potential theory, neglecting the effect of the mixing layers.

A possible example of a mixing layer observed in a geothermal region is examined. The theoretical form of the temperature distribution is calculated numerically, taking into account the large viscosity variation with temperature and also the possibility of a large permeability variation. These effects are found to have less influence upon the solution than might have been expected. Quantitative values obtained for the physical parameters are consistent with other geophysical observations.

Certain aspects of the transport of solid particles by a turbulent airstream are discussed, namely: the conveyance of particles in a horizontal pipe, including those carrying an appreciable electrostatic charge; the mechanism of deposition onto a solid wall; and the behaviour of fine particles in a shear flow, such as that in a round jet.

Rough estimates of the effect of the particles on the gaseous turbulence are made, and a primitive physical explanation is offered of the observed velocitylag and pressure drop associated with the transport of particles in a horizontal pipe, under conditions where the influence of the particles’ weight is significant.

Attention is drawn to the difficult problem of dynamically scaling a two-phase flow, and to the different types of interaction between the phases which can occur in a pipe according to its size, the gas velocity through it, and the physical characteristics of the particles.

The paper is an annotated version of a survey presented to the I.U.T.A.M. Symposium on ‘Flow of fluid-solid mixtures’ held in Cambridge during March 1969.

Receptivity of a high-speed boundary layer on a flat plate to acoustic disturbances is investigated using a combined numerical and asymptotic approach. The leading-edge receptivity problem is discussed with emphasis on physical mechanisms associated with scattering and diffraction of acoustic waves. Analytical solutions provide insight into the interplay of these mechanisms as a function of the angle of incidence of external acoustic waves. The theoretical predictions are in good agreement with the wind-tunnel experimental data of Maslov et al. obtained at free-stream Mach number 6. The leading-edge receptivity model is incorporated into the multiple-modes method to account for the inter-modal exchange downstream from the leading edge. This combined modelling resembles basic features of the direct numerical simulation of Ma & Zhong. A comparative analysis of the leading-edge receptivity and the inter-modal exchange associated with non-parallel effects is presented. The theory allows fast evaluation of the receptivity coefficients and clarifies the physics of the receptivity process. The theoretical results may guide further direct numerical simulations and experimental studies of boundary layer receptivity at supersonic and hypersonic speeds.

The spatial decay and structural evolution of grid-generated turbulence under the effect of buoyancy was studied in a ten-layer, salt-stratified water channel. The various density gradients were chosen such that the initial overturning turbulent scale was slightly smaller than any of the respective buoyancy scales. The observed general evolution of the flow from homogeneous turbulence to a composite of fossil turbulence or quasi-two-dimensional turbulence and internal wavefield is in good agreement with the predictions of Gibson (1980) and the lengthscale model of Stillinger, Helland & Van Atta (1983). The effect of the initial size of the turbulent lengthscale compared with the buoyancy scale on the decay and evolution of the turbulence is investigated and the observed influence on the rate of decay of both longitudinal and vertical velocity fluctuations pointed out by Van Atta, Helland & Itsweire (1984) is shown to be related to the magnitude of the initial internal wavefield at the grid. An attempt is made to remove the wave-component kinetic energy from the vertical-velocity-fluctuation data of Stillinger, Helland & Van Atta (1983) in order to obtain the true decay of the turbulent fluctuations. The evolution of the resulting fluctuations is similar to that of the present large-grid data and several towed-grid experiments. The rate of destruction of the density fluctuations (active-scalar dissipation rate) is estimated from the evolution equation for the potential energy, and the deduced Cox numbers are compared with those obtained from oceanic microstructure measurements. The classical Kolmogorov and Batchelor scalings appear to collapse the velocity and density spectra better than the buoyancy scaling proposed by Gargett, Osborn & Nasmyth (1984). The rise of the velocity spectra at low wavenumbers found by Stillinger, Helland & Van Atta (1983) is shown to be related to internal waves.

The solution is obtained to the problem of the steady one-dimensional flow of an incompressible, viscous, electrically conducting fluid through a circular pipe in the presence of an applied (transverse) uniform magnetic field. A no-slip condition on the velocity is assumed at the non-conducting wall. The solution is exact and thus valid for all values of the Hartmann number. Excellent agreement exists between the present theoretical results and the experimental values obtained by Hartmann & Lazarus (1937) in the low to medium Hartmann number range. The high Hartmann number case is treated by Shercliff (1962) in the following paper.

In this paper, we study the interaction of two co-rotating trailing vortices. It is well-known that vortices of like-sign ultimately merge to form a single vortex, and there has been much work on measuring and predicting the initial conditions for the onset of merger, especially concerning the critical vortex core radius. However, the physical mechanism causing this merger has received little attention. In this work, we directly measure the structure of the antisymmetric vorticity field that causes the co-rotating vortices to be pushed towards each other during merger. We discover that the form of the antisymmetric vorticity comprises two counter-rotating vortex pairs, whose induced velocity field readily pushes the two centroids together. The merging velocity computed from the antisymmetric vorticity field agrees closely with the merging velocity measured directly from the total (original) flow field.

The co-rotating vortex pair evolves through four distinct phases. The initial stage comprises a diffusive growth, which can be either viscous or turbulent. In either case, the number of turns that they rotate around one another until the vortices start to merge increases with Reynolds number (Re). If one observes the streamlines in a rotating reference frame (moving with the vortices), then one finds an inner and outer recirculating region of the flow bounded by a separatrix streamline. When the vortices grow large enough in the first stage, diffusion across the separatrix places vorticity into the outer recirculating region of the flow, and this leads to the generation of the antisymmetric vorticity, causing convective merger. This second (convective) stage corresponds to the motion of the vortex centroids towards each other, and is a process which is almost independent of viscosity. During the late part of this stage, the antisymmetric vorticity is diminished by a symmetrization process, and the final merging into one vorticity structure is achieved by a second diffusive stage. The fourth and ultimate phase is one where the merged vortex core grows by diffusion.

To study the dynamics of fluid mud with a high concentration of cohesive clay particles, we present a theory for a thin sheet of Bingham-plastic fluid flowing slowly on an inclined plane. The physics is discussed on the approximate basis of the lubrication theory. Because of the yield stress, the free surface need not be horizontal when the Bingham fluid is in static equilibrium, nor parallel to the plane bed when in steady flow. We then show that there is a variety of gravity currents that can advance at a constant speed and with the same profile. Experimental confirmation of one type is presented. By solving a nonlinear partial differential equation, transient flows due either to a steady upstream discharge or to the sudden release of a finite fluid mass on another fluid layer are studied. In the first case there is a mud front which ultimately propagates as a constant speed as a steady gravity current. In the second case, when the ambient layer is sufficiently shallow that there is no initial motion, the flow induced by the new fluid can terminate after the disturbance has travelled a finite distance. The extent of the final spread is examined. Disturbances due to an external pressure travelling parallel to the free surface are also examined. It is found in particular that a travelling localized pulse of pressure gradient not only generates a localized mud disturbance which travels along with the forcing pressure, but further leaves behind a permanent footprint.

The influence of large variations of viscosity on convection in a layer of fluid heated from below has been investigated. Solutions for the flow and temperature fields were obtained numerically assuming infinite Prandtl number, free-surface boundary conditions, and two-dimensional motion of fixed horizontal wavelength. The effects of a temperature-dependent and a depth-dependent viscosity were each studied; calculations were also carried out using a temperature- and depth-dependent viscosity model appropriate to the earth's mantle.

Laboratory experiments on barotropic vortices in a rotating fluid revealed that the instability behaviour of cyclonic and anticyclonic vortices is remarkably different. Depending on its initial vorticity distribution, the cyclonic vortex has in a number of experiments been observed to be unstable to wavenumber-2 perturbations, leading to the gradual formation of a stable tripolar vortex structure. This tripole consists of an elongated cyclonic core vortex adjoined by two anticyclonic satellite vortices.

In contrast, the anticyclonic vortex shows a rather explosive instability behaviour, in the sense that it is observed to immediately split up into two dipoles. Under somewhat different circumstances the higher-order mode-3 instability is observed, in which the anticyclonic core has a triangular shape, with three smaller cyclonic satellite vortices at its sides.

A modified version of Rayleigh's instability criterion offers a qualitative explanation for this apparent difference between unstable cyclonic and anticyclonic vortices.

Turbulent spots evolving in a laminar boundary layer on a nominally zero pressure gradient flat plate are investigated. The plate is towed through an 18 m water channel, using a carriage that rides on a continuously replenished oil film giving a vibrationless tow. Turbulent spots are initiated using a solenoid valve that ejects a small amount of fluid through a minute hole on the working surface. A novel visualization technique that utilizes fluorescent dye excited by a sheet of laser light is employed. Some new aspects of the growth and entrainment of turbulent spots, especially with regard to lateral growth, are inferred from the present experiments. To supplement the information on lateral spreading, a surbulent wedge created by placing a roughness element in the laminar boundary layer is also studied both visually and with probe measurements. The present results show that, in addition to entrainment, another mechanism is needed to explain the lateral growth characteristics of a turbulent region in a laminar boundary layer. This mechanism, termed growth by destabilization, appears to be a result of the turbulence destabilizing the unstable laminar boundary layer in its vicinity. To further understand the growth mechanisms, the turbulence in the spot is modulated using drag-reducing additives and salinity stratification.

Many processes in the Earth, such as magma migration, can be described by the flow of a low-viscosity fluid in a viscously deformable, permeable matrix. The purpose of this and a companion paper is to develop a better physical understanding of the equations governing these two-phase flows. This paper presents a series of analytic approximate solutions to the governing equations to show that the equations describe two different modes of matrix deformation. Shear deformation of the matrix is governed by Stokes equation and can lead to porosity-driven convection. Volume changes of the matrix are governed by a nonlinear dispersive wave equation for porosity. Porosity waves exist because the fluid flux is an increasing function of porosity and the matrix can expand or compact in response to variations in the fluid flux. The speed and behaviour of the waves depend on the functional relationship between permeability and porosity. If the partial derivative of the permeability with respect to porosity, ∂kϕ/∂ϕ, is also an increasing function of porosity, then the waves travel faster than the fluid in the pores and can steepen into porosity shocks. The propagation of porosity waves, however, is resisted by the viscous resistance of the matrix to volume changes. Linear analysis shows that viscous stresses cause plane waves to disperse and provide additional pressure gradients that deflect the flow of fluid around obstacles. When viscous resistance is neglected in the nonlinear equations, porosity shock waves form from obstructions in the fluid flux. Using the method of characteristics, we quantify the specific criteria for shocks to develop in one and two dimensions. A companion paper uses numerical schemes to show that in the full equations, viscous resistance to volume changes causes simple shocks to disperse into trains of nonlinear solitary waves.

An accurate method is described for integrating the Navier-Stokes equations numerically for the time-dependent flow past an impulsively started circular cylinder. Results of integrations over the range of Reynolds numbers, based on the diameter of the cylinder, from 5 to ∞ are presented and compared with previous numerical, theoretical and experimental results. In particular, the growth of the length of the separated wake behind the cylinder has been calculated for R = 40, 100 and 200 and is found to be in very good agreement with the results of recent experimental measurements. The calculated pressure distribution over the surface of the cylinder for R = 500 is also found to be in reasonable agreement with experimental measurements for the case R = 560.

For Reynolds numbers up to 100 the equations were integrated until most of the features of the flow showed a close approximation to steady-state conditions. The results obtained are in good agreement with previous calculations of the steady flow past a circular cylinder. For R > 100 the integrations were continued until the implicit method of integration broke down by reason of its failure to converge. A secondary vortex appeared on the surface of the cylinder in the case R = 500, but for higher Reynolds numbers, including the case R = ∞, the procedure broke down before the appearance of a secondary vortex. In all cases the flow was assumed to remain symmetrical.

A flow has been produced having effectively zero skin friction throughout its region of pressure rise, which extended for a distance of 3 ft. No fundamental difficulty was encountered in establishing the flow and it had, moreover, a good margin of stability. The dynamic head in the zero skin friction boundary layer was found to be linear at the wall (i.e. u ∞ y½), as predicted theoretically in the previous paper (Stratford 1959).

The flow appears to achieve any specified pressure rise in the shortest possible distance and with probably the least possible dissipation of energy for a given initial boundary layer. Thus an aerofoil which could utilize it immediately after transition from laminar flow would be expected to have a very low drag. A design pressure distribution (besides having the usual safety margin against stall) should have a slightly more gradual start to the pressure rise than in the present experiment, as small errors close to the discontinuity can cause difficulty.

The control of separated fluid flow by reduced-order models is studied using the two-dimensional incompressible flow over an open square cavity at Reynolds numbers where instabilities are present. Actuation and measurement locations are taken on the upstream and downstream edge of the cavity. A bi-orthogonal projection is introduced to arrive at reduced-order models for the compensated problem. Global modes, proper orthogonal decomposition (POD) modes and balanced modes are used as expansion bases for the model reduction. The open-loop behaviour of the full and the reduced systems is analysed by comparing the respective transfer functions. This analysis shows that global modes are inadequate to sufficiently represent the input–output behaviour whereas POD and balanced modes are capable of properly approximating the exact transfer function. Balanced modes are far more efficient in this process, but POD modes show superior robustness. The performance of the closed-loop system corroborates this finding: while reduced-order models based on POD are able to render the compensated system stable, balanced modes accomplish the same with far fewer degrees of freedom.

Through an extensive and systematic experimental investigation of two geometries of flexible plates in air, it is shown that a properly defined scaled Cauchy number allows collapsing all drag measurements of the reconfiguration number. In the asymptotic regime of large deformation, it is shown that the Vogel exponents that scale the drag with the flow velocity for different geometries of plates can be predicted with a simple dimensional analysis reasoning. These predicted Vogel exponents are in agreement with previously published models of reconfiguration. The mechanisms responsible for reconfiguration, namely area reduction and streamlining, are studied with the help of a simple model for flexible plates based on an empirical drag formulation. The model predicts well the reconfiguration observed in the experiments and shows that for a rectangular plate, the effect of streamlining is prominent at the onset of reconfiguration, but area reduction dominates in the regime of large deformation. Additionally, the model demonstrates for both geometries of plates that the reconfiguration cannot be described by a single value of the Vogel exponent. The Vogel exponent asymptotically approaches constant values for small and for very large scaled Cauchy numbers, but in between both extremes it varies significantly over a large range of scaled Cauchy number.

A numerical proceure is introduced to solve the one-dimensional equations of gasdynamics for a cylindrically or spherically symmetric flow. The method consists of a judicious combination of Glimm's method and operator splitting. The method is applied to the problem of a converging cylindrical shock.