A numerical model based on the Fokker–Planck equation for studying the orientation of paper fibers in a suspension is presented. The model is implemented in a finite volume framework. The rotational operators are discretized on a triangulated half sphere and the spatial operators are discretized on a Cartesian octree grid, where internal walls are handled with immersed boundary methods. Transient solutions of the two dimensional fiber orientation distribution are computed in the entire fluid domain. This is in contrast to earlier approaches that either only compute along streamlines or one dimensional projected orientations. The fibers are assumed to have negligible inertia and thus to follow the fluid perfectly and to be evenly distributed. A head-box geometry relevant for paper forming is used to exemplify the usefulness of the method. The influence of translational and rotational diffusion is included in the presented method.