The majority of fluid flows that are interesting from a practical point of view are turbulent flows. The problem is that turbulence is particularly difficult to model. Although simulation of turbulent flows has been the topic of important researches, it remains an open issue. We present a stabilized finite element method combine with the one equation Spalart–Allmaras model. In opposition to what is usually done, there is no splitting between the Spalart–Allmaras equation and the Navier–Stokes equations. The aim is to built high order scheme for compressible turbulents flows in order to control numerical dissipation and get accurate solutions.