A two-dimensional Navier-Stokes code has been developed for rapid numerical simulation of axisymmebic flow fields, including flow fields with an azimuthal velocity component. The azimuthal-invariant Navier-Stokes equations in a cylindrical coordinate system are mapped to a general body-fitted coordinate system, with the streamwise viscous terms then neglected by applying the thinlayer approximation. Turbulence effects are modeled using an algebraic model, typically the Baldwin-Lomax turbulence model, although a modified Cebeci-Smith model can also be used. The equations are discretized using central finite differences and solved using a multistage Runge-Kutta algorithm with a spatially-varying time step and implicit residual smoothing.Results are presented for calculations of supersonic flow over a waisted body-of-revolution, transonic flow through a normal shock wave in a straight circular duct of constant cross-sectional area, swirling supersonic (inviscid) flow through a strong shock in a straight radial duct, and swirling subsonic flow in an annular-to-circular diffuser duct. Comparisons between computed and experimental results a r e in fair to good agreement, demonstrating that the viscous code can be a useful tool for practical engineering design and analysis work.