An accurate adaptive solver for surface-tension-driven interfacial flows

A method combining an adaptive quad/octree spatial discretisation, geometrical Volume-Of-Fluid interface representation, balanced-force continuum-surface-force surface-tension formulation and height-function curvature estimation is presented. The extension of these methods to the quad/octree discretisation allows adaptive variable resolution along the interface and is described in detail. The method is shown to recover exact equilibrium (to machine accuracy) between surface-tension and pressure gradient in the case of a stationary droplet, irrespective of viscosity and spatial resolution. Accurate solutions are obtained for the classical test case of capillary wave oscillations. An application to the capillary breakup of a jet of water in air further illustrates the accuracy and efficiency of the method. The source code of the implementation is freely available as part of the Gerris flow solver.

[1]  Markus Bussmann,et al.  Adaptive VOF with curvature‐based refinement , 2007 .

[2]  David P. Schmidt,et al.  Adaptive tetrahedral meshing in free-surface flow , 2005 .

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

[4]  Heinz Pitsch,et al.  An accurate conservative level set/ghost fluid method for simulating turbulent atomization , 2008, J. Comput. Phys..

[5]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

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

[7]  Jean-François Remacle,et al.  Transient adaptivity applied to two-phase incompressible flows , 2008, J. Comput. Phys..

[8]  Dong-Yol Yang,et al.  Finite element analysis of transient fluid flow with free surface using VOF (volume-of-fluid) method and adaptive grid , 1998 .

[9]  P. Colella,et al.  A Conservative Adaptive Projection Method for the Variable Density Incompressible Navier-Stokes Equations , 1998 .

[10]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[11]  Xiaofeng Yang,et al.  An adaptive coupled level-set/volume-of-fluid interface capturing method for unstructured triangular grids , 2006, J. Comput. Phys..

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

[13]  T. Kowalewski,et al.  On the separation of droplets from a liquid jet , 1996 .

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

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

[16]  Alistair G.L. Borthwick,et al.  Finite‐volume‐type VOF method on dynamically adaptive quadtree grids , 2004 .

[17]  Sidney R. Nagel,et al.  Breakdown of scaling in droplet fission at high Reynolds number , 1997 .

[18]  Alexandre Joel Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations* , 1989 .

[19]  T. Chan,et al.  Robust multigrid methods for nonsmooth coefficient elliptic linear systems , 2000 .

[20]  J. Eggers,et al.  Universal pinching of 3D axisymmetric free-surface flow. , 1993, Physical review letters.

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

[22]  Said I. Abdel-Khalik,et al.  Accurate representation of surface tension using the level contour reconstruction method , 2005 .

[23]  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..

[24]  L. Rayleigh,et al.  XVI. On the instability of a cylinder of viscous liquid under capillary force , 1892 .

[25]  R. I. Issa,et al.  A Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes , 1999 .

[26]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[27]  S. Zaleski,et al.  Analytical relations connecting linear interfaces and volume fractions in rectangular grids , 2000 .

[28]  Alexandre J. Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations , 1969 .

[29]  Deborah Greaves,et al.  A quadtree adaptive method for simulating fluid flows with moving interfaces , 2004 .

[30]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for solving two phase incompressible flows , 2006, J. Comput. Phys..

[31]  Marcus Herrmann,et al.  A balanced force refined level set grid method for two-phase flows on unstructured flow solver grids , 2008, J. Comput. Phys..

[32]  Ronald Fedkiw,et al.  Simulating water and smoke with an octree data structure , 2004, ACM Trans. Graph..

[33]  Nikos Nikolopoulos,et al.  Three-dimensional numerical investigation of a droplet impinging normally onto a wall film , 2007, J. Comput. Phys..

[34]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

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

[36]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[37]  A. Brandt,et al.  The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous Coefficients , 1981 .

[38]  J. Strutt.  Scientific Papers: On the Instability of a Cylinder of Viscous Liquid under Capillary Force , 2009 .

[39]  Deborah Greaves,et al.  Simulation of viscous water column collapse using adapting hierarchical grids , 2006 .

[40]  L. Fuchs,et al.  High-order surface tension VOF-model for 3D bubble flows with high density ratio , 2004 .

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

[42]  G. Wittum,et al.  Two-Phase Flows on Interface Refined Grids Modeled with VOF, Staggered Finite Volumes, and Spline Interpolants , 2001 .

[43]  S. Popinet.  Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries , 2003 .

[44]  C. W. Hirt,et al.  NASA-VOF2D: A computer program for incompressible ows with free surfaces , 1985 .

[45]  S. Osher,et al.  Spatially adaptive techniques for level set methods and incompressible flow , 2006 .

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

[47]  Andreas Theodorakakos,et al.  Simulation of sharp gas–liquid interface using VOF method and adaptive grid local refinement around the interface , 2004 .

[48]  Fehmi Cirak,et al.  ZEITSCHRIFT FÜR ANGEWANDTE MATHEMATIK UND MECHANIK , 1952, Über Stammfaktoren bei ternären quadratischen Formen.

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

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

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

[52]  M. Brenner,et al.  A Cascade of Structure in a Drop Falling from a Faucet , 1994, Science.

[53]  S. Zaleski,et al.  A geometrical area-preserving volume-of-fluid advection method , 2003 .

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

[55]  Eugenio Aulisa,et al.  Interface reconstruction with least-squares fit and split advection in three-dimensional Cartesian geometry , 2007, J. Comput. Phys..

[56]  F. Durst,et al.  Comparison of volume-of-fluid methods for surface tension-dominant two-phase flows , 2006 .

[57]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[58]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[59]  Jie Li,et al.  Calcul d'Interface Affine par Morceaux , 1995 .

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

[61]  M. Davidson,et al.  An analysis of parasitic current generation in Volume of Fluid simulations , 2005 .

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

[63]  S. Zaleski,et al.  Interface reconstruction with least‐square fit and split Eulerian–Lagrangian advection , 2003 .

[64]  Vassili S. Sochnikov,et al.  Level set calculations of the evolution of boundaries on a dynamically adaptive grid , 2003 .

[65]  Mark Sussman,et al.  Improvements for calculating two-phase bubble and drop motion using an adaptive sharp interface method. , 2007 .

[66]  John B. Bell,et al.  Approximate Projection Methods: Part I. Inviscid Analysis , 2000, SIAM J. Sci. Comput..

[67]  D. Peregrine,et al.  The bifurcation of liquid bridges , 1990, Journal of Fluid Mechanics.

[68]  Andrea Prosperetti,et al.  Motion of two superposed viscous fluids , 1981 .

[69]  Takahiko Tanahashi,et al.  Numerical analysis of moving interfaces using a level set method coupled with adaptive mesh refinement , 2004 .

[70]  J. Eggers.  Nonlinear dynamics and breakup of free-surface flows , 1997 .

[71]  Elaine S. Oran,et al.  Surface tension and viscosity with Lagrangian hydrodynamics on a triangular mesh. Memorandum report , 1988 .

[72]  M. Rudman.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

[73]  C. Weber.  Zum Zerfall eines Flüssigkeitsstrahles , 1931 .

[74]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics , 2007, J. Comput. Phys..