Dynamic simulation of microstructure and rheology of fiber suspensions

A numerical simulation method has been developed to predict the motion of fiber suspensions in various flows by using the particle simulation method (PSM), in which a fiber is modeled by an array of spheres. The hydrodynamic interaction among fibers is considered by decomposing into intra-fiber and inter-fiber interactions. For the intra-fiber case, the many-body problem is solved by calculating the mobility matrix for each fiber to obtain the hydrodynamic force and torque exerted on each sphere. For the inter-fiber case, only the near-field lubrication force is considered between spheres belonging to different fibers. The methodology was applied to predicting the microstructure and rheological properties of rigid fiber suspensions in simple shear flow and fiber orientation in two-dimensional diverging flow by combining this method with the finite element method (FEM). In the former, the overshoot of suspension viscosity, due to the transient change of the microstructure, was observed in semi-dilute to concentrated systems. In the latter, the calculated fiber orientations agreed with expectations of theory.

[1]  J. Agassant,et al.  Experimental and theoretical study of short fiber orientation in diverging flows , 1985 .

[2]  S. Güçeri,et al.  On the Description of the Orientation State for Fiber Suspensions in Homogeneous Flows , 1989 .

[3]  A. B. Metzner,et al.  The rheology, fiber orientation, and processing behavior of fiber‐filled fluids , 1992 .

[4]  C. L. Tucker,et al.  Orientation Behavior of Fibers in Concentrated Suspensions , 1984 .

[5]  P. Cladis,et al.  Stability of Nematic Liquid Crystals in Couette Flow , 1975 .

[6]  E. J. Hinch,et al.  Constitutive equations in suspension mechanics. Part 2. Approximate forms for a suspension of rigid particles affected by Brownian rotations , 1976, Journal of Fluid Mechanics.

[7]  Suresh G. Advani,et al.  Closure approximations for three‐dimensional structure tensors , 1990 .

[8]  T. Matsuoka,et al.  Dynamic simulation of flow‐induced fiber fracture , 1995 .

[9]  G. B. Jeffery.  The motion of ellipsoidal particles immersed in a viscous fluid , 1922 .

[10]  Satoru Yamamoto,et al.  A method for dynamic simulation of rigid and flexible fibers in a flow field , 1993 .

[11]  E. Guyon,et al.  Two Shear-Flow Regimes in Nematic p-n-Hexyloxybenzilidene-p'- aminobenzonitrile , 1974 .

[12]  T. Matsuoka,et al.  Viscosity of dilute suspensions of rodlike particles: A numerical simulation method , 1994 .

[13]  Ronald G. Larson,et al.  Effect of molecular elasticity on out-of-plane orientations in shearing flows of liquid-crystalline polymers , 1991 .

[14]  S. Onogi,et al.  The rheology of suspensions of titanate fibers in polymer solution , 1977 .

[15]  Suresh G. Advani,et al.  The Use of Tensors to Describe and Predict Fiber Orientation in Short Fiber Composites , 1987 .