The performance of several turbulence models in computing an axisymmetric supersonic base flow is investigated. A compressible Navier-Stokes code, which incorporates k-ε, k-ω model and Reynolds stress closure with three kinds of pressure-strain correlation model, has been developed using implicit LU-SGS algorithm with second-order upwind TVD scheme. Numerical computations have been carried out for Herrin and Duttons base flow. It is observed that the two-equation models give large backward axial velocity approaching to the base and somewhat larger variation of base pressure distribution than the Reynolds stress model. It is also found that the Reynolds stress model with third order pressure-strain model in the anisotropy tensor predicts most accurate mean flow field.