The prediction of fiber orientation is critical in short-fiber-reinforced composite systems. Fiber orientation during the manufacturing process largely determines the mechanical, electrical, and thermal properties of the final product. For example, a product with aligned fibers, i.e. fibers mostly oriented in the same direction, has much better mechanical properties than that with random fiber orientation. In this paper, a three-dimensional finite element simulation tool is developed to study fiber suspensions in various viscous flows, including homogeneous flows and inhomogeneous flows (such as Poiseuille flows). Instead of using commercial finite element solver, we propose a finite element-based simulation package coupled with the analytical Jacobian matrix in the Newton–Raphson iteration in solving fiber’s velocities. The proposed method is extended to investigate the hydrodynamic interaction between fibers, and the change of fiber motion and fluid rheological properties.
JackDASmithDE. Elastic properties of short-fiber polymer composites, derivation and demonstration of analytical forms for expectation and variance from orientation tensors. J Compos Mater2008; 42: 277–308.
2.
TuckerCL IIILiangE. Stiffness predictions for unidirectional short-fiber composites: Review and evaluation. Compos Sci Technol1999; 59: 655–671.
3.
JefferyGB. The motion of ellipsoidal particles immersed in a viscous fluid. Proc Roy Soc Lond Ser A1922; 102: 161–179.
4.
SegreGSilberbergA. Radial particle displacements in Poiseuille flow of suspensions. Nature1961; 189: 209–210.
5.
GoldsmithHLMasonSG. Axial migration of particles in Poiseuille flow. Nature1961; 194: 1095–1096.
6.
OliverDR. Influence of particle rotation on radial migration in the Poiseuille flow of suspensions. Nature1962; 194: 1269–1271.
7.
KarnisAGoldsmithHLMasonSG. Axial migration of particles in Poiseuille flow. Nature1963; 200: 159–160.
8.
RepettiRVLeonardEF. Segré-Silberberg annulus formation: a possible explanation. Nature1964; 203: 1346–1348.
9.
KarnisAGoldsmithHLMasonSG. The flow of suspensions through tubes, v. inertial effects. Can J Chem Eng1966; 44: 181–193.
10.
MatasJPMorrisJFGuazzelliÉ. Inertial migration of rigid spherical particles in Poiseuille flow. J Fluid Mech2004; 515: 171–195.
11.
MiuraKItanoTSugihara-SekiM. Inertial migration of neutrally buoyant spheres in a pressure-driven flow through square channels. J Fluid Mech2014; 749: 320–330.
12.
RubinowSIKellerJB. The transverse force on a spinning sphere moving in a viscous fluid. J Fluid Mech1961; 11: 447–459.
13.
SaffmanPG. The lift on a small sphere in a slow shear flow. J Fluid Mech1965; 22: 385–400.
14.
CoxRGBrennerH. The lateral migration of solid particles in Poiseuille – i theory. Chem Eng Sci1968; 23: 147–17.
15.
HoBPLealLG. Inertial migration of rigid spheres in two-dimensional unidirectional flows. J Fluid Mech1974; 65: 365–400.
16.
BrethertonFP. The motion of rigid particles in a shear flow at low Reynolds number. J Fluid Mech1962; 14: 284–304.
17.
ChwangAT. Hydromechanics of low-Reynolds-number flow. Part 3. Motion of a spheroidal particle in quadratic flows. J Fluid Mech1975; 72: 17–34.
18.
StoverCACohenC. The motion of rodlike particles in the pressure-driven flow between two flat plates. Rheol Acta1990; 29: 192–203.
19.
Sugihara-SekiM. The motion of an elliptical cylinder in channel flow at low Reynolds numbers. J Fluid Mech1993; 257: 575–596.
20.
Sugihara-SekiM. The motion of an ellipsoid in tube flow at low Reynolds numbers. J Fluid Mech1996; 324: 287–308.
21.
FengJHuHHJosephDD. Direction simulation of initial value problems for the motion of solid bodies in a Newtonian fluid, part 2. Couette and Poiseuille flows. J Fluid Mech1994; 277: 271–301.
22.
HuHHJosephDDCrochetMJ. Direct simulation of fluid particle motions. Theor Comput Fluid Dyn1992; 3: 285–306.
23.
MortazaviSTryggvasonG. A numerical study of the motion of drops in Poiseuille flow. Part 1. Lateral migration of one drop. J Fluid Mech2000; 411: 325–350.
24.
KhaliliHMortazaviS. Numerical simulation of buoyant drops suspended in Poiseuille flow at nonzero Reynolds numbers. Acta Mech2013; 224: 269–286.
25.
ClagueDSCorneliusPJ. The hydrodynamic force and torque on a bounded sphere in Poiseuille flow. Int J Numer Methods Fluids2001; 35: 55–70.
26.
AdvaniSGTuckerCL III. The use of tensors to describe and predict fiber orientation in short fiber composites. J Rheol1987; 31: 751–784.
27.
BatchelorGK. The stress system in a suspension of force-free particles. J Fluid Mech1970; 41: 545–570.
28.
DinhSMArmstrongRC. Rheological equation of state for semi-concentrated fiber suspensions. J Rheol1984; 28: 207–227.
29.
MackaplowMBShaqfehSG. A numerical study of the rheological properties of suspensions of rigid, non-Brownian fibres. J Fluid Mech1996; 329: 155–186.
30.
RahnamaMShaqfehSG. The effect of hydrodynamic interactions on the orientation distribution in a fiber suspension subject to simple shear flow. Phys Fluids1995; 7: 487–506.
31.
YamaneYKanedaYDioM. Numerical simulation of semi-dilute suspensions of rodlike particles in shear flow. J Non-Newt Fluid Mech1994; 54: 405–421.
32.
FanXJPhan-ThienNZhengR. A direct simulation of fibre suspensions. J Non-Newt Fluid Mech1998; 74: 113–135.
33.
ClaeysILBradyJF. Suspensions of prolate spheroids in Stokes flow. Part 1. Dynamics of a finite number of particles in an unbounded fluid. J Fluid Mech1993; 251: 411–442.
34.
CoxRG. The motion of long slender bodies in a viscous fluid. Part 2. Shear flow. J Fluid Mech1971; 45: 625–657.
35.
ZhangDSmithDEJackDA. Numerical evaluation of single fiber motion for short-fiber-reinforced composite materials processing. ASME Trans J Manuf Sci Eng2011; 133: 051002–1–9.
36.
FolgarFTuckerCL III. Orientation behavior of fibers in concentrated suspensions. J Reinf Plast Compos1984; 3: 98–119.
37.
JunkMIllnerR. A new derivation of Jeffery’s equation. J Math Fluid Mech2007; 9: 455–488.
38.
DoiM. Molecular dynamics and rheological properties of concentrated solutions of rodlike polymers in isotropic and liquid crystalline phases. J Polym Sci Part B Polym Phys1981; 19: 229–243.
39.
CintraJSTuckerCL III. Orthotropic closure approximations for flow-induced fiber orientation. J Rheol1995; 39: 1095–1122.
40.
VerleyeVDupretF. Prediction of fiber orientation in complex injection molded parts. Dev Non-Newt Flows1993; 175: 139–163.
41.
DupretFVerleyeVLanguillierB. Numerical prediction of the molding of composite parts. Rheol Fluid Mech Nonlinear Mater1997; 243: 79–90.
42.
JackDASmithDE. An invariant based fitted closure of the sixth-order orientation tensor for short-fiber suspensions. J Rheol2005; 49: 1091–1115.
43.
ReddyJNGartlingDK. The finite element method in heat transfer and fluid dynamics, 2nd ed. Boca Raton, Florida: CRC Press LLC, 2001.
44.
DoneaJHuertaA. Finite element methods for flow problems, Hoboken, New Jersey: John Wiley & Sons Ltd, 2003.
45.
BaaijensFP. Mixed finite element methods for viscoelastic flow analysis: a review. J Non-Newt Fluid Mech1998; 79: 361–385.
46.
HwangWRHulsenMAMeijerHE. Direct simulations of particle suspensions in a viscoelastic fluid in sliding bi-periodic frames. J Non-Newt Fluid Mech2004; 121: 15–33.
47.
D’AvinoGHulsenMASnijkersF. Rotation of a sphere in a viscoelastic liquid subjected to shear flow. Part i: simulation results. J Rheol (1978–present)2008; 52: 1331–1346.
48.
HwangWKimWYLeeES. Direct simulations of orientational dispersion of platelet particles in a viscoelastic fluid subjected to a partial drag flow for effective coloring of polymers. Korea–Aust Rheol J2011; 23: 173–181.
49.
CherizolRSainMTjongJ. Review of non-newtonian mathematical models for rheological characteristics of viscoelastic composites. Green Sustain Chem2015; 5: 6–6.
