The present effort is the development of a multiscale modeling, simulation methodology for investigating complex phenomena arising from flowing fiber suspensions. The present approach is capable of coupling behaviors from the Kolmogorov turbulence scale through the full-scale system in which a fiber suspension is flowing. Here the key aspect is adaptive hierarchical modeling. Numerical results are presented for which focus is on fiber floc formation and destruction by hydrodynamic forces in turbulent flows. Specific consideration was given to dynamic simulations of viscoelastic fibers in which the fluid flow is predicted by a method that is a hybrid between direct numerical simulations and large eddy simulation techniques and fluid fibrous structure interactions will be taken into account. Dynamics of simple fiber networks in a shearing flow of water in a channel flow illustrate that the shear-induced bending of the fiber network is enhanced near the walls. Fiber-fiber interaction in straight ducts is also investigated and results show that deformations would be expected during the collision when the surfaces of flocs will be at contact. Smaller velocity magnitudes of the separated fibers compare to the velocity before collision implies the occurrence of an inelastic collision. In addition because of separation of vortices, interference flows around two flocs become very complicated. The results obtained may elucidate the physics behind the breakup of a fiber floc, opening the possibility for developing a meaningful numerical model of the fiber flow at the continuum level where an Eulerian multiphase flow model can be developed for industrial use.