It has been argued that the no-slip boundary condition, applicable when a viscous fluid flows over a solid surface, may be an inevitable consequence of the fact that all such surfaces are, in practice, rough on a microscopic scale: the energy lost through viscous dissipation as a fluid passes over and around these irregularities is sufficient to ensure that it is effectively brought to rest. The present paper analyses the flow over a particularly simple model of such a rough wall to support these physical ideas.

The structure of the vorticity fields in homogeneous turbulent shear flow and various irrotational straining flows is examined using results from direct numerical simulations of the unsteady, incompressible Navier-Stokes equations with up to 128 × 128 × 128 grid points. In homogeneous shear flow, the distribution of the inclination angle of the vorticity vectors and contour plots of two-point correlations of both velocity and vorticity are consistent with the existence of persistent vortical structures inclined with respect to the flow direction. Early in the development of these shear flows, the angle of inclination at which most of these structures are found is near 45°; after the flow develops, this angle lies between 35°–40°. Instantaneous vorticity-vector and vortex-line plots confirm the presence of hairpin vortices in this flow at the two Reynolds numbers simulated. These vortices are formed by the roll-up of sheets of mean spanwise vorticity. The average hairpin leg spacing decreases with increasing Reynolds number but increases relative to the Taylor microscale for developed shear flows. Examination of irrotational axisymmetric contraction, axisymmetric expansion, and plane strain flows shows, as expected, that the vorticity tends to be aligned with the direction of positive strain. For example, the axisymmetric contraction flow is dominated by coherent longitudinal vortices. Without the presence of mean shear, however, hairpin structures do not develop. The simulations strongly indicate that the vorticity occurs in coherent filaments that are stretched and strengthened by the mean strain. When compressed, these filaments appear to buckle rather than to decrease in strength.

An exact solution is found in a fairly simple form for two-dimensional progressive waves of arbitrary amplitude on a fluid of unlimited depth, when only surface tension and not gravity is taken into account as the restoring force. The calculated wave forms are exhibited graphically for various amplitudes, and the relation between wave velocity and amplitude is plotted. The wave of greatest height occurs when the vertical distance between trough and crest is 0·730 wavelengths (compared with 0·142 for gravity waves). Higher waves are prevented from appearing by the enclosing of air bubbles in the troughs.

The scattering of water waves by an array of N bottom-mounted vertical circular cylinders is solved exactly (under the assumption of linear water wave theory) using the method proposed by Spring & Monkmeyer in 1974. A major simplification to this theory has been found which makes the evaluation of quantities such as the forces on the cylinders much simpler. New formulae are given for the first and mean second-order forces together with one for the free-surface elevation in the vicinity of a particular cylinder. Comparisons are made between the exact results shown here and those generated using the approximate method of McIver & Evans (1984). The behaviour of the forces on the bodies in the long-wave limit is also examined for the special case of two cylinders with equal radii.

In this paper we describe a detailed study of the wake structures and flow dynamics associated with simulated two-dimensional flows past a circular cylinder that is either stationary or in simple harmonic cross-flow oscillation. Results are examined for Re=500 and a fixed motion amplitude of ymax/D=0.25. The study concentrates on a domain of oscillation frequencies near the natural shedding frequency of the fixed cylinder. In addition to the change in phase of vortex shedding with respect to cylinder motion observed in previous experimental studies, we note a central band of frequencies for which the wake exhibits long-time-scale relaxation oscillator behaviour. Time-periodic states with asymmetric wake structures and non-zero mean lift were also observed for oscillation frequencies near the lower edge of the relaxation oscillator band. In this regime we compute a number of bifurcations between different wake configurations and show that the flow state is not a unique function of the oscillation frequency. Results are interpreted using an analysis of vorticity generation and transport in the base region of the cylinder. We suggest that the dynamics of the change in phase of shedding arise from a competition between two different mechanisms of vorticity production.

A direct numerical simulation of a fully developed, low-Reynolds-number turbulent flow in a square duct is presented. The numerical scheme employs a time-splitting method to integrate the three-dimensional, incompressible Navier-Stokes equations using spectral/high-order finite-difference discretization on a staggered mesh; the nonlinear terms are represented by fifth-order upwind-biased finite differences. The unsteady flow field was simulated at a Reynolds number of 600 based on the mean friction velocity and the duct width, using 96 × 101 × 101 grid points. Turbulence statistics from the fully developed turbulent field are compared with existing experimental and numerical square duct data, providing good qualitative agreement. Results from the present study furnish the details of the corner effects and near-wall effects in this complex turbulent flow field; also included is a detailed description of the terms in the Reynolds-averaged streamwise momentum and vorticity equations. Mechanisms responsible for the generation of the stress-driven secondary flow are studied by quadrant analysis and by analysing the instantaneous turbulence structures. It is demonstrated that the mean secondary flow pattern, the distorted isotachs and the anisotropic Reynolds stress distribution can be explained by the preferred location of an ejection structure near the corner and the interaction between bursts from the two intersecting walls. Corner effects are also manifested in the behaviour of the pressure-strain and velocity-pressure gradient correlations.

Results of hot-wire measurements in an incompressible rectangular jet issuing into a quiet surrounding at ambient conditions are presented. The quantities measured include distributions of mean velocity and the turbulence shear stresses in the two central planes of the jet at stations up to 115 widths (small dimension of the nozzle) downstream of the nozzle exit. The flow field of the jet was found to be characterized by the presence of three distinct regions, defined by the axial mean velocity decay, which are referred to as: a potential core region, a two-dimensional-type region, and an axisymmetric-type region. The onset of the axisymmetric region occurs at a downstream location where the two shear layers from the short edges of the nozzle meet. In the central plane which contains the small dimension of the nozzle, similarity was found both in the mean velocity and shear stress profiles beyond 30 widths downstream of the nozzle exit; however, profiles of r.m.s. velocity show similarity in the second but not the third region. The mean velocity, shear stress and r.m.s. velocity profiles in the central plane containing the long dimension of the nozzle do not show geometrical similarity.

The structure of homogeneous turbulence subject to high shear rate has been investigated by using three-dimensional, time-dependent numerical simulations of the Navier–Stokes equations. The instantaneous velocity fields reveal that a high shear rate produces structures in homogeneous turbulence similar to the ‘streaks’ that are present in the sublayer of wall-bounded turbulent shear flows. Statistical quantities such as the Reynolds stresses are compared with those in the sublayer of a turbulent channel flow at a comparable shear rate made dimensionless by turbulent kinetic energy and its dissipation rate. This study indicates that high shear rate alone is sufficient for generation of the streaky structures, and that the presence of a solid boundary is not necessary.

Evolution of the statistical correlations is examined to determine the effect of high shear rate on the development of anisotropy in turbulence. It is shown that the streamwise fluctuating motions are enhanced so profoundly that a highly anisotropic turbulence state with a ‘one-component’ velocity field and ‘two-component’ vorticity field develops asymptotically as total shear increases. Because of high shear rate, rapid distortion theory predicts remarkably well the anisotropic behaviour of the structural quantities.

We consider the critical Weber number (Wec≡ ρV20D/σ) at which the transition from dripping to jetting occurs when a Newtonian liquid of density ρ and surface tension σ is injected with a velocity V0 through a tube of diameter D downward into stagnant air, under gravity g. We extend Taylor's (1959) model for the recession speed of a free edge, and obtain in the inviscid limit an exact solution which includes gravity and inertia effects. This solution provides a criterion for the transition which is shown to occur at a critical Weber number

formula here

where Bo and Boo are the Bond numbers (Bo≡[ρgD2/(2σ)]1/2), respectively based on the inside and outside diameter of the tube, and K is a constant equal to 0.37 for the case of water injected in air. This critical Weber number is shown to be in good agreement with existing experimental values as well as with new measurements performed over a wide range of Bond numbers.

In this paper, we perform direct numerical simulations (DNS) of turbulent boundary layers with nominal free-stream Mach number ranging from 0.3 to 12. The main objective is to assess the scalings with respect to the mean and turbulence behaviours as well as the possible breakdown of the weak compressibility hypothesis for turbulent boundary layers at high Mach numbers (M > 5). We find that many of the scaling relations, such as the van Driest transformation for mean velocity, Walz's relation, Morkovin's scaling and the strong Reynolds analogy, which are derived based on the weak compressibility hypothesis, remain valid for the range of free-stream Mach numbers considered. The explicit dilatation terms such as pressure dilatation and dilatational dissipation remain small for the present Mach number range, and the pressure–strain correlation and the anisotropy of the Reynolds stress tensor are insensitive to the free-stream Mach number. The possible effects of intrinsic compressibility are reflected by the increase in the fluctuations of thermodynamic quantities (p′rms/pw, ρ′rms/ρ, T′rms/T) and turbulence Mach numbers (Mt, M′rms), the existence of shocklets, the modification of turbulence structures (near-wall streaks and large-scale motions) and the variation in the onset of intermittency.

A direct numerical simulation of a supersonic turbulent boundary layer has been performed. We take advantage of a technique developed by Spalart for incompressible flow. In this technique, it is assumed that the boundary layer grows so slowly in the streamwise direction that the turbulence can be treated as approximately homogeneous in this direction. The slow growth is accounted for by a coordinate transformation and a multiple-scale analysis. The result is a modified system of equations, in which the flow is homogeneous in both the streamwise and spanwise directions, and which represents the state of the boundary layer at a given streamwise location. The equations are solved using a mixed Fourier and B-spline Galerkin method.

Results are presented for a case having an adiabatic wall, a Mach number of M = 2.5, and a Reynolds number, based on momentum integral thickness and wall viscosity, of Reθ′ = 849. The Reynolds number based on momentum integral thickness and free-stream viscosity is Reθ = 1577. The results indicate that the Van Driest transformed velocity satisfies the incompressible scalings and a small logarithmic region is obtained. Both turbulence intensities and the Reynolds shear stress compare well with the incompressible simulations of Spalart when scaled by mean density. Pressure fluctuations are higher than in incompressible flow. Morkovin's prediction that streamwise velocity and temperature fluctuations should be anti-correlated, which happens to be supported by compressible experiments, does not hold in the simulation. Instead, a relationship is found between the rates of turbulent heat and momentum transfer. The turbulent kinetic energy budget is computed and compared with the budgets from Spalart's incompressible simulations.

In premixed turbulent combustion, the modelling of the turbulent flux of the mean reaction progress variable remains somewhat controversial. Classical gradient transport assumptions based on the eddy viscosity concept are often used while both experimental data and theoretical analysis have pointed out the existence of countergradient turbulent diffusion. Direct numerical simulation (DNS) is used in this paper to provide basic information on the turbulent flux of and study the occurrence of counter-gradient transport. The numerical configuration corresponds to twoor three-dimensional premixed flames in isotropic turbulent flow. The simulations correspond to various flame and flow conditions that are representative of flamelet combustion. They reveal that different flames will feature different turbulent transport properties and that these differences can be related to basic dynamical differences in the flame-flow interactions: counter-gradient diffusion occurs when the flow field near the flame is dominated by thermal dilatation due to chemical reaction, whereas gradient diffusion occurs when the flow field near the flame is dominated by the turbulent motions. The DNS-based analysis leads to a simple expression to describe the turbulent flux of , which in turn leads to a simple criterion to delineate between the gradient and counter-gradient turbulent diffusion regimes. This criterion suggests that the occurrence of one regime or the other is determined primarily by the ratio of turbulence intensity divided by the laminar flame speed, and by the flame heat release factor, τ ≡ (Tb — Tu)/Tu, where Tu and Tb are respectively the temperature within unburnt and burnt gas. Consistent with the Bray-Moss-Libby theory, counter-gradient (gradient) diffusion is promoted by low (high) values and high (low) values of τ. DNS also shows that these results are not restricted to the turbulent transport of . Similar results are found for the turbulent transport of flame surface density, Σ. The turbulent fluxes of and Σ are strongly correlated in the simulated flames and counter-gradient (gradient) diffusion of always coincides with counter-gradient (gradient) diffusion of Σ.

There is an upper bound to the amount of power that can be generated by turbines in tidal channels as too many turbines merely block the flow. One condition for achievement of the upper bound is that the turbines are deployed uniformly across the channel, with all the flow through them, but this may interfere with other uses of the channel. An isolated turbine is more effective in a channel than in an unbounded flow, but the current downstream is non-uniform between the wake of the turbines and the free stream. Hence some energy is lost when these streams merge, as may occur in a long channel. We show here, for ideal turbine models, that the fractional power loss increases from 1/3 to 2/3 as the fraction of the channel cross-section spanned by the turbines increases from 0 to close to 1. In another scenario, possibly appropriate for a short channel, the speed of the free stream outside the turbine wake is controlled by separation at the channel exit. In this case, the maximum power obtainable is slightly less than proportional to the fraction of the channel cross-section occupied by turbines.

The problem of swimming at low Reynolds number is formulated in terms of a gauge field on the space of shapes. Effective methods for computing this field, by solving a linear boundary-value problem, are described. We employ conformal-mapping techniques to calculate swimming motions for cylinders with a variety of crosssections. We also determine the net translationl motion due to arbitrary infinitesimal deformations of a sphere.

The turbulence generated by a vertically oscillating grid in a water tank and the entrainment across a salinity interface caused by this turbulence have been investigated experimentally. Measurements were carried out in a homogeneous layer of fluid as well as a two-layered fluid, which permitted us to determine the decay law of this turbulence and the way in which the structure of the turbulence depends on the mesh size and on the frequency and amplitude of the grid oscillation. It was found that the turbulent kinetic energy decays with distance from the grid according to a power law $\overline{q^2}\propto z^{-n}$, with n close to 2, and that the turbulent Reynolds number remains approximately constant during decay. The linear dependence of the r.m.s. turbulent velocity on the grid oscillation frequency found by Thompson & Turner (1975) in the case of a square-bar grid has been confirmed. It is shown here that this linear relation remains valid when an interface is present and consequently the dependence of the entrainment velocity on the local Richardson number is of the form $u_e/u \propto Ri^{-\frac{3}{2}}$, the Péclet number being high. While the bearing of these results on the problem of the thermocline or an inversion is clear we wish to emphasize that the spatial decay of turbulence is interesting in itself.

The present work consists of two parts. Here in Part 1, a scaling law (incomplete similarity with respect to local Reynolds number based on distance from the wall) is proposed for the mean velocity distribution in developed turbulent shear flow. The proposed scaling law involves a special dependence of the power exponent and multiplicative factor on the flow Reynolds number. It emerges that the universal logarithmic law is closely related to the envelope of a family of power-type curves, each corresponding to a fixed Reynolds number. A skin-friction law, corresponding to the proposed scaling law for the mean velocity distribution, is derived.

In Part 2 (Barenblatt & Prostokishin 1993), both the scaling law for the velocity distribution and the corresponding friction law are compared with experimental data.

This paper describes an experimental programme carried out in a laboratory channel with rough and smooth beds, to investigate the interaction between gravity waves and a turbulent current. In particular, changes induced in the mean-velocity profiles, turbulent fluctuations, bed shear stresses and wave attenuation rates are considered for a range of wave heights, keeping the wave period constant. The smooth-boundary tests were carried out as a necessary preliminary to the more-realistic rough-boundary condition.

A directionally sensitive laser anemometer was used to measure horizontal, vertical, and 45° velocity components in the oscillating fluid, and an on-line minicomputer was programmed to produce ensemble averages of velocities, Reynolds stresses and wave-elevation data. The cycle was sampled at 200 separate phase positions, with 180 observations at each position. Measurements were made at up to 30 points in the vertical.

Preliminary tests were carried out on the unidirectional current and on the waves alone. These show that mean-velocity profiles and turbulence parameters of the current agree satisfactorily with previous experiments, and that the waves are approximated closely by Stokes’ second-order theory.

For combined wave and current tests, mean-velocity profiles are generally found to differ from those suggested by a linear superposition of wave and current velocities, a change in boundary-layer thickness being indicated. However, shear stresses at the smooth boundary are found to be described by such a linear addition.

The well-known analogy between the Euler equations for steady flow of an inviscid incompressible fluid and the equations of magnetostatic equilibrium in a perfectly conducting fluid is exploited in a discussion of the existence and structure of solutions to both problems that have arbitrarily prescribed topology. A method of magnetic relaxation which conserves the magnetic-field topology is used to demonstrate the existence of magnetostatic equilibria in a domain [Dscr ] that are topologically accessible from a given field B0(x) and hence the existence of analogous steady Euler flows. The magnetostatic equilibria generally contain tangential discontinuities (i.e. current sheets) distributed in some way in the domain, even although the initial field B0(x) may be infinitely differentiable, and particular attention is paid to the manner in which these current sheets can arise. The corresponding Euler flow contains vortex sheets which must be located on streamsurfaces in regions where such surfaces exist. The magnetostatic equilibria are in general stable, and the analogous Euler flows are (probably) in general unstable.

The structure of these unstable Euler flows (regarded as fixed points in the function space in which solutions of the unsteady Euler equations evolve) may have some bearing on the problem of the spatial structure of turbulent flow. It is shown that the Euler flow contains blobs of maximal helicity (positive or negative) which may be interpreted as ‘coherent structures’, separated by regular surfaces on which vortex sheets, the site of strong viscous dissipation, may be located.

In the first part of the paper, a mixing layer of tanh y form is considered, and twodimensional solutions of the non-linear inviscid equations are found representing periodic perturbations from the neutral wave of linearized stability theory. To second order in amplitude the solutions are equivalent to the equilibrium state calculated by Schade (1964), who studied the development of perturbations in time and found an evolution towards that equilibrium state. The present calculation, however, is taken to fourth-order in amplitude. It is noted that the solutions presented in this paper are regular, even though viscosity is ignored; and the relationships to the singular (if inviscid) time-dependent solutions of Schade are explained. Such regular, inviscid solutions have been found only for odd velocity profiles, such as tanh y.

Although the details of the velocity distributions are not of the form observed experimentally, it is shown that the amplitude ratios of fundamental and first harmonic, for a given absolute amplitude, are comparable with those observed.

In part 2 some exact non-linear solutions are presented of the inviscid, incompressible equations of fluid flow in two or three spatial dimensions. They illustrate the flows of part 1, since they are periodic in one co-ordinate (x), have a shear in another (y) and are independent of the third. Included, as two-dimensional cases, are (i) the tanh y velocity distribution for a flow wholly in the x-direction, (ii) the well-known solution for the flow due to a set of point vortices equi-spaced on the axis, and (iii) an example of linearized hydrodynamic (Orr-Sommerfeld) stability theory. The flows may involve concentrations of vorticity. In three-dimensional cases the z component of velocity is even in y, whereas the x component is odd. Consequently, the class of flows represents, in general, small or large periodic perturbations from a skewed shear layer. Time-dependent solutions, representing waves travelling in the x direction may be obtained by translation of axes.

The motion of free surfaces in incompressible, irrotational, inviscid layered flows is studied by evolution equations for the position of the free surfaces and appropriate dipole (vortex) and source strengths. The resulting Fredholm integral equations of the second kind may be solved efficiently in both storage and work by iteration in both two and three dimensions. Applications to breaking water waves over finite-bottom topography and interacting triads of surface and interfacial waves are given.