The current study relates to the development of a multiphase model of flexible fibre suspensions. An understanding of the rheology and dynamics of the deformation of such suspension is desirable in order to be able to fully disclose the flow behaviour from very low to very high shear rates. We present an approach for numerically simulating the dynamics of flexible fibres employing a particle-level method. This is performed by investigating the fibre dynamics against several orbit classes - i.e. rigid, springy, flexible and complex rotation of the fibres [1-3] enabling the model to have all degrees of freedom (translation, rotation, bending and twisting). The three-dimensional Navier-Stokes equations which describes the fluid motion are employed while the fibrous phase of the fluid is modeled as chains of fiber segments interacting with the fluid through viscous- and drag forces. The simulations are performed using OpenFOAM and the numerical outcomes are validated against experimental data.The purpose of the modelling framework applied in this work is to enable the numerical model to be extended to a 4-way coupling model, capturing shear thinning, shear thickening and the yield stress properties of a fibrous fluid suspension.