The compressible flow Navier-Stokes equations, in conservation form, have been solved using a modified MacCormack implicit predictor-corrector scheme. For turbulent flow the solution is modified by incorporating an eddy viscosity, to solve the Reynolds-averaged Navier-Stokes equations. The algorithm is accurate to second order in time and space and is unconditionally stable. It requires the solution of block bidiagonal matrices and the programme can run on commonly available minicomputers. It also has the advantage that a steady state solution is approached, if it exists, independent of the time step employed.
Numerical results are presented for the two-dimensional shock-boundary layer interaction in both laminar and turbulent flow and the results are compared with experimental values and the numerical results of other workers. The method is calibrated against the analytical solution for the boundary layer on a flat plate and for the time dependent Couette flow problem.