InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Materials Science » Fluids » "Viscoelasticity - From Theory to Biological Applications", book edited by Juan de Vicente , ISBN 978-953-51-0841-2, Published: November 7, 2012 under CC BY 3.0 license. © The Author(s).

Chapter 4

Die Swell of Complex Polymeric Systems

By Kejian Wang
DOI: 10.5772/50137

Article top


Plots of lnη versus lnφ1 and ln(η/φ1dln((η/φ1))/dγ˙) versus lnγ˙ from the experimental data for LDPE[77] of η~γ˙ andN1~γ˙.
Figure 1. Plots of lnη versus lnφ1 and ln(η/φ1dln((η/φ1))/dγ˙) versus lnγ˙ from the experimental data for LDPE[77] of η~γ˙ andN1~γ˙.
Predicted curves and experimental data of extrudate swell ratio for LDPE with different capillary L/D ratio and shear rates.
Figure 2. Predicted curves and experimental data of extrudate swell ratio for LDPE with different capillary L/D ratio and shear rates.
Calculated plots of B versus γ˙ from the experimental data for a semi-dilute polymer solution
Figure 3. Calculated plots of B versus γ˙ from the experimental data for a semi-dilute polymer solution
Curves of ln(t∞γ˙/Z) and Z versus γ˙ calculated from the experimental data.
Figure 4. Curves of ln(t∞γ˙/Z) and Z versus γ˙ calculated from the experimental data.
Plots of fitted swell ratio versus shear stress compared with experimental data.
Figure 5. Plots of fitted swell ratio versus shear stress compared with experimental data.
Experimental data for variation of viscosity with shear rate and the simulated values for the composites.
Figure 6. Experimental data for variation of viscosity with shear rate and the simulated values for the composites.
Comparison of predicted swell ratio with experimental data for PP/glass bead composites.
Figure 7. Comparison of predicted swell ratio with experimental data for PP/glass bead composites.
Comparison of plots ofB(γ˙,ϕ)/B(γ˙,ϕ=0) versus shear rate with experimental data.
Figure 8. Comparison of plots ofB(γ˙,ϕ)/B(γ˙,ϕ=0) versus shear rate with experimental data.
Extrudate swell ratio of PP/glass bead composites with different filler contents.
Figure 9. Extrudate swell ratio of PP/glass bead composites with different filler contents.
Plots ofB(τ,ϕ)/F(ϕ)versus shear stress from experimental data for PP/glass bead composites.
Figure 10. Plots ofB(τ,ϕ)/F(ϕ)versus shear stress from experimental data for PP/glass bead composites.

Die Swell of Complex Polymeric Systems

Kejian Wang1

1. Introduction

Die swell is a common phenomenon in polymer extrusion. When a viscoelastic fluid flows out of a die, the extrudate diameter is usually greater than the channel size. This is called die-swell, extrudate swell or the Barus effect. The degree of extrudate swell is usually expressed by the die-swell ratio (B) of extrudate diameter versus die diameter. A better understanding of such flow behavior will be beneficial for the optimization of processing parameters and the design of extrusion equipment, both of which affect product quality and production cost.

Innumerable valuable experimental data of die swell have been published, each focusing on different aspects affecting extrudate swell. There have been various interpretations of the extrusion swell of viscoelastic fluids from the macroscopic view of polymer rheology, such as a normal stress effect, an elastic energy effect, an entropy enlargement effect, an orientation effect and a memory effect. In fact, these interpretations are all related to each other[1,2]. It is generally believed that die swell is an important characteristic of the fluid elasticity during flow. The most common technique used to study rheological properties of polymer melts is capillary rheometry. In a capillary, polymer flows from the reservoir through a die and finally swells out of the exit. Under the action of extension, shear and compression, some molecular chains become disentangled, uncoiled or oriented in the convergent region resulting in an entry effect. During die flow, the resultant stress and strain cannot be relaxed completely. Simultaneously, some chains continue to be sheared and elongated during extrusion. On emerging from the die exit, the molecules are relaxed in elastic deformation by reentanglement and recoiling. The extrudate tends to contract in the flow direction and to grow in the normal direction, leading to extrudate swell[3]. As it does inside the capillary, the swelling evolves with time outside the capillary to reach a maximum at a certain extrusion distance. Graessley et al.[4] have reported that die swell occurs in two steps: (1) a very rapid swelling with a relatively large swell ratio very close to the die exit, which is known as running die swell; (2) a subsequent slow expansion to give an equilibrium die swell. Thus, the extent of the swelling will depend on both external factors as well as the intrinsic characteristics of the polymer[5]. The former include the geometry (contraction ratio and angle, die length and diameter as well as their ratio L/D) of the extrusion system (reservior, entry and die) [6,7] and the capillary operating conditions (the applied shear rate, stress and temperature) [8,9] and the external environment medium[10]. Extrudate swell can be used to assess the polymer viscoelasticity upon melt extrusion. Through rheological characterization, extrudate swell can also be correlated with the molecular structure of the polymer (molecular weight, and extents of crosslinking and long chain branching)[11].

Precise polymer processing requires an adequate quantitative description of flow. The viscous behavior determines the throughput, while the elastic properties are important for dimensional stability. Many theoretical studies have been conducted to improve our understanding and prediction of extrudate swell[12], in which a very critical aspect is the choice of an appropriate constitutive model[13,14]. In simulations, specific values of the model parameters—especially material parameters—must be defined. Some of these parameters can be measured by rheological experiments. Thus, quantitative analyses are able to relate extrudate swell to the viscoelasticity of polymer melts.

Extrudate swell was initially exploited to study the effects of die swell on processability and its correlation with the relevant rheological properties of polymer melts and the early work has been reviewed by Kontos and White[15] and Graessley[4]. In these early studies[16–18], a wide range of experiments on polymer melt swelling were performed. Some workers have reported that die swell varies with temperature[19], pressure during processing[20], molecular weight and its distribution[21–23], as well as molecular structure[24–27] and compounding ingredients[28]. Studies of die swell as a function of length to diameter ratio of die (L/D), entry speed and shape of the die have been reported by Han [29], Cotten[30], Lenk [31] and Vinogradov and Malkin[32]. Recently, the die swell of particle-filled polymers has attracted more attention because of their wider engineering applications[33,34]. These experiments have generated various valuable data of extrudate swell which need to be rationalized by theoretical studies.

Previous examinations of the existing data have identified some semi-empirical correlations relating the swell ratio to rheological parameters. In the 1970s, Bagley and Duffey[35], Graessley et al.[4], Han[36] and Tanner[37] proposed expressions for the relationship between the swell ratio B and the first normal stress difference or shear stress, on the supposition that polymer melt shear flow obeyed a simple law. One of the most famous approaches is that of Tanner[38] based on elastic recovery theory. In this model, the maximum diameter of the extrudate is related to the recoverable shear strain at the capillary wall.

Such continuum mechanics studies have been regularly reported from the early 1980s onwards. Investigations of the swelling phenomenon have been carried out using numerical simulations, especially by using the finite element method when considering complicated boundary conditions[39,40]. Most of the analyses have been conducted for laminar flow or at low Reynolds number, in which the swell ratio is greater than one. However, at higher Reynolds numbers, the dimensions of the extrudate are possibly smaller than those of the die. Few data have been published concerning the effect of Reynolds number on the die swell[41,42]. Han[29] suggested that the equilibrium die swell measurements are independent of rheometer geometry for common extrusion at relatively low flow rate through a capillary where the ratio of reservoir diameter to capillary die diameter is less than 10, and length to diameter ratio of the capillary die is greater than 20. Such a capillary is termed a long capillary, in which the effects of the entry flow from the reservoir on the die exit region flow can be neglected. Most theoretical models have been derived for a long capillary[43–47]. However, there have been few quantitative predictions of die swell in actual polymer processing operations in which the ratio of width to length in the die flow channel is relatively short[48-58]. The die swell behavior of polymeric materials should be very important in short capillary flow. However, entrance effects in the short tube flow of polymer melts are relatively complicated. Thus, it is necessary to extend the theoretical models of the die swell in a long capillary to that in a short capillary. This is one purpose of this paper.

On the other hand, particulate-filled polymer composites are widely used in engineering to improve polymer properties. Their flow properties are important when molding various products. Some swell experiments have been conducted on filled polymers in which it was noticed that the die swell was related to the filler shape, size, dispersion, concentration and interfacial modification of the inclusions in the matrix[59,60]. However, there are few quantitative theories of the swelling in composites, in contrast to the situation for pure polymer melts[37,61]. Thus, it is necessary to explore how the swell behaviors of composites differ from those of a pure polymer. A second aim of this paper is to develop a single unified quantitative extrudate swell model for both pure polymers and their composites.

The following sections first present a unified theoretical equation to describe extrusion swell in both long capillaries and short capillaries. The evolution of the finite distance of extrudate swelling and the particular features of the swelling of filled composites are also detailed. Finally, comparisons are made between the predictions of the model and experimental data.

2. Theoretical model

2.1. The extrudate swell theory for a long capillary

Up to now, most of the extrudate swell models have been established on the basis of the analysis of long capillary flow. The swell ratio B has been related to the recoverable shear strain (SR) or elastic strain energy[62]. One most pertinent systematic theories of extrudate swell of entangled polymeric liquids is that of Tanner, who first based his model on the K-BKZ constitutive equation and the free recovery from Poiseuille flow with the aspect ratio of length L to diameter D being very large[37]. The correlation was later confirmed for a wide class of constitutive equations, including PTT, pom–pom and general network type models for fully-developed tube flow[63]. It was found that the extrudate swelling ratio(B) may be expressed as a function of the first normal difference N1 and the shear stressτ12. Regardless of the fact that this clearly shows how the swell is related to the elasticity of viscoelastic polymeric fluids, the viscous heating and the time-dependent nature of swell are not considered.

More recently, Wang et al.[64,65] extensively exploited the variations of extrudate swell with both polymer characteristics (M¯n,M¯n/M¯W andM¯e) and the operational processing conditions. A double molecular mechanism of disentanglement-reentanglement and decoil-recoil was proposed to rationalize variations of polymer elasticity during flow with chain conformation. The die swell behavior essentially results from the molecular dynamics of the system. Thus, it is desirable to establish a single primary theoretical framework to relate the extrudate swell to the intrinsic viscoelasticity and external conditions.

There are three kinds of polymer segments and chains in the extrudates, which have been defined by Song[66] as extending chain, coil chain and entangled polymeric chain. Besides the change in chain conformation and the degree of reorganization of the constituent chains, the extension and flow can also induce the dynamic and reversible disentanglement and reentanglement between polymeric chains, such that the polymeric melts then undergo a partial stress relaxation leading to extrudate swelling at the die exit. The swell ratio is affected by the length to diameter ratio and the residence time. Based on such a viewpoint of dynamics, Song developed a novel molecular theory of such multiple entangled constituent chains in order to analyze non-linear viscoelasticity for the polymeric melts. His derived constitutive equations and material functions in a multiple-flow field were subsequently quantitatively verified[67]. Based on the O-W-F constitutive equation and the multiple transient-network model as well as the double relaxation dynamics of the reentanglement-disentanglement transition (RE-DT) and recoil-uncoil transition (RC-UCT) from the Poiseuille flow, Song proposed that swell evolves in three stages (instantaneous swelling, delayed swelling and ultimate extrudate swelling)[68]. A new set of swelling equations incorporating molecular parameters, operational parameters and growth time in both the steady state and dynamic state were developed. Song’s model successfully described the die swell through a long capillary of both linear polyethylene (HDPE) and linear polybutadiene (PBD) with the different molecular weights and different processing variables[69]. In this paper, Song’s extrudate swell theory will be extended in order to analyze both the die swell out of a short capillary and the swell of polymer composites.

In the steady shear flow, the shear viscosity can be written as


The above expression is verified by experimental data, and can also be deduced from the O-W-F constitutive equation together with molecular dynamics[69]. The coefficient of the first normal-stress difference Ψ1 in the steady shear flow is


From experimental data of η andψ1(γ˙), the molecular parameters of the zero shear viscosityη0, the equilibrium modulusG0, and the exponents n and a can be determined. The approximate value of the ultimate extrudate swell B was obtained as Eq. (3) by Song[69]:

Bdepends on both the molecular parameters and the operational variables [the capillary length to diameter ratioL/D and the shear rate γ˙ or shear stressτ]. (1W)is the fraction of the recoverable conformation of the entangled polymeric chain in the flow, which can be calculated from the experimental swell ratio data by Eq. (4).

Details of the derivation of the above equations (3) and (4) and their application to the extrudate swell of HDPE and PDB are given by Song[69]. Eq.3 may also be expressed in the form as shown in Eq. (5), where Bequ is the die swell ratio atL/D=, and mγis the shear rate-dependent constant.


2.2. Finite extrudate swell distance

As discussed above, the swelling evolves with time after the polymer is extruded from the capillary, during which time the molecules continue to exhibit similar disentanglement-reentanglement transitions and uncoil-recoil transitions to those occurring whilst it is still inside the capillary. To describe this phase, Eq. (3) can be reformulated as Eq. (6) by introducing two parameters kt andfw:


Here fw is a variable which replaces the fixed value of 0.5 used by Song[69] in Eq. (3) and the coefficient ktdescribes the evolution of the extrudate swell with time. ktmay be written as


As the swell approaches a maximum att=t=,kt approaches unity corresponding to the model used by Song[69]. However in contrast to Song’s model, in practice tis not infinite and actually has a finite value which can be determined experimentally[70], i.e., the maximum ultimate swelling will be realized at time t along the extrusion distance. If the extrusion distance is expressed asZ×D, where D is the capillary diameter and Z is a numerical factor, tand Zcan be correlated as follows from the shear rate,

where knis a constant for a given material, whose approximate value is

where n is the constitute exponent of the fluid. From Eqs. (7)(9), it can be shown that


Eq. (10) can be rewritten as Eq. (11).


Thus, either Eq. (10) or Eq. (11) can be used to predict the time of maximum swelling and the corresponding swell distance.

2.3. Extrudate swell theory for a short capillary

The above model is based on the assumption that the chain elongation incurred on reservoir entry is fully relaxed in the capillary. This is only approximately true for extrusion in a long capillary and a very poor assumption for a short capillary where the the entry flow is more complicated and the entry effect is more prominent[71]. Liang[72] reported the results of many swelling experiments using a short capillary and some quantitative empirical relations between the swell ratio and the material characteristics and operational parameters have been obtained[37,70,73]. In this paper, we attempt to modify Song’s model in order to describe the swelling behavior in short capillary extrusion.

Liang[74] described the die swell ratio as follows,

whereSR is the recoverable shear strain; λlshows the elastic strain induced by the stored energy in the capillary reservoir, which is related to the capillary geometry and the fluid viscoelasticity; λlis a function of the entry converging flow parameter; α0is the half-entry convergence angle of a viscoelastic fluid, which is given by a function of the constitute exponent n of the fluid, the entry pressure drop ΔP and the ratio of the capillary diameter D to its reservoir diameterDr. The pressure loss in the entry region can be approximately expressed in terms of the Bagley entrance correction factor (e).


The total die swell ratio in Eq. (3) can be approximated as Eq. (17), which has the same form as Eq. (12).


Eq. (17) can be modified to describe the extrudate swell out of a short capillary, as shown in Eq. (18),


where kt is the time-dependent coefficient reflecting the swelling evolution as defined by Eq. (7) and kLγ˙=[1+(τ0γ˙)a1+(L/D)a](1w)shows the effects of capillary length, i.e., the degree of relaxation. M/4is the recoverable effect from the stored energy in the capillary reservoir. It may be ca. 1 for a sufficiently long capillary while it may be λl for a very short capillary as defined by Eq. (14). Thus, Eq. (18) may be used for both long and short capillaries since it includes sufficient variables, whereas Eq. (3) is only appropriate for long capillaries.

2.4. Extrudate swell theory of filled composites

In particle-filled composites, the filler is distributed in an entangled matrix network and the filling affects the network relaxation[74]. Regardless of the variation in viscoelastic properties, it has been found experimentally that the values of ηand ψ1(γ˙) for filled composites are similar to those of the pure polymer matrix, i.e., the above extrudate swell theory may be modified for use with filled composites.

However, experiments have shown that the die swell is usually weakened with increasing filler concentration [75]. This suggests that the die swell of a composite can be expressed by Eq. (19),


where B(γ˙,ϕ) and B(γ˙,ϕ=0)=Bmare the die swell of the composite with a filler content of ϕ and the die swell of the pure matrix, respectively, andF(ϕ) is the filling effect. Eq. (19) implies that the viscoelastic behavior of the filled composites is dominated by the elasticity of the composite matrix in the high shear rate range. Liang[75] found thatB=λτw(1ϕς), whereτwwas the wall shear stress, andλ and ς were constants related to the elasticity of the matrix melt and the geometry of the filler particles, respectively. The function F(ϕ) is analogous to those used for the viscosity of a suspension of spheres, and may be called a ‘concentration shift factor’[76]. There are several forms forF(ϕ). Here we use the form shown in Eq. (20),


wherek,p and q are constants which depend on the filler-polymer system and the flow field; ϕcis the percentage decrease of the network elasticity caused by filling with the particles, which is directly related to the shear rate γ˙ as described by Eq. (21); γ˙cis the limiting shear rate when the network is almost completely destroyed, in accordance with percolation theory. An approximation to Eq. (21) may be rewritten in stress form as Eq. (22).

3. Comparison of die swell given by the model with experiment

3.1. Experimental data

For comparison of the predicted swell ratios using the model described above with experimental data, several sets of raw swell data were taken from the literature and replotted in appropriate formats.

In testing the validity of the model for a system with large L/D, the rheometric and swell data of a pure IUPAC-LDPE standard was used as shown in Figures 1 and 2 [77,78]. The capillary had die diameter D=3.00mm and experiments were carried out at 150℃.

To show whether the model is valid in analyzing the maximum extrusion distance during swelling, we used the data for a semi-dilute polymer solution of a partially hydrolyzed polyacrylamide manufactured by Rhone Poulenc (Rhodoflood AD37, Mw= 6 x 106, degree of hydrolysis = 24%)[70]. Purified water containing 20 g/l of NaCl was employed as a solvent for the polymer. The polymer concentration was 3000 ppm. The swell tests were conducted in a stainless steel capillary with length of 51 mm and inner diameter of 0.60 mm.

The swell data for a rubber compound[73] were used to verify the effectiveness of the suggested equations for polymer extrusion in a short capillary. Sample SI was a calcium carbonate (CaCO3) filled natural rubber (NR) compound, in which the content of CaCO3 was 20 phr. Sample SII was a carbon-black (CB)-filled NR/styrene-butadiene rubber SBR) /cis-1,4-butadiene rubber (CBR) compound. The blending ratio of NR/SBR/CBR was 45/10/45, and the content of CB was 56 phr. Both the rubber compounds included some additives. An Instron capillary rheometer (model 3211) was used in the tests. Two capillary dies with different lengths were selected in order to measure the rheological properties of the sample materials. The length of the long die was 40 mm, the length of the short die was 0.2 mm, and the diameter of both dies was 1 mm.

Another set of data was for a polypropylene/glass bead composite[75]. The matrix resin was a common polypropylene (Pro-fax 6331). The melt flow rate (2.16 kg, 230 oC) and density were 12 g/10 min and 0.9 g/cm3, respectively. The glass bead filler had a mean diameter of 219 nm and density of 2.5 g/cm3. Rheological experiments with these samples were carried out using a Rosand capillary rheometer with twin cylinders. Two dies of different length with the same diameter (1 mm) were selected in order to measure the entry pressure losses. The length/diameter of the short die was << 1 and the other was 16. The test temperature was 190 oC, and apparent shear rates were varied from 50 to10–4s–1.


Figure 1.

Plots of lnη versus lnφ1 and ln(η/φ1dln((η/φ1))/dγ˙) versus lnγ˙ from the experimental data for LDPE[77] of η~γ˙ andN1~γ˙.

3.2. Comparison of the predicted values of swell ratio with the experimental data

3.2.1. Swelling simulations for pure LDPE extruded from capillaries with different L/D ratios

The dashed curves in Figure 2 present the experimental data[77] of extrudate swell ratio for LDPE for different capillary L/D ratio and shear rates from which (1W) can be calculated assuming kt=1 using Eq. (4), as listed in Table 1. fwis determined by fitting B against L/Dunder different shear rates.

In addition, the swell ratio was also fitted by B=B0kB(L/D) as suggested by Liang[75] (as shown by the dotted curves in Figure 2) and by Eq. (5) (as shown by the solid curves in Figure 2). The resulting values of key parameters for the three cases are compared in Table 1.

Dashed curves
By Eq. (3) (kt=1)
Dotted curves
Solid curves

Table 1.

Values of parameters for LDPE obtained from the plots in Figures 1 and 2


Figure 2.

Predicted curves and experimental data of extrudate swell ratio for LDPE with different capillary L/D ratio and shear rates.

The comparisons of swell ratio predicted by Eq. (3), B=B0kB(L/D)and Eq. (5) are also compared with experimental data (points) in Figure 2. All three models predict that the swell ratio decreases with increasing capillary L/D ratio and increases with increasing shear rate. It is also found that all three equations describe well the approximate linear relationship between the experimental swell ratio and L/D at larger values of L/D, i.e., in long capillaries. However, below a certain critical value of L/D, B increases rapidly with decreasing L/D. This trend is well predicted by Eq. 5 only. Liang’s model[77] B=B0kB(L/D)gives a completely different picture to experiment. Song’s model Eq. (3) underestimates the swell ratio compared with experiment, showing that Eqs. (3) and (4) are only appropriate for swelling on extrusion from a long capillary, because of the assumption in Song’s theory[69] that the entry effect is almost relaxed in capillary flow. In a short capillary, the chain elongation on entry cannot be neglected. The critical L/D ratio, which determines whether a capillary behaves as ‘long’ or ‘short’, increases with shear rate. This may be explained by the more incomplete relaxation of the elongated chains in the entry region at higher shear rate. The degree of relaxation depends on both the capillary length and the strength of the flow field as described by Eq. (5).

3.2.2. Extrusion distance giving the maximum swell ratio for semi-dilute polymer solutions

For a capillary with a L/D ratio of 85, the swelling of polyacrylamide solution was observed to reach a maximum at a given time tat a given distance (expressed as a multiple Zof capillary diameterD) from the capillary[70], i.e., the swelling is finite. In the simulations, n, aand τ0 were obtained by fitting the experimental η~γ˙ curve[70] using Eqs. (1) and (2). Ktand (1W)can be fitted using Eq. (3). tand Z were subsequently calculated from Eq. (11) and compared with the experimental data as shown in Figure 3. The calculated swell ratio B is close to the experimental data for Kt=3.5468 and(1W)=0.5964. The simulated value of tdecreases as a power law function of shear rate as observed from the experimental data in Figure 4a, as does Z as shown in Figure 4b.


Figure 3.

Calculated plots of B versus γ˙ from the experimental data for a semi-dilute polymer solution


Figure 4.

Curves of ln(tγ˙/Z) and Z versus γ˙ calculated from the experimental data.

3.2.3. Verification of the swelling equation for rubber in a short capillary

The recoverable strain can be calculated by Eq. (23) due to Tanner[37].


Eq. (24) can be approximately fitted to give values ofMkt/4, aand 1Wwhenfw=0.5.


The experimental swell ratios[73] are shown in Figure 5, from which the values of the parameters in Eq. (24) can be obtained by fitting and they are listed in Table 2. Figure 5 also illustrates that the swell ratios predicted by Eq. (24) are in good agreement with the experimental data. It should be noted that the swell ratios reported in the literature[73] were measured when the extrudate was naturally cooled to ambient temperature. The swell ratio is considerably larger for extrusion from a very short capillary than from a long capillary. A larger fraction of the relaxation occurs in the melt state at higher temperature in the longer capillary than in the short capillary. The unrelaxed energy is not easily released from the capillary at low temperature. The fitted value of (1W)is larger for swelling on extrusion from the short capillary than from the long capillary, which reflects a greater amount of retained elasticity in the former case.


Figure 5.

Plots of fitted swell ratio versus shear stress compared with experimental data.

The filler content of sample SI with a matrix of NR rubber is lower than that of sample SII with a matrix of NR/SBR/CBR blend[73]. As seen from Table 2, the average relaxation time τ is shorter, and the values of a and (1W)are lower for sample SI than for sample SII containing a higher filler content.


Table 2.

Values of the parameters in Eq. (24) for rubber compounds

3.2.4. Application of the swelling equations for a PP/glass bead composite

As shown in Figure 6, the viscosity of the PP/glass bead composite can be calculated from the experimental data of stress versus shear rate in the literature[75]. Eq. (1) accurately predicts the variation of viscosity with shear rate for composites with different filler percentages. The fitted values ofη0,τ0 , aand n are listed in Table 3; all of these values increase with increasing filler content in the composite while shear thinning becomes more marked when increasing the filler content.


Figure 6.

Experimental data for variation of viscosity with shear rate and the simulated values for the composites.

Filler contentnτ0/saη0/kPa.s

Table 3.

Structural parameters for PP/glass bead composites


Figure 7.

Comparison of predicted swell ratio with experimental data for PP/glass bead composites.

Eq. (6) can be used to predict the extrudate swell ratio of the composites with the fitted structural parameters listed in Table 3. Comparison of the predicted data with the experimental values are shown in Figure 7, assumingkt=1. The modified theory shows better agreement at higher shear rate which indicates its validity for particulate-filled composites. Table 4 shows that(1W), which is the residual fraction of the recoverable entanglement in polymeric chains after flow shear, decreases with increasing filler content. Smaller values of(1W) correspond to lower elastic recovery, i.e., swell ratios are lower for composites with higher filler content. Furthermore, it can be seen that fw decreases (from fw=0.51370.5 for the pure matrix melt) to lower values with increasing filler concentration. This shows that recovery is harder for an entangled network with higher filler content than for the corresponding matrix network.

Dashed curves given by Eq. (3)Dotted curves given by B=(γ˙/γ˙cb)mb
Filler content1Wfwγ˙cbmB

Table 4.

Values of model parameters in the equations for the composite in Figure 7

In addition, it is interesting that the swell is well depicted by one simple equationB=(γ˙/γ˙cb)mb. γ˙cbmay be related to the composite rigidity or the density while mB is sensitive to the shear, reflecting the recovery capability of the composites. When increasing the filler content in the composite, γ˙cbbecomes larger while mB becomes smaller, showing that the elasticity is weakened while the rigidity becomes stronger, which corresponds to the observed decrease in swell ratio.


Figure 8.

Comparison of plots ofB(γ˙,ϕ)/B(γ˙,ϕ=0) versus shear rate with experimental data.

Plots ofB(γ˙,ϕ)/B(γ˙,ϕ=0) against shear rate based on the experimental data in Figures 6 and 7 are shown in Figure 8. The concentration shift factor F(ϕ)=[1ϕ/ϕc] (here k and p in Eq. (20) are both unity) illustrates the reinforcement effect while reflecting the weaker swelling capability. In our opinion, there are two levels of network structure. The microscopic one involves molecular entanglement, whilst the other involves mesoscopic particles dispersed in the compound melt. The concentration shift factor shows that the swelling functions of the two levels may be separated. The mesoscopic network probably shrinks, counteracting part of the microscopic melt swelling. The offset effect can be described by the concentration shift factor. However, the filler is far larger than the polymer chains and does not vary the basic molecular relaxation dynamics of the melt matrix except at the interfacial region around the particles. The dispersed particles perturb the flow of the melt. It is interesting that this is consistent with previous publications showing that the dimensions of the extrudate were possibly smaller than those of the die at higher Reynolds numbers[41,42].

In addition, Eq. (21) is also successful in fitting the plots ofB(γ˙,ϕ)/B(γ˙,ϕ=0) against γ˙ as Figure 8. γ˙cincreases while q decreases when adding more glass beads to the polypropylene. To a certain degree, γ˙crepresents a critical shear rate when the melt swelling is completely offset by the shrinkage of the mesoscopic network.

    Plots of die swell versus filler concentration[75] are shown in Figure 9. Eq. (20) well demonstrates this correlation. In fitting the data with Eq. (20), it is found that with increasing shear rate the critical content ϕc decreases, while the die swell ratio at ϕ=0 increases. Larger die swell ratios are observed for pure PP at higher shear rates, since more elastic energy is stored.


    Figure 9.

    Extrudate swell ratio of PP/glass bead composites with different filler contents.

    B(τ,ϕ=0)=B(τ,ϕ)/F(ϕ)is also appropriate for correlating shear stress and filler content for the composite as shown in Figure 10. In this case, in the expression ϕc=(τ/τc)q'τc increases while q' decreases with increasing filler content. More importantly, B(τ,ϕ)/F(ϕ), i.e., B(τ,ϕ=0)is almost a linear function of shear stress. It can be written as B(τ,ϕ)/F(ϕ)=(B/F)τ=0+kττ where Bτ=0 and kτ decrease with increasing amount of filler in the PP melt.

    Figure 10.

    Plots ofB(τ,ϕ)/F(ϕ)versus shear stress from experimental data for PP/glass bead composites.

    4. Conclusions

    Two limitations of Song’s polymer extrudate swell theory have been identified for the first time. Song’s model has been modified in order to predict the finite distance at which the swelling reaches a maximum. Furthermore, the model was extended to describe the die swell on extrusion from a short capillary by considering the entry effect in Song’s molecular dynamics model and incorporating Liang’s expression. The resulting modified model can be applied to extrusion swelling for both long and short capillaries, whereas Song’s model is only appropriate for long capillaries. More importantly, the modified model is also suitable for analysis of the swelling of particle-filled composites which cannot be treated by Song’s model. The composite swell ratio can be separated into the product of the matrix swell ratio and the concentration shift factor. The excellent agreements between the values predicted by the modified model and experimental data reported in the literature for a variety of different systems demonstrate its viability for a wide range of materials and experimental conditions.


    1 - Liang J.Z.Plast Rubber Compos Process Appl 1991157579
    2 - J. Z. J. Liang, Process. Mater, Technol, 199552207215
    3 - Lodge A.S. Elastic Liquids, Academic Press, NewYork, 1964
    4 - W. W. Graessley, S. D. Glasscock, R. L. Crawley, Soc. Trans, Rheol, 197014519544
    5 - A. C. Y. J. Wong, Proc. Mater, Tech. 199879163169
    6 - Koopmans R.J. J. Appl. Polym. Sci. 1992; 32:1755-1762.
    7 - J. Z. Liang, Rubber. Plast, Appl. Proc, 1995239396
    8 - Seo Y. J. Appl. Polym. Sci. 1990;30:235-240.
    9 - Richardson S. Rheol. Acta 1969;9 :193-199.
    10 - R. J. J. Koopmans, Polym. Appl, Sci. 19923217501767
    11 - R. J. J. Koopmans, Polym. Appl, Sci. 19923217411750
    12 - P. L. Swan, J. M. Dealy, Eng. Polym, Sci. 199131705710
    13 - Larson R.G., Constitutive Equations for Polymer Melts and Solutions, Butterworths, Boston, 1988.
    14 - H. J. Giesekus, Fluid. Non-Newtonian, Mech. 19821169109
    15 - White J.L. Rubber Chem. Technol.197750163185
    16 - R. S. Spencer, R. E. J. Dillon, Sci. Colloid, 19483163171
    17 - Schreiber H.P., Rudin A., Bagley E.B. J. Appl. Polym. Sci. 1965; 9:887-892.
    18 - A. B. Metzner, Soc. Trans, Rheol. 196913467470
    19 - Brydson J.A.Flow Properties of Polymer melts. London: The Plastic Institute. 1981
    20 - F. N. Cogswell, Polym. Plast, 1973413943
    21 - M. T. Shaw, Eng. Polym, Sci, 197717266268
    22 - Rogers M.G. J. Appl. Polym. Sci. 1970;14 :1679-1689.
    23 - Racin R., D.C. Bogue, J. Rheol. 1979;23:263-280.
    24 - Mendelson R.A., F.L. Finger, J. Appl. Polym. Sci. 1975;19:1061-1078.
    25 - Rokudai M. J. Appl. Polym. Sci. 1981;26:1427-1429.
    26 - Liang J. Z., Nes J. N. s,Polymer Testing 1999;18:37-46,
    27 - L. Dufrancatel-Veiller, M. F. Lacrampe, J. J. Pababiot, Polym. Appl, Sci. 20018017101724
    28 - N. Minagawa, J. L. J. White, Poly. Appl, Sci, 197620501523
    29 - C. D. Han, Rheology in Polymer Processing; Academic: New York, 1976
    30 - G. R. Cotten, Chem. Rubber, Technol, 197952187198
    31 - Lenk R.S. Polymer Rheology. London: Applied Science Publishers Ltd.Ch.10, 1011978
    32 - Vinogradov G.V. and Malkin A.Y. Rheology of Polymers: Viscoelasticity and Flow of Polymers, Moscow: Mir Publishers.Ch. 5, 3553791980
    33 - D. Rapeephun, Y. Jimmy, S. Pitt, Polymer, Testing 200524211
    34 - M. Nattaya, N. Manit, P. Brian, Rathanawan. M. Gradyb, Polymer, Testing 200827470479
    35 - E. B. Bagley, H. J. Duffey, Soc. Trans, Rheol, 197014454459
    36 - C. D. Han, M. Charles, Soc. Trans, Rheol. 197014213218
    37 - R. I. J. Tanner, Sci. Polym, A-. Part, Phys. Polym, 1970820672071
    38 - R. I. J. Tanner, Fluid. Non-Newtonian, Mech. 19806289302
    39 - Crochet M.J., Keunings R. J. Non-Newtonian Fluid Mech. 1980; 7:199-212.
    40 - G. Barakos, E. J. Mitsoulis, Rheol. 199539193209
    41 - S. Middleman, Fundamental of Polymer Processing. New York: Magraw-Hill. 1977
    42 - Georgiou G.C., Papanastasiou T.C. and Wilkes, J.O. AIChE. J. 1988; 34:1559-1562.
    43 - L. A. Utracki, Z. Bakerdjian, M. R. J. Kamal, Polym. Appl, Sci, 197519481501
    44 - J. Meissner, Appl. Pure, Chem, 197542551612
    45 - C. W. Macosko, Principles. Rheology, Measurements, New. Application, Wiley. York-Vch, . Ch, pp, 2362601990
    46 - Treloar L.R.G.The Physics of Rubber Elasticity. 3rd Edn., Oxford: Clarendon Press. 1975
    47 - J. R. Barone, N. Plucktaveesak, S. Q. J. Wang, Rheology199842813832
    48 - D. Huang, J. L. White, Polym Eng Sci 198020182186
    49 - Doelder. C. F. J.den, R. J. J. Koopmans, Fluid. Non-Newt, Mec, 2008
    50 - N. Muksing, M. Nithitanakul, B. P. Grady, R. Magaraphan, Testing. Polym, 2008274470479
    51 - K. Kumari, P. P. Kundu, of. Bulletin, Mater, Sci. 2008312159167
    52 - S. Ayutthaya, N. Isarankura, J. Wootthikanokkhan, J. Appl, Polym, Sci. 2008107638533863
    53 - S. Mitra, S. Chattopadhyay, A. K. Bhowmick, J. Appl, Sci. Polym, 2008107527552767
    54 - R. M. Rowell, J. Polym, Environ, 2007154229235
    55 - G. S. de Paulo, M. F. Tome, S. J. Mc Kee, Fluid. Non-Newtonian, Mech, 20071473149174
    56 - E. Mitsoulis, J. Fluids, Trans.of. Eng, A. S. M. the, ASME20071291113841393
    57 - S. Mohanty, S. K. Nayak, Eng. Polym, Sci, 2007471016341642
    58 - B. Peng, H. Wu, S. Y. Guo, S. Y. Lai, J. Jow, J. Appl, Polym, Polym.Sci. 2007106317251732
    59 - J. Z. Liang, R. K. Y. Li, C. Y. Tang, S. W. J. Cheung, Polym. Appl, Sci, 200076419424
    60 - D. L. Faulkner, L. R. Schmidt, Eng. Polym, Sci. 197717657667
    61 - M. B. Bush, Eng. Polym, Sci, 199333950958
    62 - K. K. Kamal, U. O. J. Joshua, Elast, Plas. 200133297336
    63 - R. I. J. Tanner, Fluid. Non-Newtonian, Mech. 20051298587
    64 - Tapadia P., Wang S.Q. Macromolecules 2004;37:9083-9099.
    65 - Zhu Z Y, Wang S Q, JRheol. 2004483571578
    66 - M. S. Song, G. X. Hu, Z. H. Yang, Q. Xu, S. Z. J. Wu, Sci. Mater, Tech, 2006
    67 - M. S. Song, Q. Xu, G. X. Hu, S. Z. J. Wu, Sci. Mater, Tech. 2006225664686
    68 - C. W. Zhu, M. S. Song, G. X. Hu, J. Zhao, D. M. Wu, J. Chin, Phys. Chem, 2007205563581
    69 - J. Zhao, M. S. Song, C. W. Zhu, G. X. Hu, K. J. Wang, D. M. Wu, J. Chin, Phys. Chem, 20082115588
    70 - C. Allain, M. Cloitre, P. J. Perrot, Fluid. Non-Newtonian, Mech, 1997735166
    71 - R. Zheng, G. J. Tang, Polymer, communication 19863161172
    72 - J. Z. J. Liang, Process. Mater, Technol, 199659268276
    73 - J. Z. Liang, J. Appl, Sci. Polym, 20071047074
    74 - Stabik J., Inter Polym Process2004194350355
    75 - Liang J.Z. Polymer Testing200221927931
    76 - W. Gleissle, B. J. Hochstein, Rheol. 200347897910
    77 - S. X. Huang, C. J. J. Lu, Fluid. Non-Newtonian, Mech, 2006136147156
    78 - J. Meissner, Appl. Pure, Chem, 197542551612