The numerical simulation of interfacial and free surface ows is a vast and interesting topic in the areas of engineering and fundamental physics, such as the study of liquid-gas interfaces, formation of droplets, bubbles and sprays, combustion problems with liquid and gas reagents, study of wave motion and others. One of the most powerful and robust methods for interface simulation in xed grids is the Volume-of-Fluid (VOF). In this method, the fluids are represented by a scalar fi eld Ck, known as volume fraction, that represents the portion of volume lled with fluid k. Given a velocity fi eld, interfaces are then tracked by evolving fluid volumes in time. At any time in the solution, an exact interface location is not known. Interface geometry is instead inferred (based on assumptions of the particular algorithm) and its location is reconstructed from local volume fraction data (Interface Reconstruction). The reconstructed interface is then used to compute the volume fluxes necessary for the convective terms in the volume evolution equation (Advection). The objective of this work is to implement a fast, accurate and parallelizable VOF method well suited to 3D unstructured meshes. The selected interface reconstruction algorithms will be the Youngs' (fi rst order) and the LVIRA (second order). In the other hand, the advection step will be computed by the means of an unsplit-advection volume tracking algorithm. In the paper, the VOF method will be tested for different test problems. First, a study of reconstruction accuracy, it is most easily assessed by analyzing the reconstruction of known geometries, such as a hollowed sphere. Second, a rotation test, where a velocity field is imposed and the advection algorithm is tested.

Peer Reviewed