A novel boundary-integral algorithm for viscous interaction of deformable drops

Alexander Z. Zinchenko, Michael A. Rother, Robert H. Davis

Research output: Contribution to journalArticle

155 Citations (Scopus)

Abstract

A new three-dimensional boundary-integral algorithm for deformable drops moving in a viscous medium at low Reynolds numbers is developed, which overcomes some familiar difficulties with boundary-integral calculations. The algorithm is used to simulate different modes of interaction between drops or bubbles, primarily for buoyancy-driven motion. The present iterative method for mean curvature calculation is found to be more robust and accurate than contour integration schemes. A novel iterative strategy based on combining biconjugate gradient and simple iterations overcomes the poor convergence of "successive substitutions" for drops in very close approach with extreme viscosity ratio. A substantially new variational method of global mesh stabilization solves the problem of mesh degradation with advantageous, soft stability constraints. A curvatureless boundary-integral formulation is also derived and shown to provide, in principle, a more accurate description of the drop breakup than the conventional formulation. The efficiency of these techniques is demonstrated by numerical examples for two drops in gravity-induced motion with high surface resolutions. The present code successfully simulates mutual approach of slightly deformable drops to extremely small separations, as well as their rotation when in "apparent contact," thus bridging the gap between finite deformation calculations and a recent asymptotic theory for small capillary numbers. Also provided is a 3D simulation of the experimental phenomenon of enhanced bubble coalescence, discovered by Manga and Stone [J. Fluid Mech. 256, 647 (1993); 300, 231 (1995)]. For drops of viscosity comparable to that of the surrounding fluid, it is shown in contrast that breakup is a typical result of hydrodynamic interaction in gravity-induced motion for large and even moderate capillary numbers. The code is readily applicable to any type of an ambient flow and may be adapted to more than two drops.

Original languageEnglish (US)
Pages (from-to)1493-1511
Number of pages19
JournalPhysics of Fluids
Volume9
Issue number6
DOIs
StatePublished - Jan 1 1997

Fingerprint

mesh
bubbles
interactions
viscosity
gravitation
formulations
fluids
low Reynolds number
buoyancy
coalescing
iteration
stabilization
hydrodynamics
curvature
rocks
substitutes
degradation
gradients
simulation

Cite this

A novel boundary-integral algorithm for viscous interaction of deformable drops. / Zinchenko, Alexander Z.; Rother, Michael A.; Davis, Robert H.

In: Physics of Fluids, Vol. 9, No. 6, 01.01.1997, p. 1493-1511.

Research output: Contribution to journalArticle

Zinchenko, Alexander Z. ; Rother, Michael A. ; Davis, Robert H. / A novel boundary-integral algorithm for viscous interaction of deformable drops. In: Physics of Fluids. 1997 ; Vol. 9, No. 6. pp. 1493-1511.
@article{f0692adc817f4d7fbd9d9d02c87aa460,
title = "A novel boundary-integral algorithm for viscous interaction of deformable drops",
abstract = "A new three-dimensional boundary-integral algorithm for deformable drops moving in a viscous medium at low Reynolds numbers is developed, which overcomes some familiar difficulties with boundary-integral calculations. The algorithm is used to simulate different modes of interaction between drops or bubbles, primarily for buoyancy-driven motion. The present iterative method for mean curvature calculation is found to be more robust and accurate than contour integration schemes. A novel iterative strategy based on combining biconjugate gradient and simple iterations overcomes the poor convergence of {"}successive substitutions{"} for drops in very close approach with extreme viscosity ratio. A substantially new variational method of global mesh stabilization solves the problem of mesh degradation with advantageous, soft stability constraints. A curvatureless boundary-integral formulation is also derived and shown to provide, in principle, a more accurate description of the drop breakup than the conventional formulation. The efficiency of these techniques is demonstrated by numerical examples for two drops in gravity-induced motion with high surface resolutions. The present code successfully simulates mutual approach of slightly deformable drops to extremely small separations, as well as their rotation when in {"}apparent contact,{"} thus bridging the gap between finite deformation calculations and a recent asymptotic theory for small capillary numbers. Also provided is a 3D simulation of the experimental phenomenon of enhanced bubble coalescence, discovered by Manga and Stone [J. Fluid Mech. 256, 647 (1993); 300, 231 (1995)]. For drops of viscosity comparable to that of the surrounding fluid, it is shown in contrast that breakup is a typical result of hydrodynamic interaction in gravity-induced motion for large and even moderate capillary numbers. The code is readily applicable to any type of an ambient flow and may be adapted to more than two drops.",
author = "Zinchenko, {Alexander Z.} and Rother, {Michael A.} and Davis, {Robert H.}",
year = "1997",
month = "1",
day = "1",
doi = "10.1063/1.869275",
language = "English (US)",
volume = "9",
pages = "1493--1511",
journal = "Physics of Fluids",
issn = "1070-6631",
publisher = "American Institute of Physics",
number = "6",

}

TY - JOUR

T1 - A novel boundary-integral algorithm for viscous interaction of deformable drops

AU - Zinchenko, Alexander Z.

AU - Rother, Michael A.

AU - Davis, Robert H.

PY - 1997/1/1

Y1 - 1997/1/1

N2 - A new three-dimensional boundary-integral algorithm for deformable drops moving in a viscous medium at low Reynolds numbers is developed, which overcomes some familiar difficulties with boundary-integral calculations. The algorithm is used to simulate different modes of interaction between drops or bubbles, primarily for buoyancy-driven motion. The present iterative method for mean curvature calculation is found to be more robust and accurate than contour integration schemes. A novel iterative strategy based on combining biconjugate gradient and simple iterations overcomes the poor convergence of "successive substitutions" for drops in very close approach with extreme viscosity ratio. A substantially new variational method of global mesh stabilization solves the problem of mesh degradation with advantageous, soft stability constraints. A curvatureless boundary-integral formulation is also derived and shown to provide, in principle, a more accurate description of the drop breakup than the conventional formulation. The efficiency of these techniques is demonstrated by numerical examples for two drops in gravity-induced motion with high surface resolutions. The present code successfully simulates mutual approach of slightly deformable drops to extremely small separations, as well as their rotation when in "apparent contact," thus bridging the gap between finite deformation calculations and a recent asymptotic theory for small capillary numbers. Also provided is a 3D simulation of the experimental phenomenon of enhanced bubble coalescence, discovered by Manga and Stone [J. Fluid Mech. 256, 647 (1993); 300, 231 (1995)]. For drops of viscosity comparable to that of the surrounding fluid, it is shown in contrast that breakup is a typical result of hydrodynamic interaction in gravity-induced motion for large and even moderate capillary numbers. The code is readily applicable to any type of an ambient flow and may be adapted to more than two drops.

AB - A new three-dimensional boundary-integral algorithm for deformable drops moving in a viscous medium at low Reynolds numbers is developed, which overcomes some familiar difficulties with boundary-integral calculations. The algorithm is used to simulate different modes of interaction between drops or bubbles, primarily for buoyancy-driven motion. The present iterative method for mean curvature calculation is found to be more robust and accurate than contour integration schemes. A novel iterative strategy based on combining biconjugate gradient and simple iterations overcomes the poor convergence of "successive substitutions" for drops in very close approach with extreme viscosity ratio. A substantially new variational method of global mesh stabilization solves the problem of mesh degradation with advantageous, soft stability constraints. A curvatureless boundary-integral formulation is also derived and shown to provide, in principle, a more accurate description of the drop breakup than the conventional formulation. The efficiency of these techniques is demonstrated by numerical examples for two drops in gravity-induced motion with high surface resolutions. The present code successfully simulates mutual approach of slightly deformable drops to extremely small separations, as well as their rotation when in "apparent contact," thus bridging the gap between finite deformation calculations and a recent asymptotic theory for small capillary numbers. Also provided is a 3D simulation of the experimental phenomenon of enhanced bubble coalescence, discovered by Manga and Stone [J. Fluid Mech. 256, 647 (1993); 300, 231 (1995)]. For drops of viscosity comparable to that of the surrounding fluid, it is shown in contrast that breakup is a typical result of hydrodynamic interaction in gravity-induced motion for large and even moderate capillary numbers. The code is readily applicable to any type of an ambient flow and may be adapted to more than two drops.

UR - http://www.scopus.com/inward/record.url?scp=0030833424&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0030833424&partnerID=8YFLogxK

U2 - 10.1063/1.869275

DO - 10.1063/1.869275

M3 - Article

AN - SCOPUS:0030833424

VL - 9

SP - 1493

EP - 1511

JO - Physics of Fluids

JF - Physics of Fluids

SN - 1070-6631

IS - 6

ER -