Effect of a controlled pre-deformation history on extensional viscosity of dilute polymer solutions

We use a modified filament stretching rheometer to quantify the influence of a known controlled pre-shear history on the transient extensional viscosity of a dilute polymer solution. Two different types of pre-deformation are explored; both influence the subsequent stretching significantly, albeit in opposite ways. Small-amplitude oscillatory straining parallel to the direction of stretching enhances strain hardening and accelerates the tensile stress growth toward the steady-state value. Conversely, steady torsional shearing orthogonal to the direction of stretching retards strain hardening and results in a delayed approach to steady-state elongational flow. In both cases, the final steady-state extensional viscosity is the same as that observed with no pre-shearing. Calculations using a finitely extensible nonlinear elastic Peterlin dumbbell model qualitatively capture the trends observed in experiments, enabling interpretation of these observations in terms of the degree of polymer chain stretching imposed by the flow before extensional stretching.     

High-rate start-up behaviour of unvulcanised rubber compounds

The recently developed machines, which are shown to be comparable with “established” devices, are used to perform the start-up of uniaxial extension and simple shear at high constant deformation rates. Henky strains in excess of 8 at constant strain rates of over 100 s^?1 are achieved for two highly filled unvulcanised rubber compounds of Butyl and Chloroprene rubbers. In general no steady state is attained in uniaxial extension; the stress reaches a maximum and then declines. It is shown that the time—temperature superposition principle holds for start-up flows resulting in a family of master curves. It is established that the extensional stresses are one or two orders of magnitude greater than the shear stresses at corresponding deformation rates over a wide rate range.     

Dehnviskosität von Kautschukmischungen

It is shown that the extensional viscosity measured by uniaxial stretch may give useful information as to the processability of rubber compounds. A simple apparatus was constructed to measure this rheological property at constant rate of strain. The influence of rubber type, filler, aging, mastication and degree of crosslinking on the reduced stress is represented by diagrams. No steady-state extensional viscosity was obtained for highly filled rubber compounds. It was found that the flow behaviour of rubber compounds can be better characterized with an extensional rheometer than with a high-pressure capillary rheometer.

     

Assessment of nonlinear strain measures for extensional and shearing flows of polymer melts

From stress-strain experiments in extensional and shearing flows, nonlinear strain measures and effective damping functions are derived for a polyisobutylene melt. The strain measures determined in planar extensional flow and in simple shear flow coincide. Experimental results are compared with predictions of two molecular theories, the Doi-Edwards model and the molecular stress function approach of Wagner and Schaeffer. Discrepancies between theories and experiment lead to a reconsideration of the classification of extensional flows. The symmetry of the flow field is identified and quantified as an important parameter influencing the strain measure, and a unifying strain measure for general extensional and shearing flows of polymer melts is presented.     

Computational analysis of techniques to determine extensional viscosity from entrance flows

Recent computational analysis of entrance flows (Mitsoulis et al. 1998) suggests that the entrance pressure drop is insensitive to large changes in steady extensional viscosity-a result that directly contradicts a large body of experimental work in this area. A re-examination of entrance flows using numerical simulations is presented in this work which shows that entrance pressure drops do depend on the steady extensional viscosity, provided the extension rate in the entrance flow is large enough. Numerical simulations are presented using both the strain thinning and thickening versions of the Phan-Thien–Tanner (PTT) constitutive model. Several techniques for extracting extensional viscosity from entrance pressure are applied to the results of these simulations. The resulting predictions of extensional viscosity are compared to the steady extensional viscosity curves predicted by the PTT constitutive model used to generate the simulated pressure drop curves. The analytical techniques examined here are shown to provide reasonably accurate estimates of the steady extensional viscosity. This work also clearly demonstrates the advantage of using variable power-law coefficients for the rheological properties, used as inputs to the analyses, to capture the extensional behavior at deformation rates below the power law region more accurately. Received: 23 July 1999/Accepted: 24 November 1999     

The triple jet: influence of shear history on the stretching of polymer solutions

A kerosene-based aircraft safety fuel and aqueous solutions of poly (ethylene oxide) and polyacrylamide are examined using the “triple jet” system. This device allows the solution to be stretched as it flows from a capillary tube and the axial stress, strain and strain rate in the liquid are measured.The shear history of the solution is altered by placing cylindrical inserts in the capillary tube. This is shown to have a large effect on the extensional behaviour of aircraft safety fuel, a moderate effect on the extensional behaviour of poly (ethylene oxide) solution and little effect on the behaviour of polyacrylamide solution. The extensional viscosity of the aircraft fuel is raised by an order of magnitude when a long period of high shear is used; the effects last for periods of up to one second, though traditional methods suggest a relaxation time of the order of 10^?3 seconds. A liquid of shear viscosity 4 centipoise may have an extensional viscosity of over 100 poise.Plots of the extensional modulus of the jet as a function of distance along the jet emphasize the importance of shear history for the first two types of solution and suggest that the latter stages of the stretching process are elastic in character. Typical extensional moduli for the solutions tested are in the range 1.3–5.0 × 10⁴ dyn.cm^?2.The relevance of the interplay between shearing and stretching flow to the phenomena of lubrication and turbulence suppression is mentioned.     

The reptation model with segmental stretch

The Doi-Edwards model with segmental stretch and a non-linear finitely extensible spring law is described and examined in simple flow situations where analytic results are derivable; namely oscillatory flow and steady state flow at high deformation rates. The model is shown to be consistent with the Bueche-Ferry hypothesis in fast large strain unidirectional flows but to violate this rule in small strain reversing flows. The discrepancy is identified with a preaveraging approximation used to describe the relative tube-chain velocity. Experimentally verifiable scaling rule for the birefringence as a universal function of a planar flow-type parameter and deformation rate are identified. Sensitivity to the extensional flow character, absent in the original tube model, manifests itself with the introduction of segmental stretch. Although the model generates a non-separable memory function kernel the deformation dependence of the memory function is quantitatively shown to have negligible impact on the predicted theological properties relative to the original Doi-Edwards model. With this simplification, relatively uncomplicated approximations to the segmental stretch model can be deduced.     

Birefringence and stress growth in uniaxial extension of polymer solutions

Simultaneous measurements of extensional stresses and birefringence are rare, especially for polymer solutions. This paper reports such measurements using the filament stretch rheometer and a phase modulated birefringence system. Both the extensional viscosity and the birefringence increase monotonically with strain and reach a plateau. Estimates of this saturation value for birefringence, using Peterlin’s formula for birefringence of a fully extended polymer chain are in agreement with the experimental results. However, estimates of the saturation value of the extensional viscosity using Batchelor’s formula for suspensions of elongated fibres are much higher than observed. Reasons for the inability of the flow field to fully unravel the polymer chain are examined using published Brownian dynamics simulations. It is tentatively concluded that the polymer chain forms a folded structure. Such folded chains can exhibit saturation in birefringence even though the stress is less than that expected for a fully extended molecule.Simultaneous measurements of stress and birefringence during relaxation indicate that the birefringence decays much more slowly than the stress. The stress-birefringence data show a pronounced hysteresis as predicted by bead-rod models. The failure of the stress optic coefficient in strong flows is noted.Experiments were also performed wherein the strain was increased linearly with time, then held constant for a short period before being increased again. The response of the stress and birefringence in such experiments is dramatically different and can be traced to the different configurations obtained during stretching and relaxation. The results cast doubt on the appropriateness of pre-averaging the non-linear terms in constitutive equations.     

Review of the entry flow problem: Experimental and numerical

This paper is concerned with a review of both experimental and numerical studies of axisymmetric and planar entry flows which have been considered as test problems for the numerical simulation of viscoelastic fluids. The test of the method is usually based upon whether the numerical model predicts vortices in the entry corners. However, it is not clear as to whether one should observe vortices for all viscoelastic fluids. Polyacrylamide solutions and Boger fluids exhibit vortices in axisymmetric flow and the size of the vortex does increase with fluid elasticity. However, the vortex is nearly suppressed in planar entry flow. On the other hand, not all polymer melts are found to exhibit vortices in either axisymmetric or planar entry flow. It is our belief that the origin of vortices is not related to the elasticity based on shear flow propertes but to the behavior of the transient extensional viscosity. Certain polymer melts such as low density polyethylene exhibit vortices in both planar and axisymmetric flow along with unbounded stress growth at the start up of extensional flow. It is believed that the constitutive equations used in the numerical simulation must reflect this extensional behavior if vortices are to be predicted. A review of the numerical simulations concerned with entry flow shows that there is considerable doubt about the accuracy of the predictions for most of the studies. Even for those where the numerical solution is thought to be accurate, the magnitude of the stream function associated with the vortices is usually very low. None of the differential models used to date predicts strain hardening extensional viscosity, but those which are thought to predict vortices do rise more rapidly to the steady-state extensional viscosity values with time. It is recommended that the search of test fluids be widened beyond polymer solutions as there may already exist a number of polymer melts which behave similarly to the predictions of existing constitutive equations.     

A new look at extensional rheology of low-density polyethylene

The nonlinear rheology of three selected commercial low-density polyethylenes (LDPE) is measured in uniaxial extensional flow. The measurements are performed using three different devices including an extensional viscosity fixture (EVF), a homemade filament stretching rheometer (DTU-FSR) and a commercial filament stretching rheometer (VADER-1000). We show that the measurements from the EVF are limited by a maximum Hencky strain of 4, while the two filament stretching rheometers are able to probe the nonlinear behavior at larger Hencky strain values where the steady state is reached. With the capability of the filament stretching rheometers, we show that LDPEs with quite different linear viscoelastic properties can have very similar steady extensional viscosity. This points to the potential for independently controlling shear and extensional rheology in certain rate ranges.     

Suspensions of monodisperse spheres in polymer melts: particle size effects in extensional flow

Suspensions in polymeric, viscoelastic liquids have been studied in uniaxial extensional flow. The fibre wind-up technique has been used for this purpose. The effects of particle size and particle volume fraction have been investigated, using monodisperse, spherical particles. The results have been compared with shear flow data on the same materials. The values of the relative extensional viscosities at low stretching rates are in agreement with the relative shear viscosities and relative moduli. This indicates that hydrodynamic forces are stronger than the particle interaction forces. At larger strain rates strain hardening occurs; it is suppressed when particles are added. Small aggregating particles reduce the strain hardening more strongly than larger particles; strain hardening can even be totally eliminated. When further increasing the stretching rate, hydrodynamic effects dominate again and the effect of particle size effect on strain hardening disappears.     

Mesoscopic dynamics of polymer chains in high strain rate extensional flows

We take a step towards accessing the physics of viscoelastic liquid breakup in high speed, high strain rate flows by performing Brownian dynamics computations of dilute uniaxial, equibiaxial, and ellipsoidal polymeric extensional flows. Our computational implementation of the bead-spring model, when tailored to the DNA molecule, consistently with recent works of Larson and co-workers, is shown: (a) to predict a coil-stretch transition at Deborah number De≈0.5, and (b) to reproduce the experimental longest relaxation time. Furthermore, after adapting the model parameters to represent the polyethylene oxide (PEO) chain (for M=10⁶ Da), we find it possible to reproduce our own experimental data of the longest relaxation time, the transient extensional viscosity of dilute solutions at small Deborah numbers, and a coil-stretch transition at Deborah number De≈0.5. Extended to large Deborah numbers, the model predicts that polymer stretching is controlled by: (a) the randomness of the initial conditions that, in combination with rapid kinematically imposed compression, lead to the formation of initially frozen chain-folds, and (b) the speed with which thermo-kinematic processes relax these folds. The slowest fold relaxation occurs during uniaxial extension. As expected, the introduction of stretching along a second direction enhances the efficiency of fold relaxation mechanisms. Even for Deborah numbers (based on the chain longest relaxation time) of the order of one thousand, there is a large variation in the time a polymer needs in order to extend fully, and the effects of Brownian motion cannot be ignored. The computed Trouton ratios and polymer contributions to the total stress as functions of Hencky strain provide information about the relative importance of elastic effects during polymeric liquid stretching. At high strain rates, the steady state elastic stresses increase linearly with the Deborah number, resembling at that stage an anisotropic Newtonian fluid (constant extensional viscosity).     

The effect of pre-shear on the extensional rheology of wormlike micelle solutions

The effect of initial microstructural deformation, alignment, and morphology on the response of wormlike micelle solutions in transient uniaxial extensional flows is investigated using a pre-shear device attached to a filament stretching rheometer. In filament stretching experiments, increasing the strength and the duration of the pre-shear just before stretch is found to delay the onset of strain hardening. In these experiments, the wormlike micelle solution filaments fail through a rupture near the axial midplane. The value of the elastic tensile stress at rupture is found to decrease with increasing pre-shear rate and duration. The most dramatic effects are observed at shear rates for which shear banding has been independently observed. The reduction in the strain hardening suggests that pre-shear before filament stretching might break down the wormlike micelles reducing their size before stretch. Strain hardening is also observed in capillary breakup rheometry experiments; however, the pre-sheared wormlike micelle solutions strain harden faster, achieve larger steady-state extensional viscosities and an increase in the extensional relaxation time with increasing shear rate and duration. The difference between the response of the wormlike micelles in filament stretching and capillary breakup experiments demonstrates the sensitivity of these self-assembling micelle networks to pre-conditioning.     

Dynamics of worm-like micelles: the Cox–Merz rule

The viscoelastic behaviour of worm-like micelles in small-amplitude oscillatory, steady simple shear and uniaxial extensional flows are analyzed with a model that couples the Oldroyd-B constitutive equation with a kinetic equation that accounts for the structural changes induced by the flow. In some cases, the constitutive equation predicts a viscoelastic behaviour that is consistent with the Cox–Merz rule. Departures from this rule are also predicted. Experimental data obtained for two worm-like micellar systems indicate that in these solutions, the Cox–Merz rule is not usually followed, in agreement with the predictions of our model. In uniaxial extensional flow, the model predicts a strain hardening in the extensional viscosity at low extensional rates and a strain-thinning at high extensional rates.     

The relevance of entry flow measurements for the estimation of extensional viscosity of polymer melts

The extensional viscosity of some flexible chain polymers and a thermotropic liquid crystalline polymer was measured in uniaxial extensional flow at constant extension rate. Power law functions were found for the dependence of the extensional viscosity at constant accumulated strain on strain rate. The stress growth curves were compared with measurements in axisymmetric entry flow, where both elongation and shear occur. The comparison showed that the values of the extensional viscosity calculated from the measurements in the entry flow correspond to the ones calculated from the viscosity growth measured in uniaxial elongation and averaged over extensional strain equal to what is accumulated on the fluid as it flows from the barrel into the capillary.     

Elongational flow and rheology of monodisperse polymers in solution

We describe the utilization of idealized stagnation point extensional flows, produced by opposed jets, for birefringence visualization of induced molecular strain and flow resistance measurements. We identify rheological changes associated with the coil---stretch transition which occurs beyond a critical strain-rate in elongational flow-fields. In dilute solutions of monodisperse atactic polystyrene, increases in extensional viscosity are observed as isolated molecules become stretched. The largest increases in extensional viscosity, however, are found only beyond a critical concentration and strain rate, and are associated with the stretching of transient networks of interacting molecules. These results parallel similar effects seen in porous media flow and capillary entrance experiments. We determine the molecular weight dependence of the critical concentration which scales as M^−0.55 in agreement with pairwise interaction of coils, but is much lower than conventional values of the critical polymer concentration, c^*. We believed that polydispersity may play an important role in the development of such transient networks, and in controlling the degradation behaviour during flow.     

Equibiaxial extensional flow of polymer melts via lubricated squeezing flow. II. Flow modeling

A model for lubricated squeezing flow of a viscoelastic fluid is developed in order to study the viability of this flow as a rheological technique for generating equibiaxial extensional deformations in polymer melts. In this simple flow model, the melt, described by an upper-convected Maxwell fluid, is squeezed between thin films of a Newtonian fluid. Comparisons of the model predictions for constant strain rate and constant stress flows are made with experimental results presented in the first paper. Predictions from the model are able to describe the effects of lubricant viscosity and experimental configuration and indicate the technique fails for these flows at Hencky strains of approximately one. The cause for this failure is lubricant thinning, which leads to significant errors in both the measured stress difference and the strain. Received: 31 January 2000 Accepted: 31 May 2000     

Evaluation by means of stress relaxation (after a step strain) experiments of the viscoelastic behavior of polymer melts in uniaxial extension

 As is widely acknowledged, morphology in most materials is far more sensitive to extensional than to shear deformations but, unfortunately, due to the experimental difficulties involved, there are no non-destructive, morphology probing techniques in such flows, i.e., the equivalent of stress relaxation and oscillatory experiments in shear flows. This paper tries to overcome some of those drawbacks by proposing an experimental technique that allows stress relaxation experiments after a step strain in uniaxial extension to be performed. The benefits of this technique are twofold: (a) while the deformation is small enough for the response to be in the linear viscoelastic regime it constitutes a probe of the microstructure of the material and (b) it allows the departure to the non-linear regime to be studied, useful, for example, for the definition of the damping function in uniaxial extensional flow or for the study of the response of materials to fast transient flows with a strong extensional component, such as contraction flows. In this work the proposed technique, which requires a correction to the apparent (theoretical) strain rate in order to allow the calculation of the true Hencky strains attained during the strain step, is tested and validated for two polyisobutylene melts. Received: 9 April 2001 Accepted: 26 July 2001     

Consistent hybrid finite volume/element formulations: model and complex viscoelastic flows

The accuracy and consistency of a new cell‐vertex hybrid finite element/volume scheme are investigated for viscoelastic flows. Finite element (FE) discretization is employed for the momentum and continuity equation, with finite volume (FV) applied to the constitutive law for stress. Here, the interest is to explore the consequences of utilizing conventional cell‐vertex methodology for an Oldroyd‐B model and to demonstrate resulting drawbacks in the presence of complex source terms on structured and unstructured grids. Alternative strategies worthy of consideration are presented. It is demonstrated how high‐order accuracy may be achieved in steady state by respecting consistency in the formulation. Both FE and FV spatial discretizations are embedded in the scheme, with FV triangular sub‐cells referenced within parent triangular finite elements. Both model and complex flow problems are selected to quantify and assess accuracy, appealing to analysis and experimental validation. The test problem is that of steady sink flow, a pure extensional flow, which reflects some of the numerical difficulties involved in solving more generalized viscoelastic flows, where both source and flux terms may contribute equally to stress propagation. In addition, a complex transient filament‐stretching flow is chosen to compute the evolution of stress fields within liquid bridges. Shortcomings of the various stress upwinding schemes are discussed in this context, whilst dealing with such free‐surface type problems. Here, stress fluctuation distribution alone is advocated, and a Lax‐scheme is found to deliver accuracy and stability to the computational results, comparing well with the literature. Copyright © 2004 John Wiley & Sons, Ltd.