Numerical simulation of buffet flows have been performed using various turbulence models, including linear and non-linear low-Re eddy-viscosity models. The accuracy of the models has been assessed against experimental data for transonic flows around the NACA-0012 aerofoil. Discretisation of the mean flow and turbulence transport equations is obtained by a Riemann solver up to third-order accurate. The time integration is obtained by an implicit unfactored method. The study shows that two-equation models associated with functional cμ coefficient in the calculation of the eddy-viscosity (henceforth labelled NL-cμ), provide better results. The Spalart-Allmaras one-equation model provides comparable results with the above models, while larger inaccuracies are presented in the case of linear and non-linear models with constant cμ coefficient. The buffet onset boundaries are similarly predicted by the one-equation and NL-cμ models.