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

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.