An Efficient Algorithm for Hydrodynamical Interaction of Many Deformable Drops

An efficient and accurate 3D algorithm for dynamical simulations of many deformable drops with strong hydrodynamical interactions at zero Reynolds numbers is developed. The drop-to-medium viscosity ratio, λ, and the Bond number are arbitrary, and the drops are subject to gravity with stationary trip...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics Vol. 157; no. 2; pp. 539 - 587
Main Authors: Zinchenko, Alexander Z., Davis, Robert H.
Format: Journal Article
Language:English
Published: Elsevier Inc 20.01.2000
Subjects:
multipole methods
deformable drops
boundary integrals
sedimentation
hydrodynamical interactions
Stokes flow
ISSN:0021-9991, 1090-2716
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:An efficient and accurate 3D algorithm for dynamical simulations of many deformable drops with strong hydrodynamical interactions at zero Reynolds numbers is developed. The drop-to-medium viscosity ratio, λ, and the Bond number are arbitrary, and the drops are subject to gravity with stationary triply periodic boundary conditions. The algorithm, at each step, is a hybrid of boundary-integral and economical multipole techniques, with extensive use of rotational transformations and economical truncation of multipole expansions to optimize near-field interactions. A significant part of the code is the new, “best paraboloid-spline” technique for calculating the normal vectors and curvatures on drop surfaces, which greatly improves the quality of long-time simulations. Examples show the phenomenon of clustering in a concentrated sedimenting emulsion for λ=0.25 and 1, which leads to an increase in the average sedimentation velocity with time. A high efficiency of the method is demonstrated, with two orders-of-magnitude gains over the standard O(N2N2▵) boundary-integral technique for N∼102 drops in a periodic cell with N▵∼103 triangular boundary elements per drop, so that typical long-time dynamical simulations can be performed in a few days or weeks on a standard workstation (as compared to the several years which would be required for the same simulations using standard boundary-integral techniques). The effects of drop triangulation and truncation of multipole expansions on dynamical simulations are assessed.
ISSN:0021-9991
1090-2716
DOI:10.1006/jcph.1999.6384