The impetus of this paper is the comparative applications of two numerical schemes for supersonic flows using computational algorithms tailored for a supercomputer. The mathematical model is the conservation form of Navier-Stokes equations with the effect of turbulence being modeled algebraically. The first scheme is an implicit, unfactored, upwind-biased, line-Gauss-Seidel relaxation scheme based on finite-volume discretization. The second scheme is the explicit-implicit MacCormack scheme based on finite-difference discretization. The best overall efficiences are obtained using the upwind relaxation scheme. The integrity of the solutions obtained for the example cases is shown by comparisons with experimental and other computational results.