A novel coupled level set and volume of fluid method for sharp interface capturing on 3D tetrahedral grids

We present a new three-dimensional hybrid level set (LS) and volume of fluid (VOF) method for free surface flow simulations on tetrahedral grids. At each time step, we evolve both the level set function and the volume fraction. The level set function is evolved by solving the level set advection equation using a second-order characteristic based finite volume method. The volume fraction advection is performed using a bounded compressive normalized variable diagram (NVD) based scheme. The interface is reconstructed based on both the level set and the volume fraction information. The novelty of the method lies in that we use an analytic method for finding the intercepts on tetrahedral grids, which makes interface reconstruction efficient and conserves volume of fluid exactly. Furthermore, the advection of volume fraction makes use of the NVD concept and switches between different high resolution differencing schemes to yield a bounded scalar field, and to preserve both smoothness and sharp definition of the interface. The method is coupled to a well validated finite volume based Navier-Stokes incompressible flow solver. The code validation shows that our method can be employed to resolve complex interface changes efficiently and accurately. In addition, the centroid and intercept data available as a by-product of the proposed interface reconstruction scheme can be used directly in near-interface sub-grid models in large eddy simulation.

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

[2]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[3]  Stuart E. Rogers,et al.  Steady and unsteady solutions of the incompressible Navier-Stokes equations , 1991 .

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

[5]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

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

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

[8]  Timothy J. Pedley,et al.  Flow along a channel with a time-dependent indentation in one wall: the generation of vorticity waves , 1985, Journal of Fluid Mechanics.

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

[10]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

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

[12]  Frederick Stern,et al.  A Sharp Interface Method for Two-Phase Flows Interacting with Moving Bodies , 2007 .

[13]  Boo Cheong Khoo,et al.  The ghost fluid method for compressible gas-water simulation , 2005 .

[14]  Stéphane Popinet,et al.  Bubble collapse near a solid boundary: a numerical study of the influence of viscosity , 2002, Journal of Fluid Mechanics.

[15]  C. W. Hirt,et al.  SOLA-VOF: a solution algorithm for transient fluid flow with multiple free boundaries , 1980 .

[16]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

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

[18]  F. Sotiropoulos,et al.  A hybrid Cartesian/immersed boundary method for simulating flows with 3D, geometrically complex, moving bodies , 2005 .

[19]  Johannes Vassiliou Oulis.  COMPUTATION OF TWO-DIMENSIONAL DAM-BREAK FLOOD FLOWS , 1992 .

[20]  M. Kawahara,et al.  LAGRANGIAN FINITE ELEMENT ANALYSIS APPLIED TO VISCOUS FREE SURFACE FLUID FLOW , 1987 .

[21]  S. Cummins,et al.  Estimating curvature from volume fractions , 2005 .

[22]  J. C. Martin.  An experimental study of the collapse of liquid column on a rigid horizontal plane , 1952 .

[23]  W. Rider,et al.  Stretching and tearing interface tracking methods , 1995 .

[24]  Pengzhi Lin,et al.  A numerical study of breaking waves in the surf zone , 1998, Journal of Fluid Mechanics.

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

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

[27]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[28]  Dimitris E. Papantonis,et al.  A characteristic-based method for incompressible flows , 1994 .

[29]  S. Koshizuka.  A particle method for incompressible viscous flow with fluid fragmentation , 1995 .

[30]  Nasser Ashgriz,et al.  FLAIR: fluz line-segment model for advection and interface reconstruction , 1991 .

[31]  Yong Zhao,et al.  Higher-order characteristics-based method for incompressible flow computation on unstructured grids , 2001 .

[32]  R. Franke.  Scattered data interpolation: tests of some methods , 1982 .

[33]  C. W. Hirt,et al.  Improved free surface boundary conditions for numerical incompressible-flow calculations , 1971 .

[34]  Yong Zhao,et al.  An efficient parallel/unstructured-multigrid preconditioned implicit method for simulating 3D unsteady compressible flows with moving objects , 2006, J. Comput. Phys..

[35]  D. Durran.  Numerical Methods for Fluid Dynamics , 2010 .

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

[37]  A. Jameson.  Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings , 1991 .

[38]  Yoshiaki Oka,et al.  A particle method for calculating splashing of incompressible viscous fluid , 1995 .

[39]  S. Afkhami,et al.  Height functions for applying contact angles to 2D VOF simulations , 2008 .

[40]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[41]  Xiaofeng Yang,et al.  Analytic relations for reconstructing piecewise linear interfaces in triangular and tetrahedral grids , 2006, J. Comput. Phys..

[42]  B. P. Leonard,et al.  The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection , 1991 .

[43]  Yong Zhao,et al.  Parallel unsteady incompressible viscous flow computations using an unstructured multigrid method , 2003 .

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

[45]  Rainald Löhner,et al.  On the simulation of flows with violent free surface motion , 2006 .

[46]  Marius Paraschivoiu,et al.  Second order accurate volume tracking based on remapping for triangular meshes , 2003 .

[47]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[48]  Yong Zhao,et al.  A high-resolution characteristics-based implicit dual time-stepping VOF method for free surface flow simulation on unstructured grids , 2002 .

[49]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

[50]  R. I. Issa,et al.  A Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes , 1999 .

[51]  Yong Zhao,et al.  A matrix-free implicit unstructured multigrid finite volume method for simulating structural dynamics and fluid-structure interaction , 2007, J. Comput. Phys..

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

[53]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[54]  P. Gaskell,et al.  Curvature‐compensated convective transport: SMART, A new boundedness‐ preserving transport algorithm , 1988 .

[55]  Djamel Lakehal,et al.  Multi-physics treatment in the vicinity of arbitrarily deformable gas-liquid interfaces , 2007, J. Comput. Phys..

[56]  Rainald Löhner,et al.  An accurate, fast, matrix-free implicit method for computing unsteady flows on unstructured grids , 2001 .

[57]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

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

[59]  Rainald Löhner,et al.  An accurate, fast, matrix-free implicit method for computing unsteady flows on unstructured grids , 1999 .

[60]  A. Jameson,et al.  Finite volume solution of the two-dimensional Euler equations on a regular triangular mesh , 1985 .