In the present work, an algorithm for the numerical solution of the two-dimensional mass-weighted, ensembleaveraged, compressible Navier-Stokes equations is presented. Turbulence closure for the Reynolds stresses is obtained using a low-Reynolds-number two-equation k-ε model, which includes near-wall effects. The same algorithm is used for the integration of both the Navier-Stokes and the turbulence-transport equations. The algorithm consists of a basic explicit finite-volume scheme, whose convergence is accelerated using local time-stepping and multiple-grid techniques. Convergence is further enhanced by the use of an implicit-residual-smoothing technique, applied both in the fine and the coarser grids. The application of the numerical scheme to the solution of the k-ε equations is examined in detail. Comparisons of computations with experimental data are presented for three shock-wave/turbulent boundary-layer interaction flows. Algorithm convergence rate and CPU-time requirements are discussed. © 1990 American Institute of Aeronautics and Astronautics, Inc., All rights reserved.