We report experiments on the shape and motion of millimetre-sized drops sliding down a plane in a situation of partial wetting. When the Bond number based on the component of gravity parallel to the plane $\Bo_{\alpha}$ exceeds a threshold, the drops start moving at a steady velocity which increases linearly with $\Bo_{\alpha}$. When this velocity is increased by tilting the plate, the drops change their aspect ratio: they become longer and thinner, but maintain a constant, millimetre-scale height. As their aspect ratio changes, a threshold is reached at which the drops are no longer rounded but develop a ‘corner’ at their rear: the contact line breaks into two straight segments meeting at a singular point or at least in a region of high contact line curvature. This structure then evolves such that the velocity normal to the contact line remains equal to the critical value at which the corner appears, i.e. to a maximal speed of dewetting. At even higher velocities new shape changes occur in which the corner changes into a ‘cusp’, and later a tail breaks into smaller drops (pearling transition). Accurate visualizations show four main results. (i) The corner appears when a critical non-zero value of the receding contact angle is reached. (ii) The interface then has a conical structure in the corner regime, the in-plane and out-of-plane angles obeying a simple relationship dictated by a lubrication analysis. (iii) The corner tip has a finite non-zero radius of curvature at the transition to a corner, and its curvature diverges at a finite capillary number, just before the cusp appears. (iv) The cusp transition occurs when the corner opening in-plane half-angle reaches a critical value of about 45°.

A cloud of solute injected into a pipe or channel is known to spread out by a dispersion process based on cross-sectional diffusion across a velocity shear. The original description of the process is due to Taylor (1953, 1954), and an important subsequent contribution was by Aris (1956), who framed and partially solved equations for the integral moments of the cloud of contaminant. The present work resolves some technical difficulties that occur when Aris’ solution method (separation of variables) is pursued in depth. In particular, it is shown that Aris’ technique has to be modified to give the moments at short and moderate times after the injection of solute into the flow. The paper is concerned with dispersion in those parallel flows for which an associated eigenvalue problem has a discrete spectrum of eigenvalues; fortunately, this case appears to be the rule rather than the exception. Expressions are obtained for the second and third moment about the mean, and the theory is applied to three cases of interest.

The rough, hairy surfaces of many insects and spiders serve to render them water-repellent; consequently, when submerged, many are able to survive by virtue of a thin air layer trapped along their exteriors. The diffusion of dissolved oxygen from the ambient water may allow this layer to function as a respiratory bubble or ‘plastron’, and so enable certain species to remain underwater indefinitely. Maintenance of the plastron requires that the curvature pressure balance the pressure difference between the plastron and ambient. Moreover, viable plastrons must be of sufficient area to accommodate the interfacial exchange of O2 and CO2 necessary to meet metabolic demands. By coupling the bubble mechanics, surface and gas-phase chemistry, we enumerate criteria for plastron viability and thereby deduce the range of environmental conditions and dive depths over which plastron breathers can survive. The influence of an external flow on plastron breathing is also examined. Dynamic pressure may become significant for respiration in fast-flowing, shallow and well-aerated streams. Moreover, flow effects are generally significant because they sharpen chemical gradients and so enhance mass transfer across the plastron interface. Modelling this process provides a rationale for the ventilation movements documented in the biology literature, whereby arthropods enhance plastron respiration by flapping their limbs or antennae. Biomimetic implications of our results are discussed.

A mathematical formulation is presented for the transient development of the fluid-flow field and the temperature field in a liquid pool, generated by a spatially variable heat flux falling on an initially solid metal block. This physical situation is an idealized representation of a TIG (tungsten-inert-gas) welding process. In the formulation allowance is made for electromagnetic, buoyancy and surface forces and the resultant equations are solved numerically.

It is found that both the convective flow field and the temperature field are markedly affected by the nature of the heat flux and the flux of electric current falling on the free surface.

In the absence of surface-tension effects a broadly distributed heat flux and corresponding current distribution cause a situation where both electromagnetic and buoyancy forces are important in determining the fluid-flow field; however, in these systems the fluid-flow field does not play a significant role in defining the heat-transfer process. In contrast, a sharply focused heat flux and current density on the free surface give rise to strong electromagnetically driven flows, which play an important role in determining the shape of the weld pool.

Calculations are also done exploring the effect of surface-tension-driven flows. It is found that surface-tension gradients may produce quite high surface velocities and can have a profound effect on determining the weld-pool shape.

Turbulent mixing and combustion are investigated in a gaseous shear layer formed between two streams: one containing a low concentration of hydrogen in nitrogen and the other containing a low concentration of fluorine in nitrogen. The resulting temperature field is measured simultaneously at eight points across the width of the layer using fast-response cold-wire thermometry. The results show the presence of large, hot structures separated by tongues of cool fluid that enter the layer from either side. The usual bell-shaped mean-temperature profiles therefore result from a duty cycle whereby a fixed probe sees alternating hot and cool fluid, which results in the local mean. The adiabatic flame temperature is not achieved in the mean, at any location across the layer. For fixed velocities, it is found that, in general, two different mean-temperature profiles result from a given pair of reactant compositions if the sides of the layer on which they are carried are exchanged (‘flipped’). This finding is a direct consequence of the asymmetric entrainment of fluid into the layer. Results are compared with the predictions of Konrad and discussed in the context of the Broadwell–Breidenthal model. By comparison with the liquid result of Breidenthal, the amount of product formed in the layer at high Reynolds number is found to be dependent upon the Schmidt number. Results for a helium–nitrogen layer are discussed briefly.

The stability of three axisymmetric jet profiles is reviewed. These profiles represent the development of an incompressible jet from a nearly top-hat profile to a fully developed jet profile. The disturbance equations for arbitrary mode number in a region of zero shear, which provide the boundary conditions for the numerical solution, are solved analytically through use of the disturbance vorticity equations. Numerical solutions for the spatial stability for the axisymmetric (n = 0) disturbance and the asymmetric n = l disturbance are presented. Previously published calculations of least stable modes are shown to be incorrectly interpreted and their actual mode types are given. The critical Reynolds number is found to increase as the profile varies from a top-hat to a fully developed jet form. Closed contours of constant amplification, which are unusual in free shear flows, are shown to exist for the n = 1 disturbance in the fully developed jet region. A fluctuation energy balance is used to justify the occurrence of this destabilizing effect of decreasing Reynolds number.

A technique for the quantitative measurement of fluid velocities in the range 0–5 cm/sec is described. The technique uses a pH indicator, is applicable in aqueous solutions and permits visualization and measurement of three-dimensional flow fields.

We present a study of the rheological and optical behaviour of Kramers bead–rod chains in dilute solution using stochastic computer simulations. We consider two model linear flows, steady shear and uniaxial extensional flow, in which we calculate the non-Newtonian Brownian and viscous stress contribution of the polymers, their birefringence and a stress-optic coefficient. We have developed a computer algorithm to differentiate the Brownian from the viscous stress contributions which also avoids the order (δt)−1/2 noise associated with the Brownian forces. The strain or shear rate is made dimensionless with a molecular relaxation time determined by simulated relaxation of the birefringence and stress after a strong flow is applied. The characteristic long relaxation time obtained from the birefringence and stress are equivalent and shown to scale with N2 where N is the number of beads in the chain. For small shear or extension rates the viscous contribution to the effective viscosity is constant and scales as N. We obtain an analytic expression which explains the scaling and magnitude of this viscous contribution. In uniaxial extensional flow we find an increase in the extensional viscosity with the dimensionless flow strength up to a plateau value. Moreover, the Brownian stress also reaches a plateau and we develop an analytic expression which shows that the Brownian stress in an aligned bead–rod chain scales as N3. Using scaling arguments we show that the Brownian stress dominates in steady uniaxial extensional flow until large Wi, Wi ≈ 0.06N2, where Wi is the chain Weissenberg number. In shear flow the viscosity decays as Wi−1/2 and the first normal stress as Wi−4/3 at moderate Wi. We demonstrate that these scalings can be understood through a quasi-steady balance of shear forces with Brownian forces. For small Wi (in shear and uniaxial extensional flow) and for long times (in stress relaxation) the stress-optic law is found to be valid. We show that the rheology of the bead–rod chain can be qualitatively described by an equivalent FENE dumbbell for small Wi.

Vorticity is deposited baroclinically by shock waves on density inhomogeneities. In two dimenslons, the circulation deposited on a planar interface may be derived analytically using shock polar analysis provided the shock refraction is regular. We present analytical expressions for Γ′, the circulation deposited per unit length of the unshocked planar interface, within and beyond the regular refraction regime. To lowest order, Γ′ scales as \[ \Gamma^\prime\propto (1-\eta^{-\frac{1}{2}})(\sin\alpha)(1+M^{-1}+2M^{-2})(M-1)(\gamma^{\frac{1}{2}}/\gamma + 1)\] where M is the Mach number of the incident shock, η is the density ratio of the gases across the interface, α is the angle between the shock and the interface and γ is the ratio of specific heats for both gases. For α ≤ 30°, the error in this approximation is less than 10% for 1.0 < M ≤ 1.32 for all η > 1, and 5.8 ≤ η ≤ 32.6 for all M. We validate our results by quantification of direct numerical simulations of the compressible Euler equations with a second-order Godunov code.

We generalize the results for total circulation on non-planar (sinusoidal and circular) interfaces. For the circular bubble case, we introduce a ‘near-normality’ ansatz and obtain a model for total circulation on the bubble surface that agrees well with results of direct numerical simulations. A comparison with other models in the literature is presented.

The deformation of a viscous drop, driven by buoyancy towards a solid surface or a deformable interface, is analysed in the asymptotic limit of small Bond number, for which the deformation becomes important only when the drop is close to the solid surface or interface. Lubrication theory is used to describe the flow in the thin gap between the drop and the solid surface or interface, and boundary-integral theory is used in the fluid phases on either side of the gap. The evolution of the drop shape is traced from a relatively undeformed state until a dimple is formed and a long-time quasi-steady-state pattern is established. A wide range of drop to suspending phase viscosity ratios is examined. It is shown that a dimple is always formed, independently of the viscosity ratio, and that the long-time thinning rates take simple forms as inverse fractional powers of time.

The three-dimensional structure of the flow behind a heaving and pitching finite-span wing is investigated using dye flow visualization at a Reynolds number of 164. Phase-locked image sequences, which are obtained from two orthogonal views, are combined to create a set of composite images that give an overall sense of the three-dimensional structure of the flow. A model of the vortex system behind the wing is constructed from the image sequences. Variations of the Strouhal number, pitch amplitude and heave/pitch phase angle are qualitatively shown to affect the structure of the wake.

The classical lubrication-theory model of steady peristaltic transport of periodic sinusoidal waves in infinite-length tubes (Shapiro et al. 1969) is generalized to arbitrary wave shape and wavenumber in tubes of finite length. Whereas the classical model is steady in a frame of reference moving with the peristaltic waves, peristaltic transport in a finite-length tube is inherently non-steady. It may be shown, however, that pumping performance is independent of tube length if there exists an integral number of peristaltic waves in the tube. Three particularly interesting characteristics of non-steady peristalsis are described: (i) fluctuations in pressure and shear stress arise due to a non-integral number of waves in the finite-length tube; (ii) retrograde motion of fluid particles during peristaltic transport (reflux) has inherently different behaviour with single peristaltic waves as compared to multiple ‘train waves’, and (iii) finite tube length, the number of peristaltic waves and the degree of tube occlusion affect global pumping performance. We find that, whereas significant increases in pressure and shear stress result from the tube-to-wave length ratio being non-integral, global pumping performance is only slightly degraded by the existence of a non-integral number of waves in the tube during peristaltic transport. Furthermore, the extent of retrograde motion of fluid particles is much greater with single waves than with train waves. These results suggest that in the design and analysis of peristaltic pumps attention should be paid to the unsteady effects of finite tube length and to the differences between single and multiple peristaltic waves.

The impact of a drop onto a liquid film with a relatively high impact velocity, leading to the formation of a crown-like ejection, is studied theoretically. The motion of a kinematic discontinuity in the liquid film on the wall due to the drop impact, the formation of the upward jet at this kinematic discontinuity and its elevation are analysed. Four main regions of the drop and film are considered: the perturbed liquid film on the wall inside the crown, the unperturbed liquid film on the wall outside the crown, the upward jet forming a crown, and the free rim bounding this jet. The theory of Yarin & Weiss (1995) for the propagation of the kinematic discontinuity is generalized here for the case of arbitrary velocity vectors in the inner and outer liquid films on the wall. Next, the mass, momentum balance and Bernoulli equations at the base of the crown are considered in order to obtain the velocity and the thickness of the jet on the wall. Furthermore, the dynamic equations of motion of the crown are developed in the Lagrangian form. An analytical solution for the crown shape is obtained in the asymptotic case of such high impact velocities that the surface tension and the viscosity effects can be neglected in comparison to inertial effects. The edge of the crown is described by the motion of a rim, formed due to the surface tension.

Three different cases of impact are considered: normal axisymmetric impact of a single drop, oblique impact of a single drop, and impact and interaction of two drops. The theoretical predictions of the height of the crown in the axisymmetric case are compared with experiments. The agreement is quite good in spite of the fact that no adjustable parameters are used.

The effect of applied magnetic fields on the deformation of a biocompatible hydrophobic ferrofluid drop suspended in a viscous medium is investigated numerically and compared with experimental data. A numerical formulation for the time-dependent simulation of magnetohydrodynamics of two immiscible non-conducting fluids is used with a volume-of-fluid scheme for fully deformable interfaces. Analytical formulae for ellipsoidal drops and near-spheroidal drops are reviewed and developed for code validation. At low magnetic fields, both the experimental and numerical results follow the asymptotic small deformation theory. The value of interfacial tension is deduced from an optimal fit of a numerically simulated shape with the experimentally obtained drop shape, and appears to be a constant for low applied magnetic fields. At high magnetic fields, on the other hand, experimental measurements deviate from numerical results if a constant interfacial tension is implemented. The difference can be represented as a dependence of apparent interfacial tension on the magnetic field. This idea is investigated computationally by varying the interfacial tension as a function of the applied magnetic field and by comparing the drop shapes with experimental data until a perfect match is found. This estimation method provides a consistent correlation for the variation in interfacial tension at high magnetic fields. A conclusion section provides a discussion of physical effects which may influence the microstructure and contribute to the reported observations.

This article is aimed at relating a certain substantial body of established material concerning wave loading on offshore structures to fundamental principles of mechanics of solids and of fluids and to important results by G. I. Taylor (1928a, b). The object is to make some key parts within a rather specialised field accessible to the general fluid-mechanics reader.

The article is concerned primarily to develop the ideas which validate a separation of hydrodynamic loadings into vortex-flow forces and potential-flow forces; and to clarify, as Taylor (1928b) first did, the major role played by components of the potential-flow forces which are of the second order in the amplitude of ambient velocity fluctuations. Recent methods for calculating these forces have proved increasingly important for modes of motion of structures (such as tension-leg platforms) of very low natural frequency.

The time-dependent deformation and burst of a viscous drop in an arbitrary shear flow at zero Reynolds number is studied. The viscosities of the drop and the suspending fluid are assumed to be equal. A numerical scheme to track the (non-axisymmetric) drop shape in time is presented, and used to investigate the deformation induced by two-dimensional shear and orthogonal rheometer flows. Steady deformations, critical flow rates and burst modes are determined, and compared with asymptotic (small and large) deformation theories, and with experiment.

Non-intrusive measurements of the streamwise velocity in turbulent round jets in air are performed by recording short-time displacements and distorsions of very thin tagging lines written spanwise into the flow. The lines are written by Raman-exciting oxygen molecules and are interrogated by laser-induced electronic fluorescence (relief). This gives access to the spatial structure of transverse velocity increments without recourse to the Taylor hypothesis. The resolution is around 25 μm, less than twice the Kolmogorov scale η for the experiments performed (with Rλ≈360–600).

The technique is validated by comparison with results obtained from other techniques for longitudinal or transverse structure functions up to order 8. The agreement is consistent with the estimated errors – a few percent on exponents determined by extended-self-similarity – and indicates significant departures from Kolmogorov (1941) scaling.

Probability distribution functions of transverse velocity increments Δu over separations down to 1:8η are reported for the first time. Violent events, with Δu comparable to the r.m.s. turbulent velocity fluctuation, are found to take place with statistically significant probabilities. The shapes of the corresponding lines suggest the effect of intense slender vortex filaments.

The location of the asymptotic virtual origin of positively buoyant turbulent plumes with a deficit of initial momentum flux when compared with equivalent pure plumes is investigated. These lazy plumes are generated by continuous steady releases of momentum, buoyancy and volume into a quiescent uniform environment from horizontal sources (at z = 0) of finite area, and are shown to be equivalent to the far-field flow above point source pure plumes, of buoyancy only, rising from the asymptotic virtual source located below the actual source at z = −zavs.

An analytical expression for the location of the asymptotic virtual source relative to the actual source of the lazy plume is developed. The plume conservation equations are solved for the volume flow rate, and the position of the asymptotic virtual origin is deduced from the scaling for the volume flow rate at large distances from the source.

The displacement zavs of the asymptotic virtual origin from the actual origin scales on the source diameter and is a function of the source parameter Γ ∝ Qˆ20Fˆ0/Mˆ5/20 which is a measure of the relative importance of the initial fluxes of buoyancy Fˆ0, momentum Mˆ0, and volume Qˆ0 in the plume. The virtual origin correction developed is valid for Γ > 1/2 and is therefore applicable to lazy plumes for which Γ > 1, pure plumes for which Γ = 1, and forced plumes in the range 1/2 < Γ < 1. The dimensionless correction z*avs decreases as Γ increases, and for Γ [Gt ] 1, z*avs → 0.853Γ−1/5. Comparisons made between the predicted location of the asymptotic virtual origin and the location inferred from measurements of lazy saline plumes in the laboratory show close agreement.

General dimensional and similarity arguments are applied to derive a heat and mass transfer law for fully turbulent flow along a rough wall. The derivation is quite analogous to Millikan's (1939) derivation of a skin-friction law for smooth-and rough-wall flows and to the derivation of the heat and mass transfer law for smooth-wall flows by Fortier (1968a, b) and Kader & Yaglom (1970, 1972).

The equations derived for the heat or mass transfer coefficient (Stanton number) Ch and Nusselt number Nu include the constant term β of the logarithmic equation for the mean temperature or concentration of a diffusing substance. This term is a function of the Prandtl number, the dimensionless height of wall protrusions and of the parameters describing the shapes and spatial distribution of the protrusions. The general form of the function β is roughly estimated by a simplified analysis of the eddy-diffusivity behaviour in the proximity of the wall (in the gaps between the wall protrusions). Approximate values of the numerical coefficients of the equation for β are found from measurements of the mean velocity and temperature (or concentration) above rough walls. The equation agrees satisfactorily with all the available experimental data. It is noted that the results obtained indicate that roughness affects heat and mass transfer in two ways: it produces the additional disturbances augmenting the heat and mass transfer and simultaneously retards the fluid flow in the proximity of the wall. This second effect leads in some cases to deterioration of heat and mass transfer from a rough wall as compared with the case of a smooth wall at the same values of the Reynolds and Prandtl numbers.

Water wave diffraction by an array of bottom-mounted circular cylinders is analysed under the assumptions of linear theory. The cylinders are identical, and equally spaced along the array. When the number of cylinders is large, but finite, near-resonant modes occur between adjacent cylinders at critical wavenumbers, and cause unusually large loads on each element of the array. These modes are associated with the existence of homogeneous solutions for the diffraction by an array which extends to infinity in both directions. This phenomenon is related to the existence of trapped waves in a channel. A second trapped wave is established, corresponding to Dirichlet boundary conditions on the channel walls, as well as a sequence of higher wavenumbers where ‘nearly trapped’ modes exist.