The present paper describes an evaluation of three Navier-Stokes computational fluid dynamics codes simulating supersonic retropropulsion flow fields from a Ames 9'7' Unitary Tunnel entry. Three codes-DPLR, FUN3D, and OVERFLOW-have been exercised (using lessons learned from a previous round of simulations) for both single-and multi-nozzle configurations for two Mach numbers and three thrust coefficients, all at zero degrees angle of attack. The focus of the present work is on high thrust coefficients and low supersonic Mach numbers, which were not considered by a previous supersonic retropropulsion wind tunnel test performed in the Langley Unitary Plan Wind Tunnel. Surface pressure measurements and shadowgraphs have been used to evaluate the flow prediction tools. All codes predict pressure measurements to within ±1.8% for the smooth (nonozzle) configuration. The effect of the sting mount was also simulated to access its influence on predicted flow fields. The single-and three-nozzle configurations show periodic oscillation in triple point regions but little variation in the bow shock and plume termination shock features. The three-nozzle configuration becomes steadier with increasing thrust coefficient while the four-nozzle configurations become more chaotic. The predictions of unsteadiness and periodicity, however, do not seem to effect the predictions of time-averaged surface pressure coefficients.