A numerical study has been conducted to evaluate the performance of the k-e turbulence model for axisymmetric, unconfined, swirling and nonswirling c o a x i a l jet flows. Two lower order schemes, hybrid and power-law, and two higher order schemes, flux-spline and hounded skew upwind differencing, were employed in this investigation. The predicted results indicate that the higher order numerical schemes have greater potential for future model improvements and complex flow calculations in terms of storage, accuracy, and execution time. For the nonswirling flow, computations using an algebraic stress model (ASM) have also been presented.Nomrnclature anisotropy tensor coefficients in the turbulence a i j cl, cp, cEl, cE2. . moaei antisymmetric part of the production rate of Reynolds stress tenturbulent kinetic energy characteristic turbulence length scale static pressure Peclet number symmetric Part of the production rate of Reynolds stress tensor production rete of turbulent kinetic energy Reynolds number strain rate tensor fluctuating velocity component (i = 1, 2, 3 ) time-averaged velocity component (i = 1, 2, 3) sor coefficient in the turbulence model coefficient in the turbulence model coefficient in the turbulence model Kronecker delta dissipation rate of turbulent kinetic energy dynamic ViSCosfty Of fluid turbulent viscosity density of fluid turbulent Prandtl number for k 1 turbulent Prandtl number for E pressure-strain term