A FULLY ADAPTIVE, CONSERVATIVE FRONT TRACKING METHOD FOR THE SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOWS

This paper presents a fully adaptive formulation of the Front Tracking method for the simulation of incompressible, multiphase, bubbly flows, based on the Tryggvason formulation. The Navier-Stokes equations are discretized using a finite difference scheme, and domain discretization is carried out with Berger & Collela’s structured adaptive mesh refinement (SAMR) algorithm. Time discretization is based on SBDF scheme, with adaptive time stepping. The lagrangian interface is represented using the GTS library, which provides a volumeand shapepreserving remeshing algorithm, therefore minimizing the volume change due to non-conservative interpolation of the eulerian velocity field. Nevertheless, a simple volume recovery algorithm is also provided, along with a subgrid undulation removal algorithm based on the TSUR-3D algorithm [6]. Rising bubble flows were simulated under several regimes, showing small errors when comparing to experimental results.

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

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

[3]  Héctor D. Ceniceros,et al.  Short Note: A multi-phase flow method with a fast, geometry-based fluid indicator , 2005 .

[4]  James Glimm,et al.  Front tracking for hyperbolic systems , 1981 .

[5]  Kensuke Yokoi,et al.  A Numerical Method for Free-Surface Flows and Its Application to Droplet Impact on a Thin Liquid Layer , 2008, J. Sci. Comput..

[6]  D. Meiron,et al.  Efficient algorithms for solving static hamilton-jacobi equations , 2003 .

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

[8]  Peter Lindstrom,et al.  Fast and memory efficient polygonal simplification , 1998, Proceedings Visualization '98 (Cat. No.98CB36276).

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

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

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

[12]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[13]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[14]  Shuangshuang Jin,et al.  A comparison of algorithms for vertex normal computation , 2005, The Visual Computer.

[15]  R. LeVeque,et al.  Adaptive Mesh Refinement Using Wave-Propagation Algorithms for Hyperbolic Systems , 1998 .

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

[17]  M. Berger,et al.  Automatic adaptive grid refinement for the Euler equations , 1985 .

[18]  M. Berger,et al.  An Adaptive Version of the Immersed Boundary Method , 1999 .

[19]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .