There are rich vortex phenomena presented in swirling circulation flow, such as secondary flow and vortex breakdown. For numerical simulation of such complicated recirculating flow, higher requirements are needed for grid density and calculation precision. With the increase of grid points, both convergence rate and calculation precision of traditional iterative numerical methods decrease. On the basis of finite control volume scheme and pressure correction method, the FAS multi staggered grid scheme is presented here for Navier Stokes equation in cylinder coordination. The axisymmetric swirling flow confined in cylindrical container is used as test problem. The results indicate that both convergence rate and calculation precision are highly improved.