Fiber orientation kinetics of a concentrated short glass fiber suspension in startup of simple shear flow

The common approach for simulating the evolution of fiber orientation during flow in concentrated suspensions is to use an empirically modified form of Jeffery's equation referred to as the Folgar–Tucker (F-T) model. Direct measurements of fiber orientation were performed in the startup of shear flow for a 30 wt% short glass fiber-filled polybutylene terephthalate (PBT-30); a matrix that behaves similar to a Newtonian fluid. Comparison between predictions based on the F-T model and the experimental fiber orientation show that the model over predicts the rate of fiber reorientation. Rheological measurements of the stress growth functions show that the stress overshoot phenomenon approaches a steady state at a similar strain as the fiber microstructure, at roughly 50 units. However, fiber orientation measurements suggest that a steady state is not reached as the fiber orientation continues to slowly evolve, even up to 200 strain units. The addition of a “slip” parameter to the F-T model improved the model predictions of the fiber orientation and rheological stress growth functions.     

Finite-length effects on dynamical behavior of rod-like particles in wall-bounded turbulent flow

Combined Particle Image Velocimetry (PIV) and Particle Tracking Velocimetry (PTV) measurements have been performed in dilute suspensions of rod-like particles in wall turbulence. PIV results for the turbulence field in the water table flow apparatus compared favorably with data from Direct Numerical Simulations (DNS) of channel flow turbulence and the universality of near-wall turbulence justified comparisons with DNS of fiber-laden channel flow. In order to examine any shape effects on the dynamical behavior of elongated particles in wall-bounded turbulent flow, fibers with three different lengths but the same diameter were used. In the logarithmic part of the wall-layer, the translational fiber velocity was practically unaffected by the fiber length l. In the buffer layer, however, the fiber dynamics turned out to be severely constrained by the distance z to the wall. The short fibers accumulated preferentially in low-speed areas and adhered to the local fluid speed. The longer fibers (l/z > 1) exhibited a bi-modal probability distribution for the fiber velocity, which reflected an almost equal likelihood for a long fiber to reside in an ejection or in a sweep. It was also observed that in the buffer region, high-speed long fibers were almost randomly oriented whereas for all size cases the slowly moving fibers preferentially oriented in the streamwise direction. These phenomena have not been observed in DNS studies of fiber suspension flows and suggested l/z to be an essential parameter in a new generation of wall-collision models to be used in numerical studies.     

Optimization of the rod chain model to simulate the motions of a long flexible fiber in simple shear flows

A rod-chain model was used to study the dynamics of a long flexible fiber which is suspended in a Newtonian fluid. It was shown that the fiber exhibited different apparent flexibility, which depended on the stiffness of the fiber, the strength of the flow field and the initial configuration of the fiber. Different combinations of rod length and rod number were used to simulate the motions of a long fiber in simple shear flow. Long rods would greatly save the computational resources with the loss of accuracy, while short rods would generate accurate results with the loss of efficiency. An optimum rod length was suggested to balance the requirement of accuracy and efficiency, and it was the function of the relative strength of flow field and the initial angle between fiber and vortex axis.     

Highly spatially resolved velocity measurements of a turbulent channel flow by a fiber-optic heterodyne laser-Doppler velocity-profile sensor

Velocity measurements with a high spatial resolution are important in turbulent flow research. In this paper, we report on the development of a new fiber-optic laser-Doppler velocity-profile sensor exhibiting a spatial resolution of up to 5 μm and its application to turbulent boundary layers. The sensor developed in the present work employs a frequency-division-multiplexing technique in order to separate two measurement signals from the two fringe systems. Velocity measurements close to zero at the solid wall were realized using heterodyne technique. The use of fiber optics improved a robustness of the sensor. The measurement accuracy of the sensor was experimentally investigated with respect to the spatial resolution and velocity. Universal velocity profile of a turbulent flow was obtained in a fully developed channel flow. Mean and fluctuating velocity are presented with a high spatial resolution.     

Green water void fraction due to breaking wave impinging and overtopping

The present study uses laboratory measurements to investigate the void fraction of an overtopping flow on a structure. The overtopping flow, also called green water, was generated by the impingement of a plunging breaking wave on the structure following the Froude similarity of an extreme hurricane wave and a simplified offshore structure. The flow is multi-phased and turbulent with significant aeration. A fiber optic reflectometer (FOR) and bubble image velocimetry (BIV) were employed to measure the void fraction and velocity in the flow, respectively, and to determine the water level on the deck. Mean properties of void fraction and velocity were obtained by ensemble-averaging and time-averaging the repeated instantaneous measurements. The temporal and spatial distributions of void fraction reveal that the flow is very highly aerated near the front of green water and has relatively low aeration near the deck surface. The mean void fraction and velocity distributions were also depth-averaged for simplicity and potential use in engineering applications. Using the measured data, similarity profiles for depth-averaged void fraction, depth-averaged velocity, and water level were found. The study suggests that using only the velocity data is insufficient if the flow momentum or the flow rate is to be determined. The accuracy of the void fraction measurements was validated by comparing the directly measured water volume of the overtopping flow with the calculated water volume based on the measured velocity and void fraction.     

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.     

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     

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.     

Micro-interferometry for measurement of thermal displacements at fiber/matrix interfaces

A micro-interferometric technique for measuring out-of-plane thermal displacements on a scale commensurate with the dimensions of the fiber/matrix unit cell is described. A scanning micro-interferometer is used to image surface displacements of samples containing a single-pitch-based carbon fiber embedded in an epoxy matrix. The interferometer design gives the necessary resolution to detect small changes in thermal displacements in the fiber/matrix interface region. The samples were heated electrically through the fiber to create radially symmetric temperature and displacement fields. Repeatable displacement measurements were obtained on a radial line across the interface region with an accuracy of ±25 Å. A sharp expansion of the matrix surrounding the fiber was observed with each heating. Overall, the experiments demonstrate the utility of micro-interferometry for measuring submicron displacements.     

Simultaneous PIV and concentration measurements in a gas-turbine combustor model

Simultaneous velocity and concentration measurements have been performed in a gas-turbine combustor model. Particle image velocimetry (PIV) was used to acquire planar velocity information and to identify coherent flow structures. The Mie scattering technique, based on a slightly modified experimental setup, was used for concentration measurements in this mixing flow. The degree of mixing was assessed by examining local concentration measurements while inhomogeneously seeding the primary and secondary stream of the mixing layer. Connections between flow field and concentration distribution were highlighted using the proper orthogonal decomposition algorithm (POD). Uncertainties and systematic errors for the PIV measurements due to the suboptimal seeding are discussed using a comparison with a second test series at optimal seeding conditions. Results are presented for several flow parameters and at various lateral planes.     

Gravity Effects on Fiber Dynamics in Wall Turbulence

The present study investigated gravity effects on the dynamical behavior of inertial fibers suspended in a vertical channel flow. Direct numerical simulations were performed to obtain the turbulent flow field and the fibers were modelled as prolate spheroidal point particles. For each of the four fiber classes, three different gravity configurations were considered: upward flow with gravity opposing, downward flow with aiding gravity, and channel flow in absence of gravity. Results for the fiber distribution and the translational and rotational fiber motion were reported. In the near-wall region, the presence of gravity resulted in an increased fiber density in the downward flow but a nearly uniform distribution of fibers in upward flow. However, the preferential clustering of fibers in near-wall low-speed streaks was unaffected by gravity. The mean wall-normal or drift velocity of the fibers was higher in the downward flow and lower in the upward flow as compared to the case with no gravity. The suppressed drift velocity in the upward flow resulted in a more uniform fiber distribution throughout the channel in contrast to the near-wall accumulation of fibers in the two other cases. Overall gravity turned out to have negligible effects on some of the statistics of the least inertial fibers whereas the inclusion of gravity had a strong impact for heavier fibers.     

Champ dynamique dans une cellule d'échangeur à vortex : 1. Régime d'entrée laminaire

Numerical simulation and LDV measurements are performed for laminar inlet flow condition in a vortex exchange chamber with isothermal conditions. Newtonian and non-Newtonian fluids were used. Experimental measurements are compared to results of numerical computations with good agreement. Both show the existence of a secondary flow generated by hydrodynamical instabilities due to streamline curvature. The main vortex flow appears then to be restricted by this secondary flow.     

Experimental study of vortex emission behind bluff obstacles in a gas liquid vertical two-phase flow

Vortex emission behind cylinders with trapezoidal cross section was experimentally studied in air-water vertical two-phase flows (liquid velocities vary from 45 cm/s to 2 m/s inside a 15 cm ID pipe); the void fraction ranged from 0 to 25%. The measurements were performed at room pressure and temperature. Two flow regimes were observed. For void fraction smaller than 10% vortex emission remained stable and its frequency sharply defined. However, the rms amplitude of the associated pressure fluctuations strongly decreased. These results were explained by bubble trapping inside the vortex cores. This effect was verified experimentally and analyzed using optical fiber probe measurements. Above a 10% void fraction, vortex emission became erratic. Its spectrum became broader but could be identified up to 25% void fraction.     

The flow behavior of fiber suspensions in Newtonian fluids and polymer solutions.

The main purpose of this study was to examine the viscous and elastic properties and capillary flow of fiber suspensions in Newtonian fluids as well as in polymer solutions. The fillers used were glass, carbon, nylon and vinylon fibers. Glycerin was used as a Newtonian suspending medium and HEC and Separan solutions as viscoelastic suspending media. The viscosity and the first normal-stress difference were measured using a coaxial cylindrical rotating viscometer and a parallel-plate rheogoniometer respectively. The influence of the concentration, aspect ratio, diameter and flexibility of the fibers on the viscous and elastic properties of the fiber suspensions was investigated. Empirical equations were obtained for the relative viscosity and first normal-stress difference for the fiber suspensions in glycerin. The capillary flow of these suspensions is discussed in part II.     

A time dependent micro-mechanical fiber composite model for inelastic zone growth in viscoelastic matrices

In this work, a fiber composite model is developed to predict the time dependent stress transfer behavior due to fiber fractures, as driven by the viscoelastic behavior of the polymer matrix, and the initiation and propagation of inelastic zones. We validate this model using in situ, room temperature, micro-Raman spectroscopy fiber strain measurements. Multifiber composites were placed under constant load creep tests and the fiber strains were evaluated with time after one fiber break occurred. These composite specimens ranged in fiber volume fraction and strain level. Comparison between prediction and MRS measurements allows us to characterize key in situ material parameters, the critical matrix shear strain for inelastic zones and interfacial frictional slip shear stress. We find that the inelastic zone is predominately either shear yielding or interfacial slipping, and the type depends on the local fiber spacing.     

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.     

EFFECTS OF TENSOR CLOSURE MODELS AND 3-D ORIENTATION ON THE STABILITY OF FIBER SUSPENSIONS IN A CHANNEL FLOW

IntroductionFlowoffibresuspensionshasbeenveryfamiliarinmanyindustrialfields.Fibreadditivesplayanimportantroleindragreductioninmanytypesofflow[1- 3].Inthesuspensions,somebehavioroftheflowmaybealteredbythefibres.Oneoftheimportantexamplesisthehydrodynamicsta…     

Fiber spinning of very dilute solutions of polyacrylamide in water

The extensional viscometer developed earlier by the authors was refined and used to extend very dilute (50 ppm) solutions of polyacrylamide in distilled water. A slender liquid filament was stretched by the use of a suction device, and this resulted in the spinning of the fiber. By varying the volumetric flow rate and the filament length, stretch rates in the 100–1000 s⁻¹ range were easily obtained. The corresponding tensile stresses were very large, and these gave apparent extensional viscosities of the order of 200 P (20 Pa s). In contrast to this, the material functions in shear were difficult to measure, except for the shear viscosity which showed pronounced shear thinning. It was found that all the measurements, in shear as well as extension, could be explained based on the four constant Johnson-Segalman constitutive equation.     

Fiber-optic probe measurements of void fraction and bubble size distributions beneath breaking waves

Experiments were performed to determine the accuracy of single-tip fiber-optic probes for making simultaneous measurements of the void fraction and bubble size distributions under breaking waves. Tests in a vertical bubble column showed that the normalized RMS error in the void fraction measurements was ∼10%. The relationship between bubble rise time and bubble velocity was investigated in a unidirectional flow cell. Similar to previous studies the rise time and bubble velocity were found to be related by a power law equation. The probes can provide mean bubble velocities accurate to ±10% when a minimum of ∼15 individual bubble velocities are averaged. The fiber-optic probes were deployed beneath a plunging breaking wave in a laboratory wave channel. The slope and shape of the bubble cord length size distribution measured with the probes was found to agree closely with the size distribution measured from digital video recordings. The probes were then positioned in the splash-up zone of a plunging breaker and the resulting cord length distribution had a shape and slope that was in agreement with previous measurements. These results demonstrate that single-tip fiber optic probes can provide accurate simultaneous measurements of the void fraction and bubble sizes inside the dense bubble clouds entrained by breaking waves.     

Non-linear viscoelastic response of magnetic fiber suspensions in oscillatory shear

This paper reports the first study on the large amplitude oscillatory shear flow for magnetic fiber suspensions subject to a magnetic field perpendicular to the flow. The suspensions used in our experiments consisted of cobalt microfibers of the average length of 37 μm and diameter of 4.9 μm, dispersed in a silicon oil. Rheological measurements have been carried out at imposed stress using a controlled stress magnetorheometer. The stress dependence of the shear moduli presented a staircase-like decrease with, at least, two viscoelastic quasi-plateaus corresponding to the onset of microscopic and macroscopic scale rearrangement of the suspension structure, respectively. The frequency behavior of the shear moduli followed a power-law trend at low frequencies and the storage modulus showed a high-frequency plateau, typical for Maxwell behavior. Our simple single relaxation time model fitted reasonably well the rheological data. To explain a relatively high viscous response of the fiber suspension, we supposed a coexistence of percolating and pivoting aggregates. Our simulations revealed that the former became unstable beyond some critical stress and broke in their middle part. At high stresses, the free aggregates were progressively destroyed by shear forces that contributed to a drastic decrease of the moduli. We have also measured and predicted the output strain waveforms and stress–strain hysteresis loops. With the growing stress, the shape of the stress–strain loops changed progressively from near-ellipsoidal one to the rounded-end rectangular one due to a progressive transition from a linear viscoelastic to a viscoplastic Bingham-like behavior.