Three-dimensional upwinding Navier-Stokes code with k-epsilon model for supersonic flows