Instabilities in multi-hole converging flow of viscoelastic fluids

Studies of the onset of instabilities were conducted on single hole and multi-hole contractions using laser speckle visualization. A well characterized elastic fluid was used with constant viscosity of 13.1 Pa · s and elasticity characterized by a longest relaxation time constant of 2.233 s. The onset of instabilities was characterized in terms of the Deborah number and the contraction ratio. Three types of instabilities were observed: pulsing vortices, azimuthally rotating vortices, and swirling vortices. For the single hole contractions the critical Deborah number for instability increased from 4.4 to 5.07 to 5.25 as the contraction ratio increased from 4: 1 to 8: 1 to 12: 1. The magnitude of the instabilities was much greater for the 4: 1 contraction than for the other two contraction ratios. For the multi-hole contraction a square array of nine holes was used and the ratio of the hole diameter to hole spacing was varied. The height of the vortices is very similar for the single hole and multi-hole contractions at low Deborah numbers. At high Deborah numbers the effect of adjacent holes is to reduce the height of the vortices by a factor of three. For the 4: 1 spacing no secondary vortex was observed below a Deborah number of De = 3.7. Secondary vortices occurred for the 8:1 and 10:1 spacing at all Deborah numbers. Unstable pulsing vortices appeared for all spacings at a critical Deborah number around 5.5. Adjacent holes decreased the strength of the unsteady vortex motions. The centerline velocities were measured for the multi-hole contraction at shear rates of 5, 30, and 300 s^–1. The elongational strain rates are similar at a low shear rate of 5 s^–1. As shear rate is increased the onset of stretching occurs closer to the plane of the contraction for the smaller contraction ratios.     

Analysis of Taylor vortex flow by means of laser light scattering

Rheological measurements and light-scattering experiments were performed on dilute solutions of high molecular polystyrene. We are able to describe the orientation behavior of chain molecules under shear flow by means of light-scattering. Beyond that these investigations of light-scattering of flowing polymer solutions are an useful and suitable tool for detection and characterization of Taylor vortex formation. We can estimate the appearance of these hydrodynamic instabilities, which overlay the laminar main flow and we can observe a typical influence of the solvent power on it.Presented in part at the meeting of the Deutsche Rheologische Gesellschaft, Berlin, 13–15 May, 1991.     

The evolution of linear viscoelasticity during the vulcanization of polyethylene

The evolution of linear viscoelasticity during the vulcanization of polyethylene is studied through the gel point. The material in the vicinity of the gel point is described by two scaling laws: one characterizes the viscoelasticity at the critical point and a second characterizes the evolution of viscoelasticity near the gel point. Time Resolved Mechanical Spectroscopy is used to observe both scaling phenomena. The material at the gel point displays power law relaxation over five decades of time with a power-law relaxation exponent equal to 0.32. This study conforms with previous findings that this exponent is composition-dependent. The longest relaxation time diverges in the vicinity of the gel point as _max |p _c - p| ^–1/, and we find = 0.2. This result conforms with previous reports that this exponent may be independent of composition. The Arrhenius flow activation energy for this material undergoes three-fold changes during crosslinking up to the gel point. A single-adjustable-parameter stretched-exponential-power law relaxation function is an inadequate representation of crosslinked materials over any significant range of extent of the reaction up to the gel point.     

Numerical simulation of delayed die swell

In a recent paper, Joseph et al. showed that, for a number of viscoelastic fluids, one can observe the phenomenon of delayed die swell beyond a critical extrusion velocity, or beyond a critical value of the viscoelastic Mach number. Giesekus had also observed that delayed die swell is a critical phenomenon.In the present paper, we find a set of material and flow parameters under which it is possible to simulate delayed die swell. For the viscoelastic flow calculation, we use the finite element algorithm with sub-elements for the stresses and streamline upwinding in the discretized constitutive equations. For the free surface, we use an implicit technique which allows us to implement Newton's method for solving the non-linear system of equations. The fluid is Oldroyd-B which, in the present problem, is a singular perturbation of the Maxwell fluid. The results show very little sensitivity to the size of the retardation time. We also show delayed die swell for a Giesekus fluid.This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

The effect of dissolved polymers on the rheological properties of coating colours

Coating colours used for the coating of paper and board consist mainly of a mineral pigment, which is very often clay, a synthetic binder such as a styrenebutadiene latex, dispersion agents and water retention aids. The latter are often water soluble polymers. These polymers have a very strong influence on the rheological properties of the coating colours, both on the strain rate dependence of the apparent viscosity and on the viscoelasticity. The effects of two different grades of carboxymethylcellulose (CMC) and one grade of hydroxyethylcellulose (HEC), on the rheological properties at room temperature of a clay-based coating colour at pH 8, have been investigated. It is concluded that the high values of the dynamic modulus of the colours are due to interactions between the cellulose derivatives and the solid particles, i.e. mainly the clay particles. For HEC this interaction is associated with adsorption of the polymeric molecules on the clay particles. In the case of CMC, the adsorption is strongly retarded by the presence of the dispersant (a polyacrylate salt). It is suggested that the marked elasticity of the CMC-containing colour in addition to a possible polymer adsorption may be due to charge interactions and/or depletion flocculation. The effect of CMC and HEC on the water-retention properties of the colour is also discussed.     

Measurement of rheological and ultrasonic properties of food and synthetic polymer melts

A slit die viscometer has been used in conjunction with a co-rotating twin screw extruder to study the rheological behaviour of maize grits, potato powder and low density polyethylene, as a function of feed rate, screw speed and temperature. The shear viscosity of both maize and potato decreased with increasing feed rate. Increasing the temperature or screw speed at any given feed rate also reduced the viscosity. The ultrasonic velocity through the material has also been shown to be sensitive to the extruder operating conditions. Overall, the ultrasonic velocity decreased as screw speed and temperature increased and feed rate decreased.     

Characterization of yield stress fluids with concentric cylinder viscometers

The methods normally employed for shear rate calculations from concentric cylinder viscometer data generally are not applicable for fluids with a yield stress. In cylindrical systems with large radius ratios, as usually is the case with suspensions, the yield stress induces two possible flow regimes in the annulus. Unless the yield value is exceeded everywhere in the gap only part of the fluid can be sheared while the remaining region behaves like a solid plug. A correct calculation of the shear rate must take into account the presence of a variable effective gap width determined by the extent of the sheared layer. For time-independent yield stress fluids, a two-step procedure, which does not require any specific flow model, is proposed for analysing the experimental torque-speed data. Under the partially sheared condition, the shear rate can be computed exactly, whereas for the fully sheared flow the Krieger and Elrod approximation is satisfactory. The method is assessed by examining both semi-ideal data generated with a Casson fluid with known properties, and experimental data with an industrial suspension. A more complicated problem associated with characterization of time-dependent yield stress fluids is also identified and discussed. An approximate procedure is used to illustrate the dependence of the shear rate on time of shear in constant-speed experiments.     

Dynamic-mechanical properties of interfacially modified glass sphere filled polyethylene

Model composites of spherical glass particles dispersed in a matrix of high density polyethylene were prepared both with and without interfacial modification by an azidofunctional trialkoxysilane. Dynamic mechanical measurements of the composites in the melt state were recorded. The unmodified composites behave as theoretically predicted and the effect of particle—particle interaction at high volume fractions can be measured. The composites with a modified interfacial region have greater shear moduli due to the effect of a region surrounding the particle modified by the silane. The material in this region is largely bound to the glass surface and was examined by Fourier transform infrared spectroscopy after extraction of the bulk matrix. Theoretical calculations are shown to be useful in calculating the mechanical properties and volume fraction of the interfacial region.     

Rheological characteristics of sewage sludge: A granuloviscous material

Problems have been experienced in the screw conveying of a centrifuge dewatered sewage sludge, at about 66% by weight water content. The rheology of the sludge was therefore tested in a viscometer in which the normal stress between the platens could be controlled. Rough platens were used to measure internal shear strengths while smooth steel platens were used to study wall slip.The sludge displayed a number of the characteristics of a granuloviscous material. It was found to possess high initial shear strength, but rapidly softened with shear to lower steady-state strengths. The initial shear strength increased with increasing normal stress, indicating that the sludge was frictional. The initial strength was lower against the smooth platen.The steady-state shear strengths were both rate and stress dependent using the rough platen, and this internal shear strength could be described by a frictional pseudoplastic equation. Using the smooth platen the material was nearly rate independent, and was best described by a Coulomb friction equation. At the higher rates the wall shear strengths were lower than the internal shear strengths.The strengths (both initial and steady-state) were characterised by a high degree of variability from sample to sample, and also showed marked hysteresis with the rough platen. They did not show any consistent variation with moisture content.In terms of screw conveying, the major conclusion is that smooth walls to the flights of the screw and barrel of the conveyor are desirable.     

In situ technique for monitoring the gelation of UV curable polymers

Ultraviolet-radiation-cured cross-linked systems are used extensively in optoelectronic applications. We describe a new in situ technique for doing controlled curing of photosensitive materials and monitoring the gelation through dynamic theological measurements. In this method, a sample placed between quartz windows of a modified parallel plate fixture of a Rheometrics Mechanical Spectrometer is radiated from a collimated UV source using a liquid light guide and a mirror. By varying the sample exposure to UV radiation, the degree of cross-linking can be precisely controlled to observe the material behavior in the pre, post, and critical gel state. Time-dependent measurements on urethane based materials show the gelation behavior to have an induction period with no change in dynamic moduli (G, G), followed by their sharp increase as the materials develop a network structure. The critical gel point is characterized by a power-law dependence of the dynamic moduli on frequency. The gelation kinetics scales with sample thickness and radiation intensity, the scaling factor in both cases being the critical gelation time. The rheological measurements correlate with differential photocalorimetric studies.     

On the cone and plate deformation of a rubber-like material: Instability

In this report it is shown that the place boundary-value problem, for a small deformation superimposed on the large cone-and-plate deformation of a Mooney-Rivlin material, has no unique solution at some displacement angles of the cone. The kinematics of the small deformation are chosen such that the problem is reduced to a set of ordinary differential equations. With respect to this restricted class of kinematics the cone-and-plate deformation is unstable.     

The flow past a sphere in a cylindrical tube: effects of intertia,shear-thinning and elasticity

Numerical solutions are presented for the flow past a sphere placed at the centreline of a cylindrical tube for Reynolds numbers ranging from 0 to 150, using a boundary element method. Fluids are modelled by a variety of constitutive equations including the Newtonian, the Carreau and the Phan-Thien-Tanner models. The influence of inertia, shear-thinning and fluid elasticity on the flow field, drag and the pressure drop force-drag ratio is examined. Some results are compared with available experimental data.     

Fragile networks and rheology of concentrated suspensions

Suspensions consisting of particles of colloidal dimensions have been reported to form connected structures. When attractive forces act between particles in suspension they may flocculate and, depending on particle concentration, shear history and other parameters, flocs may build-up in a three-dimensional network which spans the suspension sample. In this paper a floc network model is introduced to interpret the elastic behavior of flocculated suspensions at small deformations. Elastic percolation concepts are used to explain the variation of the elastic modulus with concentration. Data taken from the suspension rheology literature, and new results with suspensions of magnetic -Fe₂O₃ and non-magnetic -Fe₂O₃ particles in mineral oil are interpreted with the model proposed.Non-zero elastic modulus appeared at threshold particle concentrations of about 0.7 vol.% and 0.4 vol.% of the magnetic and non-magnetic suspensions, respectively. The difference is attributed to the denser flocs formed by magnetic suspensions. The volume fraction of particles in the flocs was estimated from the threshold particle concentration by transforming this concentration into a critical volume concentration of flocs, and identifying this critical concentration with the theoretical percolation threshold of three-dimensional networks of different coordination numbers. The results obtained indicate that the flocs are low-density structures, in agreement with cryo-scanning electron micrographs. Above the critical concentration the dynamic elastic modulus G was found to follow a scaling law of the type G ( _f - _f ^c ) ^f , where _f is the volume fraction of flocs in suspension, and _f ^c is its threshold value. For magnetic suspensions the exponent f was found to rise from a low value of about 1.0 to a value of 2.26 as particle concentration was increased. For the non-magnetic a similar change in f was observed; f changed from 0.95 to 3.6. Two other flocculated suspension systems taken from the literature showed a similar change in exponent. This suggests the possibility of a change in the mechanism of stress transport in the suspension as concentration increases, i.e., from a floc-floc bond-bending force mechanism to a rigidity percolation mechanism.     

Direct measurement of static yield properties of cohesive suspensions

A technique of yield stress investigation based upon the combined use of two devices (an applied stress rheometer and an instrument for measuring the propagation velocity of small amplitude, torsional shear waves) is described. Investigations into the low shear rate rheological properties of illitic suspensions are reported for shear rates, typically, in the range 10^–4— 10^–1 s^–1 under applied stresses in the range 0.01 — 10 Nm^–2 and involving shear strains between 10^–1 and 10^–4. Results are presented which demonstrate that the technique does not invoke the excessive structural disruption of material associated with applied shear rate based methods (direct and otherwise) and the widely encountered problem of wall slip at the surface of rotational measuring devices is avoided using miniature vane geometries. Results are compared with those obtained using smooth-walled cyclindrical measuring devices in both applied stress and applied shear rate instruments.Yield measurements are considered in relation to the structural properties of the undisturbed material state and shear moduli obtained by studying the propagation of small amplitude (10^–5 rad), high frequency (~ 300 Hz) torsional shear waves through the test materials are reported. Experimental techniques and instrument modifications to permit these measurements are described.     

Der elektroviskose Effekt als Folge elektrostatischer Kraft

Zusammenfassung Die Viskosität vieler Flüssigkeiten ändert sich als Folge elektrischer Felder. Die größten elektroviskosen Effekte findet man unter äußerem Feld bei Suspensionen aus dielektrischen Flüssigkeiten und feinen Festkörperteilchen.Die Flüssigkeiten zeigen ausgeprägtes Bingham-Verhalten. Solange die auf die Wand ausgeübte Scherspannung kleiner als die von der Feldstärke abhängige Haftspannung ist, verhalten sich die Materialien wie Festkörper. Man beobachtet, daß Elektroviskosität proportional zur Feldstärke im Quadrat und umgekehrt proportional zur Schergeschwindigkeit ist, und daß sich die Teilchen der Flüssigkeiten in Richtung der Feldstärke orientieren und faserige Strukturen bilden. In dieser Arbeit erklären wir diese Eigenschaften durch ein einfaches Modell. In diesem Modell resultiert die Scherkraftzunahme aus der elektrostatischen Kraft, die sich als Gradient der im Feld abgespeicherten Energie ergibt.The electroviscosity effect following electrostatic force

---

Apparent viscosity of many fluids changes under an electric field. The largest effects are found with dispersions of fine particles in dielectric liquids subjected to an externally applied field. The large effect is apparently due to induced fibration when the particles align themselves as fibers in field direction. The dispersions exhibit Bingham-behavior, i.e., as long as the shear stress is smaller than the field-dependent yield stress, the material behaves as a solid, a behavior which can be made plausible as the electrodes are mechanically linked together by the fibers (Winslow 1949). In this paper arguments are put forward that the induced fibration is the consequence of electrostatic forces acting on the particles and pulling them into regions of higher field strength. The resulting configuration is an equilibrium configuration and it is postulated that even under shear load, the electric field tends to maintain this configuration such that only the particles close to the electrodes are sheared off, which then form a suspension near the electrode surfaces. A phenomenological theory based on the assumption that the fibers are columns of rectangular particles gives quantitative results for the electroviscosity and the yield stress which agree with experimental observations.

     

Bimodal model of slurry viscosity with application to coal-slurries. Part 1. Theory and experiment

The rheological behavior of stable slurries is shown to be characterized by a bimodal model that represents a slurry as made up of a coarse fraction and a colloidal size fine fraction. According to the model, the two fractions behave independently of each other, and the non-Newtonian behavior of the viscosity is solely caused by the colloidal fraction, while the coarse fraction increases the viscosity level through hydrodynamic interactions. Data from experiments run with colloidal coal particles of about 2–3 µm average size dispersed in water show the viscosity of these colloidal suspensions to exhibit a highly shearrate-dependent behavior and, in the high shear limit, to match very closely the viscosity of suspensions of uniform size rigid spheres although the coal volume fraction must be determined semi-empirically. Different amounts of coarse coal particles are added to the colloidal suspension and the viscosity of the truly bimodal slurries measured as a function of shear rate. In agreement with the bimodal model, the measured shear viscosities show the coarse fraction to behave independently of the colloidal fraction and its contribution to the viscosity rise to be independent of the shear rate. It is shown that the shear rate exerted on the colloidal fraction is higher than that applied by the viscometer as a result of hydrodynamic interactions between the coarse particles, and that it is this effective higher shear rate which is necessary to apply in the correlations. For determining the coal volume fraction a relatively simple and quite accurate measurement technique is developed for determining the density and void fraction of coarse porous particles; the technique directly relates volume fraction to mass fraction.     

The rotation of a suspended axisymmetric ellipsoid in a magnetic field

Under the influence of a uniform and parallel magnetic field, a ferromagnetic fiber suspended in a Newtonian fluid rotates to align with the field direction. This study examines the field-induced rotation process for an individual non-Brownian axisymmetric ellipsoid suspended in a stagnant Newtonian fluid. Theoretical predictions are derived by a perturbation analysis for the limiting case where the strength of the applied magnetic field far exceeds the saturation magnetization of the ellipsoid. Numerical calculations are performed for the more general problem of an ellipsoid with known isotropic, non-hysteretic magnetic properties, using nickel and a stainless steel as examples. The analysis encompasses materials with field-induced, nonlinear magnetic properties, distinguishing these results from the simpler cases where the particle magnetization is either independent of, or linearly dependent on, the strength of the applied external field. In this study, predictions indicate that when the ellipsoid is magnetically saturated, the particle rotation is governed by the magnetoviscous time constant, _MV = _s/₀ M _s ² . It is found that the rotation rate depends strongly on the aspect ratio,a/b, of the ellipsoid, but only weakly on the dimensionless magnetization,M _s/H ₀. A geometric parameter for an ellipsoid, defined in eq. (2.5) - a, b major, minor semi-axes of an axisymmetric ellipsoid - D demagnetization tensor for an ellipsoid - D ^M magnetometric demagnetization tensor, the volume-average ofD ^P (r) - D ^P (r) position dependent demagnetization tensor, implicitly defined in eq. (2.12) - D _xx,D _yy,D _zz demagnetization factors, the diagonal elements ofD. Values for ellipsoids are defined in eq. (2.15) - F ^(m) magnetic force exerted on a body in a magnetic field - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - h _i ;h _ix,h _iy Cartesian components of dimensionless internal magnetic field,h _i =H _i /H ₀ - I moment of inertia tensor - k geometric parameter for hydrodynamic resistance of a body rotating in a Newtonian fluid given in eq. (2.3) - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component of the hydrodynamic torque - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.10); thez-component of the magnetic torque - M;M the magnetization, or dipole moment density, of a magnetic material; the magnitude ofM - M _s the saturation value ofM, approached by all ferromagnetic materials asH _i becomes large (figure 3) - m _s the dimensionless saturation magnetization,M _s/H ₀ - r position vector of a point within a ferromagnetic body - s dummy integration variable in eq. (2.5) - t time - U magnetoquasistatic potential energy of a magnetic body in a magnetic field, given in eq. (2.8) - u curve-fitting variable in eq. (4.1);u = logH _i - V volume of a magnetic particle; for an axisymmetric ellipsoid,V = (4/3) ab ² - x, y, z rectangular coordinate axes fixed in the ellipsoid (figure 1) - angle of inclination of the major axis of the ellipsoid with respect toH ₀ - _s viscosity of the Newtonian suspending medium - µ ₀ the magnetic permeability of free space,µ ₀ =4 · 10^–7H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - the magnetic susceptibility of a magnetic material, = M/H _i - ; _z angular velocity of a rotating body; angular velocity about thez-axis of an ellipsoid, _z=–d/dt     

A novel geometry for rheological characterization of viscoelastic materials

A novel geometry for generating a viscometric flow presents advantages of both cone and plate and parallel plate geometries, regarding uniform shear field and adjustable range of measurement. Kinematics and dynamics of the generated flow have been described mathematically utilizing an orthogonal curvilinear coordinate system based on the shapes of the shearing surfaces which are similar to the surface that generates the flow. Simple equations that allow the calculation of quantities of experimental interest in the rheological characterization of liquid materials, namely, shear rate, shear stress and two normal stress differences, have also been derived.The geometry, called pseudosphere, was tested with two types of fluids (Newtonian and pseudoplastic). Results show that the geometry can be used with low viscosity liquids (Newtonian liquids) by only adjusting the gapH. The behavior of pseudoplastic fluids for both low and moderately high viscosity could also be studied with this geometry. Very reproducible results were obtained when compared with those obtained with cone and plate geometry. Regions of lower shear rate could be studied using only the pseudosphere geometry.     

Anisotropic conduction of heat in a flowing polymeric material

In this paper a theory is presented which relates the thermal conductivity tensor of an amorphous polymeric material to the history of deformation of the material. The basis of the theory is formed by the network theory for polymeric materials. It will be shown that the results obtained here are in good agreement with experimental results on rubber. The effect of anisotropic heat conduction on the flow of a polymeric material will be demonstrated by the simple example of viscous heating in shear flow.Presented at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.     

Using regularization methods for the determination of relaxation and retardation spectra of polymeric liquids

The application of Tikhonov regularization to the determination of relaxation and retardation spectra of viscoelastic fluids is discussed using simulated and real experimental data. It is thereby shown that with this method consistent results for relaxation and retardation spectra can be obtained from experimental data for different material functions. Furthermore, the differences between the most frequently used variants of Tikhonov regularization and the maximum entropy method are discussed. For most calculations the program FTIKREG has been used, which is an especially fast and reliable implementation of Tikhonov's regularization method.