The development of a two-dimensional viscous incompressible flow generated from a circular cylinder impulsively started into rectilinear motion is studied computationally. An adaptative numerical scheme, based on vortex methods, is used to integrate the vorticity/velocity formulation of the Navier–Stokes equations for a wide range of Reynolds numbers (Re = 40 to 9500). A novel technique is implemented to resolve diffusion effects and enforce the no-slip boundary condition. The Biot–Savart law is employed to compute the velocities, thus eliminating the need for imposing the far-field boundary conditions. An efficient fast summation algorithm was implemented that allows a large number of computational elements, thus producing unprecedented high-resolution simulations. Results are compared to those from other theoretical, experimental and computational works and the relation between the unsteady vorticity field and the forces experienced by the body is discussed.

Approximate constitutive equations are derived for a dilute suspension of rigid spheroidal particles with Brownian rotations, and the behaviour of the approximations is explored in various flows. Following the suggestion made in the general formulation in part 1, the approximations take the form of Hand's (1962) fluid model, in which the anisotropic microstructure is described by a single second-order tensor. Limiting forms of the exact constitutive equations are derived for weak flows and for a class of strong flows. In both limits the microstructure is shown to be entirely described by a second-order tensor. The proposed approximations are simple interpolations between the limiting forms of the exact equations. Predictions from the exact and approximate constitutive equations are compared for a variety of flows, including some which are not in the class of strong flows analysed.

The three-dimensional resonant interaction of a plane Tollmien-Schlichting wave, having a frequency f1, with a pair of oblique waves having frequencies ½ f1, was observed and studied experimentally. In the initial stages, the interaction proved to be a parametric resonance, resulting in the amplification of small random priming (background) oscillations of frequency ½ f1, and of a packet of low-frequency oscillations. The resonant interaction of waves in a boundary layer was investigated also by introducing a priming oscillation with frequency f’ = ½ f1 + Δf for different values of the frequency detuning Δf. The importance of the discovered wave interaction in boundary-layer transition is demonstrated. Causes of realization of different types of laminar-flow breakdown are discussed.

Direct numerical simulations have been carried out for a fully developed turbulent channel flow with a smooth upper wall and a lower wall consisting of square bars separated by a rectangular cavity. A wide range of $w/k$, the cavity width to roughness height ratio, was considered. For $w/k\,{\ge}\,7$, recirculation zones occur immediately upstream and downstream of each element while mean streamlines and spatial distributions of the skin frictional drag indicate that each element is virtually isolated. The maximum form drag occurs at $w/k\,{=}\,7$ and coincides with the minimum skin frictional drag. The dependence on $w/k$ of the Clauser roughness function reflects that of the form drag.

The behaviour of heavy particles in isotropic, homogeneous, decaying turbulence has been experimentally studied. The settling velocity of the particles has been found to be much larger than in a quiescent fluid. It has been determined that the enhancement of the settling velocity depends on the particle loading, increasing as the volume fraction of particles in the flow increases. The spatial and temporal distribution of the particle concentration field is shown to exhibit large inhomogeneities. As the particles interact with the underlying turbulence they concentrate preferentially in certain regions of the flow. A characteristic dimension of these particle clusters is found to be related to the viscous scales of the flow. Measurements of the settling velocity conditioned on the local concentration of particles in the flow have shown that there is a monotonic increase in the settling velocity with the local concentration (the relation being quasi-linear). A simple phenomenological model is proposed to explain this behaviour.

In an experimental study of the vortex-induced oscillations of a long flexible circular cylinder, the observed stationary amplitudes describe an hysteresis loop partially different from earlier studies. Each branch of the loop is associated with a vortex shedding mode and, as a jump from one branch to the other occurs, the phase difference between the cylinder displacement and the vortex shedding undergoes an abrupt change. The critical flow velocities at which the jump occurs concur with the flow visualization observations of Williamson & Roshko (1988) on the vortex shedding modes near the fundamental synchronization region. Impulsive regimes, obtained at a given flow velocity with the cylinder initially at rest or pre-excited, and progressive regimes resulting from a variation of the flow velocity, are examined. The occurrence of bifurcations is detected for a flow velocity range in the case of the impulsive regimes. The coordinates of the bifurcations define a boundary between two vortex shedding modes, a boundary that verifies the critical curve obtained by Williamson & Roshko (1988). The experimental set-up of this study simulates half the wavelength of a vibrating cable, eliminates the end effects present in oscillating rigid cylinder set-up and has one of the lowest damping ratios reported for the study of this phenomenon.

The flow around a circular cylinder has been examined over the Reynolds number range 105 to 7·5 × 105, Reynolds number being based on cylinder diameter. Narrow-band vortex shedding has been observed up to a Reynolds number of 5·5 × 105, i.e. well into the critical régime. At this Reynolds number the Strouhal number reached the unusually high value of 0·46. Spectra of the velocity fluctuations measured in the wake are presented for several values of Reynolds number.

The development of a small three-dimensional disturbance in plane Poiseuille flow is considered. Its kinetic energy is expressed in terms of the velocity and vorticity components normal to the wall. The normal vorticity develops according to the mechanism of vortex stretching and is described by an inhomogeneous equation, where the spanwise variation of the normal velocity acts as forcing. To study specifically the effect of the forcing, the initial normal vorticity is set to zero and the energy density in the wavenumber plane, induced by the normal velocity, is determined. In particular, the response from individual (and damped) Orr–Sommerfeld modes is calculated, on the basis of a formal solution to the initial-value problem. The relevant timescale for the development of the perturbation is identified as a viscous one. Even so, the induced energy density can greatly exceed that associated with the initial normal velocity, before decay sets in. Initial conditions corresponding to the least-damped Orr–Sommerfeld mode induce the largest energy density and a maximum is obtained for structures infinitely elongated in the streamwise direction. In this limit, the asymptotic solution is derived and it shows that the spanwise wavenumbers at which the largest amplification occurs are 2.60 and 1.98, for symmetric and antisymmetric normal vorticity, respectively. The asymptotic analysis also shows that the propagation speed for induced symmetric vorticity is confined to a narrower range than that for antisymmetric vorticity. From a consideration of the neglected nonlinear terms it is found that the normal velocity component cannot be nonlinearly affected by the normal vorticity growth for structures with no streamwise dependence.

The purpose of this article is to show that the Navier–Stokes equations can be rewritten as a set of linearized inhomogeneous Euler equations (in convective form) with source terms that are exactly the same as those that would result from externally imposed shear-stress and energy-flux perturbations. These results are used to develop a mathematical basis for some existing and potential new jet noise models by appropriately choosing the ‘base flow’ about which the linearization is carried out.

Experiments were conducted in a two-stream planar mixing layer at convective Mach numbers, Mc, of 0.28, 0.42, 0.50, 0.62 and 0.79. Planar laser Mie scattering (PLMS) from a condensed alcohol fog and planar laser-induced fluorescence (PLIF) of nitric oxide were used for flow visualization in the side, plan and end views. The PLIF signals were also used to characterize the turbulent mixture fraction fluctuations.

Visualizations using PLMS indicate a transition in the turbulent structure from quasi-two-dimensionality at low convective Mach number, to more random three-dimensionality for $M_c\geqslant 0.62$. A transition is also observed in the core and braid regions of the spanwise rollers as the convective Mach number increases from 0.28 to 0.62. A change in the entrainment mechanism with increasing compressibility is also indicated by signal intensity profiles and perspective views of the PLMS and PLIF images. These show that at Mc = 0.28 the instantaneous mixture fraction field typically exhibits a gradient in the streamwise direction, but is more uniform in the cross-stream direction. At Mc = 0.62 and 0.79, however, the mixture fraction field is more streamwise uniform and with a gradient in the cross-stream direction. This change in the composition of the structures is indicative of different entrainment motions at the different compressibility conditions. The statistical results are consistent with the qualitative observations and suggest that compressibility acts to reduce the magnitude of the mixture fraction fluctuations, particularly on the high-speed edge of the layer.

A porous-wavemaker theory is developed to analyse small-amplitude surface waves on water of finite depth, produced by horizontal oscillations of a porous vertical plate. Analytical solutions in closed forms are obtained for the surface-wave profile, the hydrodynamic-pressure distribution and the total force on the wavemaker. The influence of the wave-effect parameter C and the porous-effect parameter G, both being dimensionless, on the surface waves and on the hydrodynamic pressures is discussed in detail.

A turbulent mixing layer in a water channel was observed at Reynolds numbers up to 3 × 106. Flow visualization with dyes revealed (once more) large coherent structures and showed their role in the entrainment process; observation of the reaction of a base and an acid indicator injected on the two sides of the layer, respectively, gave some indication of where molecular mixing occurs. Autocorrelations of streamwise velocity fluctuations, using a laser-Doppler velocimeter (LDV) revealed a fundamental periodicity associated with the large structures. The surprisingly long correlation times suggest time scales much longer than had been supposed; it is argued that the mixing-layer dynamics at any point are coupled to the large structure further downstream, and some possible consequences regarding the effects of initial conditions and of the influence of apparatus geometry are discussed.

Diffuse-interface models may be used to compute moving contact lines because the Cahn–Hilliard diffusion regularizes the singularity at the contact line. This paper investigates the basic questions underlying this approach. Through scaling arguments and numerical computations, we demonstrate that the Cahn–Hilliard model approaches a sharp-interface limit when the interfacial thickness is reduced below a threshold while other parameters are fixed. In this limit, the contact line has a diffusion length that is related to the slip length in sharp-interface models. Based on the numerical results, we propose a criterion for attaining the sharp-interface limit in computing moving contact lines.

These exact solutions correspond to time-dependent motions in parabolic basins. A characteristic feature is that the shoreline is not fixed. It is free to move and must be determined as part of the solution. In general, the motion is oscillatory and has the appropriate small-amplitude limit. For the case in which the parabolic basin reduces to a flat plane, there is a solution for a flood wave. These solutions provide a valuable test for numerical models of inundating storm tides.

The extent or existence of similarities between fully developed turbulent pipes and channels, and in zero-pressure-gradient turbulent boundary layers has come into question in recent years. This is in contrast to the traditionally accepted view that, upon appropriate normalization, all three flows can be regarded as the same in the near-wall region. In this paper, the authors aim to provide clarification of this issue through streamwise velocity measurements in these three flows with carefully matched Reynolds number and measurement resolution. Results show that mean statistics in the near-wall region collapse well. However, the premultiplied energy spectra of streamwise velocity fluctuations show marked structural differences that cannot be explained by scaling arguments. It is concluded that, while similarities exist at these Reynolds numbers, one should exercise caution when drawing comparisons between the three shear flows, even near the wall.

Velocity measurements in a 51 mm diameter turbulent jet are presented. The measurement programme is conducted in two parts. The first part is devoted to the validation of laser velocimeter (LV) data. This consists of comparative measurements with the LV and a hot-wire anemometer. The second part consists of a survey of the jet flow field at Mach 0·28, 0·90, and 1·37 under ambient temperature conditions. Radial and centre-line distributions of the axial and radial, mean and fluctuating velocities are obtained. The distributions indicate a decrease in the spreading rate of the mixing layer with increasing Mach number and a corresponding lengthening of the potential core. The results further indicate that these two parameters vary with the square of the jet Mach number. Radial distributions collapse when plotted in terms of ση*, where σ = 10.7/(1 - 0.273 MJ2) and η* = (r − r0·5)/x. This is true for distributions in planes located as far downstream as two potential core lengths. The collapsed data of mean velocity can be approximated by a Görtler error function profile: \[ U/U_J = 0.5[1-{\rm erf}(\sigma\eta^{*})]. \] Centre-line distributions at various Mach numbers also collapse if plotted in terms of x/xc, xc being the potential core length. A general equation for the collapsed data of mean velocity is given by: U/UJ = 1 - exp{1.35/(1 - x/xc)}, for the range of Mach numbers 0·3-1·4, where xc = 4.2 + 1.1 MJ2.

From an analysis of small oscillations of a viscous fluid droplet immersed in another viscous fluid a general dispersion equation is derived by which frequency and rate of damping of oscillations can be calculated for arbitrary values of droplet size, physical properties of the fluids, and interfacial viscosity and elasticity coefficients. The equation is studied for two distinct extremes of interfacial characteristics: (i) a free interface between the two fluids in which only a constant, uniform interfacial tension acts; (ii) an ‘inextensible’ interface between the two fluids, that is, a highly condensed film or membrane which, to first order, cannot be locally expanded or contracted. Results obtained are compared with those previously published for various special cases.

When the viscosities of both fluids are low, the primary contribution to the rate of damping of oscillations is generally the viscous dissipation in a boundary layer near the interface, in both the free and inextensible interface situations. For this reason inviscid velocity profiles, which do not account for the boundarylayer flow, do not lead to good approximations to the damping rate. The two exceptions in which the approximation based on inviscid profiles is adequate occur when the interface is free and either the interior or exterior fluid is a gas of negligible density and viscosity.

An experimental study has been made of the motion of long bubbles in closed tubes. The influence of viscosity and surface tension on the bubble velocity is clarified. A correlation of bubble velocities in vertical tubes is suggested and is shown to be useful for the whole range of parameters investigated. In addition, the effect of tube inclination angle on bubble velocity is presented, and certain features of the flow are described qualitatively.

We study the temporal stability of the Orr–Sommerfeld and Squire equations in channels with turbulent mean velocity profiles and turbulent eddy viscosities. Friction Reynolds numbers up to $Re_\tau \,{=}\, 2\,{\times}\, 10^4$ are considered. All the eigensolutions of the problem are damped, but initial perturbations with wavelengths $\lambda_x \,{>}\, \lambda_z$ can grow temporarily before decaying. The most amplified solutions reproduce the organization of turbulent structures in actual channels, including their self-similar spreading in the logarithmic region. The typical widths of the near-wall streaks and of the large-scale structures of the outer layer, $\lambda_z^+ \,{=}\, 100$ and $\lambda_z/h \,{=}\, 3$, are predicted well. The dynamics of the most amplified solutions is roughly the same regardless of the wavelength of the perturbations and of the Reynolds number. They start with a wall-normal $v$ event which does not grow but which forces streamwise velocity fluctuations by stirring the mean shear ($uv\,{<}\,0$). The resulting $u$ fluctuations grow significantly and last longer than the $v$ ones, and contain nearly all the kinetic energy at the instant of maximum amplification.

This paper describes an experimental investigation of the form of the large scale motions in turbulent flow. These motions have been found to be more ordered than has usually been supposed and their origin and dynamics are discussed in terms of physical models of typical eddies.

Nine components of the double velocity correlation tensor have been measured at a number of positions in the wake of a circular cylinder and in a ‘flat plate’ boundary layer. These have been supplemented by measurements of correlations with separations in directions other than the axial ones. In the wake, the correlations at large values of the separation are explained in terms of two types of large scale motion. One of these is a pair of vortices, side by side and rotating in opposite directions with axes aligned approximately normal to the plane of the wake. The other typical motion is a series of jets in which turbulent fluid is projected outward from the core of the wake. It is suggested that these are the result of an instability of the turbulent shear stress. A qualitative explanation of the apparent structural equilibrium of the wake is given in terms of this instability. The vortex pair eddies were not found in the boundary layer but there is evidence of jets much like those in the wake.

Correlations measured in grid turbulence have been found to be highly anisotropic and consistent with the presence of vortex pair eddies. When a plane strain was applied to grid turbulence, the effect on the correlations suggested the presence of a stress instability similar to that postulated for the wake.