Effects of bending and torsion rigidity on deformation and breakage of flexible fibers: A direct simulation study J. Chem. Phys. 136, 074903 (2012) Clouds of particles in a periodic shear flow Phys. Fluids 24, 021703 (2012) Steady-state hydrodynamics of a viscous incompressible fluid with spinning particles J. Chem. Phys. 135, 234901 (2011) Orientational order in concentrated suspensions of spherical microswimmers Phys. The orientation of fibers in simple two-dimensional flows is investigated. According to different ranges of the Péclet number, Pe, defined as the ratio of a characteristic rotational speed of fibers and the orientational diffusivity, three methods are developed: characteristic method for Pe= ϱ, regular perturbation method for Peӷ 1, and spectral method for everything else. All the methods subtly utilize the evolving solution of the rotational dynamics of fibers, which is also given in this paper. Especially, the adoption of spherical harmonics in the spectral method eliminates the singularity of the Fokker-Planck equation in spherical coordinates, and provides high precision and efficiency. The evolving solution of orientation distribution with Pe= ϱ is obtained through the solution of rotational dynamics. Using a regular perturbation method, the solution of orientation distribution with Pe= ϱ is extended for the condition of Peӷ 1. This paper provides systematical and high efficient techniques to deal with the fiber orientation.