摘要

A successive iteration method is proposed to numerically simulate fiber suspensions. The computational field is discretized with the collocated finite volume method, and an ergodic hypothesis is adopted to greatly accelerate the solution to the Fokker-Planck equation. The method is employed in channel flows with different fiber volume fractions and aspect ratios, and its effectiveness is proved. The numerical results show that the existence of fibers significantly changes the pressure distribution, and the fiber aspect ratio has a greater effect on the velocity profile than on the volume faction. At the center of the channel, the velocity increases along the streamwise direction, while the velocity near the wall decreases slightly. The uncoupling and coupling solutions of the additional stress of the fiber suspensions are quite different.