A normalization scheme for the numerical solution of the moment approximation equation in fiber suspension flows is presented. Here, normalization refers to rescaling the trace of the second moment tensor to unity at each time step. The equivalence between the normalization scheme and the quadratic closure model is analytically proved. The performance of the scheme is investigated in simple shear flow with respect to the quadratic and hybrid closures, and a stochastic Monte-Carlo simulator that provides exact solution. The proposed scheme is a computationally efficient alternative to the quadratic closure: it performs equally well and is more efficient regarding computational time. (c) 2010 Elsevier B.V. All rights reserved.