Today, Vortex Generators (VGs) are becoming an integral part of a Wind Turbine blade design. However, the challenges involved in the computation of the flow around VGs are yet to be dealt with in a satisfactory manner. A large number of VG models for Reynolds Averaged Navier Stokes (RANS) solvers has been proposed and, among them, the Bender–Anderson–Yagle (BAY) model (ASME Pap. FEDSM99-6919) is one of the most popular, due to its ease of use and relatively low requirements for user input. In the present paper a thorough investigation on the performance and application of the BAY model for aerodynamic VG flows is presented. A fully resolved RANS simulation is validated against experiments and then used as a benchmark for the BAY model simulations. A case relevant to wind turbines is examined, which deals with the flow past a wind turbine airfoil at Reynolds number 0.87e6. When the grid related errors are excluded, it is found that the generated vortices are weaker in the BAY model simulations than in the fully resolved computation. The latter finding is linked to an inherent deficiency of the model, which is first found in this study and which is explained in detail. •Most detailed analysis of the BAY model to date.•First time the BAY model was applied to a Fully Resolved grid.•Model inherent deficiency highlighted.•Identification of areas and ways to improve the model.