The study focuses on the fundamental nonorthogonal grid coverage of an axisymmetric flow domain with irregular boundaries, and involves incorporating swirl and turbulence effects into a stream function-vorticity simulation. Turbulence is simulated by way of a two-equation k-e turbulence model, together with wall functions for the treatment of near-wall flows. The finite difference equations are derived at points of a general boundary-fitted nonorthogonal grid using displaced, linear, and quadratic interpolation. The solution procedure is implemented in a new computer code which features generalized boundaries, rectangular and/or nonrectangular grids, swirl, turbulence, time-dependent calculations, and adaptive stability. Model evaluation is accomplished through a comparative study between predictions and available experimental data.