Use of high accuracy schemes to handle free surfaces in computing unsteady two-phase flows