Numerical method for calculation of surface tension flows in arbitrary grids

We outline a scheme to define and track interfaces in an arbitrary quadrilateral grid by making use of a single fractional volume if fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid-particles flowfield. The present methodology, however, represents a substantial refinement over the original VOF programs in that nonorthogonal, body-fitted, and even deforming grids can be used, and the gas and liquid phased can take on totally different properties, including compressibility and incompressibility.