Hostname: page-component-76fb5796d-2lccl Total loading time: 0 Render date: 2024-04-28T05:37:44.927Z Has data issue: false hasContentIssue false
  • English
  • Français

Numerical and theoretical solutions for a drop spreading below a free fluid surface

Published online by Cambridge University Press:  26 April 2006

Dorothy M. Koch  and
Donald L. Koch
Dorothy M. Koch
Affiliation:
Department of Geology and Geophysics, Yale University, New Haven, CT 065120–8109, USA
Donald L. Koch
Affiliation:
School of Chemical Engineering, Cornell University, Ithaca, NY 14853, USA
Get access Rights & Permissions [Opens in a new window]

Abstract

Numerical solutions are derived for a viscous, buoyant drop spreading below a free fluid surface. The drop has zero interfacial tension, and we consider viscosity contrasts 0.1 < λ < 10 with the surrounding fluid half-space. The density contrast between the drop and outer fluid is assumed to be small compared with the density contrast at the fluid surface. The numerical solutions for the approach and initial spread of the drop below the fluid surface are obtained using the boundary integral method. To facilitate an investigation over a larger range of viscosity contrasts and for longer time periods, we solve for the motion of gravity currents at the fluid surface. For this geometry we also solve the boundary integral equations for the cases λ = 0 and 1/λ = 0.

For extensive drop spreading, the motion is described by asymptotic solutions. Three asymptotic solutions are derived, which apply for different values of the viscosity contrast relative to the aspect ratio ((radial extent R)/(drop thickness a)). For very low-viscosity drops (λ [Lt ] a/R[ln(R/a)]-1), the greatest resistance to spreading occurs at the drop rim, and the asymptotic solution is found using slender body theory. Drops with intermediate viscosity contrast (a/R [Lt ] λ [Lt ] R/a) are slowed primarily by shear stresses at the lower drop surface, and a lubrication solution (Lister & Kerr 1989) applies. The greatest resistance to the spread of very viscous drops (λ [Gt ] R/a) comes from the radial stresses within the drop, and the asymptotic solution is independent of the outer fluid. All drops having 0 [Lt ] λ [Lt ] ∞ will eventually spread according to lubrication theory, when their aspect ratio becomes sufficiently large relative to viscosity contrast.

Theoretical results are compared with numerical and experimental results for drops and gravity currents spreading at a fluid surface. The solutions can be applied to aspects of planetary mantle flow where temperature variations cause significant viscosity contrasts. The low-viscosity solution has been applied to study the encounter of a hot, low-viscosity upwelling plume with a planet surface (Koch 1994). Here we apply the high-viscosity asymptotic solution to study how cold downwelling slabs spread at a depth of neutral buoyancy in the Earth's mantle.

Type
Research Article
Information
Journal of Fluid Mechanics , Volume 287 , 25 March 1995 , pp. 251 - 278
Copyright
© 1995 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Ascoli, E. P., Dandy, D. S. & Leal, L. G. 1990 Buoyancy-driven motion of a deformable drop toward a planar wall at low Reynolds number. J. Fluid Mech. 213, 287–311.Google Scholar
Batchelor, G. K. 1970 Slender-body theory for particles of arbitrary cross-section in Stokes flow. J. Fluid Mech. 44, 419–440.Google Scholar
Bercovici, D., Schubert, G. & Tackley, P. J. 1993 On the penetration of the 660 km phase change by mantle downflows. Geophys. Res. Lett. 20, 2599–2602.Google Scholar
Campbell, I. H. & Griffiths, R. W. 1992 The changing nature of mantle hotspots through time: Implications for the chemical evolution of the mantle. J. Geology 92, 497–523.Google Scholar
Canright, D. R. 1987 A finite-amplitude analysis of the buoyant instability of a highly viscous film over a less viscous half-space. PhD thesis, University of California, Berkeley.
Chi, B. K. & Leal, L. G. 1989 A theoretical study of the motion of a viscous drop toward a fluid interface at low Reynolds number. J. Fluid Mech. 201, 123–146.Google Scholar
Davis, R. H., Schonberg, J. A. & Rallison, J. M. 1989 The lubrication force between two viscous drops. Phys. Fluids A 1, 77–81.Google Scholar
Geller, A. S., Lee, S. H. & Leal, L. G. 1986 The creeping motion of a spherical partical normal to a deformable interface. J. Fluid Mech. 169, 27–69.Google Scholar
Griffiths, R. W. & Campbell, I. H. 1990 Stirring and structure in mantle starting plumes. Earth Planet. Sci. Lett. 99, 66–78.Google Scholar
Griffiths, R. W. & Campbell, I. H. 1991 Interaction of mantle plume heads with the Earth's surface and onset of small-scale convection. J. Geophys. Res. 96, 18295–18310.Google Scholar
Griffiths, R. W. & Turner, J. S. 1988 Folding of viscous plumes impinging on a density or viscosity interface. Geophys. J. 95, 397–419.Google Scholar
Honda, S., Balachandar, S., Yuen, D. A. & Reuteler, D. 1993 Three-dimensional mantle dynamics with an endothermic phase transition. Geophys. Res. Lett. 20, 221–224.Google Scholar
Huppert, H. E. 1982 The propagation of two-dimensional and axisymmetric viscous gravity currents over a rigid horizontal surface. J. Fluid Mech. 121, 43–58.Google Scholar
Kerr, R. C. & Lister, J. R. 1987 The spread of subducted lithospheric material along the mid-mantle boundary. Earth Planet. Sci. Lett. 85, 241–247.Google Scholar
Kincaid, C. & Olson, P. 1987 An experimental study of subduction and slab migration. J. Geophys. Res. 92, 13832–13840.Google Scholar
Koch, D. M. 1993 A spreading drop model for plumes on Venus. PhD thesis, Yale University.
Koch, D. M. 1994 A spreading drop model for mantle plumes and volcanic features on Venus. J. Geophys. Res. 99, 2035–2052.Google Scholar
Koch, D. M. & Ribe, N. M. 1989 The effect of lateral viscosity variations on surface observables. Geophys. Res. Lett. 16, 535–538.Google Scholar
Koh, C. J. & Leal, L. G. 1989 The stability of drop shapes for translation at zero Reynolds number through a quiescent fluid. Phys. Fluids A 8, 1309–1313.Google Scholar
Kojima, M., Hinch, E. J. & Acrivos, A. 1984 The formation and expansion of a toroidal drop moving in a viscous fluid. Phys. Fluids A 27, 19–32.Google Scholar
Lee, S. H. & Leal, L. G. 1982 The motion of a sphere in the presence of a deformable interface. J. Fluid Mech. 87, 81–106.Google Scholar
Lister, J. R. & Kerr, R. C. 1989 The propagation of two-dimensional and axisymmetric viscous gravity currents at a fluid interface. J. Fluid Mech. 203, 215–249.Google Scholar
Machetel, P. & Weber, P. 1991 Intermittent layered convection in a model with an endothermic phase change at 670 km. Nature 350, 55–57.Google Scholar
Manga, M., Stone, H. A. & O'Connell, R. J. 1993 The interaction of plume heads with compositional discontinuities in the Earth's mantle. J. Geophys. Res. 98, 19979–19990.Google Scholar
Pozrikidis, C. 1990 a The instability of a moving viscous drop. J. Fluid Mech. 210, 1–21.Google Scholar
Pozrikidis, C. 1990 b The deformation of a liquid drop moving normal to a plane wall. J. Fluid Mech. 215, 331–363.Google Scholar
Rallison, J. M. & Acrivos, A. 1978 A numerical study of the deformation and burst of a viscous drop in an extensional flow. J. Fluid Mech. 1978, 191–200.Google Scholar
Richards, M. A. & Hager, B. H. 1984 Geoid anomalies in a dynamic earth. J. Geophys. Res. 89, 5987–6002.Google Scholar
Stone, H. A. & Leal, L. G. 1989 Relaxation and breakup of an initially extended drop in an otherwise quiescent fluid. J. Fluid Mech. 198, 399–427.Google Scholar
Tackley, P. J., Stevenson, D. J., Glatzmaier, G. A. & Schubert, G. 1993 Effects of an endothermic phase transition at 670 km depth in a spherical model of convection in the Earth's mantle. Nature 361, 699–704.Google Scholar
Tanzosh, J., Manga, M. & Stone, H. A. 1992 Boundary integral methods for viscous free-boundary problems: Deformation of single and multiple fluid-fluid interfaces. In Boundary Element Technology VI (ed. C. A. Brebbia & M. S. Ingber), pp. 19–39. Elsevier Applied Science.
Whttehead, J. A. & Luther, D. S. 1975 Dynamics of laboratory diapir and plume models. J. Geophys. Res. 80, 705–717.Google Scholar
Yiantsios, S. G. & Davis, R. H. 1990 On the buoyancy-driven motion of a drop towards a rigid surface or a deformable interface. J. Fluid Mech. 217, 547–573.Google Scholar