The Local Front Reconstruction Method for direct simulation of two- and three-dimensional multiphase flows

We present a new interface reconstruction technique, the Local Front Reconstruction Method (LFRM), for incompressible multiphase flows. This new method falls in the category of Front Tracking methods but it shares automatic topology handling characteristics of the previously proposed Level Contour Reconstruction Method (LCRM). The LFRM tracks the phase interface explicitly as in Front Tracking but there is no logical connectivity between interface elements thus greatly easing the algorithmic complexity. Topological changes such as interfacial merging or pinch off are dealt with automatically and naturally as in the Level Contour Reconstruction Method. Here the method is described for both two- and three-dimensional flow geometries. The interfacial reconstruction technique in the LFRM differs from that in the LCRM formulation by foregoing using an Eulerian distance field function. Instead, the LFRM uses information from the original interface elements directly to generate the new interface in a mass conservative way thus showing significantly improved local mass conservation. Because the reconstruction procedure is independently carried out in each individual reconstruction cell after an initial localization process, an adaptive reconstruction procedure can be easily implemented to increase the accuracy while at the same time significantly decreasing the computational time required to perform the reconstruction. Several benchmarking tests are performed to validate the improved accuracy and computational efficiency as compared to the LCRM. The results demonstrate superior performance of the LFRM in maintaining detailed interfacial shapes and good local mass conservation especially when using low-resolution Eulerian grids.

[1]  Héctor D. Ceniceros,et al.  Study of the long-time dynamics of a viscous vortex sheet with a fully adaptive nonstiff method , 2004 .

[2]  A. Chorin Numerical Solution of the Navier-Stokes Equations* , 1989 .

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

[4]  Hong Yi,et al.  A survey of the marching cubes algorithm , 2006, Comput. Graph..

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

[6]  Héctor D. Ceniceros,et al.  Three-dimensional, fully adaptive simulations of phase-field fluid models , 2010, J. Comput. Phys..

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

[8]  Wei Shyy,et al.  Three-dimensional adaptive Cartesian grid method with conservative interface restructuring and reconstruction , 2007, J. Comput. Phys..

[9]  Gihun Son,et al.  EFFICIENT IMPLEMENTATION OF A COUPLED LEVEL-SET AND VOLUME-OF-FLUID METHOD FOR THREE-DIMENSIONAL INCOMPRESSIBLE TWO-PHASE FLOWS , 2003 .

[10]  Eugenio Aulisa,et al.  A novel representation of the surface tension force for two-phase flow with reduced spurious currents , 2006 .

[11]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

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

[13]  I. Yoon,et al.  Tetra-marching procedure for high order Level Contour Reconstruction Method , 2010 .

[14]  Frank Losasso,et al.  A fast and accurate semi-Lagrangian particle level set method , 2005 .

[15]  Gretar Tryggvason,et al.  Direct numerical simulations of bubbly flows. Part 1. Low Reynolds number arrays , 1998, Journal of Fluid Mechanics.

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

[17]  Gretar Tryggvason,et al.  Direct numerical simulations of three-dimensional bubbly flows , 1999 .

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

[19]  Mikhail J. Shashkov,et al.  Reconstruction of multi-material interfaces from moment data , 2008, J. Comput. Phys..

[20]  Emil Coyajee,et al.  Numerical simulation of drop impact on a liquid-liquid interface with a multiple marker front-capturing method , 2009, J. Comput. Phys..

[21]  Hender López,et al.  Oscillation of viscous drops with smoothed particle hydrodynamics. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  Gretar Tryggvason,et al.  Numerical simulations of three-dimensional drop collisions , 1996 .

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

[24]  R. Clift,et al.  Bubbles, Drops, and Particles , 1978 .

[25]  Seungwon Shin,et al.  Modeling three-dimensional multiphase flow using a level contour reconstruction method for front tracking without connectivity , 2002 .

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

[27]  David P. Schmidt,et al.  A moving mesh interface tracking method for 3D incompressible two-phase flows , 2007, J. Comput. Phys..

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

[29]  David J. Torres,et al.  The point-set method: front-tracking without connectivity , 2000 .

[30]  O. Lebaigue,et al.  The second gradient method for the direct numerical simulation of liquid—vapor flows with phase change , 2001 .

[31]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[32]  D. Juric,et al.  A hybrid interface method for three‐dimensional multiphase flows based on front tracking and level set techniques , 2009 .

[33]  B. Engquist,et al.  Discretization of Dirac delta functions in level set methods , 2005 .

[34]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

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

[36]  A. Chorin Numerical solution of the Navier-Stokes equations , 1968 .

[37]  L. G. Leal,et al.  Numerical solution of free-boundary problems in fluid mechanics. Part 1. The finite-difference technique , 1984, Journal of Fluid Mechanics.

[38]  Tadashi Watanabe,et al.  Numerical simulation of oscillations and rotations of a free liquid droplet using the level set method , 2008 .

[39]  Damir Juric,et al.  High order level contour reconstruction method , 2007 .

[40]  R. Scardovelli,et al.  A mixed markers and volume-of-fluid method for the reconstruction and advection of interfaces in two-phase and free-boundary flows , 2003 .

[41]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[42]  Seungwon Shin,et al.  Computation of the curvature field in numerical simulation of multiphase flow , 2007, J. Comput. Phys..

[43]  J. Sethian Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[44]  Qiang Zhang,et al.  Three-Dimensional Front Tracking , 1998, SIAM J. Sci. Comput..

[45]  R. Scardovelli,et al.  A surface marker algorithm coupled to an area-preserving marker redistribution method for three-dimensional interface tracking , 2004 .

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

[47]  Frederick Stern,et al.  An improved particle correction procedure for the particle level set method , 2009, J. Comput. Phys..

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

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

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

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

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