A new model for turbulent fibre suspension flow is proposed by introducing a model for the fibre orientation distribution function (ODF). The coupling between suspended fibres and the fluid momentum is then introduced through the second and fourth order fibre orientation tensors, respectively. From the modelled ODF, a method to construct explicit expressions for the components of the orientation tensors as functions of the flow field is derived. The implementation of the method provides a fibre model that includes the anisotropic detail of the stresses introduced due to presence of the fibres, while being significantly cheaper than solving the transport of the ODF and computing the orientation tensors from numerical integration in each iteration. The model was validated and trimmed using experimental data from flow over a backwards facing step. The model was then further validated with experimental data from a turbulent fibre suspension channel flow. Simulations were also carried out using a Bingham viscoplastic fluid model for comparison. The ODF model and the Bingham model performed reasonably well for the turbulent flow areas, and the latter model showed to be slightly better given the parameter settings tested in the present study. The ODF model may have good potential, but more rigorous study is needed to fully evaluate the model.