“…The aforementioned drawbacks of Lagrangian methods have made Eulerian methods, such as the Volumeof-Fluid (VOF) method [20,32] or the Level-Set (LS) method [30,33,37,27,29,34] (or combinations thereof [36,39]) to be preferred in the modeling of complex interfaces undergoing topological changes such as bubble coalescence, wave breaking, etc [40,38,47,8,9]. Eulerian methods are based on the advection of a scalar field φ defined on the whole flow domain.…”