Connectivity-free front tracking method for multiphase flows with free surfaces

In this study, a connectivity-free front tracking method is developed to simulate multiphase flows with free surfaces. This method is based on the point-set method which does not require any connectivities between interfacial points to represent the interface. The main advantage of the connectivity-free approach is the easiness in re-constructing the interface when large topology change occurs. It requires an indicator field to be constructed first based on the existing interface and the surface curvature and normal are then computed using the indicator field. Here, we adopt the reproducing kernel particle method (RKPM) interpolation function that provides the ability to deal with free-surface flows and the flexibility of using non-uniform meshes when local fine resolution is needed. A points regeneration scheme is developed to construct smooth interfaces and to automatically handle topology changes. The mass conservation is verified by performing a single vortex advection test. Several 2-D and 3-D numerical tests including an oscillating droplet, dam-breaking, two droplet impacting and multi-bubble merging are presented to show the accuracy and the robustness of the method.

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

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

[3]  Zhilin Li,et al.  A level-set method for interfacial flows with surfactant , 2006, J. Comput. Phys..

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

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

[6]  David P. Schmidt,et al.  Modeling merging and breakup in the moving mesh interface tracking method for multiphase flow simulations , 2009, J. Comput. Phys..

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

[8]  Axel Klar,et al.  A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method , 2009, J. Comput. Phys..

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

[10]  Metin Muradoglu,et al.  A front-tracking method for computation of interfacial flows with soluble surfactants , 2008, J. Comput. Phys..

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

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

[13]  Olivier Desjardins,et al.  A ghost fluid, level set methodology for simulating multiphase electrohydrodynamic flows with application to liquid fuel injection , 2010, J. Comput. Phys..

[14]  Tayfun E. Tezduyar,et al.  Stabilized-finite-element/interface-capturing technique for parallel computation of unsteady flows with interfaces , 2000 .

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

[16]  Jeroen A. S. Witteveen,et al.  An improved front tracking method for the Euler equations , 2005, J. Comput. Phys..

[17]  Sandro Manservisi,et al.  Simulation of axisymmetric jets with a finite element Navier-Stokes solver and a multilevel VOF approach , 2010, J. Comput. Phys..

[18]  Lucy T. Zhang,et al.  Immersed finite element method , 2004 .

[19]  Seungwon Shin,et al.  The Local Front Reconstruction Method for direct simulation of two- and three-dimensional multiphase flows , 2011, J. Comput. Phys..

[20]  Thomas-Peter Fries,et al.  The extended finite element method for two-phase and free-surface flows: A systematic study , 2011, J. Comput. Phys..

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

[22]  Wing Kam Liu,et al.  Reproducing kernel particle methods for structural dynamics , 1995 .

[23]  Ping Lin,et al.  Numerical simulation of 3D bubbles rising in viscous liquids using a front tracking method , 2008, J. Comput. Phys..

[24]  C. K. Thornhill,et al.  Part IV. An experimental study of the collapse of liquid columns on a rigid horizontal plane , 1952, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

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

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

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

[28]  Pao-Hsiung Chiu,et al.  A conservative phase field method for solving incompressible two-phase flows , 2011, J. Comput. Phys..

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