Implementation of three different Navier–Stokes solvers in an isogeometric finite element framework is presented in this paper. The three solvers Chorin projection method and Coupled formulation , both with the Spalart–Allmaras turbulence model, and Variational Multiscale (VMS) method have been applied to simulate flow past a two-dimensional NACA0012 airfoil at a high Reynolds number (Re=3×106Re=3×106) for four different angles of attack. The predicted flow characteristics have been compared and the effects of increasing the order of the spline element on the accuracy of prediction and computational efficiency are evaluated. In this study it turns out that flow separation does not take place up to an angle of attack of 16°. Up to this angle of attack all three solvers predict similar results in good agreement with each other and with available experimental results. However, a big spread in lift and drag coefficients is observed in the stall regime. Our study also shows that for linear spline elements all three solvers are computationally similar. For quadratic spline elements the Chorin solver compares favorably to the others based on the results presented here.