A two-dimensional analysis, based on linear surface-wave theory, is developed for an oscillating-water-column wave-energy device in water of arbitrary constant depth. The immersed part of the structure is assumed of shallow draught except for a submerged vertical reflecting wall. Both the cases of linear and nonlinear power take-off are considered. The results show that air compressibility can be important in practice, and its effects may in general be satisfactorily represented by linearization. The analysis indicates that using a turbine whose characteristic exhibits a phase difference between pressure and flow rate may be a method of strongly reducing the chamber length and turbine size with little change in the capability of energy extraction from regular waves. It was found in two examples of devices with strongly nonlinear power take-off that the maximum efficiency is only marginally inferior to what can be achieved in the linear case.

The modulation of isotropic turbulence by particles has been investigated using direct numerical simulation (DNS). The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Gravitational settling is also neglected and particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian time scale of the turbulence and particle mass loadings up to 1. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 963 collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier–Stokes equations using the point-force approximation in which 963 particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviours. The spectrum of the fluid–particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high-wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid–particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of Squires & Eaton (1990) and Elghobashi & Truesdell (1993).

This article deals with the first Hopf bifurcation of a cylinder flow, and more particularly with the properties of the unsteady periodic Kármán vortex street regime that sets in for supercritical Reynolds numbers Re > 46. Barkley (Europhys. Lett. vol.75, 2006, p. 750) has recently studied the linear properties of the associated mean flow, i.e. the flow which is obtained by a time average of this unsteady periodic flow. He observed, thanks to a global mode analysis, that the mean flow is marginally stable and that the eigenfrequencies associated with the global modes of the mean flow fit the Strouhal to Reynolds experimental function well in the range 46 < Re < 180. The aim of this article is to give a theoretical proof of this result near the bifurcation. For this, we do a global weakly nonlinear analysis valid in the vicinity of the critical Reynolds number Rec based on the small parameter ε = Rec−1 − Re−1 ≪ 1. We compute numerically the complex constants λ and μ′ which appear in the Stuart-Landau amplitude equation: dA/dt = ε λA − εμ′ A|A|2. Here A is the scalar complex amplitude of the critical global mode. By analysing carefully the nonlinear interactions yielding the term μ′, we show for the cylinder flow that the mean flow is approximately marginally stable and that the linear dynamics of the mean flow yields the frequency of the saturated Stuart-Landau limit cycle. We will finally show that these results are not general, by studying the case of the bifurcation of an open cavity flow. In particular, we show that the mean flow in this case remains strongly unstable and that the frequencies associated with the eigenmodes do not exactly match those of the nonlinear unsteady periodic cavity flow. It will be demonstrated that two precise conditions must hold for a linear stability analysis of a mean flow to be relevant and useful.

Experiments on transition to turbulence in a purely oscillatory pipe flow were performed for values of the Reynolds number Rδ, defined using the Stokes-layer thickness δ = (2ν/ω)½ and the cross-sectional mean velocity amplitude Û, from 19 to 1530 (or for values of the Reynolds number Re, defined using the pipe diameter d and Û, from 105 to 5830) and for values of the Stokes parameter λ = ½d(ω/2ν)½ (ν = kinematic viscosity and ω = angular frequency) from 1·35 to 6·19. Three types of turbulent flow regime have been detected: weakly turbulent flow, conditionally turbulent flow and fully turbulent flow. Demarcation of the flow regimes is possible on Rλ, λ or Re, λ diagrams. The critical Reynolds number of the first transition decreases as the Stokes parameter increases. In the conditionally turbulent flow, turbulence is generated suddenly in the decelerating phase and the profile of the velocity distribution changes drastically. In the accelerating phase, the flow recovers to laminar. This type of partially turbulent flow persists even at Reynolds numbers as high as Re = 5830 if the value of the Stokes parameter is high.

The effect of Mach number on the evolution of instabilities in the compressible mixing layer is investigated. The full time-dependent compressible Navier–Stokes equations are solved numerically for a temporally evolving mixing layer using a mixed spectral and high-order finite difference method. The convective Mach number Mc (the ratio of the velocity difference to the sum of the free-stream sound speeds) is used as the compressibility parameter. Simulations with random initial conditions confirm the prediction of linear stability theory that at high Mach numbers (Mc > 0.6) oblique waves grow more rapidly than two-dimensional waves. Simulations are then presented of the nonlinear temporal evolution of the most rapidly amplified linear instability waves. A change in the developed large-scale structure is observed as the Mach number is increased, with vortical regions oriented in a more oblique manner at the higher Mach numbers. At convective Mach numbers above unity the two-dimensional instability is found to have little effect on the flow development, which is dominated by the oblique instability waves. The nonlinear structure which develops from a pair of equal and opposite oblique instability waves is found to resemble a pair of inclined A-vortices which are staggered in the streamwise direction. A fully nonlinear computation with a random initial condition shows the development of large-scale structure similar to the simulations with forcing. It is concluded that there are strong compressibility effects on the structure of the mixing layer and that highly three-dimensional structures develop from the primary inflexional instability of the flow at high Mach numbers.

Upper bounds for the heat flux through a horizontally infinite layer of fluid heated from below are obtained by maximizing the heat flux subject to (a) two integral constraints, the ‘power integrals’, derived from the equations of motion, and (b) the continuity equations. This variational problem is solved completely, for all values of the Rayleigh number R, when only the constrains (a) are imposed, and it is thus shown that the Nusselt number N for any statistically steady convective motion cannot exceed a certain value N1(R), which for large R is approximately (3R/64)½. When (b) is included as a constraint, the variational problem is solved for large R, under the additional hypothesis that the solution has a single horizontal wave number; the associated upper bound on the Nusselt number is (R–248)⅜. The mean properties of this maximizing ‘flow’, in particular the mean temperature and mean square temperature deviation fields, are found to resemble the mean properties of the real flow observed by Townsend; the results thus tend to support Malkus's hypothesis that turbulent convection maximizes heat flux.

The relationship between the upstream boundary layer and the low-frequency, large-scale unsteadiness of the separated flow in a Mach 2 compression ramp interaction is investigated by performing wide-field particle image velocimetry (PIV) and planar laser scattering (PLS) measurements in streamwise–spanwise planes. Planar laser scattering measurements in the upstream boundary layer indicate the presence of spanwise strips of elongated regions of uniform momentum with lengths greater than 40δ. These long coherent structures have been observed in a Mach 2 supersonic boundary layer (Ganapathisubramani, Clemens & Dolling 2006) and they exhibit strong similarities to those that have been found in incompressible boundary layers (Tomkins & Adrian 2003; Ganapathisubramani, Longmire & Marusic 2003). At a wall-normal location of y/δ=0.2, the inferred instantaneous separation line of the separation region is found to oscillate between x/δ=−3 and −1 (where x/δ=0 is the ramp corner). The instantaneous spanwise separation line is found to respond to the elongated regions of uniform momentum. It is shown that high- and low-momentum regions are correlated with smaller and larger size of the separation region, respectively. Furthermore, the instantaneous separation line exhibits large-scale undulations that conform to the low- and high-speed regions in the upstream boundary layer. The low-frequency unsteadiness of the separation region/shock foot observed in numerous previous studies can be explained by a turbulent mechanism that includes these elongated regions of uniform momentum.

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.

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.

Previous measurements of the decay rate of the fluctuation intensity of passive scalars in grid-generated turbulence show large variation. New results presented here show that the decay rate of passive temperature fluctuations produced by heating the grid is a function of the initial temperature fluctuation intensity. Although a full reason for this is wanting, spectra of the temperature fluctuations show that, by varying the heat applied to the grid, the wavenumber of the maximum in the temperature spectrum changes, indicating that the geometry of the thermal fluctuations is being altered in some way. In these experiments the one-dimensional temperature spectrum shows an anomalous $-\frac{5}{3}$ slope. In order to eliminate the dependence of the decay rate of the temperature fluctuations on their intensity, we describe a new way of generating temperature fluctuations by means of placing a heated parallel array of fine wires (a mandoline) downstream from the unheated grid. Results of this experiment show that the decay rate of passive thermal fluctuations is uniquely determined by the wave-number of the initial temperature fluctuations. In this type of flow there appears to be no equilibrium value for the thermal fluctuation decay rate and hence for the mechanical/thermal time-scale ratio since the thermal fluctuation decay rate does not change within the tunnel length, which is the equivalent of nearly one turbulence decay time.

The stability of non-Newtonian jets was investigated. A linearized stability analysis shows that a liquid column of a viscoelastic fluid exhibits more rapid growth of axisymmetric wave disturbances than a Newtonian fluid of the same zero shear viscosity. This result is independent of the form of constitutive equation chosen. Experiments in weakly elastic fluids confirm this expectation, whereas data on fluids with more pronounced elastic properties indicate that non-linear phenomena are dominating. The disturbances appear as a series of droplets connected by random lengths of threads, which thin with distance and eventually lead to jet breakup. Even in dilute viscoelastic solutions, jet breakup does not occur by the growth of clearly defined waves.

Experimental and theoretical studies were undertaken to investigate the particulatefree secondary flow interaction in the wire-type electrostatic precipitator. The secondary flow generated by the corona discharge is not negligible, and strong flow interactions take place owing to the induced circulatory cells. The calculated numerical results demonstrate close agreement with experiment.

Weis-Fogh (1973) proposed a new mechanism of lift generation of fundamental interest. Surprisingly, it could work even in inviscid two-dimensional motions starting from rest, when Kelvin's theorem states that the total circulation round a body must vanish, but does not exclude the possibility that if the body breaks into two pieces then there may be equal and opposite circulations round them, each suitable for generating the lift required in the pieces’ subsequent motions! The ‘fling’ of two insect wings of chord c (figure 1) turning with angular velocity Ω generates irrotational motions associated with the sucking of air into the opening gap which are calculated in § 2 as involving circulations −0·69Ωc2 and + 0.69Ωc2 around the wings when their trailing edges, which are stagnation points of those irrotational motions, break apart (position (f)). Viscous modifications to this irrotational flow pattern by shedding of vorticity at the boundary generate (§ 3) a leading-edge separation bubble, and tend to increase slightly the total bound vorticity. Its role in a three-dimensional picture of the Weis-Fogh mechanism of lift generation, involving formation of trailing vortices at the wing tips, and including the case of a hovering insect like Encarsia formosa moving those tips in circular paths, is investigated in § 4. The paper ends with the comment that the far flow field of such very small hovering insects should take the form of the exact solution (Landau 1944; Squire 1951) of the Navier-Stokes equations for the effect of a concentrated force (the weight mg of the animal) acting on a fluid of kinematic viscosity v and density p, whenever the ratio mg/pv2 is small enough for that jet-type induced motion to be stable.

This paper describes the flow structure observed over a bell-shaped hill with height h (the profile of which is the reciprocal of a fourth-order polynomial) when it was placed first in a large towing tank containing stratified saline solutions with uniform stable density gradients and second in an unstratified wind tunnel. (A similarly shaped model hill was also studied in a small towing tank.) Observations were made at values of the Froude number F (≃ U/Nh) in the range 0·1 to 1·7 and at F = ∞, where U is the towing speed and N is the Brunt-Väisälä frequency, and at values of the Reynolds number from 400 to 275000. For F ≲ 0·4, the observations verify Drazin's (1961) theory for low-Froude-number flow over three-dimensional obstacles and establish limits of applicability. For Froude numbers of the order of unity, it is found that a classification of the lee-wave patterns and separated-flow regions observed in two-dimensional flows also appears to apply to three-dimensional hills.

Flow-visualization techniques were used extensively in obtaining both qualitative and quantitative information on the flow structure around the hill. Representative photographs of dye tracers, potassium permanganate dye streaks, shadowgraphs, surface dye smears, and hydrogen-bubble patterns are included here. While emphasis is centred on obtaining a basic understanding of the flow around three-dimensional hills, the results are applicable to the estimation of air pollutant dispersion around hills.

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.

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 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.

The suggestion made by the authors in a previous paper (Hurle & Jakeman 1969) that the Soret effect could give rise to overstable solutions of the thermosolutal Rayleigh–Jeffreys problem is investigated theoretically and experimentally.

Oscillatory instability is shown to occur in initially homogeneous layers of water-methanol mixtures when they are heated from below. This instability triggers a finite-amplitude steady mode. The magnitude and sign of the Soret coefficient was changed by varying the composition of the mixture; as predicted, overstable modes were observed when the sign of the coefficient was such as to produce a stabilizing contribution to the density gradient. The observed critical Rayleigh numbers and temporal frequencies are consistent with theory.

A chemical explosive mode analysis (CEMA) was developed as a new diagnostic to identify flame and ignition structure in complex flows. CEMA was then used to analyse the near-field structure of the stabilization region of a turbulent lifted hydrogen–air slot jet flame in a heated air coflow computed with three-dimensional direct numerical simulation. The simulation was performed with a detailed hydrogen–air mechanism and mixture-averaged transport properties at a jet Reynolds number of 11000 with over 900 million grid points. Explosive chemical modes and their characteristic time scales, as well as the species involved, were identified from the Jacobian matrix of the chemical source terms for species and temperature. An explosion index was defined for explosive modes, indicating the contribution of species and temperature in the explosion process. Radical and thermal runaway can consequently be distinguished. CEMA of the lifted flame shows the existence of two premixed flame fronts, which are difficult to detect with conventional methods. The upstream fork preceding the two flame fronts thereby identifies the stabilization point. A Damköhler number was defined based on the time scale of the chemical explosive mode and the local instantaneous scalar dissipation rate to highlight the role of auto-ignition in affecting the stabilization points in the lifted jet flame.

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.