Improvements for calculating two-phase bubble and drop motion using an adaptive sharp interface method.

In this paper, we describe new tech- niques for numerically approximating two-phase flows. Specifically, we present new techniques fortreatingtheviscosityandsurfacetensionterms that appear in the Navier-Stokes equations for in- compressible two-phase flow. Our resulting nu- merical method has the property that results com- putedusingourtwo-phasealgorithmapproachthe corresponding "one-phase" algorithm in the limit of zero gas density/viscosity; i.e. the two-phase results approach the one-phase free-boundary re- sultsinthelimitthatthegas isassumed to become a uniform pressure void. By grid convergence checks and comparison with previous experimen- tal data, we shall demonstrate the advantages of our new proposed discretizations. Keyword: volume-of-fluid method, sharp inter- face, multiphase flow, surface tension, viscosity

[1]  Wei Shyy,et al.  A fixed-grid, sharp-interface method for bubble dynamics and phase change , 2001 .

[2]  Mark Sussman,et al.  A Computational Study of Bubble Motion in Newtonian and Viscoelastic Fluids , 2005 .

[3]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[4]  S. Zaleski,et al.  Volume-of-Fluid Interface Tracking with Smoothed Surface Stress Methods for Three-Dimensional Flows , 1999 .

[5]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[6]  Martin E. Weber,et al.  Bubbles in viscous liquids: shapes, wakes and velocities , 1981, Journal of Fluid Mechanics.

[7]  Gretar Tryggvason,et al.  Head-on collision of drops—A numerical investigation , 1993 .

[8]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[9]  P. Colella,et al.  Non-convex profile evolution in two dimensions using volume of fluids , 1997 .

[10]  D. B. Kothe,et al.  Convergence and accuracy of kernel-based continuum surface tension models , 1998 .

[11]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[12]  Mark Sussman,et al.  A Discontinuous Spectral Element Method for the Level Set Equation , 2003, J. Sci. Comput..

[13]  Stéphane Popinet,et al.  A front-tracking algorithm for accurate representation of surface tension , 1999 .

[14]  Mark Sussman,et al.  A parallelized, adaptive algorithm for multiphase flows in general geometries , 2005 .

[15]  J. Li,et al.  Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method , 2000 .

[16]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[17]  S. Cummins,et al.  Estimating curvature from volume fractions , 2005 .

[18]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[19]  Osamu Tatebe,et al.  The multigrid preconditioned conjugate gradient method , 1993 .

[20]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

[21]  V. Cristini,et al.  Scalings for Droplet Sizes in Shear-Driven Breakup: Non-Microfluidic Ways to Monodisperse Emulsions , 2006 .

[22]  J. G. Hnat,et al.  Spherical cap bubbles and skirt formation , 1976 .

[23]  H. Udaykumar,et al.  Sharp interface Cartesian grid method II: A technique for simulating droplet interactions with surfaces of arbitrary shape , 2005 .

[24]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

[25]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[26]  L. G. Leal,et al.  Numerical solution of free-boundary problems in fluid mechanics. Part 2. Buoyancy-driven motion of a gas bubble through a quiescent liquid , 1984, Journal of Fluid Mechanics.

[27]  W. Shyy,et al.  Elafint: a Mixed Eulerian-Lagrangian Method for Fluid Flows with Complex and Moving Boundaries , 1996 .

[28]  Nasser Ashgriz,et al.  A computational method for determining curvatures , 1989 .

[29]  Ronald Fedkiw,et al.  Animation and rendering of complex water surfaces , 2002, ACM Trans. Graph..

[30]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

[31]  V. Cristini,et al.  Adaptive unstructured volume remeshing - II: Application to two- and three-dimensional level-set simulations of multiphase flow , 2005 .

[32]  M. Sussman.  A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles , 2003 .

[33]  P. Colella,et al.  An Adaptive Level Set Approach for Incompressible Two-Phase Flows , 1997 .

[34]  Alexandre J. Chorin,et al.  Curvature and solidification , 1985 .

[35]  M. Renardy,et al.  PROST: a parabolic reconstruction of surface tension for the volume-of-fluid method , 2002 .

[36]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .