We consider a head-on collision between two solitary waves on the surface of an inviscid homogeneous fluid. A perturbation method which in principle can generate an asymptotic series of all orders, is used to calculate the effects of the collision. We find that the waves emerging from (i.e. long after) the collision preserve their original identities to the third order of accuracy we have calculated. However a collision does leave imprints on the colliding waves with phase shifts and shedding of secondary waves. Each secondary wave group trails behind its primary, a solitary wave. The amplitude of the wave group diminishes in time because of dispersion. We have also calculated the maximum run-up amplitude of two colliding waves. The result checks with existing experiments.

The collisional dynamics of equal-sized water and normal-alkane droplets, in the 150 μm radius range, have been experimentally studied for situations involving 0(1) droplet Weber numbers and head-on to grazing impact parameters. Results show that in the parametric range investigated the behaviour of hydrocarbon droplets is significantly more complex than that of water droplets. For head-on collisions, while permanent coalescence always results for water droplets, the outcome is quite nonmonotonic for the hydrocarbon droplets in that, with increasing droplet Weber number, the collision can result in permanent coalescence, bouncing, permanent coalescence again, and coalescence followed by separation with or without production of satellite droplets. Similar complexities exist for off-centre collisions. Phenomenological explanations are offered for these observations based on the material properties of the fluids, the relative influences of the normal and shearing aspects of the collision, and the nature and extent of energy dissipation due to droplet deformation during collision.

Sound waves can be generated in a colloid by the application of an alternating electric field. In this paper we describe the method for calculating this and the related electro-acoustic phenomenon of electric fields generated by sound waves. As an illustration of the procedure, we obtain formulae for these two effects for a suspension of spherical particles with thin double layers, in a parallel plate geometry.

We propose a diagnostic technique to detect instability waves in a subsonic round jet using a phased microphone array. The detection algorithm is analogous to the beam-forming technique, which is typically used with a far-field microphone array to localize noise sources. By replacing the reference solutions used in the conventional beam-forming with eigenfunctions from linear stability analysis, the amplitudes of instability waves in the axisymmetric and first two azimuthal modes are inferred. Experimental measurements with particle image velocimetry and a database from direct numerical simulation are incorporated to design a conical array that is placed just outside the mixing layer near the nozzle exit. The proposed diagnostic technique is tested in experiments by checking for consistency of the radial decay, streamwise evolution and phase correlation of hydrodynamic pressure. The results demonstrate that in a statistical sense, the pressure field is consistent with instability waves evolving in the turbulent mean flow from the nozzle exit to the end of the potential core, particularly near the most amplified frequency of each azimuthal mode. We apply this technique to study the effects of jet Mach number and temperature ratio on the azimuthal mode balance and evolution of instability waves. We also compare the results from the beam-forming algorithm with the proper orthogonal decomposition and discuss some implications for jet noise.

Detailed understanding of flowing granular materials is severely hampered by the deficiencies of present experimental methods. To help increase the information base, a computer simulation has been developed to describe two-dimensional unidirectional flows of inelastic fully rough particles. This paper presents the results of a Couette shear-flow simulation. The results include distributions of velocity, density and granular temperature (a measure of the kinetic energy contained in the random particle motions). The effects of density and shear rate on the granular temperature are explored. Shear and normal forces on the solid walls are compared with experimental and theoretical results. The behaviour of the particles in the simulated flow is examined and assessments are made of the collision angle and velocity distributions. The development of a distinct, ‘layered’ microstructure is observed in high-density granular flows.

The drying of a drop of human blood exhibits coupled physical mechanisms, such as Marangoni flow, evaporation and wettability. The final stage of a whole blood drop evaporation reveals regular patterns with a good reproducibility for a healthy person. Other experiments on anaemic and hyperlipidaemic people were performed, and different patterns were revealed. The flow motion inside the blood drop is observed and analysed with the use of a digital camera: the influence of the red blood cells motion is revealed at the drop periphery as well as its consequences on the final stage of drying. The mechanisms which lead to the final pattern of the dried blood drops are presented and explained on the basis of fluid mechanics in conjunction with the principles of haematology. The blood drop evaporation process is evidenced to be driven only by Marangoni flow. The same axisymmetric pattern formation is observed, and can be forecast for different blood drop diameters. The evaporation mass flux can be predicted with a good agreement, assuming only the knowledge of the colloids mass concentration.

The discrete infinite-dimensional modal system describing nonlinear sloshing of an incompressible fluid with irrotational flow partially occupying a tank performing an arbitrary three-dimensional motion is derived in general form. The tank has vertical walls near the free surface and overturning waves are excluded. The derivation is based on the Bateman–Luke variational principle. The free surface motion and velocity potential are expanded in generalized Fourier series. The derived infinite-dimensional modal system couples generalized time-dependent coordinates of free surface elevation and the velocity potential. The procedure is not restricted by any order of smallness. The general multidimensional structure of the equations is approximated to analyse sloshing in a rectangular tank with finite water depth. The amplitude–frequency response is consistent with the fifth-order steady-state solutions by Waterhouse (1994). The theory is validated by new experimental results. It is shown that transients and associated nonlinear beating are important. An initial variation of excitation periods is more important than initial conditions. The theory is invalid when either the water depth is small or water impacts heavily on the tank ceiling. Alternative expressions for hydrodynamic loads are presented. The procedure facilitates simulations of a coupled vehicle–fluid system.

A drop of fluid, initially held spherical by surface tension, will deform when an electric or magnetic field is applied. The deformation will depend on the electric/magnetic properties (permittivity/permeability and conductivity) of the drop and of the surrounding fluid. The full time-dependent low-Reynolds-number problem for the drop deformation is studied by means of a numerical boundary-integral technique. Fluids with arbitrary electrical properties are considered, but the viscosities of the drop and of the surrounding fluid are assumed to be equal.

Two modes of breakup have been observed experimentally: (i) tip-streaming from drops with pointed ends, and (ii) division of the drop into two blobs connected by a thin thread. Pointed ends are predicted by the numerical scheme when the permittivity of the drop is high compared with that of the surrounding fluid. Division into blobs is predicted when the conductivity of the drop is higher than that of the surrounding fluid. Some experiments have been reported in which the drop deformation exhibits hysteresis. This behaviour has not in general been reproduced in the numerical simulations, suggesting that the viscosity ratio of the two fluids can play an important role.

Three-dimensional flows over impulsively translated low-aspect-ratio flat plates are investigated for Reynolds numbers of 300 and 500, with a focus on the unsteady vortex dynamics at post-stall angles of attack. Numerical simulations, validated by an oil tow-tank experiment, are performed to study the influence of aspect ratio, angle of attack and planform geometry on the wake vortices and the resulting forces on the plate. Immediately following the impulsive start, the separated flows create wake vortices that share the same topology for all aspect ratios. At large time, the tip vortices significantly influence the vortex dynamics and the corresponding forces on the wings. Depending on the aspect ratio, angle of attack and Reynolds number, the flow at large time reaches a stable steady state, a periodic cycle or aperiodic shedding. For cases of high angles of attack, an asymmetric wake develops in the spanwise direction at large time. The present results are compared to higher Reynolds number flows. Some non-rectangular planforms are also considered to examine the difference in the wakes and forces. After the impulsive start, the time at which maximum lift occurs is fairly constant for a wide range of flow conditions during the initial transient. Due to the influence of the tip vortices, the three-dimensional dynamics of the wake vortices are found to be quite different from the two-dimensional von Kármán vortex street in terms of stability and shedding frequency.

In this paper we intend to predict the magnitude of the contribution to the Reynolds stress of bursting events: ‘ejections’, ‘sweeps’, ‘inward interactions’ and ‘outward interactions’. We shall do this by making use of the conditional probability distribution of the Reynolds stress − uv, which can be derived by applying the cumulant-discard method to the Gram-Charlier probability distribution of the two variables u and v. The Reynolds-stress fluctuations in openchannel flows over smooth and rough beds are measured by dual-sensor hot-film anemometers, whose signals are conditionally sampled and sorted into the four quadrants of the u, v plane by using a high-speed digital data processing system.

We shall verify that even the third-order conditional probability distribution of the Reynolds stress shows fairly good agreement with the experimental results and that the sequence of events in the bursting process, i.e. ejections, sweeps and interactions, is directly related to the turbulent energy budget in the form of turbulent diffusion. Also, we shall show that the roughness effect is marked in the area from the wall to the middle of the equilibrium region, and that sweeps appear to be more important than ejections as the roughness increases and as the distance from the wall decreases.

It is well known that subgrid models such as Smagorinsky's cannot be used for the spatially growing simulation of the transition to turbulence of flat-plate boundary layers, unless large-amplitude perturbations are introduced at the upstream boundary: they are over-dissipative, and the flow simulated remains laminar. This is also the case for the structure-function model (SF) of Métais & Lesieur (1992). In the present paper we present a sequel to this model, the filtered-structure-function (FSF) model. It consists of removing the large-scale fluctuations of the field before computing its second-order structure function. Analytical arguments confirm the superiority of the FSF model over the SF model for large-eddy simulations of weakly unstable transitional flows. The FSF model is therefore used for the simulation of a quasi-incompressible (M∞ = 0.5) boundary layer developing spatially over an adiabatic flat plate, with a low level of upstream forcing. With the minimal resolution 650 × 32 × 20 grid points covering a range of streamwise Reynolds numbers Rex1 ε [3.4 × 105, 1.1 × 106], transition is obtained for 80 hours of time-processing on a CRAY 2 (whereas DNS of the whole transition takes about ten times longer). Statistics of the LES are found to be in acceptable agreement with experiments and empirical laws, in the laminar, transitional and turbulent parts of the domain. The dynamics of low-pressure and high-vorticity distributions is examined during transition, with particular emphasis on the neighbourhood of the critical layer (defined here as the height of the fluid travelling at a speed equal to the phase speed of the incoming Tollmien–Schlichting waves). Evidence is given that a subharmonic-type secondary instability grows, followed by a purely spanwise (i.e. time-independent) mode which yields peak-and-valley splitting and transition to turbulence. In the turbulent region, flow visualizations and local instantaneous profiles are provided. They confirm the presence of low- and high-speed streaks at the wall, weak hairpins stretched by the flow and bursting events. It is found that most of the vorticity is produced in the spanwise direction, at the wall, below the high-speed streaks. Isosurfaces of eddy viscosity confirm that the FSF model does not perturb transition much, and acts mostly in the vicinity of the hairpins.

Results of an experimental investigation of the evolution of a nonlinear wave train on deep water are reported. The initial stage of evolution is found to be characterized by exponential growth of a modulational instability, as was first discovered by Benjamin ' Feir. At later stages of evolution it is found that the instability does not lead to wave-train disintegration or loss of coherence. Instead, the modulation periodically increases and decreases, and the wave train exhibits the Fermi–Pasta–Ulam recurrence phenomenon. Results of an earlier study of nonlinear wave packets by Yuen ' Lake, in which solutions of the nonlinear Schrödinger equation were shown to provide quantitatively correct descriptions of the properties of nonlinear wave packets, are applied to describe the experimentally observed wave-train phenomena. A comparison between the laboratory data and numerical solutions of the nonlinear Schrödinger equation for the long-time evolution of nonlinear wave trains is given.

The sound generated by vortex pairing in a two-dimensional compressible mixing layer is investigated. Direct numerical simulations (DNS) of the Navier–Stokes equations are used to compute both the near-field region and a portion of the acoustic field. The acoustic analogy due to Lilley (1974) is also solved with acoustic sources determined from the near-field data of the DNS. It is shown that several commonly made simplifications to the acoustic sources can lead to erroneous predictions for the acoustic field. Predictions based on the quadrupole form of the source terms derived by Goldstein (1976a, 1984) are in excellent agreement with the acoustic field from the DNS. However, despite the low Mach number of the flow, the acoustic far field generated by the vortex pairings cannot be described by considering compact quadrupole sources. The acoustic sources have the form of modulated wave packets and the acoustic far field is described by a superdirective model (Crighton & Huerre 1990). The presence of flow–acoustic interactions in the computed source terms causes the acoustic field predicted by the acoustic analogy to be very sensitive to small changes in the description of the source.

Measurements of skin friction and mean-velocity profiles have been made in fully developed flows in pipes and channels in the Reynolds number range 1000 < Re < 10000. These measurements, and observations of hot-wire signals, indicate rather remarkable differences between two-dimensional and axially symmetric flows and also make it difficult to give a precise definition of the term ‘fully developed turbulent flow’.

Results of laboratory experiments are presented in which a finite volume of homogeneous fluid was released instantaneously into another fluid of slightly lower density. The experiments were performed in a channel of rectangular cross-section, and the two fluids used were salt water and fresh water. As previously reported, the resulting gravity current, if viscous effects are negligible, passes through two distinct phases: an initial adjustment phase, during which the initial conditions are important, and an eventual self-similar phase, in which the front speed decreases as t−1/3 (where t is the time measured from release). The experiments reported herein were designed to emphasize the inviscid motion. From our observations we argue that the current front moves steadily in the first phase, and that the transition to the inviscid self-similar phase occurs when a disturbance generated at the endwall (or plane of symmetry) overtakes the front. If the initial depth of the heavy fluid is equal to or slightly less than the total depth of the fluid in the channel, the disturbance has the appearance of an internal hydraulic drop. Otherwise, the disturbance is a long wave of depression. Measurements of the duration of the initial phase and of the speed and depth of the front during this phase are presented as functions of the ratio of the initial heavy fluid depth to the total fluid depth. These measurements are compared with numerical solutions of the shallow-water equations for a two-layer fluid.

The paper presents the results of measurements of fluctuating lift and drag on a long square cylinder. The measurements include the correlation of lift along the cylinder and the distribution of fluctuating pressure on a cross-section. The magnitude of the fluctuating lift was found to be considerably greater than that for a circular cross-section and the spanwise correlation much stronger.

It was found that the presence of large-scale turbulence in the stream had a marked influence on both the steady and the fluctuating forces. The most significant changes were at small angles of attack (%alpha; < 10°) and included a reduction in base suction and a decrease in fluctuating lift of about 50%.

A numerical model for two-dimensional flow around an airfoil undergoing prescribed heaving motions in a viscous flow is presented. The model is used to examine the flow characteristics and power coefficients of a symmetric airfoil heaving sinusoidally over a range of frequencies and amplitudes. Both periodic and aperiodic solutions are found. Additionally, some flows are asymmetric in that the upstroke is not a mirror image of the downstroke. For a given Strouhal number – defined as the product of dimensionless frequency and heave amplitude – the maximum efficiency occurs at an intermediate heaving frequency. This is in contrast to ideal flow models, in which efficiency increases monotonically as frequency decreases. In accordance with Wang (2000), the separation of the leading-edge vortices at low heaving frequencies leads to diminished thrust and efficiency. At high frequencies, the efficiency decreases similarly to inviscid theory. Interactions between leading- and trailing-edge vortices are categorized, and the effects of this interaction on efficiency are discussed. Additionally, the efficiency is related to the proximity of the heaving frequency to the frequency of the most spatially unstable mode of the average velocity profile of the wake; the greatest efficiency occurs when the two frequencies are nearly identical. The importance of viscous effects for low-Reynolds-number flapping flight is discussed.

In Part 1 by Stuart (1960), a study was made of the growth of an unstable infinitesimal disturbance, or the decay of a finite disturbance through a stable infinitesimal disturbance to zero, in plane Poiseuille flow, and that paper gave the most important terms in a solution of the equations of motion. The greater part of the present paper is concerned with a re-formulation of this problem which readily yields the complete solution. By the same method a solution for Couette flow is obtained. This solution is only a formal one for the present because the conditions imposed in deriving the solution may not be valid for Couette flow; this flow is believed to be stable to infinitesimal disturbances of the type considered.

Direct simulations of the turbulent shear layer are performed for subsonic to supersonic Mach numbers. Fully developed turbulence is achieved with profiles of mean velocity and turbulence intensities that agree well with laboratory experiments. The thickness growth rate of the shear layer exhibits a large reduction with increasing values of the convective Mach number, Mc. In agreement with previous investigations, it is found that the normalized pressure–strain term decreases with increasing Mc, which leads to inhibited energy transfer from the streamwise to cross-stream fluctuations, to the reduced turbulence production observed in DNS, and, finally, to reduced turbulence levels as well as reduced growth rate of the shear layer. An analysis, based on the wave equation for pressure, with supporting DNS is performed with the result that the pressure–strain term decreases monotonically with increasing Mach number. The gradient Mach number, which is the ratio of the acoustic time scale to the flow distortion time scale, is shown explicitly by the analysis to be the key quantity that determines the reduction of the pressure–strain term in compressible shear flows. The physical explanation is that the finite speed of sound in compressible flow introduces a finite time delay in the transmission of pressure signals from one point to an adjacent point and the resultant increase in decorrelation leads to a reduction in the pressure–strain correlation.

The dependence of turbulence intensities on the convective Mach number is investigated. It is found that all components decrease with increasing Mc and so does the shear stress.

DNS is also used to study the effect of different free-stream densities parameterized by the density ratio, s = ρ2/ρ1, in the high-speed case. It is found that changes in the temporal growth rate of the vorticity thickness are smaller than the changes observed in momentum thickness growth rate. The momentum thickness growth rate decreases substantially with increasing departure from the reference case, s = 1. The peak value of the shear stress, uv, shows only small changes as a function of s. The dividing streamline of the shear layer is observed to move into the low-density stream. An analysis is performed to explain this shift and the consequent reduction in momentum thickness growth rate.

The separation of a two-dimensional laminar boundary layer under the influence of a suddenly imposed external adverse pressure gradient was studied by time-accurate numerical solutions of the Navier–Stokes equations. It was found that a strong adverse pressure gradient created periodic vortex shedding from the separation. The general features of the time-averaged results were similar to experimental results for laminar separation bubbles. Comparisons were made with the ‘steady’ separation experiments of Gaster (1966). It was found that his ‘bursting’ occurs under the same conditions as our periodic shedding, suggesting that bursting is actually periodic shedding which has been time-averaged. The Strouhal number based on the shedding frequency, local free-stream velocity, and boundary-layer momentum thickness at separation was independent of the Reynolds number and the pressure gradient. A criterion for onset of shedding was established. The shedding frequency was the same as that predicted for the most amplified linear inviscid instability of the separated shear layer.