Abstract: The orientation distribution function of cylindrical particle suspensions was deduced and numerically simulated, and an application was taken in a wedge-shaped flow field. The relationship between the orientation distribution function and particle orientation angles was obtained. The results show that comparing with the most probable angle distribution which comes to being in short time, the distribution of the steady state doesn’t vary much in range; the main difference is the anti-clockwise rotation in the right and upper field, that is, particles rotate more at the points where the velocity gradients are larger. The most probable orientations are close to the direction of local streamlines. In the direction of streamlines, with poleradius decreasing, the most probable angles increase, but the angles between their orientations and the local streamlines decrease.