Numerical investigation of three-dimensional fiber suspension flow by using finite volume method

Fiber suspension flow is common in many industrial processes like papermaking and fiber-reinforcing polymer-based material forming. The investigation of the mechanism of fiber suspension flow is of significant importance, since the orientation distribution of fibers directly influences the mechanical and physical properties of the final products. A numerical methodology based on the finite volume method is presented in the study to simulate three-dimensional fiber suspension flow within complex flow field. The evolution of fiber orientation is described using different formulations including FT model and RSC model. The pressure implicit with splitting of operators algorithm is adopted to avoid oscillations in the calculation. A laminate structure of fiber orientation including the shell layer, the transition layer and the core layer along radial direction within a center-gated disk flow channel is predicted through a three-dimensional simulation, which agrees well with Mazahir’s experimental results. The evolution of fiber orientation during the filling process within the complex flow field is further discussed. The mathematical model and numerical method proposed in the study can be successfully adopted to predict fiber suspension flow patterns and hence to reveal the fiber orientation mechanism.