On the Theory and Computation of Surface Tension: The Elimination of Parasitic Currents through Energy Conservation in the Second-Gradient Method

Errors in the computation of fluid flows with surface tension are examined. These errors cause large parasitic flows when the capillary number is large and have often been attributed to truncation error in underresolved interfacial regions. A study using the second-gradient method reveals that when truncation error is eliminated in the computation of energy exchanges between surface and kinetic energies so that energy is strictly conserved, the parasitic currents are reduced to round-off. The results are based on general thermodynamic arguments and can be used to guide improvements in other methods, such as the continuum-surface-force (CSF) method, which is commonly used with the volume-of-fluid (VOF) method.