Investigation of fiber motion near solid boundaries in simple shear flow

In this paper, fiber motion near a planar wall was investigated using a planar shear flow apparatus. Fibers were placed (one at a time) perpendicular to the flow direction at various locations throughout the flow field. The location and orientation of each fiber versus time was measured, using an image processing system, until the fiber aligned with the flow direction. When the centroid of the fiber was located at distances greater than a fiber length from the wall, Jeffery's equations governing particle motion were verified. For distances less than a fiber length and greater than a fiber diameter from the wall, the fiber experienced an increased rate of rotation. In this regime, the motion of the fiber could be described by Jeffery's equations if an increased effective shear rate was used. The effective shear rate was found to increase logarithmically with decreasing separation distance. The wall effect was higher for longer aspect ratio fibers and was also a function of orientation; fibers oriented perpendicular to the wall rotated faster than those oriented parallel to the wall at the same separation distance. Once the fiber aligned with the flow direction, it ceased to rotate within the field of view. In this orientation, the wall had a stabilizing effect on the fiber. In efforts to relate the increase in shear rate to the aspect ratio of the fiber and the separation distance between the fiber and a solid wall, a translation model based on the work of De Mestre and Russel was explored. This model allows one to quantify the increase in shear rate experienced by the fiber due to the presence of a wall or obstruction in the flow field. However, the model has its limitations and care should be taken when applying this model outside its realm of validity. When compared to experimental data, the translation model provides a very good estimate of the increased shear rate experienced by the fiber when it is located less than 2/3 of a fiber length from a planar wall. Received: 20 April 2000 Accepted: 28 September 2000     

Numerical simulation of flow kinematics and fiber orientation for multi-disperse suspension

The development of flow kinematics and fiber orientation distribution from the parabolic velocity profile and isotropic orientation at the channel inlet was computed in multi-disperse suspension flow through a parallel plate channel and their predictions were compared with those of mono- and bi-disperse suspensions. A statistical scheme (orientations of a large number of fibers are evaluated from the solution of the Jeffery equation along the streamlines) was confirmed to be very useful and feasible method to analyze accurately the orientation distribution of fibers in multi-disperse fiber suspension flow as well as mono- and bi-dispersions, instead of direct solutions of the orientation distribution function of fibers or the evolution equation of the orientation tensor which involves a closure equation. It was found that the flow kinematics and the fiber orientation depend completely on both the fiber aspect-ratio and the fiber parameter for multi-disperse suspension when the fiber–fiber and fiber-wall interactions are neglected. Furthermore, the addition of large aspect-ratio fibers as well as an increase in the fiber parameter related to the large aspect-ratio fibers could suppress the complex velocity field and stress distributions which are observed in suspensions containing small aspect-ratio fibers. From a practical point of view, therefore, the mechanical and physical properties of fiber composites should be improved with an increase in the volume fraction of large aspect-ratio fibers.     

An experimental study of flow-induced fiber orientation and concentration distributions in a concentrated suspension flow through a slit channel containing a cylinder

Flow-induced fiber orientation and concentration distributions were measured in a concentrated fiber suspension (CFS) and a dilute one (DFS). The channel has a thin slit geometry containing a circular cylinder. In the previous work, many researchers have qualitatively studied fiber orientation and concentration distributions in injection-molded products of fiber-reinforced plastics. In the present work, however, they are quantitatively estimated by direct observation of fibers in the concentrated suspension flow. For the CFS, some fibers rotate in an expansion part between the channel wall and the circular cylinder, and the fiber orientation becomes almost random state. On the other hand, fibers are perfectly aligned along the flow direction owing to the elongational flow near the centerline downstream of the cylinder. The fiber concentration has a flat distribution except near the channel wall and the centerline. For the DFS a minimum in the fiber concentration distribution was clearly observed on the centerline, and two peaks beside the centerline and near the channel wall. This characteristic distribution is caused by the fiber-wall and fiber-cylinder interactions. It is found that the obstacle such as the circular cylinder in the channel significantly affects the fiber orientation downstream of the obstacle for the CFD, while it affects the fiber concentration distribution for the DFS.     

Well-posedness of the problem of fiber suspension flows

The well-posedness of the equations governing the flow of fiber suspensions is studied. The fluid is assumed to be Newtonian and incompressible, and the presence of fibers is accounted for through the use of second- and fourth-order orientation tensors, which model the effects of the orientation of fibers in an averaged sense. The fourth-order orientation tensor is expressed in terms of the second-order tensor through various closure relations. It is shown that the linear closure relation leads to anomalous behavior, in that the rest state of the fluid is unstable, in the sense of Liapounov, for certain ranges of the fiber particle number. No such anomalies arise in the case of quadratic and hybrid closure relations. For the quadratic closure relation, it is shown that a unique solution exists locally in time for small data.     

Orientation distribution of fibers and rheological property in fiber suspensions flowing in a turbulent boundary layer

A model relating the translational and rotational transport of orientation distribution function （ODF） of fibers to the gradient of mean ODF and the dispersion coefficients is proposed to derive the mean equation for the ODE Then the ODF of fibers is predicted by numerically solving the mean equation for the ODF together with the equations of turbulent boundary layer flow. Finally the shear stress and first normal stress difference of fiber suspensions are obtained. The results, some of which agree with the available relevant experimental data, show that the most fibers tend to orient to the flow direction. The fiber aspect ratio and Reynolds number have significant and negligible effects on the orientation dis- tribution of fibers, respectively. The additional normal stress due to the presence of fibers is anisotropic. The shear stress of fiber suspension is larger than that of Newtonian solvent, and the first normal stress difference is much less than the shear stress. Both the additional shear stress and the first normal stress difference increase with increasing the fiber concentration and decreasing fiber aspect ratio.     

Research on the motion of particles in the turbulent pipe flow of fiber suspensions

The motion of fibers in turbulent pipe flow was simulated by 3-D integral method based on the slender body theory and simplified model of turbulence. The orientation distribution of fibers in the computational area for different Re numbers was computed. The results which were consistent with the experimental ones show that the fluctuation velocity of turbulence cause fibers to orient randomly. The orientation distributions become broader as the Re number increases. Then the fluctuation velocity and angular velocity of fibers were obtained. Both are affected by the fluctuation velocity of turbulence. The fluctuation velocity intensity of fiber is stronger at longitudinal than at lateral, while it was opposite for the fluctuation angular velocity intensity of fibers. Finally, the spatial distribution of fiber was given. It is obvious that the fiber dispersion is strenghened with the increase of Re numbers.     

Fiber suspension flow in a tapered channel: The effect of flow/fiber coupling

A numerical model for predicting the flow and orientation state of semi-dilute, rigid fiber suspensions in a tapered channel is presented. The effect of the two-way flow/fiber coupling is investigated for low Reynolds number flow using the constitutive model of Shaqfeh and Fredrickson. An orientation distribution function is used to describe the local orientation state of the suspension and evolves according to a Fokker–Plank type equation. The planar orientation distribution function is determined along streamlines of the flow and is coupled with the fluid momentum equations through a fourth-order orientation tensor. The coupling term accounts for the two-way interaction and momentum exchange between the fluid and fiber phases. The fibers are free to interact through long range hydrodynamic fiber–fiber interactions which are modeled using a rotary diffusion coefficient, an approach outlined by Folgar and Tucker. Numerical predictions are made for two different orientation states at the inlet to the contraction, namely a fully random and a partially aligned fiber orientation state. Results from these numerical predictions show that the streamlines of the flow are altered and that velocity profiles change from Jeffery–Hamel, to something resembling a plug flow when the fiber phase is considered in the fluid momentum equations. This phenomenon was found when the suspension enters the channel in either a pre-aligned, or in a fully random orientation state. When the suspension enters the channel in an aligned orientation state, fiber orientation is shown to be only marginally changed when the two-way coupling is included. However, significant differences between coupled and uncoupled predictions of fiber orientation were found when the suspension enters the channel in a random orientation state. In this case, the suspension was shown to align much more quickly when the mutual coupling was accounted for and profiles of the orientation anisotropy were considerably different both qualitatively and quantitatively.     

Stability of axially compressed cylindrical shells made of reinforced materials with specific fiber orientation within each layer

A technique for stability analysis of anisotropic cylindrical shells is developed. It permits us to examine the cases of reinforcement where the elastic axes of layers do not coincide with the coordinate axes of the shell. The solution is obtained using the mixed equations of the Donnell-Mushtari-Vlasov theory of shells. The deflection and force functions are approximated by trigonometric series. Single-layer and multilayer cylindrical shells with fiber orientation of two types are analyzed for stability. It is revealed that when layers are few, failure to incorporate the direction of fibers in layers into the design model results in highly inaccurate values of critical loads __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 3, pp. 80–88, March 2006.     

Structure and properties of sheared fiber suspensions with mechanical contacts

The properties of fiber suspensions are highly sensitive to the suspension microstructure. In dilute or semi-dilute suspensions, nL²d≪1, the fibers' orientation distribution is controlled by hydrodynamic interactions among the fibers. However, direct mechanical contacts among the fibers play an important role in semi-concentrated suspensions, nL²d=O(1). Here, n is the number of fibers per unit volume, L is the fiber length and d is the fiber diameter. We have performed dynamic simulations of fiber suspensions including contact forces that prevent any two fibers from passing through one another. Collisions between the fibers cause them to flip more frequently in the shear flow, leading to a spread of the orientation distribution away from the flow direction. Both this increased orientational dispersion and the direct stress transmitted through the contacts enhance the shear viscosity of the suspension significantly. The contacts also give rise to normal stress differences. The results of the simulation are compared with experiments and the relative importance of contacts and hydrodynamic interactions is discussed.     

Comparative numerical study of two concentrated fiber suspension models

We consider two rheological models for concentrated fiber suspensions. In both models the equations for orientation and flow are fully coupled, i.e., the orientation influences the flow via a constitutive relation for the viscosity and the orientation of the fibers is determined by the flow field. The orientation state of the fibers is characterized by the Advani–Tucker orientation tensor. We are investigating suspensions of fibers in which the kinetic energies of the fibers are large compared to the thermal energies, i.e., the influence of Brownian motion may be neglected. The first model is the Folgar–Tucker model with backcoupling to the flow (FT model). The second model is an extension of Folgar–Tucker, which models phenomenologically the topological exclusion interaction in dense suspensions (FTMS model). As test cases for the simulation are considered channel flow, 8:1 contraction flow and flow around a cylinder.     

Numerical simulation of flow kinematics for suspensions containing continuously dispersed aspect ratio fibers

Fiber suspension flow and fiber orientation through a parallel-plate channel were numerically simulated for fiber suspensions including continuously dispersed aspect ratios from 10 to 50. In the simulations, both the fiber–fiber and fiber–wall interactions were not taken into account. A statistical scheme that proceeds by evaluating the orientation evolution of a large number of fibers from the solution of the Jeffery equation along the streamlines was confirmed to be a very useful and feasible method to accurately analyze the orientation distribution of fibers with continuously dispersed aspect ratios. For monodisperse suspensions with small-aspect-ratio fibers, flip-over or oscillation phenomenon of the orientation ellipsoid caused the wavy patterns of the velocity profile and the streamlines as well as the abrupt and complex variation of the shear stress and the normal stress difference near the channel wall as proven in one of our former works. On the other hand, continuous dispersions containing from small- to large-aspect-ratio fibers were able to induce smoother evolutions of the fiber orientation and the flow kinematics. In the processing of fiber composites, the length of suspended fibers is always continuously distributed because of fiber breakage during processing; thus, the smooth evolutions of the flow kinematics and the stress distribution can be attained.This paper was presented at the Annual Meeting of the European Society of Rheology, Grenoble, April 2005.     

Torsional Swelling of a Hyperelastic Tube with Helically Wound Reinforcement

This paper examines the combination of radial deformation with torsion for a circular cylindrical tube composed of a transversely isotropic hyperelastic material subject to finite deformation swelling. The stored energy function involves separate matrix and fiber contributions such that the fiber contribution is minimized when the fiber direction is at a natural length. This natural length is not affected by the swelling. Hence swelling preferentially expands directions that are orthogonal to the fiber. The swelling itself is described via a swelling field that prescribes the local free volume at each location in the body. Such a treatment is a relatively simple generalization of the conventional incompressible theory. The direction of transverse isotropy associated with the fiber reinforcement is described by a helical winding about the tube axis. The swelling induced preferential expansion orthogonal to this direction induces the torsional aspect of the deformation. For a specific class of strain energy functions we find that the twist increases with swelling and approaches a limiting asymptotic value as the swelling becomes large. The fibers reorient such that fibers at the inner portion of the tube assume a more circumferential orientation whereas, at least for small and moderate swelling, the fibers in the outer portion of the tube assume a more axial orientation. For large swelling the fibers in the outer portion of the tube reorient beyond the axial orientation, and so are described by helices with orientation in the opposite sense to that in the reference configuration.      

Flow-driven orientation dynamics of semiflexible fiber systems

We study the flow-induced orientation dynamics of semiflexible fibers in dilute fiber suspensions. Starting from the equations of motion for a two-rod model of flexible fibers in Stokes flow, the Smoluchowski equation for a connected monomer orientation distribution function is derived. We then obtain a set of equations for the time dependence of the first and second moments of the orientation distribution function, thus extending the Folgar Tucker equations for short rigid fiber suspensions to flexible fiber suspensions. The resulting generalized equations for the orientation dynamics of a suspension of flexible fibers are solved for simple channel flow. It is shown that all qualitative effects of bending and straightening of fibers and their influence on the orientation of flexible fibers are captured within our model. A scalar measure for the distribution of bending in a flow is introduced, which allows to detect the degree of bending of fibers. Paper was presented at the 3rd Annual Rheology Conference, AERC 2006, April 27–29, 2006, Crete, Greece.     

Investigation of the abrupt contraction flow of fiber suspensions in polymeric fluids

Numerical simulations of the flow of rigid fibres through a 4:1 planar contraction, and the predicted flow pattern and fiber orientation are presented. Entirely new is the examination of the nature of the suspending matrix which may consist of either a Newtonian fluid or a polymer melt. In the case of a polymer matrix three rheological models, the Phan-Thien–Tanner, FENE-CR, and Carreau models have been used to investigate the effects of shear-thinning and elasticity on the flow and the orientation of the fibers. The effects of inertia are neglected, and the governing equations for the flow field, polymer stress, and fiber orientation are coupled and simultaneously solved. A parametric study is used to explore the effects of different dimensionless parameters on the velocity field, the fiber orientation, the pressure drop, as well as the vortex size measured by the dimensionless reattachment length. We particularly focus on the role of the fibers aspect ratio, volume fraction, and interaction coefficient which measures the intensity of fiber interaction in the suspension. Furthermore, we evaluate and compare the results of four different closure approximations: the quadratic, linear, hybrid A and T, and natural closures.     

A new constitutive theory for fiber-reinforced incompressible nonlinearly elastic solids

We consider an incompressible nonlinearly elastic material in which a matrix is reinforced by strong fibers, for example fibers of nylon or carbon aligned in one family of curves in a rubber matrix. Rather than adopting the constraint of fiber inextensibility as has been previously assumed in the literature, here we develop a theory of fiber-reinforced materials based on the less restrictive idea of limiting fiber extensibility. The motivation for such an approach is provided by recent research on limiting chain extensibility models for rubber. Thus the basic idea of the present paper is simple: we adapt the limiting chain extensibility concept to limiting fiber extensibility so that the usual inextensibility constraint traditionally used is replaced by a unilateral constraint. We use a strain-energy density composed with two terms, the first being associated with the isotropic matrix or base material and the second reflecting the transversely isotropic character of the material due to the uniaxial reinforcement introduced by the fibers. We consider a base neo-Hookean model plus a special term that takes into account the limiting extensibility in the fiber direction. Thus our model introduces an additional parameter, namely that associated with limiting extensibility in the fiber direction, over previously investigated models. The aim of this paper is to investigate the mathematical and mechanical feasibility of this new model and to examine the role played by the extensibility parameter. We examine the response of the proposed models in some basic homogeneous deformations and compare this response to those of standard models for fiber reinforced rubber materials. The role of the strain-stiffening of the fibers in the new models is examined. The enhanced stability of the new models is then illustrated by investigation of cavitation instabilities. One of the motivations for the work is to apply the model to the biomechanics of soft tissues and the potential merits of the proposed models for this purpose are briefly discussed.     

An anisotropic rotary diffusion model for fiber orientation in short- and long-fiber thermoplastics

The Folgar–Tucker model, which is widely-used to predict fiber orientation in injection-molded composites, accounts for fiber–fiber interactions using isotropic rotary diffusion. However, this model does not match all aspects of experimental fiber orientation data, especially for composites with long discontinuous fibers. This paper develops a fiber orientation model that incorporates anisotropic rotary diffusion. From kinetic theory we derive the evolution equation for the second-order orientation tensor, correcting some errors in earlier treatments. The diffusivity is assumed to depend on a second-order space tensor, which is taken to be a function of the orientation state and the rate of deformation. Model parameters are selected by matching the experimental steady-state orientation in simple shear flow, and by requiring stable steady states and physically realizable solutions. Also, concentrated fiber suspensions align more slowly with respect to strain than models based on Jeffery's equation, and we incorporate this behavior in an objective way. The final model is suitable for use in mold filling and other flow simulations, and it gives improved predictions of fiber orientation for injection molded long-fiber composites.     

On the fiber orientation in steady recirculating flows involving short fibers suspensions

In this work we present a new numerical strategy to treat the 3D Fokker–Planck equation in steady recirculating flows. This strategy combines some ideas of the method of particles, with a more original treatment of the periodicity condition, which characterizes the steady solution of the FP equation in steady recirculating flows, as usually encountered in some rheometric devices. Using this numerical technique the fiber orientation distribution can be computed accurately in any steady recirculating flow. The simulation results can be used to identify some rheological parameters of the suspension, using an inverse technique, as well as to analyze the validity of some simplified models widely used, which require a closure relation. Thus, in this paper several closure relations of the fourth-order orientation tensor will be discussed in the context of a numerical example involving a steady recirculating flow.     

Experimental-computational study of fibrous particle transport and deposition in a bifurcating lung model

Experiments carried out using a lung model with a single horizontal bifurcation under different steady inhalation conditions explored the orientation of depositing carbon fibers, and particle deposition fractions. The orientations of deposited fibers were obtained from micrographs. Specifically, the effects of the sedimentation parameter (γ), fiber length, and flow rate on orientations were analyzed. Our results indicate that gravitational effect on deposition cannot be neglected for 0.0228 < γ < 0.247. The absolute orientation angle of depositing fibers decreased linearly with increasing γ for values 0.0228 < γ < 0.15. Correspondence between Stokes numbers and γ suggests these characteristics can be used to estimate fiber deposition in the lower airways. Computer simulations with sphere-equivalent diameter models for the fibers explored deposition efficiency vs. Stokes number. Using the volume-equivalent diameter model, our experimental data for the horizontal bifurcation were replicated. Results for particle deposition using a lung model with a vertical bifurcation indicate that body position also affects deposition.     

Investigation of closure approximations for fiber orientation distribution in contracting turbulent flow

We have investigated the orientation state of a dilute fiber suspension flow in a planar contraction at high Reynolds numbers in turbulent flow. High speed imaging is used to directly measure the orientation distribution function at different downstream positions along the contraction centerline. The results from the direct measurement of the orientation distribution are used to evaluate the existing closure models. The results show that the fitted orthotropic and natural closure approximations give almost identical results with the best agreement to the orientation distribution in the contraction flow considered here.     

Numerical prediction of fiber motion in a branching channel flow of fiber suspensions

Fiber orientation and dispersion in the dilute fiber suspension that flows through a T-shaped branching channel are simulated numerically based on the slender-body theory. The simulated results are consistent qualitatively with the experimental data available in the literature. The results show that the spatial distribution of fibers is dependent on the fiber aspect ratio, but has no relation with the volume fraction of fiber. The content ratio of fibers near the upper wall increases monotonically with an increasing Re number, and the situation is reverse for the region near the bottom wall. The orientation of fibers depends on Re number, however, the function of fiber volume fraction and aspect ratio is negligible. The fibers near the wall and in the central region of the channel align along the flow direction at all times, but the fibers in the other parts of the channel tend to align along the flow direction only in the downstream region.The project supported by the National Natural Science Foundation of China (10372090) and Doctoral Program of Higher Education in China (20030335001)The English text was polished by Ron Marshall