A 3D moving mesh Finite Element Method for two-phase flows

A 3D ALE Finite Element Method is developed to study two-phase flow phenomena using a new discretization method to compute the surface tension forces. The computational method is based on the Arbitrary Lagrangian-Eulerian formulation (ALE) and the Finite Element Method (FEM), creating a two-phase method with an improved model for the liquid-gas interface. An adaptive mesh update procedure is also proposed for effective management of the mesh to remove, add and repair elements, since the computational mesh nodes move according to the flow. The ALE description explicitly defines the two-phase interface position by a set of interconnected nodes which ensures a sharp representation of the boundary, including the role of the surface tension. The proposed methodology for computing the curvature leads to accurate results with moderate programming effort and computational cost. Static and dynamic tests have been carried out to validate the method and the results have compared well to analytical solutions and experimental results found in the literature, demonstrating that the new proposed methodology provides good accuracy to describe the interfacial forces and bubble dynamics. This paper focuses on the description of the proposed methodology, with particular emphasis on the discretization of the surface tension force, the new remeshing technique, and the validation results. Additionally, a microchannel simulation in complex geometry is presented for two elongated bubbles. (C) 2014 Elsevier Inc. All rights reserved.

[1]  Mark Meyer,et al.  Implicit fairing of irregular meshes using diffusion and curvature flow , 1999, SIGGRAPH.

[2]  Gretar Tryggvason,et al.  Bifurcation of tracked scalar waves , 1986 .

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

[4]  H. Si,et al.  Adaptive tetrahedral mesh generation by constrained Delaunay refinement , 2008 .

[5]  W. J. Gordon,et al.  Transfinite mappings and their application to grid generation , 1982 .

[6]  F. Brezzi.  On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers , 1974 .

[7]  J. Remacle,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

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

[9]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[10]  N. Mangiavacchi,et al.  A Lagrangian level‐set approach for the simulation of incompressible two‐fluid flows , 2005 .

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

[12]  Ramesh Nallapati,et al.  A moving unstructured staggered mesh method for the simulation of incompressible free-surface flows , 2003 .

[13]  Peter Wriggers,et al.  Arbitrary Lagrangian Eulerian finite element analysis of free surface flow , 2000 .

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

[15]  Jinsong Hua,et al.  Numerical simulation of bubble rising in viscous liquid , 2007, J. Comput. Phys..

[16]  Howard H. Hu,et al.  Direct numerical simulations of fluid-solid systems using the arbitrary Langrangian-Eulerian technique , 2001 .

[17]  Lingling Wu,et al.  A simple package for front tracking , 2006, J. Comput. Phys..

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

[19]  J. Thome,et al.  3D ALE Finite-Element Method for Two-Phase Flows With Phase Change , 2014 .

[20]  I. Babuska.  Error-bounds for finite element method , 1971 .

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

[22]  A. A. Amsden,et al.  The SMAC method: A numerical technique for calculating incompressible fluid flow , 1970 .

[23]  Luis Gustavo Nonato,et al.  Front tracking with moving-least-squares surfaces , 2008, J. Comput. Phys..

[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]  Sashikumaar Ganesan,et al.  Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants , 2012, J. Comput. Phys..

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

[27]  G. Tallini,et al.  ON THE EXISTENCE OF , 1996 .