The fluid particle dynamics method is employed to study the fluid-solid interaction problem, which can avoid the explicit implementation of fluid particle boundary condition and capture the hydrodynamic interaction well. We solve the incompressible Navier-Stokes equation coupled with the rigid body motion equation in polar or cylindrical coordinates. A pressure stabilization scheme is used to solve the system in polar coordinates for two dimensional case and cylindrical coordinates for three dimensional case. Our objective is to understand numerically the fluid-solid interaction in rotational extrusion flow. We numerically verify the correctness of method presented here and give comparative analysis for different parameters. We present the Jeffery orbit formulation in annular region. Numerical experiments show that the fluid particle dynamics method is reliable and efficient for numerical simulation of particulate flow in cylindrical coordinate system.