Three low-Reynolds turbulence models are implemented in an unstructuredgrid, Navier-Stokes, finite-volume solver. These are the one-equation model of Spalart- Allmaras, the Shear-Stress-Transport k - ω model and a hybrid one/two-equation k - ε model. All of them are shown to be robust and compliant with unstructured grids. They are used to numerically predict the subsonic flow around two-dimensional high-lift configurations, namely the two-element NLR-7301 airfoil and the single-element A-airfoil, where turbulent separation occurs near the trailing edge at high incidence angles. Both global coefficient distributions (Cp, Cf , C l, Cd) and local boundary layer quantities or velocity profiles are compared with measured data. Due to the very stretched grid cells used close to solid boundaries, the definition of the finite-volumes (median-dual or containment circle-dual tesselation of the domain) is crucial for the quality of the results. Besides, the effect of including or not a transition mechanism into the turbulence model is investigated. The results from all models are satisfactory; the Spalart-Allmaras model performs slightly better than the other two models in terms of both predictive capability and numerical robustness.