VOFTools - A software package of calculation tools for volume of fluid methods using general convex grids

The VOFTools library includes efficient analytical and geometrical routines for (1) area/volume computation, (2) truncation operations that typically arise in VOF (volume of fluid) methods, (3) area/volume conservation enforcement (VCE) in PLIC (piecewise linear interface calculation) reconstruction and(4) computation of the distance from a given point to the reconstructed interface. The computation of a polyhedron volume uses an efficient formula based on a quadrilateral decomposition and a 2D projection of each polyhedron face. The analytical VCE method is based on coupling an interpolation procedure to bracket the solution with an improved final calculation step based on the above volume computation formula. Although the library was originally created to help develop highly accurate advection and reconstruction schemes in the context of VOF methods, it may have more general applications. To assess the performance of the supplied routines, different tests, which are provided in FORTRAN and C, were implemented for several 2D and 3D geometries. Program summary: Program Title: VOFTools Program Files doi: http://dx.doi.org/10.17632/brrgt645bh.1 Licensing provisions: GNU General Public License, version 3. Programming language: FORTRAN and C, with C interfaces. Nature of problem: The package of routines includes simple and efficient analytical and geometrical tools for area/volume computation, truncation operations that typically arise in VOF (volume of fluid) methods, area/volume conservation enforcement (VCE) in PLIC (piecewise linear interface calculation) reconstruction and computation of the distance from a given point to the reconstructed interface. Solution method: The volume (area in 2D) computation of a polyhedron (polygon in 2D) uses an efficient formula based on a quadrilateral decomposition and a 2D projection of each polyhedron face. The analytical VCE method is based on coupling an interpolation bracketing procedure with an improved final calculation step based on the above volume computation formula. Also, the exact distance from a given point to a reconstructed polygonal interface is calculated. Restrictions: Convex 2D and 3D polytopes.

[1]  A. Oliva,et al.  A 3-D Volume-of-Fluid advection method based on cell-vertex velocities for unstructured meshes , 2014 .

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

[3]  F. Faura,et al.  A two-grid adaptive volume of fluid approach for dendritic solidification☆ , 2013 .

[4]  J. López,et al.  An improved PLIC-VOF method for tracking thin fluid structures in incompressible two-phase flows , 2005 .

[5]  F. Faura,et al.  A new volume of fluid method in three dimensions—Part II: Piecewise‐planar interface reconstruction with cubic‐Bézier fit , 2008 .

[6]  J. López,et al.  A volume of fluid method based on multidimensional advection and spline interface reconstruction , 2004 .

[7]  Julio Hernández,et al.  Analytical and geometrical tools for 3D volume of fluid methods in general grids , 2008, J. Comput. Phys..

[8]  Julio Hernández,et al.  A volume of fluid approach for crystal growth simulation , 2010, J. Comput. Phys..

[9]  S. Zaleski,et al.  Analytical relations connecting linear interfaces and volume fractions in rectangular grids , 2000 .

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

[11]  Michael Dodd,et al.  A mass-conserving volume-of-fluid method: Volume tracking and droplet surface-tension in incompressible isotropic turbulence , 2014 .

[12]  F. Faura,et al.  A new volume of fluid method in three dimensions—Part I: Multidimensional advection method with face‐matched flux polyhedra , 2008 .

[13]  Julio Hernández,et al.  A new volume conservation enforcement method for PLIC reconstruction in general convex grids , 2016, J. Comput. Phys..

[14]  Dominique Legendre,et al.  On the combined effects of surface tension force calculation and interface advection on spurious currents within Volume of Fluid and Level Set frameworks , 2015, J. Comput. Phys..

[15]  Patrice Koehl.  Fast Recursive Computation of 3D Geometric Moments from Surface Meshes , 2012, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[16]  Kokichi Sugihara,et al.  A Robust and Consistent Algorithm for Intersecting Convex Polyhedra , 1994, Comput. Graph. Forum.

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

[18]  B. Premachandran,et al.  A Coupled Level Set and Volume of Fluid method with multi-directional advection algorithms for two-phase flows with and without phase change , 2014 .

[19]  Marianne M. Francois,et al.  An interface reconstruction method based on analytical formulae for 2D planar and axisymmetric arbitrary convex cells , 2014, J. Comput. Phys..

[20]  Parviz Moin,et al.  Accurate interface normal and curvature estimates on three-dimensional unstructured non-convex polyhedral meshes , 2015, J. Comput. Phys..