“…The accurate representation of surface tension within these methods is typically more delicate than for schemes using an explicit representation of the interface [14,15,16,17] and care must be taken to avoid any imbalance between discrete surface tension and pressure gradient terms [18,19,15,20,21,22,23]. For vof methods the accurate evaluation of the geometrical properties of the interface such as curvature is also an active research topic [24,21,25,23,26]. Finally dynamic mesh-adaptive methods can deal efficiently with phenomena involving a wide range of spatial scales [27,28,29,6,7,30,8,31,17,9,32,33,34].…”