Skip to main content Accessibility help

A numerical study of the effect of fibre stiffness on the rheology of sheared flexible fibre suspensions

  • JINGSHU WU (a1) and CYRUS K. AIDUN (a1)


A recently developed particle-level numerical method is used to simulate flexible fibre suspensions in Newtonian simple shear flow. In this method, the flow is computed on a fixed regular ‘lattice’ using the lattice Boltzmann method, where each solid particle, or fibre in this case, is mapped onto a Lagrangian frame moving continuously through the domain. The motion and orientation of the fibre are obtained from Newtonian dynamics equations. The effect of fibre stiffness on the rheology of flexible fibre suspensions is investigated and a relation for the relative viscosity is obtained. We show that fibre stiffness (bending ratio, BR) has a strong impact on rheology in the range BR < 3. The relative viscosity increases significantly as BR decreases. These results show that the primary normal stress difference has a minimum value at BR ~ 1. The primary normal stress difference for slightly deformable fibres reaches a minimum and increases significantly as BR decreases below one. The results are explained based on Batchelor's relation for non-Brownian suspensions.


Corresponding author

Email address for correspondence:


Hide All
Aidun, C. K., Lu, Y. & Ding, E. 1998 Direct analysis of particulate suspensions with inertia using the discrete Boltzmann equation. J. Fluid Mech. 373, 287311.
Batchelor, G. K. 1971 The stress generated in a non-dilute suspensions of elongated particles by pure straining motion. J. Fluid Mech. 46, 813829.
Bibbo, M. A. 1987 Rheology of semi-concentrated fiber suspensions. PhD thesis, Massachusetts Institute of Technology, Cambridge, MA.
Bibbo, M. A., Dinh, S. M. & Armstrong, R. C. 1985 Shear flow properties of semi-concentrated fibre suspensions. J. Rheol. 29, 905929.
Blakeney, W. R. 1966 The viscosity of suspensions of straight, rigid rods. J. Colloid Interface Sci. 22, 324330.
Carter, L. F. 1967 A study of the rheology of suspensions of rod-shaped particles in a Navier–Stokes liquid. PhD thesis, University of Michigan, Ann Arbor, MI.
Claeys, I. L. & Brady, J. F. 1993 Suspensions of prolate spheroids in Stokes flow. Part 1. Dynamics of a finite number of particles in an unbounded fluid. J. Fluid Mech. 251, 411442.
Cox, R. G. 1971 The motion of long slender bodies in a viscous fluid. Part 2. Shear flow. J. Fluid Mech. 45, 625657.
Doi, M. & Edwards, S. F. 1978 Dynamics of concentrated polymer systems. Part I. Brownian motion in equilibrim state. J. Chem. Soc. Faraday Trans. II 74, 17891801.
Fan, X., Phan-Thien, N. & Zheng, R. 1998 A direct simulation of fibre suspensions. J. Non-Newtonian Fluid Mech. 74, 113135.
Forgacs, O. L. & Mason, S. G. 1959 Particle motions in sheared suspensions. Part 9. Spin and deformation of thread-like particles. J. Colloid Interface Sci. 14, 457472.
Goldsmith, H. L. & Mason, S. G. 1967 The microrheology of dispersions. Rheology 4, 85201.
Goto, S., Nagazono, H. & Kato, H. 1986 The flow behaviour of fiber suspensions in Newtonian fluids and polymer solutions. Part I. Mechanical properties. Rheol. Acta 25, 119129.
Jeffery, G. B. 1922 The motion of ellipsoidal particles immersed in a viscous fluid. Proc. R. Soc. Lond. A 102, 161179.
Joung, C. G., Phan-Thien, N. & Fan, X. J. 2001 Direct simulation of flexible fibers. J. Non-Newtonian Fluid Mech. 99, 136.
Keshtkar, M., Heuzey, M. C. & Carreau, P. J. 2009 Rheological behavior of fiber–filled model suspensions: Effect of fiber flexibility. J. Rheol. 53, 631650.
Kitano, T. & Kataoka, T. 1981 The rheology of suspensions of vinylon fibers in polymer liquids. Part I. Suspensions in silicone oil. Rheol. Acta 20, 390402.
Lees, A. W. & Edwards, S. F. 1972 The computer study of transport processes under extreme conditions. J. Phys. C 5, 19211928.
Lindstrom, S. B. & Uesaka, T. 2008 Simulation of semi-dilute suspensions of non-Brownian fibers in shear flow. J. Chem. Phys. 128, 024901.
MacMeccan, R. M., Clausen, J. R., Neitzel, G. P. & Aidun, C. K. 2009 Simulating deformable particle suspensions using a coupled lattice-Boltzmann and finite-element method. J. Fluid Mech. 618, 1339.
Petrich, M. P., Koch, D. L. & Cohen, C. 2000 An experimental determination of the stress-microstructure relationship in semi-concentrated fiber suspensions. J. Non-Newtonian Fluid Mech. 95, 101133.
Ross, R. F. & Klingenberg, D. J. 1997 Dynamic simulation of flexible fibers composed of linked rigid bodies. J. Chem. Phys. 106, 29492960.
Schmid, C. F., Switzer, L. H. & Klingenberg, D. J. 2000 Simulation of fiber flocculation: effects of fiber properties and interfiber friction. J. Rheol. 44, 781809.
Sepehr, M., Carreau, P. J., Moan, M. & Ausias, G. 2004 Rheological properties of short fiber model suspensions. J. Rheol. 48, 10231048.
Stover, C. A., Koch, D. L. & Cohen, C. 1992 Observations of fiber orientation in simple shear flow of semi-dilute suspensions. J. Fluid Mech. 238, 277296.
Wagner, A. J. & Pagonabarraga, I. 2002 Lees–Edwards boundary conditions for lattice Boltzmann. J. Stat. Phys. 107, 521537.
Wu, J. & Aidun, C. K. 2010 a A method for direct simulation of flexible fiber suspensions using lattice-Boltzmann equation with external boundary force. Intl J. Multiphase Flow 36, 202209.
Wu, J. & Aidun, C. K. 2010 b Simulating 3d deformable particle suspensions using lattice B oltzmann method with discrete external boundary force. Intl J. Numer. Method Fluids 62, 765783.
Yamamoto, S. & Matsuoka, T. 1992 A method for dynamic simulation of rigid and flexible fibers in a flow field. J. Chem. Phys. 98, 644650.
Yamane, Y., Kaneda, Y. & Doi, M. 1994 Numerical simulation of semi-dilute suspensions of rod-like particles in shear flow. J. Non-Newtonian Fluid Mech. 54, 405421.
Zirnsak, M. A., Hur, D. U. & Boger, D. V. 1994 Normal stresses in fibre suspensions. J. Non-Newtonian Fluid Mech. 54, 153193.
MathJax is a JavaScript display engine for mathematics. For more information see

JFM classification

A numerical study of the effect of fibre stiffness on the rheology of sheared flexible fibre suspensions

  • JINGSHU WU (a1) and CYRUS K. AIDUN (a1)


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed