L. G. Leal

Inertial migration of rigid spheres in two-dimensional unidirectional flows

The familiar Segre-Silberberg effect of inertia-induced lateral migration of a neutrally buoyant rigid sphere in a Newtonian fluid is studied theoretically for simple shear flow and for two-dimensional Poiseuille flow. It is shown that the spheres reach a stable lateral equilibrium position independent of the initial position of release. For simple shear flow, this position is midway between the walls, whereas for Poiseuille flow, it is 0·6 of the channel half-width from the centre-line. Particle trajectories are calculated in both cases and compared with available experimental data. Implications for the measurement of the rheological properties of a dilute suspension of spheres are discussed.

An experimental investigation of drop deformation and breakup in steady, two-dimensional linear flows

We consider the deformation and burst of small fluid droplets in steady linear, two-dimensional motions of a second immiscible fluid. Experiments using a computer-controlled, four-roll mill to investigate the effect of flow type are described, and the results compared with predictions of several available asymptotic deformation and burst theories, as well as numerical calculations. The comparison clarifies the range of validity of the theories, and demonstrates that they provide quite adequate predictions over a wide range of viscosity ratio, capillary number, and flow type.

The effect of Brownian motion on the rheological properties of a suspension of non-spherical particles

The effect of rotary Brownian motion on the rheology of a dilute suspension of rigid spheroids in shear flow is considered for various limiting cases of the particle aspect ratio r and dimensionless shear rate γ/ D . As a preliminary the probability distribution function is calculated for the orientation of a single, axisymmetric particle in steady shear flow, assuming small particle Reynolds number. The result for the case of weak-shear flows, γ/ D [Lt ] 1, has been known for many years. After briefly reviewing this limiting case, we present expressions for the case of strong shear where ( r 3 + r −3 ) [Lt ] γ/ D , and for an intermediate regime relevant for extreme aspect ratios where 1 [Lt ] γ/ D [Lt ] ( r 3 + r −3 ). The bulk stress is then calculated for these cases, as well as the case of nearly spherical particles r ∼ 1, which has not hitherto been discussed in detail. Various non-Newtonian features of the suspension rheology are discussed in terms of prior continuum mechanical and experimental results.

Numerical solution of free-boundary problems in fluid mechanics. Part 1. The finite-difference technique

We present here a brief description of a numerical technique suitable for solving axisymmetric (or two-dimensional) free-boundary problems of fluid mechanics. The technique is based on a finite-difference solution of the equations of motion on an orthogonal curvilinear coordinate system, which is also constructed numerically and always adjusted so as to fit the current boundary shape. The overall solution is achieved via a global iterative process, with the condition of balance between total normal stress and the capillary pressure at the free boundary being used to drive the boundary shape to its ultimate equilibrium position.

Relaxation and breakup of an initially extended drop in an otherwise quiescent fluid

In this paper we examine some general features of the time-dependent dynamics of drop deformation and breakup at low Reynolds numbers. The first aspect of our study is a detailed numerical investigation of the ‘end-pinching’ behaviour reported in a previous experimental study. The numerics illustrate the effects of viscosity ratio and initial drop shape on the relaxation and/or breakup of highly elongated droplets in an otherwise quiescent fluid. In addition, the numerical procedure is used to study the simultaneous development of capillary-wave instabilities at the fluid-fluid interface of a very long, cylindrically shaped droplet with bulbous ends. Initially small disturbances evolve to finite amplitude and produce very regular drop breakup. The formation of satellite droplets, a nonlinear phenomenon, is also observed.

Constitutive equations in suspension mechanics. Part 2. Approximate forms for a suspension of rigid particles affected by Brownian rotations

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

An experimental study of transient effects in the breakup of viscous drops

A computer-controlled four-roll mill is used to examine two transient modes of deformation of a liquid drop: elongation in a steady flow and interfacial-tension-driven motion which occurs after the flow is stopped abruptly. For modest extensions, drop breakup does not occur with the flow on, but may occur following cessation of the flow as a result of deterministic motions associated with internal pressure gradients established by capillary forces. These relaxation and breakup phenomena depend on the initial drop shape and the relative viscosities of the two fluids. Capillary-wave instabilities at the fluid-fluid interface are observed only for highly elongated drops. This study is a natural extension of G. I. Taylors original studies of the deformation and burst of droplets in well-defined flow fields.

Numerical solution of free-boundary problems in fluid mechanics. Part 2. Buoyancy-driven motion of a gas bubble through a quiescent liquid

In this paper numerical results are presented for the buoyancy-driven rise of a deformable bubble through an unbounded quiescent fluid. Complete solutions, including the bubble shape, are obtained for Reynolds numbers in the range 1 [less-than-or-equal] R [less-than-or-equal] 200 and for Weber numbers up to 20. For Reynolds numbers R [less-than-or-equal] 20 the shape of the bubble changes from nearly spherical to oblate-ellipsoidal to spherical-cap depending on Weber number; at higher Reynolds numbers ‘disk-like’ and ‘saucer-like’ shapes appear at W = O(10). The present results show clearly that flow separation may occur at a smooth free surface at intermediate Reynolds numbers; this fact suggests a qualitative explanation of the often-observed irregular (zigzag or helical) paths of rising bubbles.

The effects of surfactants on drop deformation and breakup

The effects of surface-active agents on drop deformation and breakup in extensional flows at low Reynolds numbers are described. In this free-boundary problem, determination of the interfacial velocity requires knowledge of the distribution of surfactant, which, in turn, requires knowledge of the interfacial velocity field. We account for this explicit coupling of the unknown drop shape and the evolving surfactant distribution. An analytical result valid for nearly spherical distortions is presented first. Finite drop deformation is studied numerically using the boundaryintegral method in conjunction with the time-dependent convective-diffusion equation for surfactant transport. This procedure accurately follows interfacial tension variations, produced by non-uniform surfactant distribution, on the evolving interface. The numerical method allows for an arbitrary equation of state relating interfacial tension to the local concentration of surfactant, although calculations are presented only for the common linear equation of state. Also, only the case of insoluble surfactant is studied. The analytical and numerical results indicate that at low capillary numbers the presence of surfactant causes larger deformation than would occur for a drop with a constant interfacial tension equal to the initial equilibrium value. The increased deformation occurs owing to surfactant being swept to the end of the drop where it acts to locally lower the interfacial tension, which therefore requires increased deformation to satisfy the normal stress balance. However, at larger capillary numbers and finite deformations, this convective effect competes with ‘dilution ’ of the surfactant due to interfacial area increases. These two different effects of surfaceactive material are illustrated and discussed and their influence on the critical capillary number for breakup is presented.

Natural convection in a shallow cavity with differentially heated end walls. Part 1. Asymptotic theory

The problem of natural convection in a cavity of small aspect ratio with differentially heated end walls is considered. It is shown by use of matched asymptotic expansions that the flow consists of two distinct regimes : a parallel flow in the core region and a second, non-parallel flow near the ends of the cavity. A solution valid at all orders in the aspect ratio A is found for the core region, while the first several terms of the appropriate asymptotic expansion are obtained for the end regions. Parametric limits of validity for the parallel flow structure are discussed. Asymptotic expressions for the Nusselt number and the single free parameter of the parallel flow solution, valid in the limit as A → 0, are derived.