A numerical model for the simulation of turbulent free surface flows is proposed. This model is based on large-eddy simulations (LES) and on a marker-variable based Eulerian interface capturing strategy. The LES model is introduced in the context of free surface flows problems, using the Wall-adapting local eddy-viscosity (WALE) model for the turbulent viscosity. The filtered free surface capturing model is also introduced. The proposed methodology is validated based on two well-documented free surface flow problems where turbulence is observed: the Rayleigh-Taylor instability and the two-phase transient mixing layer problem.