The goal of this paper is to implement an accurate and robust solver for compressible Navier-Stokes equations coupled with the Spalart–Allmaras model, which possesses the capability of shock-capturing and predication of boundary layer and separated flow. In a given stencil width, a WENO-Z scheme equipped with Roe flux difference split method is used to calculate the inviscid flux, and central differencing scheme for the viscous terms are employed. The explicit Runge-Kutta is adopted for the temporal discretization. The simulation results of selected cases are given to verify the validation of the solver.