Plane turbulent mixing between two streams of water which contained dilute chemical reactants was studied in a new blow-down water tunnel. In a diffusion-limited reaction, a pH indicator, phenolphthalein, in one stream mixed and reacted with a base, sodium hydroxide, in the other stream to form a visible reaction product. The product was found to exist, as expected, in concentrated regions associated with the large, span-wise-coherent structures of the turbulence. A transition in the mixing was observed in which the aqueous mixing product increased by an order of magnitude. The transition is a consequence of the appearance and development of small-scale three-dimensional motions in the flow. Downstream of the transition, the amount of mixing product was independent of Reynolds number (for an order-of-magnitude increase in the latter) and at most only weakly dependent on Schmidt number.

The energy budget and dissipation mechanisms during droplet impact on solid surfaces are studied numerically and theoretically. We find that for high impact velocities and negligible surface friction at the solid surface (i.e. free slip), approximately one-half of the initial kinetic energy is transformed into surface energy, independent of the impact parameters and the detailed energy loss mechanism(s). We argue that this seemingly universal rule is related to the deformation mode of the droplet and is reminiscent of pipe flow undergoing a sudden expansion, for which the head loss can be calculated by multiplying the kinetic energy of the incoming flow by a geometrical factor. For impacts on a no-slip surface also dissipation in the shear boundary layer at the solid surface is important. In this case the geometric head loss acts as a lower bound on the total dissipation (i.e. the spreading on a no-slip surface approaches that on a free-slip surface when the droplet viscosity is sent to zero). This new view on the impact problem allows for simple analytical estimates of the maximum spreading diameter of impacting drops as a function of the impact parameters and the properties of the solid surface. It bridges the gap between previous momentum balance approaches and energy balance approaches, which hitherto did not give consistent predictions in the low viscosity limit. Good agreement is found between our models and experiments, both for impacts on ‘slippery’ or lubricated surfaces (e.g. Leidenfrost droplet impacts and head-on droplet–droplet collisions) and for impacts on no-slip surfaces.

In this paper the normal-mode small-amplitude waves of high-speed jets are investigated analytically and computationally. Three families of instability waves, each having a distinct wave pattern and propagation characteristics, have been found. One of the families of waves is the familiar Kelvin-Helmholtz instability wave. The other two families of waves do not appear to have been clearly identified and systematically studied before. Waves of one of the new wave family propagate with supersonic phase velocities relative to the ambient gas. They are, therefore, referred to as supersonic instability waves. Waves of the other family have subsonic phase velocities. Accordingly they are called subsonic waves. The subsonic waves have the unusual property that they are unstable only for jets with finite thickness mixing layers. They are neutral waves when calculated by a vortex-sheet jet model.

Earlier Oertel (1979, 1980, 1982) using a novel optical technique observed in a series of experiments three sets of waves in high-speed jets. The origin of these waves, however, remains so far unexplained and a theory has yet to be developed. In the present study it will be shown that the computed wave patterns and propagation characteristics of the Kelvin-Helmholtz, the supersonic and the subsonic instability waves match essentially those observed by Oertel. The physical mechanisms which give rise to the three families of waves as well as some of the most salient characteristic features of each set of waves are discussed and reported here.

In this paper, we perform direct numerical simulation (DNS) of turbulent boundary layers at Mach 5 with the ratio of wall-to-edge temperature Tw/Tδ from 1.0 to 5.4 (Cases M5T1 to M5T5). The influence of wall cooling on Morkovin's scaling, Walz's equation, the standard and modified strong Reynolds analogies, turbulent kinetic energy budgets, compressibility effects and near-wall coherent structures is assessed. We find that many of the scaling relations used to express adiabatic compressible boundary-layer statistics in terms of incompressible boundary layers also hold for non-adiabatic cases. Compressibility effects are enhanced by wall cooling but remain insignificant, and the turbulence dissipation remains primarily solenoidal. Moreover, the variation of near-wall streaks, iso-surface of the swirl strength and hairpin packets with wall temperature demonstrates that cooling the wall increases the coherency of turbulent structures. We present the mechanism by which wall cooling enhances the coherence of turbulence structures, and we provide an explanation of why this mechanism does not represent an exception to the weakly compressible hypothesis.

The Chapman–Enskog expansion is generalized in order to derive constitutive relations for flows of inelastically colliding spheres in three dimensions – to Burnett order. To this end, the pertinent (nonlinear) Boltzmann equation is perturbatively solved by performing a (double) expansion in the Knudsen number and the degree of inelasticity. One of the results is that the normal stress differences and the ‘temperature anisotropy’, characterizing granular fluids, are Burnett effects. The constitutive relations derived in this work differ, both qualitatively and quantitatively, from those obtained in previous studies. In particular, the Navier–Stokes (order) terms have a different dependence on the degree of inelasticity and the number density than in previously derived constitutive relations; for instance, the expression for the heat flux contains a term which is proportional to ε∇ log n, where ε is a measure of the degree of inelasticity and n denotes the number density. This contribution to the heat flux is of zeroth order in the density; a similar term, i.e. one that is proportional to ε∇n, has been previously obtained by using the Enskog correction but this term is O(n) and it vanishes in the Boltzmann limit. These discrepancies are resolved by an analysis of the Chapman–Enskog and Grad expansions, pertaining to granular flows, which reveals that the quasi-microscopic rate of decay of the temperature, which has not been taken into account heretofore, provides an important scale that affects the constitutive relations. Some (minor) quantitative differences between our results and previous ones exist as well. These are due to the fact that we take into account an isotropic correction to the leading Maxwellian distribution, which has not been considered before, and also because we consider the full dependence of the corrections to the Maxwellian distribution on the (fluctuating) speed.

Structural characteristics of transitionally rough and fully rough turbulent boundary layers are presented. These were measured in flows at different roughness Reynolds numbers developing over uniform spheres roughness. Inner regions of the longitudinal component of normal Reynolds stress profiles and log regions of mean profiles continuously change in the transitionally rough regime, as the roughness Reynolds number, Rek, varies. These properties asymptotically approach fully rough behaviour as Rek increases, and smooth behaviour at low Rek Profiles of other Reynolds-stress tensor components, turbulence kinetic energy, turbulence-kinetic-energy production, and the turbulence-kinetic-energy dissipation are also given, along with appropriate scaling variables. Fully rough, one-dimensional spectra of longitudinal velocity fluctuations from boundary-layer inner regions are similar to smooth-wall results for k1y > 0.2 when non-dimensionalized using distance from the wall y as the lengthscale, and (τ/ρ)½ as the velocity scale, where τ is local shear stress, ρ is static density, and k1 is one-dimensional wavenumber in the flow direction.

The thermohaline stability problem previously treated by Stern, Walin and Veronis is examined in greater detail. An error in an earlier paper is corrected and some new calculations made. It is shown, for instance, that direct convection can occur for thermal Rayleigh number R much less than 100 Rs when Rs [gsim ] 0·1, where Rs is the salinity Rayleigh number. A graphical presentation is devised to show the relative importance of the different terms in the equations of motion as a function of R and Rs. The most unstable mode over all wave-numbers for each R, Rs is found and it is shown that where both unstable direct and oscillating modes are present, the most unstable mode is direct in most cases.

Measurements in the mixing region of a 1 in. diameter cold air jet are described for Mach numbers ranging from 0.2 to 0.55. The statistical characteristics of the turbulence in the first few diameters of the flow may be expressed in terms of simple kinematic similarity relationships. These are based on the jet diameter and the distance downstream from the jet orifice as length-scales, and the inverse of the local shear as a time-scale. The experiments show that the integral time scale of the turbulence in a frame convected with the maximum energy of the turbulent motion is inversely proportional to the local shear.

The most interesting result obtained is that the local intensity of the turbulence is equal to 0.2 times the shear velocity. This velocity is defined as the product of the local integral length-scale of the turbulence with the local shear. The local intensity is defined as the R.M.S. value of the local velocity fluctuations divided by the jet efflux velocity. It was found that the length-scale is proportional to the distance from the jet orifice, while the maximum shear is also related to this distance as well as to the jet efflux velocity. These two similarity relations break down close to the jet orifice and change beyond the first six or so diameters downstream. The convection velocity is not equal to the local mean velocity but varies slowly over the region of maximum shear when it is just over half the jet efflux velocity. The measurements of other observers fit the relationships obtained quite well. From these relationships it is possible to calculate the noise generated by the mixing region of a given jet directly, using expressions derived by Lilley (1958).

A two-dimensional model of flow and bed topography in sinuous channels with erodible boundaries is developed and applied in order to investigate the mechanism of meander initiation. By reexamining the problem recently tackled by Ikeda, Parker & Sawai (1981), a previously undiscovered ‘resonance’ phenomenon is detected which occurs when the values of the relevant parameters fall within a neighbourhood of certain critical values. It is suggested that the above resonance controls the bend growth, and it is shown that it is connected in some sense with bar instability. In fact, by performing a linear stability analysis of flow in straight erodible channels, resonant flow in sinuous channels is shown to occur when curvature ‘forces’ a ‘natural’ solution represented by approximately steady perturbations of the alternate bar type. A comparison with experimental observations appears to support the idea that resonance is associated with meander formation.

An experimental investigation of entrainment and mixing in reacting and non-reacting turbulent mixing layers at large Schmidt number is presented. In non-reacting cases, a passive scalar is used to measure the probability density function (p.d.f.) of the composition field. Chemically reacting experiments employ a diffusion-limited acid–base reaction to directly measure the extent of molecular mixing. The measurements make use of laser-induced fluorescence diagnostics and high-speed, real-time digital image-acquisition techniques.

Our results show that the vortical structures in the mixing layer initially roll-up with a large excess of fluid from the high-speed stream entrapped in the cores. During the mixing transition, not only does the amount of mixed fluid increase, but its composition also changes. It is found that the range of compositions of the mixed fluid, above the mixing transition and also throughout the transition region, is essentially uniform across the entire transverse extent of the layer. Our measurements indicate that the probability of finding unmixed fluid in the centre of the layer, above the mixing transition, can be as high as 0.45. In addition, the mean concentration of mixed fluid across the layer is found to be approximately constant at a value corresponding to the entrainment ratio. Comparisons with gas-phase data show that the normalized amount of chemical product formed in the liquid layer, at high Reynolds number, is 50% less than the corresponding quantity measured in the gas-phase case. We therefore conclude that Schmidt number plays a role in turbulent mixing of high-Reynolds-number flows.

So-called ‘moonpools’ are vertical openings through the deck and hull of ships or barges, used for marine and offshore operations, such as pipe laying or recovery of divers. In the present study rectangular moonpools of large horizontal dimensions are considered. The natural modes of oscillation of the inner free surfaces are determined, under the assumption of infinite water depth and infinite length and beam of the barges that contain the moonpools. The problem is treated in two and three dimensions, via linearized potential flow theory. Results are given for the natural frequencies and the associated shapes of the free surface, for wide ranges of the geometric parameters. Simple quasi-analytical approximations are derived that yield the natural frequencies. The most striking result is that the natural frequencies of the longitudinal sloshing modes increase without bounds when both the draught and the width decrease to zero, the length of the moonpool being kept constant. As a corollary the problem of waves travelling in a channel through a rigid ice sheet is addressed and their dispersion equation is derived. The same behaviour is obtained: the waves travel increasingly faster as both the draught and the width of the channel are reduced.

The drag and lift forces acting on a rotating rigid sphere in a homogeneous linear shear flow are numerically studied by means of a three-dimensional numerical simulation. The effects of both the fluid shear and rotational speed of the sphere on the drag and lift forces are estimated for particle Reynolds numbers of 1[les ]Rep[les ]500.

The results show that the drag forces both on a stationary sphere in a linear shear flow and on a rotating sphere in a uniform unsheared flow increase with increasing the fluid shear and rotational speed. The lift force on a stationary sphere in a linear shear flow acts from the low-fluid-velocity side to the high-fluid-velocity side for low particle Reynolds numbers of Rep<60, whereas it acts from the high-velocity side to the low-velocity side for high particle Reynolds numbers of Rep>60. The change of the direction of the lift force can be explained well by considering the contributions of pressure and viscous forces to the total lift in terms of flow separation. The predicted direction of the lift force for high particle Reynolds numbers is also examined through a visualization experiment of an iron particle falling in a linear shear flow of a glycerin solution. On the other hand, the lift force on a rotating sphere in a uniform unsheared flow acts in the same direction independent of particle Reynolds numbers. Approximate expressions for the drag and lift coefficients for a rotating sphere in a linear shear flow are proposed over the wide range of 1[les ]Rep[les ]500.

We compare the turbulence statistics of the canopy/roughness sublayer (RSL) and the inertial sublayer (ISL) above. In the RSL the turbulence is more coherent and more efficient at transporting momentum and scalars and in most ways resembles a turbulent mixing layer rather than a boundary layer. To understand these differences we analyse a large-eddy simulation of the flow above and within a vegetation canopy. The three-dimensional velocity and scalar structure of a characteristic eddy is educed by compositing, using local maxima of static pressure at the canopy top as a trigger. The characteristic eddy consists of an upstream head-down sweep-generating hairpin vortex superimposed on a downstream head-up ejection-generating hairpin. The conjunction of the sweep and ejection produces the pressure maximum between the hairpins, and this is also the location of a coherent scalar microfront. This eddy structure matches that observed in simulations of homogeneous-shear flows and channel flows by several workers and also fits with earlier field and wind-tunnel measurements in canopy flows. It is significantly different from the eddy structure educed over smooth walls by conditional sampling based only on ejections as a trigger. The characteristic eddy was also reconstructed by empirical orthogonal function (EOF) analysis, when only the dominant, sweep-generating head-down hairpin was recovered, prompting a re-evaluation of earlier results based on EOF analysis of wind-tunnel data. A phenomenological model is proposed to explain both the structure of the characteristic eddy and the key differences between turbulence in the canopy/RSL and the ISL above. This model suggests a new scaling length that can be used to collapse turbulence moments over vegetation canopies.

This paper describes the results of a study examining the flow and acoustic characteristics of an axisymmetric supersonic jet issuing from a sonic and a Mach 1.5 converging–diverging (C–D) nozzle and impinging on a ground plane. Emphasis is placed on the Mach 1.5 nozzle with the sonic nozzle used mainly for comparison. A large-diameter circular plate was attached at the nozzle exit to measure the forces generated on the plate owing to jet impingement. The experimental results described in this paper include lift loss, particle image velocimetry (PIV) and acoustic measurements. Suckdown forces as high as 60% of the primary jet thrust were measured when the ground plane was very close to the jet exit. The PIV measurements were used to explain the increase in suckdown forces due to high entrainment velocities. The self-sustained oscillatory frequencies of the impinging jet were predicted using a feedback loop that uses the measured convection velocities of the large-scale coherent vortical structures in the jet shear layer. Nearfield acoustic measurements indicate that the presence of the ground plane increases the overall sound pressure levels (OASPL) by approximately 8 dB relative to a corresponding free jet. For moderately underexpanded jets, the influence of the shock cells on the important flow features was found to be negligible except for close proximity of the ground plane.

The advancing front of a buoyant plume which is being established in uniform surroundings has some properties in common with the plume, while in other respects it behaves more like a ‘thermal’ released from rest. The solutions for these two cases cannot be matched directly, since the dependence of velocity on height is different. It is shown here that a similarity solution, which is consistent with the equations describing both parts of the flow, can be obtained once it is recognized that the velocity of the front may be less than that of the steady plume. The cap moves with a constant fraction of the plume velocity at the same level, and the total buoyancy in the cap is increasing, so a modification of the simple relations for thermals is required.

This prediction is verified experimentally, and numerical values for the ratio of the velocities and the rate of increase of the cap radius with height determined. The extreme front of plume cap advances at about 0·6 times the mean velocity on the axis of a steady plume, and it spreads at just over half the angle of a thermal. This implies a smaller rate of entrainment and therefore a smaller rate of dilution per unit height compared with a thermal, especially since about half of the fluid mixed into the cap comes from the plume below. The model allows one to estimate the time necessary for convection from a known steady source on the ground to lead to the formation of a cloud.

We consider the linear stability of the cocurrent flow of two fluids of different viscosity in an infinite region (the viscous analogue of the classical Kelvin-Helmholtz problem). Attention is confined to the simplest case, Couette flow, and we solve the problem using both numerical and asymptotic techniques. We find that the flow is always unstable (in the absence of surface tension). The instability arises at the interface between the two fluids and occurs for short wavelengths, when viscosity rather than inertia is the dominant physical effect.

A hybrid multiscale method is developed for simulating micro- and nano-scale fluid flows. The continuum Navier–Stokes equation is used in one flow region and atomistic molecular dynamics in another. The spatial coupling between continuum equations and molecular dynamics is achieved through constrained dynamics in an overlap region. The proposed multiscale method is used to simulate sudden-start Couette flow and channel flow with nano-scale rough walls, showing quantitative agreement with results from analytical solutions and full molecular dynamics simulations for different parameter regimes. Potential applications of the proposed multiscale method are discussed.

The sound generated by a circular cylinder in a flow at low Mach numbers is investigated by direct solution of the two-dimensional unsteady compressible Navier–Stokes equations. Results show that sound pressure waves are generated primarily by vortex shedding from the cylinder surface into its wake. When a vortex is shed from one side of the cylinder, a negative pressure pulse is generated from that side whereas a positive pressure pulse is generated from the other side; alternate vortex shedding from the upper and lower sides of the cylinder produces negative and positive pulses alternately and thus produces sound pressure waves on both sides. The dipolar nature of the generated sound is confirmed; lift dipole dominates the sound field. The Doppler effect is shown to play an important role at finite Mach numbers. The direct solutions are also compared with the solutions obtained by Curle's acoustic analogy. The results show that Curle's solution describes well not only the generation mechanism of the sound but also the propagation process if we take the Doppler effect into consideration.

We investigate the near wake of a cylinder at values of Reynolds number corresponding to the onset and development of shear-layer instabilities. By combining quantitative experimental imaging (particle image velocimetry, PIV) and direct numerical simulations at $Re \,{=}\, 3900/4000$ and 10000, we show that the flow structure is notably altered. At higher Reynolds number, the lengths of both the wake bubble and the separating shear layer decrease substantially. Corresponding patterns of velocity fluctuations and Reynolds stress contract towards the base of the cylinder. The elevated values of Reynolds stress at upstream locations in the separated layer indicate earlier onset of shear-layer transition. These features are intimately associated with the details of the shear-layer instability, which leads to small-scale vortices. The simulated signatures of the shear-layer vortices are characterized by a broadband peak at $Re \,{=}\, 3900$ and a broadband high spectral-density ‘plateau’ at $Re \,{=}\, 10\,000$ in the power spectra. The shear-layer frequencies from the present direct numerical simulations study agree well with previous experimentally measured values, and follow the power law suggested by other workers.

The rapid shearing of a mixture of cohesionless glass spheres and air or water was studied in an annular, parallel-plate shear cell designed after Savage (1978). Two types of flow were observed. In the first type of flow the entire mass of the granular material was mobilized. At high shear rates the shear and normal stresses were found to be quadratically dependent upon the mean shear rate (at constant volume concentration), in general agreement with the observations of Bagnold (1954) and Savage & Sayed (1984), and the ‘kinetic’ theory of Jenkins & Savage (1983). The stresses were found to be weakly dependent on the volume concentration up to approximately 0.5, and strongly dependent above this concentration. For flows in which water was the interstitial fluid, the ratio of the shear stress to the normal stress was slightly higher (than in air), and the stresses at lower shear rates were found to be more nearly linearly related to the shear rate. It is suggested that these effects are contributed to by the viscous dampening of grain motions by the water. The second type of flow was distinguished by the existence of an internal boundary above which the granular material deformed rapidly, but below which the granular material remained rigidly locked in place. The thickness of the shearing layer was measured to be between 5 and 15 grain diameters. The stress ratio at the bottom of the shearing layer was found to be nearly constant, suggesting the internal boundary is a consequence of the immersed weight of the shearing grains, and may be described by a Coulomb yield criterion. A scaled concentration is proposed to compare similar data obtained using different-sized materials or different apparatus. An intercomparison of the two types of flow studied, along with a comparison between the present experiments and those of Bagnold (1954) and Savage & Sayed (1984), suggests that the nature of the boundaries can have a significant effect upon the dynamics of the entire flow.