A numerical scheme for the prediction of free surface flows is presented and investigated. The method is based on an adaptive grid Eulerian finite-volume method, where non-orthogonal boundary-fitted moving grids are employed to follow the free surface. The underlying flow solver consists in a pressure-correction scheme of SIMPLE type with multigrid acceleration, which is iteratively combined with the moving grid technique. Several numerical examples are considered to illustrate the capabilities of the approach.