Transient shear flow behavior of concentrated long glass fiber suspensions in a sliding plate rheometer

In order to eventually predict the behavior of long fiber suspensions in complex flows commonly found in processing operations, it is necessary to understand their rheology and its connection to the evolution of fiber orientation and configuration in well defined flows. In this paper we report the transient behavior at the startup of shear flow of a polymer melt containing long glass fibers with a length (L) >1 mm, using a sliding plate rheometer (SPR). The operation of the SPR was confirmed by comparing the transient shear viscosity (η⁺) for a polymer melt and a melt containing short glass fibers (L < 1 mm) with measurements obtained from a cone-and-plate device, using a modified sample geometry that was designed to avoid wall effects. For the long fiber systems, measurements could only be obtained in the SPR because these systems would not stay within the gap of the rotational rheometer. Transient stress growth behavior of the long fiber systems was obtained as a function of shear rate and fiber concentration for samples prepared with three different initial orientations. Results showed that, unlike short fiber systems (with a random planar initial orientation) that usually exhibit a single overshoot peak followed by a steady state, η⁺ of the long fiber suspensions often passed through multiple transient regions, depending on the fiber concentration and applied shear rate. Additionally, η⁺ of the long fiber suspensions was found to be highly dependent on the initial orientation of the sheared samples. Finally, the initial and final fiber orientations of the long glass fiber samples were measured and used to initiate an explanation of the viscosity behavior. The results obtained in this research will be useful for future assessment of a quantitative correlation between transient rheology and the evolution of fiber orientation.     

Nonlinear vibrations of a viscoelastic plate with concentrated masses

The problem of vibrations of a viscoelastic plate with concentrated masses is studied in a geometrically nonlinear formulation. In the equation of motion of the plate, the action of the concentrated masses is taken into account using Dirac δ-functions. The problem is reduced to solving a system of Volterra type ordinary nonlinear integrodifferential equations using the Bubnov-Galerkin method. The resulting system with a singular Koltunov-Rzhanitsyn kernel is solved using a numerical method based on quadrature formulas. The effect of the viscoelastic properties of the plate material and the location and amount of concentrated masses on the vibration amplitude and frequency characteristics is studied. A comparison is made of numerical calculation results obtained using various theories. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 158–169, November–December, 2007.     

Stability of some shear flows for concentrated suspensions

In this paper we investigate the stability of some viscometric flows for a concentrated suspension model which allows for the effects of shear-induced migration, including plane and circular Couette and Poiseulle flows, and unbounded and bounded torsional flows. In the bounded torsional flow, where its radial outer boundary is assumed frictionless, an exact closeform solution is given. With the exception of torsional flows, we find that a limit point for all the steady-state solutions can exist for certain range in the parameter values. In all cases, disturbances can persist for a long time, O (H ²/a ²), where H is a dimension of the flow field, and a is the particles' radius.     

Axisymmetric bending for thick laminated circular plate under a concentrated load

Ｎｏｗａｄａｙｓ,ｔｈｅｃｕｒｒｅｎｔｔｈｅｏｒｉｅｓｏｆｐｌａｔｅｓａｎｄｓｈｅｌｌｓ,ｓｕｃｈａｓＫｉｒｃｈｈｏｆｆ’ｓｔｈｉｎｐｌａｔｅｔｈｅｏｒｙａｎｄＲｅｉｓｓｎｅｒ’ｓｍｏｄｅｒａｔｅｌｙｔｈｉｃｋｐｌａｔｅｔｈｅｏｒｙｅｔｃ．,ａｒｅｅ...     

Study of geometry effects in torsional rheometry of fibre suspensions

This work concerns the problem of measuring the viscosity of a suspension of fibres that are not short compared to the dimensions of the measurement device. We have examined various geometry effects in shear using parallel-plate and cone-and-plate configurations. Steady state viscosity, transient viscosity and first normal stress difference have been studied. Silicone oil with milled glass fibres at 8 vol.% was used as a model fibre suspension. The conventional parallel-plate geometry exhibits a significant gap dependence of the apparent viscous response. Too small gaps constrain the fibre rotations, thus lowering the initial stress peak and delaying the development of steady state. Too large gaps lead to a loss of liquid at the sample perimeter, which reduces the apparent steady state viscosity. The steady state response, however, seems to be correct for gaps in the range 1–2 mm when the maximum fibre length is 1.5 mm. The cone-and-plate geometry is less sensitive to large gaps. Too small a gap at the (truncated) cone apex leads to anomalous response due to bridging of the gap by fibres. The use of unusually large plates (100 mm diameter) and different cone angles and truncations was explored. It was shown that cone-and-plate configurations of large diameter and large truncation can effectively decrease the geometry disturbances both due to fibre bridging at the cone apex and loss of liquid at the edge.     

Dynamic stability of a viscoelastic plate with concentrated masses

The paper addresses the geometrically nonlinear problem of dynamic stability of a viscoelastic plate with concentrated masses. The Bubnov-Galerkin method based on polynomial approximation is used to reduce the problem to a system of nonlinear Volterra-type integro-differential equations with singular relaxation kernels. This system is solved by numerical method based on quadrature formulas. The critical loads are found and their dependence on the arrangement and number of concentrated masses is studied for a wide range of mechanical and geometrical parameters of the plate. The choice of a relaxation kernel for dynamic problems for viscoelastic thin-walled plate-like structures is justified. Results produced by different theories are compared __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 2, pp. 109–118, February 2008.     

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.     

Bending of rectangular plate with each edges arbitrary a point supported under a concentrated load

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｅｎｇｉｎｅｅｒｉｎｇｔｅｃｈｎｏｌｏｇｙ,ｔｈｅｂｅｎｄｉｎｇｐｒｏｂｌｅｍｏｆｅｌａｓｔｉｃｔｈｉｎｐｌａｔｅｈａｓｉｍｐｏｒｔａｎｔｐｒａｃｔｉｃａｌｓｉｇｎｉｆｉｃａｎｃｅ,ａｎｄａｌｏｔｏｆｒｅｓｕｌｔｓｈａｖｅｂｅｅｎｏｂｔａｉｎｅｄ.Ｂｅｃａｕｓｅｔｈｅｐｒｏｂｌｅｍｓｏｆｅａｃｈｅｄｇｅｓｐｏｉｎｔｓｕｐｐｏｒｔｅｄａｒｅｃｏｍｐｌｅｘｅｒ,ｔｈｅｒｅａｒｅｆｅｗｒｅｓｕｌｔｓｏｆｔｈｅｐｒｏｂｌｅｍｓ.Ｉｎｔｈｅｐａｐｅｒ[1],ｔｈｅｍｅｔｈｏｄｏｆｓｕｐ…     

Bending of cantilever rectangular plate with concentrated load

In this paper the solutions for the bending of cantilever rectangular plates with concentrated load acting at any point of the middle line perpendicular to the clamped edge are given by means of a conception named modified simply supported edges and the method of superposition. Some numerical examples are presented. The total bending moment checks very well with the value determined statically.     

NMR flow imaging of Newtonian liquids and a concentrated suspension through an axisymmetric sudden contraction

Using nuclear magnetic resonance (NMR) flow imaging to examine fluid motions at constant velocities or flows that change relatively slowly has been well-documented in the literature. Application of this technique to accelerative flows, on the other hand, has been limited. This study reports the use of an NMR flow imaging method, for which acceleration is not explicitly compensated in the NMR pulse sequence, to measure axial and radial fluid motions during flow through an axisymmetric sudden contraction. In this flow geometry, both velocity and acceleration are spatially dependent. The flow contraction ratio was 2:1. The method was first applied to examine Newtonian liquids at low and high Reynolds numbers under laminar flow conditions. The measured axial and radial velocity profiles, without accounting for acceleration effects in the data analysis, across the contraction are in excellent qualitative agreement with previous experimental data and theoretical calculations reported in the literature. Quantitative comparison of the axial and radial velocities with numerical results indicates that the maximum error from acceleration effects is about 10%. The method has also been used to examine the flow of a concentrated suspension (50% by volume of solid particles) through the contraction. The flow kinematics of the suspension at creeping flow conditions appear to mimic those of the Newtonian fluid with some slight differences. NMR images taken immediately following the cessation of flow suggest a slight degree of particle migration toward the center of the pipe downstream of the contraction.     

Nonlinear rheology of a concentrated spherical silica suspension:

Time-dependent nonlinear flow behavior was investigated for a model hard-sphere suspension, a 50 wt% suspension of spherical silica particles (radius = 40 nm; effective volume fraction = 0.53) in a 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The suspension had two stress components, the Brownian stress ^B and the hydrodynamic stress ^H After start-up of flow at various shear rates , the viscosity growth function ₊ (t, ) was measured with time t until it reached the steady state. The viscosity decay function _– (t, ) was measured after cessation of flow from the steady as well as transient states. At low where the steady state viscosity ( ) exhibited the shear-thinning, the _– (t, ) and ⁺ (t, ) data were quantitatively described with a BKZ constitutive equation utilizing data for nonlinear relaxation moduli G (t, ). This result enabled us to attribute the thinning behavior to the decrease of the Brownian contribution ^B = ^B / (considered in the BKZ equation through damping of G (t, )). On the other hand, at high where ( ) exhibited the thickening, the BKZ prediction largely deviated from the ₊ (t, ) and ₊ (t, ) data, the latter obtained after cessation of steady flow. This result suggested that the thickening was due to an enhancement of the hydrodynamic contribution ^H = ^H / (not considered in the BKZ equation). However, when the flow was stopped at the transient state and only a small strain (<0.2) was applied, ^H was hardly enhanced and the _– (t, ) data agreed with the BKZ prediction. Correspondingly, the onset of thickening of ₊ (t, ) was characterized with a -insensitive strain ( 0.2). On the basis of these results, the enhancement of ^H (thickening mechanism) was related to dynamic clustering of the particles that took place only when the strain applied through the fast flow was larger than a characteristic strain necessary for close approach/collision of the particles.     

Stability analysis in spatial mode for channel flow of fiber suspensions

Different from previous temporal evolution assumption, the spatially growing mode was employed to analyze the linear stability for the channel flow of fiber suspensions. The stability equation applicable to fiber suspensions was established and solutions for a wide range of Reynolds number and angular frequency were given numerically . The results show that, the flow instability is governed by a parameter H which represents a ratio between the axial stretching resistance of fiber and the inertial force of the fluid. An increase of H leads to a raise of the critical Reynolds number, a decrease of corresponding wave number, a slowdown of the decreasing of phase velocity , a growth of the spatial attenuation rate and a diminishment of the peak value of disturbance velocity. Although the unstable region is reduced on the whole, long wave disturbances are susceptible to fibers.     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

Use of a sliding plate rheometer to measure the first normal stress difference at high shear rates

The use of a sliding plate rheometer (SPR) to determine the first normal stress difference of molten polymers and elastomers at high shear rates is demonstrated. The simple shear flow in this instrument is not subject to the flow instabilities that limit the use of rotational rheometers to shear rates often below 1 s⁻¹. However, issues of secondary flow and wall slip must be addressed to obtain reliable data using an SPR. A highly entangled, monodisperse polybutadiene and a commercial polystyrene were the polymers studied. The inclusion of the polystyrene made it possible to compare data with those obtained by Lodge using a stressmeter, which is an instrument based on the measurement of the hole pressure. The data from the two instruments are in good agreement and are also close to the predictions of an empirical equation of Laun based on the storage and loss moduli.     

The large deflection on problem of circular plate on elastic foundation and in conjunction with linear elastic structure

This paper deals with the axisymmetrical deformation of the circular plate in large deflection, which is on elastic foundation and in conjunction with a certain linear elastic structure. The governing integral equations are established by the method of mixed boundary condition 1 and the simplified form is given. The pertrubation method is used to obtain the solutions and an example of the composite structure made up of a circular plate and a cylindrical shell is presented.     

Theoretical research on rotational dispersion coefficient of fiber in turbulent shear flow of fiber suspension

The rotational dispersion coefficient of the fiber in the turbulent shear flow of fiber suspension was studied theoretically. The function of correlation moment between the different fluctuating velocity gradients of the flow was built firstly. Then the expression, dependent on the characteristic length, time, velocity and a dimensionless parameter related to the effect of wall, of rotational dispersion coefficient is derived. The derived expression of rotational dispersion coefficient can be employed to the inhomogeneous and non-isotropic turbulent flows. Furthermore it can be expanded to three-dimensional turbulent flows and serves the theoretical basis for solving the turbulent flow of fiber suspension.     

research on rotational dispersion coefficient of fiber in turbulent shear flow of fiber suspension

The rotational dispersion coefficient of the fiber in the turbulent shear flow of fiber suspension was studied theoretically. The function of correlation moment between the different fluctuating velocity gradients of the flow was built firstly. Then the expres- sion, dependent on the characteristic length, time, velocity and a dimensionless parameter related to the effect of wall, of rotational dispersion coefficient is derived. The derived expression of rotational dispersion coefficient can be employed to the inhomogeneous and non-isotropic turbulent flows. Furthermore it can be expanded to three-dimensional turbulent flows and serves the theoretical basis for solving the turbulent flow of fiber suspension.