special care must be taken in computing , these fluxes to preserve the sharp definition of free boundaries. The method employed in SOLA-VOF uses a type of Donor-Acceptor flux approximation. The essential idea is to use information about F downstream as well as upstream of a flux boundary to establish a crude interface shape , and then to use this shape in computing the flux.