Cellwise conservative unsplit advection for the volume of fluid method

We present a cellwise conservative unsplit (CCU) advection scheme for the volume of fluid method (VOF) in 2D. Contrary to other schemes based on explicit calculations of the flux balances, the CCU advection adopts a cellwise approach where the pre-images of the control volumes are traced backwards through the flow map. The donating regions of the fluxes are calculated via the streaklines of the grid intersections, represented as polygonal chains whose vertices are determined by backward tracing of particles injected in the flow at different times. High order accuracy is obtained from the fourth-order Runge-Kutta method, where intermediate velocities along pathlines are determined with quadratic temporal and bicubic spatial interpolations. The volumes of the donating regions are corrected in order to fulfill the discrete continuity of incompressible flows. Consequently, the calculation produces non-overlapping donating regions and pre-images with conforming edges to their neighbors, resulting in the conservativeness and the boundedness (liquid volume fraction inside the interval 0 , 1 ] ) of the CCU advection scheme. Finally, the update of the liquid volume fractions is computed from the intersections of the pre-image polygons with the reconstructed interfaces.The CCU scheme is tested on several benchmark tests for the VOF advection, together with the standard piecewise linear interface calculation (PLIC). The geometrical errors of the CCU compare favorably with other unsplit VOF-PLIC schemes. Finally, potential improvements of the VOF method with the use of more precise interface representation techniques and the future extension of the CCU scheme to 3D are discussed. Cellwise VOF-update through backward tracing of the cell's pre-image.Accuracy enhanced via bicubic spatial interpolation and polygonal chain streakline.The CCU advection scheme is fully conservative and bounded.No limitation on the Courant number.Geometrical errors dominated by PLIC interface representation errors.

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

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

[3]  Qinghai Zhang.  Highly accurate Lagrangian flux calculation via algebraic quadratures on spline-approximated donating regions , 2013 .

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

[5]  Djamel Lakehal,et al.  A Newton--Krylov Solver for Remapping-Based Volume-of-Fluid Methods , 2008, SIAM J. Sci. Comput..

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

[7]  P. Colella,et al.  A conservative three-dimensional Eulerian method for coupled solid-fluid shock capturing , 2002 .

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

[9]  D. Fletcher,et al.  A New Volume of Fluid Advection Algorithm , 2000 .

[10]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[11]  S. Zaleski,et al.  Volume-of-Fluid Interface Tracking with Smoothed Surface Stress Methods for Three-Dimensional Flows , 1999 .

[12]  E. Puckett,et al.  A High-Order Projection Method for Tracking Fluid Interfaces in Variable Density Incompressible Flows , 1997 .

[13]  T. Sundararajan,et al.  A Quadratic Spline based Interface (QUASI) reconstruction algorithm for accurate tracking of two-phase flows , 2009, J. Comput. Phys..

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

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

[16]  G. Wittum,et al.  Two-Phase Flows on Interface Refined Grids Modeled with VOF, Staggered Finite Volumes, and Spline Interpolants , 2001 .

[17]  Joseph O'Rourke,et al.  Computational geometry in C (2nd ed.) , 1998 .

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

[19]  Richard R. Drake,et al.  A smoothed two- and three-dimensional interface reconstruction method , 2009 .

[20]  Aaron L. Fogelson,et al.  Fourth-Order Interface Tracking in Two Dimensions via an Improved Polygonal Area Mapping Method , 2014, SIAM J. Sci. Comput..

[21]  H. Marschall,et al.  voFoam - A geometrical Volume of Fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM , 2013, 1305.3417.

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

[23]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

[24]  Dick K. P. Yue,et al.  Conservative Volume-of-Fluid method for free-surface simulations on Cartesian-grids , 2010, J. Comput. Phys..

[25]  S. Zaleski,et al.  Interface reconstruction with least‐square fit and split Eulerian–Lagrangian advection , 2003 .

[26]  Qinghai Zhang,et al.  A new interface tracking method: The polygonal area mapping method , 2008, J. Comput. Phys..

[27]  S. Zaleski,et al.  A geometrical area-preserving volume-of-fluid advection method , 2003 .

[28]  Qinghai Zhang,et al.  On a Family of Unsplit Advection Algorithms for Volume-of-Fluid Methods , 2013, SIAM J. Numer. Anal..

[29]  D. Fletcher,et al.  A new volume of fluid advection algorithm: the defined donating region scheme , 2001 .

[30]  Dieter Bothe,et al.  On Adaptive Mesh Refinement for a 3D Geometrical Volume of Fluid Transport Algorithm on Unstructured Meshes using OpenFOAM , 2013 .

[31]  Stojan Petelin,et al.  Numerical errors of the volume‐of‐fluid interface tracking algorithm , 2002 .

[32]  Heinz Pitsch,et al.  Combination of 3D unsplit forward and backward volume-of-fluid transport and coupling to the level set method , 2013, J. Comput. Phys..

[33]  M. Renardy,et al.  PROST: a parabolic reconstruction of surface tension for the volume-of-fluid method , 2002 .

[34]  Sandro Manservisi,et al.  A geometrical predictor-corrector advection scheme and its application to the volume fraction function , 2009, J. Comput. Phys..

[35]  M. Rudman,et al.  A 3D unsplit-advection volume tracking algorithm with planarity-preserving interface reconstruction , 2006 .

[36]  J. Marsden,et al.  Tricubic interpolation in three dimensions , 2005 .

[37]  Y. L. Zhang,et al.  3D Impact and Toroidal Bubbles , 2001 .

[38]  R. LeVeque.  High-resolution conservative algorithms for advection in incompressible flow , 1996 .

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

[40]  Qinghai Zhang,et al.  On Donating Regions: Lagrangian Flux through a Fixed Curve , 2013, SIAM Rev..

[41]  Eugenio Aulisa,et al.  Interface reconstruction with least-squares fit and split advection in three-dimensional Cartesian geometry , 2007, J. Comput. Phys..

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

[43]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[44]  R. Scardovelli,et al.  A surface marker algorithm coupled to an area-preserving marker redistribution method for three-dimensional interface tracking , 2004 .

[45]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[46]  Iztok Žun,et al.  A PLIC-VOF method suited for adaptive moving grids , 2011, J. Comput. Phys..

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

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

[49]  E. Puckett,et al.  Second-Order Accurate Volume-of-Fluid Algorithms for Tracking Material Interfaces , 2013 .

[50]  J. O´Rourke,et al.  Computational Geometry in C: Arrangements , 1998 .

[51]  Olivier Desjardins,et al.  A computational framework for conservative, three-dimensional, unsplit, geometric transport with application to the volume-of-fluid (VOF) method , 2014, J. Comput. Phys..

[52]  J. Hattel,et al.  A Two-Phase Flow Solver for Incompressible Viscous Fluids, Using a Pure Streamfunction Formulation and the Volume of Fluid Technique , 2014 .

[53]  Shamit Bakshi,et al.  Noniterative interface reconstruction algorithms for volume of fluid method , 2013 .