SUMMARYThis paper presents a ÿnite di erence technique for solving incompressible turbulent free surface uid ow problems. The closure of the time-averaged Navier-Stokes equations is achieved by using the twoequation eddy-viscosity model: the high-Reynolds k-(standard) model, with a time scale proposed by Durbin; and a low-Reynolds number form of the standard k-model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non-linear terms, a second=third-order upwinding technique is adopted. The computational method is validated by applying it to the at plate boundary layer problem and to impinging jet ows. The method is then applied to a turbulent planar jet ow beneath and parallel to a free surface. Computations show that the high-Reynolds k-model yields favourable predictions both of the zero-pressure-gradient turbulent boundary layer on a at plate and jet impingement ows. However, the results using the low-Reynolds number form of the k-model are somewhat unsatisfactory.