Abstract Two nonlinear eddy viscosity models (EVMs) proposed by Shih et al. (Shih, T.H., Zhu, J., Lumley, J.L., 1995. Comput. Methods Appl. Mech. Engrg. 125, 287–302) and Craft et al. (Craft, T.J., Launder, B.E., Suga, K., 1996. Int. J. Heat Fluid Flow 17, 108–115) and one algebraic Reynolds stress model of Gatski and Speziale (Gatski, T.B., Speziale, C.G., 1993. J. Fluid Mech. 254, 59–78) are compared to the standard k–ϵ model and the Reynolds stress model of Launder et al. (Launder, B.E., Reece, G.J., Rodi, W., 1975. J. Fluid Mech. 68, 537–566) with respect to model accuracy and computational robustness. Three flow situations are considered to analyze the models' capabilities in flows with streamline curvature: A backward facing step flow (Driver, D.M., Seegmiller, H.L., 1985. AIAA J. 23, 163–171), a plane U-duct flow (Monson, D.J., Seegmiller, H.L., 1992. NASA Technical Memorandum, 103931) and a curved mixing layer (Castro, I.P., Bradshaw, P., 1976. J. Fluid Mech. 73, 265–304). Results indicate that an improved prediction of the flow fields especially for important flow features like recirculation zones can be obtained with the tested models compared to calculations with the standard k–ϵ model. The decrease of computational robustness compared to the calculations with the standard k–ϵ model is rather moderate once a quasi-linear form of the models is employed.