One basic effect of turbulence in turbulent premixed combustion is for the fluctuating velocity field to wrinkle the flame and greatly increase its surface area. In the flamelet theory, this effect is described by the flame surface density. An exact evolution equation for the flame surface density, called the Σ-equation, may be written, where basic physical mechanisms like production by hydrodynamic straining and destruction by propagation effects are described explicitly. Direct numerical simulation (DNS) is used in this paper to estimate the different terms appearing in the Σ-equation. The numerical configuration corresponds to three-dimensional premixed flames in isotropic turbulent flow. The simulations are performed for various mixture Lewis numbers in order to modify the strength and nature of the flame-flow coupling. The DNS-based analysis provides much information relevant to flamelet models. In particular, the flame surface density, and the source and sink terms for the flame surface density, are resolved spatially across the turbulent flame brush. The geometry as well as the dynamics of the flame differ quite significantly from one end of the reaction zone to the other. For instance, contrary to the intuitive idea that flame propagation effects merely counteract the wrinkling due to the turbulence, the role of flame propagation is not constant across the turbulent brush and switches from flame surface production at the front to flame surface dissipation at the back. Direct comparisons with flamelet models are also performed. The Bray-Moss-Libby assumption that the flame surface density is proportional to the flamelet crossing frequency, a quantity that can be measured in experiments, is found to be valid. Major uncertainties remain, however, over an appropriate description of the flamelet crossing frequency. In comparison, the coherent flame model of Marble & Broadwell achieves closure at the level of the Σ-equation and provides a more promising physically based description of the flame surface dynamics. Some areas where the model needs improvement are identified.

The interaction between internal waves (IW) and wind waves (WW) is studied. Three types of interaction are considered: spontaneous IW generation by a random field of WWs, and two feedback mechanisms - modulation and friction.

The latter mechanism has not been studied before. Its influence on the IW-WW coupling is of primary importance. The modulation and friction mechanisms result in exponential attenuation of the IWs. Attenuation of IWs propagating against wind is the strongest. The IW attenuation has a dimensionless decrement of order 10-3, whereas for storm winds it attains the value of 10-2. Joint action of the spontaneous generation of IWs and their attenuation due to feedback mechanisms permits a stationary ‘wind-IW’ spectrum to exist. For strong winds the ‘wind-IW’ energy is of order 105 erg cm-2. The effect of IWs on currents in the ocean's upper layer is considered. Momentum and energy lost by IWs due to their interaction with WWs generates inertial oscillations. Under the attenuation of intensive IWs, the amplitude of inertial oscillations may be compared with the background Ekman current.

It is a significant feature of most gas-fluidized beds that they contain rising ‘bubbles’ of almost clear gas. The purpose of this paper is to account plausibly for this remarkable property first by supposing that primary and secondary instabilities of the fluidized bed generate compact regions of above-average or below-average particle concentration, and second by invoking a mechanism for the expulsion of particles from a buoyant compact blob of smaller particle concentration. We postulate that the rising of such an incipient bubble generates a toroidal circulation of the gas in the bubble, roughly like that in a drop of liquid rising through a second liquid of larger density, and that particles in the blob carried round by the fluid move on trajectories which ultimately cross the bubble boundary. Numerical calculations of particle trajectories for practical values of the relevant parameters show that a large percentage of particles, of such small concentration that they move independently, are expelled from a bubble in the time taken by it to rise through a distance of several bubble diameters.

Similar calculations for a liquid-fluidized bed show that the expulsion mechanism is much weaker, as a consequence of the larger density and viscosity of a liquid, which is consistent with the absence of observations of relatively empty bubbles in liquid-fluidized beds.

It is found to be possible, with the help of the Richardson-Zaki correlation, to adjust the results of these calculations so as to allow approximately for the effect of interaction of particles in a bubble in either a gas- or a liquid-fluidized bed. The interaction of particles at volume fractions of 20 or 30 % lengthens the expulsion times, although without changing the qualitative conclusions.

In this paper we calculate how a pendant drop evolves at the end of a nozzle when the volume of the drop increases steadily with time. We find that the character of the evolution is strongly dependent on the growth rate of the drop and the radius of the nozzle. Typically we find that once the drop has become unstable, two bifurcations occur shortly after each other when the growth rate of the drop is slow. For large growth rates the bifurcations are well-separated in time. We are able to calculate the volumes of the drops after the bifurcations. A comparison with experimental data shows a satisfactory agreement.

We report the results of an investigation of the weakly nonlinear evolution of a triad of waves, each slightly amplified on a linear basis, that are superimposed on a tanh y mixing layer. The triad consists of a plane wave and a pair of oblique modes that act as a subharmonic of order 1/2. The oblique modes are inclined at approximately ±60°. to the mean flow direction and because the resonance conditions are satisfied exactly the analysis is entirely self-consistent as an asymptotic theory. The nonlinearity first occurs within the critical layer and the initial interaction is of the parametric resonance type. This produces faster than exponential growth of the oblique waves, behaviour observed recently in the experiments of Corke & Kusek (1993). The critical-layer dynamics lead subsequently to coupled integro-differential equations governing the amplitude evolution and, as first shown in related work by Goldstein & Lee (1992) on boundary layers in an adverse pressure gradient, these equations develop singularities in a finite time.

Parametric excitation of surface waves via forced vertical oscillation of a container filled with fluid (the Faraday instability) is investigated experimentally in a small-depth large-aspect-ratio system, with a viscous fluid and with two simultaneous forcing frequencies. The asymptotic pattern observed just above the threshold for the first instability of the flat surface is found to depend strongly on the frequency ratio and the amplitudes and phases of the two sinusoidal components of the driving acceleration. Parallel lines, squares, and hexagons are observed. With viscosity 100 cS, these stable standing-wave patterns do not exhibit strong sidewall effects, and are found in containers of various shapes including an irregular shape. A ‘quasi-pattern’ of twelvefold symmetry, analogous to a two-dimensional quasi-crystal, is observed for some even/odd frequency ratios. Many of the experimental phenomena can be modelled via cubic-order amplitude equations derived from symmetry arguments.

The oscillatory motion of an electrically charged non-spherical colloidal particle in an oscillating electric field is investigated. The particle is immersed in an incompressible viscous fluid and assumed to have a thin electric double layer. For moderate-aspect-ration spheroids and cylinders, a simple algebraic expression is derived that accurately describes oscillatory electrophoretic particle motion in terms of the steady Stokes resistance, added mass, and Basset force. The effects of double-layer conduction and displacement currents within dielectric particles are included. The results indicate that electroacoustic measurements may be able to determine the ζ-potential, dielectric constant, surface conductivity (and microstructural information contained therein), size, density, volume fraction, and possibly shape of non-spherical particles in a dilute suspension. A simple formula is obtained for the high-frequency electrical conductivity of a dilute suspension of colloidal spheroids with arbitrary charge and dielectric constant; only the added mass and Basset force are required and the requisite parameters are given. The result is needed for electroacoustic measurements but it may also be independently useful for determining the dielectric constant, surface conductivity, volume fraction, and possibly the shape of non-spherical particles in a dilute suspension. Electroacoustic energy dissipation is described for a dilute colloidal suspension. It is shown that resistive electrical heating and viscous dissipation occur independently. Electrical and viscous dissipation coefficients that characterize the order volume fraction contributions of the suspended particles are calculated; the electrical dissipation coefficient is O(1) for all oscillation frequencies, whereas the latter vanishes at low- and high-frequencies. The fluid motion is shown to be a superposition of unsteady, viscous and potential flows past an oscillating particle with no applied electric field. The electro-osmotic flow field is insensitive to particle geometry and qualitatively different from the flow past an oscillating particle with no applied field.

Rotation strongly affects the stability of turbulent flows in the presence of large eddies. In this paper, we examine the applicability of the classic Bradshaw-Richardson criterion to flows more general than a simple combination of rotation and pure shear. Two approaches are used. Firstly the linearized theory is applied to a class of rotating two-dimensional flows having arbitrary rates of strain and vorticity and streamfunctions that are quadratic. This class includes simple shear and elliptic flows as special cases. Secondly, we describe a large-eddy simulation of initially quasi-homogeneous three-dimensional turbulence superimposed on a periodic array of two-dimensional Taylor-Green vortices in a rotating frame.

The results of both approaches indicate that, for a large structure of vorticity W and subject to rotation Ω, maximum destabilization is obtained for zero tilting vorticity (½W + 2Ω = 0) whereas stability occurs for zero absolute vorticity (2Ω = 0) These results are consistent with the Bradshaw-Richardson criterion; however the numerical results show that in other cases the Bradshaw-Richardson number $B=2\Omega(W+2\Omega)/W^2$ is not always a good indicator of the flow stability.

We find that the RMS wave height (square root of the potential energy) rather than peak-to-peak wave height is a better experimental and analytic criterion for determining when a regular, two-dimensional deep-water wave will break. A spectral algorithm for two-dimensional potential flow is developed and used to compare breaking onset criteria for energy input from (i) converging sidewalls, (ii) a submerged disturbance, and (iii) wave focusing. We also find that wave-breaking criteria (potential energy or the more classical peak-to-peak wave height) are a function of the rate of energy input. Large plunging waves occur when energy input rates are large. As energy input rates become smaller there is a smooth transition to smaller spilling waves. The various energy input methods show similar breaking trends in the limit as the energy input rate becomes small - waves break when the potential energy becomes approximately 52 % of the energy for the most energetic Stokes wave, with the formation of a singularity immediately before the crest. The effects of wave modulation and reflection are briefly discussed and shown not to affect the potential energy breaking criterion significantly. The experimental scatter of the RMS wave height is shown to be half that of wave steepness during incipient breaking in wave packets.

The scattering of water waves by a varying bottom topography is considered using two-dimensional linear water-wave theory. A new approach is adopted in which the problem is first transformed into a uniform strip resulting in a variable free-surface boundary condition. This is then approximated by a finite number of sections on which the free-surface boundary condition is assumed to be constant. A transition matrix theory is developed which is used to relate the wave amplitudes at ±∞. The method is checked against examples for which the solution is known, or which can be computed by alternative means. Results show that the method provides a simple accurate technique for scattering problems of this type.

The boundary-layer flow over a cooled horizontal plate is considered. It is shown that the real part of the spectrum of the evolution operator of the linearized equations is not bounded uniformly from above which explains the difficulties encounterd by a numerical solution. Furthermore it is shown that near the leading edge an asymptotic expansion of the solution is not unique. A one-parametric family of asymptotic expansions of solutions can be constructed.

An exact result is calculated numerically for the dilute limiting, zero shear viscosity of bimodal suspensions of hard spheres. The required hydrodynamic functions are calculated from recent results for the resistivities of unequal spheres. Both the hydrodynamic and Brownian contributions to the Huggins coefficient exhibit a minimum that is symmetric in mixing volume fraction. The resultant minimum deepens with increasing size ratio. The results are discussed in the light of published measurements of the viscosity for bimodal suspensions and previous phenomenological theories. The reduction of viscosity upon mixing is seen to be a result of near-field hydrodynamic shielding of asymmetric particle pairs. It is also shown that the use of far-field hydrodynamic interactions yields qualitatively incorrect results for the viscosity of binary mixtures. A parametrization of the bimodal results allows an estimation of the effects of suspension polydispersity on the Huggins coefficient. For polydispersities of ten percent or less, the Huggins coefficient is essentially unchanged from the value calculated for an equivalent, monodisperse suspension at equal volume fraction. A parametrization of these results is provided for relating the reduction in Huggins coefficient to the polydispersity index.

By means of direct numerical solution of the kinetic equation for surface gravity waves, it is shown that under certain conditions the constant flux spectra of nonlinear waves, first predicted by Zakharov & Filonenko (1966) for an infinite frequency domain, can be formed in a finite frequency interval. For the case of angular isotropic spectra the conditions and timescales of this flux spectra formation are evaluated.

A linear-stability analysis is performed on a vertical falling film with a surface-active solute. It is assumed in the present model that the surfactant is soluble and volatile. In addition to the surface wave mode and the ‘wall wave’ mode which originate from the gravity-driven flow of the falling film itself, a new mode of instability related to the Marangoni effect induced by surface tension gradients is found for low Reynolds numbers and for moderate- or short-wavelength disturbances. The new mode is thought to be analogous to the thermocapillary instability examined first by Pearson (1958). The Marangoni instability of large-wavelength disturbances, revealed by Goussis & Kelly (1990) in a study of a liquid layer heated from below, may be completely suppressed in the present system by the effect of surface-excess concentration of the surfactant. The influence of the desorption of the solute and of its adsorption at the gas-liquid interface is determined for both the surface wave mode and the new wave mode. Desorption of the surfactant is shown to be responsible for the Marangoni instability of the new mode.

The motion of two equal spherical bubbles moving along their line of centres in a viscous liquid is studied numerically in bispherical coordinates. The unsteady Navier-Stokes equations are solved using a mixed spectral/finite-difference scheme for Reynolds numbers up to 200. Free-slip conditions at the bubble surfaces are imposed, while the normal stress condition is replaced by the sphericity constraint under the assumption of small Weber number. The vorticity shed by the upstream bubble affects the drag on the trailing bubble in a very complex fashion that appears to be quite beyond the power of existing asymptotic analyses. The separation between two equal bubbles rising in line under the action of buoyancy is predicted to reach an equilibrium value dependent on the Reynolds number. This result is at variance with experiment. The explanation offered of this difference casts further doubt on the feasibility of a simplified simulation of bubbly liquid dynamics.

The turbulent stress tensor in large-eddy simulation is examined from a theoretical point of view. Realizability conditions for the components of this tensor are derived, which hold if and only if the filter function is positive. The spectral cut-off, one of the filters frequently used in large-eddy simulation, is not positive. Consequently, the turbulent stress tensor based on spectrally filtered fields does not satisfy the realizability conditions, which leads to negative values of the generalized turbulent kinetic energy k. Positive filters, e. g. Gaussian or top-hat, always give rise to a positive k. For this reason, subgrid models which require positive values for k should be used in conjunction with e. g. the Gaussian or top-hat filter rather than with the spectral cutoff filter. If the turbulent stress tensor satisfies the realizability conditions, it is natural to require that the subgrid model for this tensor also satisfies these conditions. With respect to this point of view several subgrid models are discussed. For eddy-viscosity models a lower bound for the generalized turbulent kinetic energy follows as a necessary condition. This result provides an inequality for the model constants appearing in a ‘Smagorinsky-type’ subgrid model for compressible flows.

General representations are derived for both the velocity potential and the surface pressure fluctuations induced by an arbitrary distribution of vorticity near a manoeuvring cylinder. The cylinder is inextensible and in unsteady motion. Its axis may be slightly curved, with radius of curvature large in comparison with the cylinder radius.

Two model problems are considered in detail to investigate the effect of lateral displacements of a cylinder with an established boundary layer. The boundary layer on the flexible cylinder is found to be shed once the lateral displacement of the cylinder axis exceeds the boundary-layer thickness. The unsteady pressures generated by this vortex shedding are investigated.

Based on the linearized long-wave equation, two new analytical solutions are obtained respectively for the propagation of long surface gravity waves around a conical island and over a paraboloidal shoal. Having been intensively studied during the last two decades, these two problems have practical significance and are physically revealing for wave propagation over variable water depth. The newly derived analytical solutions are compared with several previously obtained numerical solutions and the accuracy of those numerical solutions is discussed. The analytical method has the potential to be used to find solutions for wave propagation over more natural bottom topographies.