Computing three-dimensional two-phase flows with a mass-conserving level set method

A method is described to compute three- dimensional two-phase flow, allowing large density ratios and coalescence and break-up of bubbles. The level set method is used to describe interfaces, and the volume-of-fluid method is used to ensure mass conservation. Efficiency in computing the interface dynamics is achieved by using a functional relation between the level set and volume-of-fluid functions. Difficulties and remedies in re-initialization of the level set function and inaccurate compution of surface tension are discussed. Test cases for validation are described, and demanding two-bubble computations to show the generality and the versatility of the method are presented.

[1]  Zhilin Li,et al.  The immersed interface method for the Navier-Stokes equations with singular forces , 2001 .

[2]  S. Osher,et al.  Computing interface motion in compressible gas dynamics , 1992 .

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

[4]  S. Osher,et al.  An improved level set method for incompressible two-phase flows , 1998 .

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

[6]  P. Wesseling.  Principles of Computational Fluid Dynamics , 2000 .

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

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

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

[10]  A. Segal,et al.  Modeling of Multi-Phase Flows with a Level-Set Method , 2004 .

[11]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[12]  Karol Mikula,et al.  High-Resolution Flux-Based Level Set Method , 2007, SIAM J. Sci. Comput..

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

[14]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

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

[16]  Cornelis Vuik,et al.  A mass conserving level set (MCLS) method for modeling of multi-phase flows , 2003 .

[17]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[18]  R. Fedkiw,et al.  USING THE PARTICLE LEVEL SET METHOD AND A SECOND ORDER ACCURATE PRESSURE BOUNDARY CONDITION FOR FREE SURFACE FLOWS , 2003 .

[19]  Arnold Reusken,et al.  An extended pressure finite element space for two-phase incompressible flows with surface tension , 2007, J. Comput. Phys..

[20]  P. Wesseling,et al.  A mass‐conserving Level‐Set method for modelling of multi‐phase flows , 2005 .

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

[22]  J. Kan.  A second-order accurate pressure correction scheme for viscous incompressible flow , 1986 .

[23]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

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

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

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

[27]  M. Gunsing.  Modelling of bubbly flows using volume of fluid, front tracking and discrete bubble models , 2004 .

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

[29]  Luis Gustavo Nonato,et al.  A front-tracking/front-capturing method for the simulation of 3D multi-fluid flows with free surfaces , 2004 .

[30]  Ng Niels Deen,et al.  Numerical simulation of gas bubbles behaviour using a three-dimensional volume of fluid method , 2005 .

[31]  D. Fletcher,et al.  A new volume of fluid advection algorithm: the defined donating region scheme , 2001 .

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

[33]  M. Rudman.  A Volume-Tracking Method for Incompressible Multifluid Flows with Large Density Variations , 1998 .