ALBERT

Description: A boundary integral method is presented for analysing particle motion in a rotating fluid for flows where the Taylor number 2T is arbitrary and the Reynolds number is small. The method determines the surface traction and drag on a particle, and also the velocity field at any location in the fluid. Numerical results show that the dimensionless drag on a spherical particle translating along the rotation axis of an unbounded fluid is determined by the empirical formula D/6n = 1 +(4/7) ^”1/2 +(8/9TI)2T, which incorporates known results for the low and high Taylor number limits. Streamline portraits show that a critical Taylor number c « 50 exists at which the character of the flow changes. For 3 “ 〈 2Tcthe flow field appears as a perturbation of a Stokes flow with a superimposed swirling motion. For T 〈 2TCthe flow field develops two detached recirculating regions of trapped fluid located fore and aft of the particle. The recirculating regions grow in size and move farther from the particle with increasing Taylor number. This recirculation functions to deflect fluid away from the translating particle, thereby generating a columnar flow structure. The flow between the recirculating regions and the particle has a plug-like velocity profile, moving slightly slower than the particle and undergoing a uniform swirling motion. The flow in this region is matched to the particle velocity in a thin Ekman layer adjacent to the particle surface. A further study examines the translation of spheroidal particles. For large Taylor numbers, the drag is determined by the equatorial radius; details of the body shape are less important. © 1994, Cambridge University Press. All rights reserved.

Description: Time-dependent interactions between two buoyancy-driven deformable drops are studied in the low Reynolds number flow limit for sufficiently large Bond numbers that the drops become significantly deformed. The first part of this paper considers the interaction and deformation of drops in axisymmetric configurations. Boundary integral calculations are presented for Bond numbers ℬ = ∆ρga2/σ in the range 0.25 ≤ ℬ 〈 ∞ and viscosity ratios λ in the range 0.2 ∆ λ ∆ 20. Specifically, the case of a large drop following a smaller drop is considered, which typically leads to the smaller drop coating the larger drop for ℬ ≫ 1. Three distinct drainage modes of the thin film of fluid between the drops characterize axisymmetric two-drop interactions: (i) rapid drainage for which the thinnest region of the film is on the axis of symmetry, (ii) uniform drainage for which the film has a nearly constant thickness, and (iii) dimple formation. The initial mode of film drainage is always rapid drainage. As the separation distance decreases, film flow may change to uniform drainage and eventually to dimpled drainage. Moderate Bond numbers, typically ℬ = 0(10) for λ = 0(1), enhance dimple formation compared to either much larger or smaller Bond numbers. The numerical calculations also illustrate the extent to which lubrication theory and analytical solutions in bipolar coordinates (which assume spherical drop shapes) are applicable to deformable drops. The second part of this investigation considers the ‘stability’ of axisymmetric drop configurations. Laboratory experiments and two-dimensional boundary integral simulations are used to study the interactions between two horizontally offset drops. For sufficiently deformable drops, alignment occurs so that the small drop may still coat the large drop, whereas for large enough drop viscosities or high enough interfacial tension, the small drop will be swept around the larger drop. If the large drop is sufficiently deformable, the small drop may then be ‘sucked’ into the larger drop as it is being swept around the larger drop. In order to explain the alignment process, the shape and translation velocities of widely separated, nearly spherical drops are calculated using the method of reflections and a perturbation analysis for the deformed shapes. The perturbation analysis demonstrates explicitly that drops will tend to be aligned for ℬ 〉 0(d/a) where d is the separation distance between the drops. © 1993, Cambridge University Press. All rights reserved.

Description: When a small air bubble bursts from an equilibrium position at an air/water interface, a complex motion ensues resulting in the production of a high-speed liquid jet. This free-surface motion following the burst is modelled numerically using a boundary integral method. Jet formation and liquid entrainment rates from jet breakup into drops are calculated and compared with existing experimental evidence. In order to investigate viscous effects, a boundary layer is included in the calculations by employing a time-stepping technique which allows the boundary mesh to remain orthogonal to the surface. This allows an approximation of the vorticity development in the region of boundary-layer separation during jet formation. Calculated values of pressure and energy dissipation rates in the fluid indicate a violent motion, particularly for smaller bubbles. This has important implications for the biological industry where animal cells in bioreactors have been found to be killed by the presence of small bubbles. © 1993, Cambridge University Press. All rights reserved.

Description: The behaviour of concentric double emulsion droplets in linear flows is examined analytically, for the case when both fluid-fluid interfaces remain nearly spherical, and numerically, for the effect of finite interface deformation. The theoretical analysis is used to calculate the velocity fields interior and exterior to the particle, the first effects of flow-induced deformation, and the effective viscosity of a dilute emulsion of compound droplets. The numerical simulations allow for a complete investigation of the finite deformation of both the outer drop and the encapsulated particle. For concentric multiphase particles, there appear to be two distinct mechanisms of globule breakup: (i) continuous extension of the globule corresponding to non-existence of a steady particle shape or (ii) contact of the two interfaces at the globule centre, owing to incompatibility of the steady inner and outer interface shapes, even though the globule is only modestly deformed. Finally, the effect of different flow-types, i.e, uniaxial or biaxial extensional flows, is shown, in one example, to suggest breakup of the inner droplet even though the outer droplet maintains a steady shape. © 1990, Cambridge University Press. All rights reserved.

Description: Motivated by the recent work of Bajer & MOFFATT (1990), we investigate the kinematics of bounded steady Stokes flows. Specifically, we consider the streamlines inside a neutrally buoyant spherical drop immersed in a general linear flow. The Eulerian velocity field internal to the drop, known analytically, is a cubic function of position. For a wide range of parameters the internal streamlines, hence the fluid particle paths, may wander chaotically. Typical Poincare sections show both ordered and chaotic regions. The extent and existence of chaotic wandering is related to (i) the orientation of the vorticity vector relative to the principal axes of strain of the undisturbed flow and (ii) the magnitude of the vorticity relative to the magnitude of the rate-of-strain tensor. In the limit of small vorticity, we use the method of averaging to predict the size of the dominant island region. This yields the critical orientation of the vorticity vector at which this dominant island disappears so that particle paths fill almost the entire Poincare section. The problem studied here appears to be one of the simplest, physically realizable, bounded steady Stokes flows which produces chaotic streamlines. © 1991, Cambridge University Press. All rights reserved.

Description: 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 boundary-integral 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 ‘dilution5of the surfactant due to interfacial area increases. These two different effects of surface-active material are illustrated and discussed and their influence on the critical capillary number for breakup is presented. © 1990, Cambridge University Press. All rights reserved.

Description: An investigation of the interfacial-tension-driven fragmentation of a very long fluid filament in a quiescent viscous fluid is presented. Experiments covering almost three orders of magnitude in viscosity ratio reveal as many as 19 satellite droplets in between the largest droplets; complementary boundary-integral calculations are used to study numerically the evolution of the filament as a function of the viscosity ratio of the fluids and the initial wavenumber of the interface perturbation. Satellite drops are generated owing to multiple breakup sequences around the neck region of a highly deformed filament. In low-viscosity ratio systems, p 〈 0(0.1), the breakup mechanism is self-repeating in the sense that every pinch-off is always associated with the formation of a neck, the neck undergoes pinch-off, and the process repeats. In general the agreement between computations and experiments is excellent; both indicate that the initial wavenumber of the disturbance is important in the quantitative details of the generated drop size distributions. However, these details are insignificant when compared with the large variations produced in the drop size distributions owing to variation in the viscosity ratio. © 1992, Cambridge University Press. All rights reserved.