The Reynolds averaged Navier-Stokes equation was solved numerically with the Reynolds stress model to get the mean fluid velocity and the turbulent kinetic energy in a round turbulent jet of fiber suspension. The fluctuating fluid velocity was described as a Fourier series with random coefficients. Then the slender-body theory was used to calculate the fiber orientation distribution and orientation tensor. Numerical results of mean axial velocity and turbulent shear stress along the lateral direction were validated by comparing with the experimental ones. The results show that most fibers are aligned with the flow direction as they go downstream, and fibers are more aligned with the flow direction within the region near the jet core. The fibers with high aspect ratio tend much easier to align with the flow direction, and the fiber orientation distribution is not sensitive to fiber aspect ratio when fiber aspect ratio is larger than 5. Fiber density has no obvious influence on the fiber orientation distribution and fiber orientation tensor. The randomizing effect of turbulence is insignificant in the regions near outside and jet core, and becomes significant in the region between outside and jet core. The randomizing effect increases with the increasing of the distance from the jet exit. Different components of fiber orientation tensor show a similar distribution pattern. (C) 2011 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved.