Superhydrophobic surfaces generate very high contact angles as a result of their microstructure. The impact of a water drop on such a surface shows unusual features, such as total rebound at low impact speed. We report experimental and numerical investigations of the impact of approximately spherical water drops. The axisymmetric free surface problem, governed by the Navier–Stokes equations, is solved numerically with a front-tracking marker-chain method on a square grid. Experimental observations at moderate velocities and capillary wavelength much less than the initial drop radius show that the drop evolves to a staircase pyramid and eventually to a torus. Our numerical simulations reproduce this effect. The maximal radius obtained in numerical simulations precisely matches the experimental value. However, the large velocity limit has not been reached experimentally or numerically. We discuss several complications that arise at large velocity: swirling motions observed in the cross-section of the toroidal drop and the appearance of a thin film in the centre of the toroidal drop. The numerical results predict the dry-out of this film for sufficiently high Reynolds and Weber numbers. When the drop rebounds, it has a top-heavy shape. In this final stage, the kinetic energy is a small fraction of its initial value.

The effectiveness of a cylindrical perforated liner with mean bias flow in its absorption of planar acoustic waves in a duct is investigated. The liner converts acoustic energy into flow energy through the excitation of vorticity fluctuations at the rims of the liner apertures. A one-dimensional model that embodies this absorption mechanism is developed. It utilizes a homogeneous liner compliance adapted from the Rayleigh conductivity of a single aperture with mean flow. The model is evaluated by comparing with experimental results, with excellent agreement. We show that such a system can absorb a large fraction of incoming energy, and can prevent all of the energy produced by an upstream source in certain frequency ranges from reflecting back. Moreover, the bandwidth of this strong absorption can be increased by appropriate placement of the liner system in the duct. An analysis of the acoustic energy flux is performed, revealing that local differences in fluctuating stagnation enthalpy, distributed over a finite length of duct, are responsible for absorption, and that both liners in a double-liner system are absorbant. A reduction of the model equations in the limit of long wavelength compared to liner length reveals an important parameter grouping, enabling the optimal design of liner systems.

The object of this paper is to derive the added mass and damping coefficients associated with the periodic motions of a floating hemisphere. Two physically distinct cases are considered; namely those of heave and surge, where these nautical terms refer respectively to a vertical or horizontal oscillation of the body. Computations have been done and the values found for the various force coefficients are presented in tabulated form. A brief derivation of the long- and short-wave asymptotics of these coefficients has also been included.

Combined thermocapillary–buoyancy convection in a thin rectangular geometry is investigated experimentally, with an emphasis on the generation of hydrothermal-wave instabilities. For sufficiently thin layers, pure hydrothermal waves are observed, and are found to be oblique as predicted by a previous linear-stability analysis (Smith & Davis 1983). For thicker layers, both a steady multicell state and an oscillatory state are found to exist, but the latter is not in the form of a pure hydrothermal wave.

To investigate the phenomena of skin-friction drag reduction in a turbulent boundary layer (TBL) at large scales and high Reynolds numbers, a set of experiments has been conducted at the US Navy's William B. Morgan Large Cavitation Channel (LCC). Drag reduction was achieved by injecting gas (air) from a line source through the wall of a nearly zero-pressure-gradient TBL that formed on a flat-plate test model that was either hydraulically smooth or fully rough. Two distinct drag-reduction phenomena were investigated; bubble drag reduction (BDR) and air-layer drag reduction (ALDR).

The streamwise distribution of skin-friction drag reduction was monitored with six skin-friction balances at downstream-distance-based Reynolds numbers to 220 million and at test speeds to 20.0ms−1. Near-wall bulk void fraction was measured at twelve streamwise locations with impedance probes, and near-wall (0 < Y < 5mm) bubble populations were estimated with a bubble imaging system. The instrument suite was used to investigate the scaling of BDR and the requirements necessary to achieve ALDR.

Results from the BDR experiments indicate that: significant drag reduction (>25%) is limited to the first few metres downstream of injection; marginal improvement was possible with a porous-plate versus an open-slot injector design; BDR has negligible sensitivity to surface tension; bubble size is independent of surface tension downstream of injection; BDR is insensitive to boundary-layer thickness at the injection location; and no synergetic effect is observed with compound injection. Using these data, previous BDR scaling methods are investigated, but data collapse is observed only with the ‘initial zone’ scaling, which provides little information on downstream persistence of BDR.

ALDR was investigated with a series of experiments that included a slow increase in the volumetric flux of air injected at free-stream speeds to 15.3ms−1. These results indicated that there are three distinct regions associated with drag reduction with air injection: Region I, BDR; Region II, transition between BDR and ALDR; and Region III, ALDR. In addition, once ALDR was established: friction drag reduction in excess of 80% was observed over the entire smooth model for speeds to 15.3ms−1; the critical volumetric flux of air required to achieve ALDR was observed to be approximately proportional to the square of the free-stream speed; slightly higher injection rates were required for ALDR if the surface tension was decreased; stable air layers were formed at free-stream speeds to 12.5ms−1 with the surface fully roughened (though approximately 50% greater volumetric air flux was required); and ALDR was sensitive to the inflow conditions. The sensitivity to the inflow conditions can be mitigated by employing a small faired step (10mm height in the experiment) that helps to create a fixed separation line.

We present experimental observations of the disk of air caught under a drop impacting onto a solid surface. By imaging the impact through an acrylic plate with an ultra-high-speed video camera, we can follow the evolution of the air disk as it contracts into a bubble under the centre of the drop. The initial size and contraction speed of the disk were measured for a range of impact Weber and Reynolds numbers. The size of the initial disk is related to the bottom curvature of the drop at the initial contact, as measured in free-fall. The initial contact often leaves behind a ring of micro-bubbles, marking its location. The air disk contracts at a speed comparable to the corresponding air disks caught under a drop impacting onto a liquid surface. This speed also seems independent of the wettability of the liquid, which only affects the azimuthal shape of the contact line. For some impact conditions, the dynamics of the contraction leaves a small droplet at the centre of the bubble. This arises from a capillary wave propagating from the edges of the contracting disk towards the centre. As the wave converges its amplitude grows until it touches the solid substrate, thereby pinching off the micro-droplet at the plate, in the centre of the bubble. The effect of increasing liquid viscosity is to slow down the contraction speed and to produce a more irregular contact line leaving more micro-bubbles along the initial ring.

The averaged equations of slow flow in random arrays of fixed spheres are developed as a hierarchy of integro-differential equations, and an iteration procedure is described for obtaining the mean drag in the case of small volume concentration c. The leading approximation is that given by Brinkman's model of flow past a single fixed sphere, in which the effects of all other spheres are treated as a Darcy resistance. The higher approximations take account of the modification to the mean flow, particularly in the near field, due to the localized nature of the actual resistance. Thus the second approximation finds the change due to a second sphere, and averages over all its possible positions. The result for the mean drag confirms Childress’ terms in clogc and c (apart from an arithmetical correction to the latter), but indicates that for practical values of c numerical evaluation of integrals is needed, rather than expansion in powers of c and log c. The last section of the paper develops the corresponding results for flow through random arrays of fixed parallel circular cylinders.

A study has been made of the wake of a cylinder vibrating in line with an incident steady flow. The Reynolds number for the experiments was 190, and the vortex shedding was at all times synchronized with the vibrations of the cylinder, which were in a range of frequencies near twice the Strouhal shedding frequency for the stationary cylinder. Two distinct vortex wake patterns were encountered. The first is a complex regime in which two vortices are shed during each cycle of the vibration and form an alternating pattern of vortex pairs downstream. The second pattern is an alternating street which results from the shedding of a single vortex during each cycle of the cylinder's motion. The street geometry in the latter case shares many basic characteristics with the wake of a cylinder vibrating in cross-flow. These include the effects of vibration amplitude and frequency on the longitudinal and transverse spacing of the vortices. The results obtained from these experiments in air are in agreement with previous findings from free- and forced-vibration experiments in water at both higher and lower Reynolds numbers.

Wave-induced oscillations in harbours of constant depth but arbitrary shape in the horizontal plane connected to the open-sea are investigated both theoretically and experimentally. A theory termed the ‘arbitrary-shape harbour’ theory is developed. The solution of the Helmholtz equation is formulated as an integral equation which is then approximated by a matrix equation. The final solution is obtained by equating, at the harbour entrance, the wave amplitudes and their normal derivatives obtained from the solutions for the regions outside and inside the harbour. Special solutions using the method of separation of variables for the region inside circular and rectangular harbours are also obtained. Experiments were conducted to verify the theories. Four specific harbours were investigated: two circular harbours with 10° and 60° openings respectively, a rectangular harbour, and a model of the East and West Basins of Long Beach Harbour, California. In each case, the theoretical results agreed well with the experimental data.

The Rayleigh number R, in a horizontal layer with temperature-dependent viscosity can be based on the viscosity at T0, the mean of the boundary temperatures. The critical Rayleigh number Roc for fluids with exponential and super-exponential viscosity variation is nearly constant at low values of the ratio of the viscosities at the top and bottom boundaries; increases at moderate values of the viscosity ratio, reaching a maximum at a ratio of about 3000, and then decreases. This behaviour is explained by a simple physical argument based on the idea that convection begins first in the sublayer with maximum Rayleigh number. The prediction of Palm (1960) that certain types of temperature-dependent viscosity always decrease Roc is confirmed by numerical results but is not relevant to the viscosity variations typical of real liquids. The infinitesimal-amplitude state assumed by linear theory in calculating Roc does not exist because the convection jumps immediately to a finite amplitude at R0c. We observe a heat-flux jump at R0c exceeding 10% when the viscosity ratio exceeds 150. However, experimental measurements of R0c for glycerol up to a viscosity ratio of 3400 are in good agreement with the numerical predictions when the effects of a temperature-dependent expansion coefficient and thermal diffusivity are included.

The problem of electric charge convection in a dielectric liquid layer of high ionic purity, when subjected to unipolar injection, is in many ways analogous to that of thermal convection in a horizontal fluid layer heated from below, although no formal analogy can be established. The problem treated is intrinsically more nonlinear than the thermal problem. We consider two asymptotic states of convection: one where the whole motion is dominated by viscosity, and one where inertial effects dominate. In each state, two or three spatial regions are distinguished. From the approximate equations that hold in the different regions, information about the variation of the different quantities with distance from the injector is obtained, and further approximations permit us to establish the dependence of the current density ratio I/I0 (called the electric Nusselt number) on the stability parameter T = M2R = εϕ0/Kρν, and on 1/R = ν/Kϕ0, which is an equivalent Prandtl number (ε is the permittivity, ρ the fluid density, K the mobility, ν the kinematic viscosity, and ϕ0 the applied voltage). In the viscous state, the analysis gives I/I0 ∞ T½; in the inertial state the law I/I0 ∞ (T/R)1/4 = M½ is obtained. Since M is independent of the applied voltage, the latter law shows the saturation in the electric Nusselt number observed in earlier experiments. The transition in the states is associated with a transition number (MR)T [gap ] 30, which is an electric Reynolds number, related to an ordinary Reynolds number of about 10.

The experimental results, obtained in liquids of very different viscosities and dielectric constants, verify these theoretical predictions; further, they yield more precise numerical coefficients. As for the transition criteria, the experiments confirm that the viscous and inertial effects are of the same order when Re [gap ] 10. It was also possible to determine roughly the limits of the viscous and inertial states. The viscous analysis remains valid up to a Reynolds number of about 1; the inertial state can be considered valid down to a Reynolds number of 60. Schlieren observations show that the motion has the structure of very stable hexagonal cells at applied voltages just above the critical voltage, which are transformed into unstable filaments when the voltage is increased further. At even higher voltages, the motion finally breaks down into turbulence. It may be of interest to point out that, when M < 3, the electric Nusselt number approaches 1, which is equivalent to the situation in thermal convection at low Prandtl numbers.

The scenario of transition to chaos for a sphere falling or ascending under the action of gravity in a Newtonian fluid is investigated by numerical simulation. The mathematical formulation is parameterized using two non-dimensional parameters: the solid/fluid density ratio and the generalized Galileo number expressing the ratio between the gravity–buoyancy and viscosity effects. The study is carried out fully in this two-parameter space. The results show that for all density ratios the vertical fall or ascension becomes unstable via a regular axisymmetry breaking bifurcation. This bifurcation sets in slightly earlier for light spheres than for dense ones. A steady oblique fall or ascension follows before losing stability and giving way to an oscillating oblique movement. The secondary Hopf bifurcation is shown not to correspond to that of a fixed sphere wake for density ratios lower than 2.5, for which the oscillations have a significantly lower frequency. Trajectories of falling spheres become chaotic directly from the oblique oscillating regime. Ascending spheres present a specific behaviour before reaching a chaotic regime. The periodically oscillating oblique regime undergoes a subharmonic transition yielding a low-frequency oscillating ascension which is vertical in the mean (zigzagging regime). In all these stages of transition, the trajectories are planar with a plane selected randomly during the axisymmetry breaking. The chaotic regime appears to result from an interplay of a regular and of an additional Hopf bifurcation and the onset of the chaotic regime is accompanied by the loss of the remaining planar symmetry. The asymptotic chaotic states present an intermittent character, the relaminarization phases letting the subcritical plane and periodic trajectories reappear.

Direct numerical simulations of a turbulent fluid laden with finite-sized particles are performed. The computations, on a 1283 grid along with a maximum of 262 144 particles, incorporated both direct particle interactions via hard-sphere collisions and particle feedback. The ‘reverse’ coupling is accomplished in a manner ensuring correct discrete energy conservation (Sundaram & Collins 1996). A novel two-field formalism (Sundaram & Collins 1994a) is employed to calculate two-point correlations and equivalent spectral densities. An important consideration in these simulations is the initial state of fluid and particles. That is, the initial conditions must be chosen so as to allow a meaningful comparison of the different runs. Using such a carefully initialized set of runs, particle inertia was observed to increase both the viscous and drag dissipations; however, simultaneously, it also caused particle velocities to correlate for longer distances. The combination of effects suggests a mechanism for turbulence enhancement or suppression that depends on the parameter values. Like previous investigators, ‘pivoting’ or crossover of the fluid energy spectra was observed. A possible new scaling for this phenomenon is suggested. Furthermore, investigations of the influence of particle mass and number densities on turbulence modulation are also carried out.

This paper is concerned with the dynamics of large bubbles subject to various strengths of buoyancy effects, which are associated with applications for underwater explosion. The bubble is produced by electric discharge in a low-pressure tank to enhance the buoyancy effects. Experiments are carried out for a bubble in an infinite field, below a free surface and above a rigid boundary. The effects of buoyancy are reflected by the dimensionless parameter ${\it\delta}=\sqrt{{\it\rho}gR_{m}/(p_{amb}-p_{v})}$, where $R_{m}$, $p_{amb}$, $p_{v}$, ${\it\rho}$ and $g$ are the maximum bubble radius, ambient pressure, saturated vapour pressure, density of water and the acceleration of gravity respectively. A systematic study of buoyancy effects is carried out for a wide range of ${\it\delta}$ from 0.034 to 0.95. A series of new phenomena and new features is observed. The bubbles recorded are transparent, and thus we are able to display and study the jet formation, development and impact on the opposite bubble surface as well as the subsequent collapsing and rebounding of the ring bubble. Qualitative analyses are carried out for the bubble migration, jet velocity and jet initiation time, etc. for different values of ${\it\delta}$. When a bubble oscillates below a free surface or above a rigid boundary, the Bjerknes force due to the free surface (or rigid boundary) and the buoyancy are in opposite directions. Three situations are studied for each of the two configurations: (i) the Bjerknes force being dominant, (ii) the buoyancy force being dominant and (iii) the two forces being approximately balanced. For case (iii), we further consider two subcases, where both the balanced Bjerknes and buoyancy forces are weak or strong. When the Bjerknes and buoyancy forces are approximately balanced over the pulsation, some representative bubble behaviours are observed: the bubble near free surface is found to split into two parts jetting away from each other for small ${\it\delta}$, or involutes from both top and bottom for large ${\it\delta}$. A bubble above a rigid wall is found to be subject to contraction from the lateral part leading to bubble splitting. New criteria are established based on experimental results for neutral collapses where there is no dominant jetting along one direction, which correlate well with the criteria of Blake et al. (J. Fluid Mech., vol. 170, 1986, pp. 479–497; J. Fluid Mech., vol. 181, 1987, pp. 197–212) but agree better with the experimental and computational results.

In addition to the random processes described as horizontal turbulence, there are certain more regular processes which may contribute to horizontal mixing. One of these occurs in a shearing current, where the vertical gradient of velocity combined with vertical turbulent mixing leads to an effective diffusion in the horizontal direction. It is shown that this effect occurs in an alternating flow, such as a tidal current, as well as in a steady flow. In estuaries and coastal waters horizontal mixing by the shear effect may be associated with tidal currents, density currents or wind-driven currents. In each case the effective coefficient of horizontal diffusion, Kx or Ky, is inversely proportional to the coefficient of vertical eddy diffusion Kz. The occurrence of a stable gradient of density therefore increases the effective horizontal mixing very considerably. Results obtained from observations in the Mersey estuary and Irish Sea are compared with theoretical estimates of these effects.

The noise produced by mean flow-turbulence interaction of a circular subsonic jet is investigated theoretically, and expanded in azimuthal constituents of the turbulent pressure fluctuations. It is found that the low-order azimuthal constituents are the most efficient sound sources. On the basis of pressure correlation measurements, the azimuthal constituents are determined in a low Mach number jet. It is found that, in a range of Strouhal numbers between 0·2 and 1, the first three to four azimuthal constituents clearly dominate over the rest of the turbulent source quantity. A strictly axisymmetric ring vortex model for the coherent structure of the turbulence is, however, shown to be inappropriate.

In Part 1 a study is made of the internal solitary wave on the pycnocline of a continuously stratified fluid. A Korteweg–de Vries (KdV) equation for the ‘interfacial’ displacement is developed following Benney's method for long nonlinear waves. Experiments were conducted in a long wave tank with the pycnocline at several different depths below the free surface, while keeping the total depth approximately constant. A step-like pool of light water, trapped behind a sliding gate, served as the initial disturbance condition. The number of solitons generated was verified to satisfy the prediction of inverse-scattering theory. The fully developed soliton was found to satisfy the KdV theory for all ratios of upper-layer thickness to total depth.

In Part 2 of this study we investigate experimentally the evolution and breaking of an internal solitary wave as it shoals on a sloping bottom connecting the deeper region where the waves were generated to a shallower shelf region. It is found through quantitative measurements that the onset of wave-breaking was governed by shear instability, which was initiated when the local gradient Richardson number became less than ¼. The internal solitary wave of depression was found to steepen at the back of the wave before breaking, in contrast with waves of elevation. Two slopes were used, with ratios 1:16 and 1:9, and the fluid was a Boussinesq fluid with weak stratification using brine solutions.

Conditionally sampled hot-wire measurements were taken in a pipe at low Reynolds numbers (2700 > Re > 2000) corresponding to the onset of turbulence as a result of a large perturbation in the flow. This type of transition gives rise to a turbulent puff which maintains itself indefinitely at around Re = 2200. The structure of puffs was investigated in some detail and was found to be very different from the structure of fully developed turbulent pipe flow. Nevertheless, it is independent of the character of the disturbance which created it. The purpose of the study was to gain some insight into the mechanism of transition in a pipe.

When buoyant fluid is released into the base of a crack in an elastic medjura the crack will propagate upwards, driven by the buoyancy of the fluid. Viscous fluid flow in such a fissure is described by the equations of lubrication theory with the pressure given by the sum of the hydrostatic pressure of the fluid and the elastic pressures exerted by the walls of the crack. The elastic pressure and the width of the crack are further coupled by an integro-differential equation derived from the theory of infinitesimal dislocations in an elastic medium. The steady buoyancy-driven propagation of a two-dimensional fluid-filled crack through an elastic medium is analysed and the governing equations for the pressure distribution and the shape of the crack are solved numerically using a collocation technique. The fluid pressure in the tip of an opening crack is shown to be very low. Accordingly, a region of relatively inviscid vapour or exsolved volatiles in the crack tip is predicted and allowed for in the formulation of the problem. The solutions show that the asymptotic width of the crack, its rate of ascent and the general features of the flow are determined primarily by the fluid mechanics; the strength of the medium and the vapour pressure in the crack tip affect only the local structure near the advancing tip of the crack. When applied to the transport of molten rock through the Earth's lithosphere by magma-fracture, this conclusion is of fundamental importance and challenges the geophysicist's usual emphasis on the controlling influence of fracture mechanics rather than that of fluid mechanics.

A new experimental technique is described for the study of the interactions between the large-scale vortical features in the two-dimensional mixing layer. Detector probes above and below the mixing layer are used to monitor the large-scale structure. Conditional sampling is performed in a moderate Reynolds number developing flow, by using phase and amplitude information from these detector probes. It is shown that significant Reynolds-stress production is associated with the pairing interaction in which two vortical structures combine to form a single, larger vortical structure.