Wavelets and fast summations for particle simulations of gravitational flows of miscible drops

Abstract For sedimentation of miscible drops through quiescent liquid of the same viscosity, a recent paper [J. Fluid. Mech. 447 (2001) 299] has shown the effectiveness of computer simulations based upon a swarm of point forces in tracking coalescence, mixing and rupture. Robustness of the approach was offset by the slow O(Np2) summations needed to calculate the mutual viscous interactions between all Np particles. Motivated by applications of wavelets to linear operators, this paper develops a conceptually simple scheme for dramatically accelerating the simulations. After lumping the particles together into N local clusters (three-dimensional “pixels”), a Haar discrete wavelet transform (DWT) is used to “compress” the “bitmapped” six-dimensional image of pixel–pixel hydrodynamic interactions. Depending upon the criterion of accuracy, the numerically observed scaling of the operation count seems to be either O(N) or O[N( log N) α ] . The DWT also works without modification for hydrodynamic wall effects, where the kernel is not purely translational and therefore fast convolutions (FFT) cannot be used.

[1]  Alexander Z. Zinchenko,et al.  Cusping, capture, and breakup of interacting drops by a curvatureless boundary-integral algorithm , 1999, Journal of Fluid Mechanics.

[2]  A. Acrivos,et al.  A numerical study of the deformation and burst of a viscous drop in an extensional flow , 1978, Journal of Fluid Mechanics.

[3]  F. Arecchi,et al.  An Experimental Investigation of the Break-up of a Liquid Drop Falling in a Miscible Fluid , 1989 .

[4]  J. CARRIERt,et al.  A FAST ADAPTIVE MULTIPOLE ALGORITHM FOR PARTICLE SIMULATIONS * , 2022 .

[5]  Moti Lal,et al.  Dynamics of a drop at a liquid/solid interface in simple shear fields: A mesoscopic simulation study , 1999 .

[6]  William H. Press,et al.  Numerical Recipes in Fortran 77 , 1992 .

[7]  A swarm of Stokeslets with interfacial tension , 2001 .

[8]  William H. Press,et al.  Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing , 1992 .

[9]  Johannes M. Nitsche,et al.  Break-up of a falling drop containing dispersed particles , 1997, Journal of Fluid Mechanics.

[10]  Vittorio Cristini,et al.  Drop breakup in three-dimensional viscous flows , 1998 .

[11]  C. Pozrikidis.  Boundary Integral and Singularity Methods for Linearized Viscous Flow: Index , 1992 .

[12]  U. Schaflinger,et al.  Interfacial phenomena in suspensions , 1999 .

[13]  Michael Manga,et al.  Buoyancy-driven interactions between two deformable viscous drops , 1993, Journal of Fluid Mechanics.

[14]  Walter Meile,et al.  Coalescence, torus formation and breakup of sedimenting drops: experiments and computer simulations , 2001, Journal of Fluid Mechanics.

[15]  Daniel D. Joseph,et al.  Fundamentals of two-fluid dynamics , 1993 .

[16]  C. Pozrikidis,et al.  The instability of a moving viscous drop , 1990, Journal of Fluid Mechanics.

[17]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[18]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[19]  J. Koelman,et al.  Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics , 1992 .

[20]  A. Sangani,et al.  An O(N) algorithm for Stokes and Laplace interactions of particles , 1996 .

[21]  R. Phillips.  A singularity method for calculating time-dependent viscoelastic flows with integral constitutive equations , 2003, Journal of Fluid Mechanics.

[22]  J. Blake,et al.  A note on the image system for a stokeslet in a no-slip boundary , 1971, Mathematical Proceedings of the Cambridge Philosophical Society.

[23]  C. Law,et al.  Molecular simulation of droplet collision in the presence of ambient gas , 1999 .

[24]  V. Rokhlin.  Rapid solution of integral equations of classical potential theory , 1985 .

[25]  H. Stone,et al.  The interaction of plume heads with compositional discontinuities in the Earth's mantle , 1993 .

[26]  R. Coifman,et al.  Fast wavelet transforms and numerical algorithms I , 1991 .

[27]  Andreas Acrivos,et al.  The formation and expansion of a toroidal drop moving in a viscous fluid , 1984 .

[28]  M. Manga.  Interactions between mantle diapirs , 1997 .

[29]  L. Nitsche,et al.  Atomistic SPH and a Link between Diffusion and Interfacial Tension , 2002 .

[30]  L Greengard,et al.  Fast Algorithms for Classical Physics , 1994, Science.

[31]  Gerald Kaiser,et al.  A Friendly Guide to Wavelets , 1994 .

[32]  T. Dąbroś,et al.  A singularity method for calculating hydrodynamic forces and particle velocities in low-Reynolds-number flows , 1985, Journal of Fluid Mechanics.

[33]  Jian-Ming Jin,et al.  Fast and Efficient Algorithms in Computational Electromagnetics , 2001 .

[34]  Robert H. Davis.  Buoyancy-driven viscous interaction of a rising drop with a smaller trailing drop , 1999 .

[35]  Vittorio Cristini,et al.  An adaptive mesh algorithm for evolving surfaces: simulation of drop breakup and coalescence , 2001 .

[36]  Joseph John Thomson,et al.  V. On the formation of vortex rings by drops falling into liquids, and some allied phenomena , 1886, Proceedings of the Royal Society of London.

[37]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[38]  U. Landman,et al.  Formation, stability, and breakup of nanojets , 2000, Science.

[39]  Patrick B. Warren,et al.  Mesoscopic Simulation of Drops in Gravitational and Shear Fields , 2000 .

[40]  Ami Harten,et al.  Fast multiresolution algorithms for matrix-vector multiplication , 1994 .

[41]  L. G. Leal,et al.  The stability of drop shapes for translation at zero Reynolds number through a quiescent fluid , 1989 .

[42]  James S. Walker,et al.  A Primer on Wavelets and Their Scientific Applications , 1999 .

[43]  Naoya Yoshioka,et al.  The behavior of a swarm of particles moving in a viscous fluid , 1978 .

[44]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[45]  Quasimolecular simulation of large liquid drops , 1990 .