Computer simulations are important to increase the knowledge of the processes involved in paper making. Because of the complexity of the paper processes the development of simulation tools requires models and numerical methods which capture as much as possible of the real physical phenomena. In this paper a framework for simulation of fiber suspension flow is presented. The framework is founded on a fluid solver for the Navier-Stokes equations and a fiber model based on finite strain beam theory including shearing. Moreover, the coupling between fibers and fluid is resolved with the immersed boundary method and the interaction between the fibers is calculated using a model based on DLVO-theory.