A comparison between the experimental visualization and numerical simulations of the occurrence of vortex breakdown in laminar swirling flows produced by a rotating endwall is presented. The experimental visualizations of Escudier (1984) were the first to detect the presence of multiple recirculation zones and the numerical model presented here, consisting of a numerical solution of the unsteady axisymmetric Navier-Stokes equations, faithfully reproduces these phenomena and all other observed characteristics of the flow. Further, the numerical calculations elucidate the onset of oscillatory flow, an aspect of the flow that was not clearly resolved by the flow visualization experiments. Part 2 of the paper examines the underlying physics of these vortex flows.

The physical mechanisms for vortex breakdown which, it is proposed here, rely on the production of a negative azimuthal component of vorticity, are elucidated with the aid of a simple, steady, inviscid, axisymmetric equation of motion. Most studies of vortex breakdown use as a starting point an equation for the azimuthal vorticity (Squire 1960), but a departure in the present study is that it is explored directly and not through perturbations of an initial stream function. The inviscid equation of motion that is derived leads to a criterion for vortex breakdown based on the generation of negative azimuthal vorticity on some stream surfaces. Inviscid predictions are tested against results from numerical calculations of the Navier-Stokes equations for which breakdown occurs.

An X-ray attenuation technique is used to obtain the local concentration of spherical particles in a polydisperse suspension as a function of vertical position and time. From these experimental data, the average velocity of sedimentation in the homogeneous part of the suspension is derived by considering the variation with time of the total volume of particles located above a given fixed horizontal plane. Measurements have been performed in suspensions of particles which differ from each other in size with a total volume concentration in particles between 0.13% and 2.5%, and also in suspensions of particles which differ from each other both in size and in density, the total volume concentration being 2%. For the first kind of suspension, the experimental hindered settling factor is plotted versus the concentration and a linear regression analysis provides the slope with its 90% confidence limits: Se = −5.3 ± 1.1. This experimental average coefficient of sedimentation is in good agreement with the theoretical average coefficient St = −5.60 obtained from the results of Batchelor & Wen (1982). The second kind of suspension, for which permanent doublets of spheres may theoretically exist, is not in the range of validity of Batchelor & Wen's results. The experimental average coefficient of sedimentation for this case is found to be much larger than the prediction obtained by extrapolating Batchelor & Wen's results out of their range of validity. This increased velocity may be experimental evidence of the existence of permanent doublets.

We study the relationship between dynamical structure and shape for vortex pairs, now usually named ‘modons’. When the boundary between the exterior irrotational flow and the inner core of non-zero vorticity is a circle, an analytical solution is known. Here, we generalize the circular modons to solitary vortex pairs whose vorticity boundary is an ellipse. We find that as the eccentricity of the ellipse increases, the vorticity becomes concentrated in narrow ridges which run just inside the elliptical vorticity boundary and continue just inside the line of zero vorticity which divides the two vortices. Each vortex becomes increasingly ‘hollow’ in the sense that each contains a broad valley of low vorticity which is completely enclosed by the ridge of high vorticity already described. The relationship between vorticity ζ and streak function Ψ, which is linear for the circular modons, becomes strongly nonlinear for highly eccentric modons, qualitatively resembling ζ ∝ Ψe−λΨ for some constant λ. In this study, we neglect the Earth's rotation, but our method is directly applicable to quasi-geostrophic modons, too. An efficient and simple spectral method for modon problems is provided.

An experimental investigation of the sedimentation of monodisperse colloidal silica spheres with grafted octadecyl chains with three different interaction potentials is presented. Small particles (0.27 μm) behaved as hard spheres in cyclohexane, but larger ones (0.60 and 0.94 μm) are weakly flocculated by van der Waals attractions. The smallest particles (0.08 μm) in hexadecane are strongly flocculated by attractions between the octadecyl layers. A medical computer tomography (CT) scanner provided an accurate and absolute density measurement without disrupting the process. For the hard spheres and the weakly flocculated systems, the kinetics of sedimentation for the dispersed phase could readily be predicted utilizing the flux curve. For flocculated networks, we found a power-law relationship between compressive yield stresses and solids fractions comparable with other experimental systems.

The characteristics of near-wall turbulence are examined and the result is used to assess the behaviour of the various terms in the Reynolds-stress transport equations. It is found that all components of the velocity-pressure-gradient correlation vanish at the wall. Conventional splitting of this second-order tensor into a pressure diffusion part and a pressure redistribution part and subsequent neglect of the pressure diffusion term in the modelled Reynolds-stress equations leads to finite near-wall values for two components of the redistribution tensor. This, therefore, suggests that, in near-wall turbulent flow modelling, the velocity-pressure-gradient correlation rather than pressure redistribution should be modelled. Based on this understanding, a methodology to derive an asymptotically correct model for the velocity-pressure-gradient correlation is proposed. A model that has the property of approaching the high-Reynolds-number model for pressure redistribution far away from the wall is derived. A similar analysis is carried out on the viscous dissipation term and asymptotically correct near-wall modifications are proposed. The near-wall closure based on the Reynolds-stress equations and a conventional low-Reynolds-number dissipation-rate equation is used to calculate fully-developed turbulent channel and pipe flows at different Reynolds numbers. A careful parametric study of the model constants introduced by the near-wall closure reveals that one constant in the dissipation-rate equation is Reynolds-number dependent, and a preliminary expression is proposed for this constant. With this modification, excellent agreement with near-wall turbulence statistics, measured and simulated, is obtained, especially the anisotropic behaviour of the normal stresses. On the other hand, it is found that the dissipation-rate equation has a significant effect on the calculated Reynolds-stress budgets. Possible improvements could be obtained by using available direct simulation data to help formulate a more realistic dissipation-rate equation. When such an equation is available, the present approach can again be used to derive a near-wall closure for the Reynolds-stress equations. The resultant closure could give improved predictions of the turbulence statistics and the Reynolds-stress budgets.

The problem is solved of the flow due to a rapidly oscillating source situated near a plane surface with a hemispherical indentation. The strength of the far-field flow is evaluated, and also the natural frequency of oscillation of a small bubble centred at the same position as the source.

Cocke (1969) showed that, on average, infinitesimal material lines and surfaces are stretched in incompressible isotropic turbulence. We have extended those results to obtain upper and lower bounds for the stretching of such infinitesimal elements in terms of the eigenvalues of the Green deformation tensor. These bounds are in turn used to find bounds for the stretching of finite material lines and surfaces.