SUMMARYThis study extends the reduced Navier-Stokes (RNS) global pressure relaxation procedure developed by Rubin and co-workers for external flow to internal flow applications. The streamwise pressure gradient is split into a backward-differenced or initial value component, as in boundary layer marching, and a forwarddifferenced or boundary value component that represents the elliptic downstream effects. The streamwise convection terms are upwind-differenced and all other streamwise derivatives are backward-differenced. We thus obtain a standard boundary layer marching technique imbedded in a conventional line relaxation technique. For compressible flow the pressure iteration determines the interior flow interation as well as the inlet mass flux that is consistent with the outflow pressure boundary condition. Results have been computed for incompressible flow in both rectangular and curved channels, and for subsonic compressible flow in the simulation of an aerofoil in a wind tunnel. Converged solutions were obtained over a range of Reynolds numbers generating small to moderately large separation bubbles.