A SLIC–VOF Method Based on Unstructured Grid

Abstract

In this paper, an interface reconstruction technology similar with SLIC method on structured grid is developed for the unstructured grid. For a cell containing the interface, a straight line segment parallel with one edge of the triangle cell is set to the interface according to the VOF function values in the three neighbor cells and volume flux of the particular fluid is determined by the geometry relation between the interface line and volume flux on each edge. Finally, four classical cases are used to test this unstructured SLIC-VOF interface reconstruction method including advection test of a right-angled triangle and a hollowed circle, Zalesak slotted disk rotation test as a test of scalar advection, and single-vortex shearing flow test, which are effective to test interface capturing method. It was found that the interface transported back by the reversed velocity field is comparative with SLIC-VOF on structure grid. Through the comparison between two triangle grids generated by Delaunay Triangulation Method and Bubble Packing Method, it is found that the grid quality has great effect on the proposed unstructured SLIC-VOF method. More regular triangle cells will lead to better interface capturing result.

A numerical investigation of various phase change models on simulation of saturated film boiling heat transfer

“…These methods keep the interface sharp and accurate. Simple line interface calculation (SLIC) method estimates the interface as straight horizontal and vertical lines in each cell. Piecewise linear interface construction (PLIC) and least squares volume‐of‐fluid interface reconstruction algorithm (LVIRA) by Pilliod and Puckett are other reconstruction methods which separate the two phases in each cell on interface region by an inclined piecewise linear line.…”

A mass‐conserving level‐set method for simulation of multiphase flow in geometrically complicated domains

“…The interface location is only approximately defined: it will intersect those cells that have intermediate values of the VoF fraction between 0 and 1, but its exact location is unknown. In each time step, the position of the interface is reconstructed within the control volumes that are intersected by the interface, under the assumption that it is either aligned with one of the Cartesian coordinate directions [5] or simply planar [6], before it is advanced in time in a Lagrangian manner. Computation of the fluxes across the cell faces is performed using the concept of donor regions, leading to an algorithm that is mass conserving to machine precision.…”