摘要

A theoretical model of turbulent fiber suspension is developed by deriving the equations of Reynolds averaged Navier-Stokes, turbulence kinetic energy and turbulence dissipation rate with the additional term of fibers. In order to close the above equations, the equation of probability distribution function for mean fiber orientation is also derived. The theoretical model is applied to the turbulent channel flow and the corresponding equations are solved numerically. The numerical results are verified by comparisons with the experimental ones. The effects of Reynolds number, fiber concentration and fiber aspect-ratio on the velocity profile, turbulent kinetic energy and turbulent dissipation rate are analyzed. Based on the numerical data, the expression for the velocity profile in the turbulent fiber suspension channel flow, which includes the effect of Reynolds number, fiber concentration and aspect-ratio, is proposed.