Modeling surface tension phenomena is computationally difficult because it requires the evaluation of second derivatives.
This is particulary true in the FLOW-3D program where the capability to represent highly complicated and multiple free surfaces difficulties are further compounded in three-dimensional calculations because one is often forced, for reasons of economy, to use marginal numerical resolution.