Measurement and modeling of uniaxial and planar extensional viscosities for linear isotactic polypropylenes

In this work, novel rectangular and circular oriﬁce (zero-length) dies were used to measure planar and uniaxial extensional viscosities as a function of strain rate for various linear isotactic polypropylene melts by using Cogswell methodology. The obtained experimental data were combined with shear and uniaxial extensional viscosity data determined at very high strain rates. The ability of the molecularized generalized Newtonian ﬂuid (mGNF) [M. Zatloukal and J. Drabek, “Generalized Newtonian ﬂuid constitutive equation for polymer liquids considering chain stretch and monomeric friction reduction for very fast ﬂows modeling,” Phys. Fluids 33 (8), 083106 (2021)], Giesekus, and explicit Yao constitutive equations to describe the measured data was tested. It has been shown that including the effect of the chemical environment (i


I. INTRODUCTION
Planar extension flow is when the polymer melt stretches in one direction, shrinks in the thickness direction, and has no shrinkage in the width direction.6][7][8] While uniaxial extensional viscosity is a material property that can be relatively easily measured using a variety of experimental techniques, [9][10][11][12] determining planar extensional viscosity is one of the most challenging rheological tasks because of the difficulty of generating and controlling planar extensional flow.It is therefore not surprising that the planar extensional viscosity is not correctly handled when modeling polymer processing, and the flow phenomena associated with this material property are not yet fully understood.0][21][22] Second, we assess the ability of the recently proposed molecularized generalized Newtonian fluid (mGNF), 23,24 Giesekus, [25][26][27] and recently proposed explicit Yao 28 constitutive equations to describe the shear and both extensional viscosities with a single set of parameters.Finally, we generalize Newton's law for the oligomeric solvent contribution to the stress tensor (i.e., to include the role of chemical environment) in the used viscoelastic constitutive equations to improve their behavior at very high strain rates.

II. EXPERIMENTAL
Three linear isotactic polypropylenes (iPPs) Borflow HL504FB (76 k), HL508FB (64 k), and HL512FB (56 k) were used in this work.The molecular characteristics of the linear isotactic polypropylene melts used are summarized in Table I.Polypropylenes were characterized at 230 C using a Rosand RH7-2 twin-bore capillary rheometer equipped with polyether ether ketone (PEEK) piston tips.Two circular dies with the same diameter of D d ¼ 0.5 mm but different length-todiameter ratios, L/D, (0.12 14,29  and 16), and also two rectangular dies of width, w ¼ 10 mm, gap size h ¼ 0.15 mm and different L/h ratios (0.12 14,29 and 16) were used to determine the uniaxial and planar extensional viscosities using the Cogswell model 14,17,18 (see Table II).In this table, Q is the volumetric flow rate; R is the radius of the capillary die; w and h are the width and gap size of the rectangular die, respectively; L is the length of the long die; P L,U and P L,P are the pressure drop through the long die, which is circular and rectangular in shape, respectively; and P 0,U and P 0,P represent the entrance pressure drop measured on the circular and rectangular orifice die, respectively.0][21][22] The maximum Hencky strain, e max , achieved when measuring uniaxial and planar extensional viscosity during the abrupt contraction flow is given as where A b and A d are the cross-sectional area of the barrel and the die, respectively. 30The maximum Hencky strain for circular (e max;U ) and rectangular (e max;P ) orifice dies is then given as where D b is the barrel diameter (15 mm), D d is the orifice circular die diameter (0.5 mm), and w and h are the width (10 mm) and gap size (0.15 mm) of the rectangular orifice die.According to Eqs. ( 2) and ( 3), e max;U ¼ 6:8 and e max;P ¼ 4:8.Note that e max;U is the same for uniaxial extensional viscosity measured at low strain rates (this work) and high strain rates (experimental data taken from our previous work [19][20][21][22] ).It should also be mentioned that for the same volume flow rate and circular and rectangular die dimensions used, the Hencky strain, extensional strain rate (including wall shear rate and wall shear stress used in the Cogswell analysis) is smaller for planar extensional flow compared to uniaxial extensional flow.While maintaining the same volume flow rate in circular and rectangular geometries, the pressure measured on rectangular dies is much less compared to circle dies.Thus, given the limited resolution of pressure transducers, experimental data from rectangular dies are typically quantifiable at much higher volume flow rates compared to circular geometries.

III. THEORETICAL
The following constitutive equations were chosen to describe the measured rheological data.

A. Molecularized generalized Newtonian fluid model (mGNF)
A recently proposed molecularized generalized Newtonian fluid model was used in this work 23 where s stands for the extra stress tensor, D represents the strain rate tensor, L is the objective velocity gradient (defined as where L is velocity gradient and X is the tensor that gives the rate of rotation of the eigenvectors of D [31][32][33][34][35] ) and gðII D ; II L Þ means the viscosity, which can vary with the second as well as on the second invariant of the objective velocity gradient II L ¼ 2trð L 2 Þ according to Eq. ( 5) (note that L and the velocity gradient tensor L are same in the steady-state flows 34,35 ) Here, gðII D Þ is given by the well-known Carreau-Yasuda model, Eq. ( 6), and f ðII L Þ varies between 1 (II L ! 0) and 2 (II L ! 1) and is given by Eq. ( 7), Here, the parameters g 0 (zero-shear-rate viscosity), g 1 (infinite-shearrate viscosity), k 1 (reptation-mode relaxation time), a (parameter characterizing transition from Newtonian to power-law regime), and n (power-law index) are related to the shear viscosity curve while b (extensional strain hardening parameter), k 2 (Rouse mode relaxation time), and A (parameter related to the asymptotic value for extensional viscosity, i.e., for _ e ! 1) are obtained by fitting the steady-state extensional viscosity.
In pure shear flow and and thus, f ðII L Þ defined by Eq. ( 7) becomes equal to 1 because II L ¼ 0; that is, the shear viscosity becomes dependent on the second invariant of the strain rate tensor only with II D ¼ _ c 2 as follows: In general steady-state extensional flow For simple uniaxial extensional flow m ¼ À0.5, for a planar extensional flow m ¼ 0, and for equibiaxial extensional flow, m ¼ 1. 9 A combination of Eqs. ( 4)-( 7), (11), and (12) leads to the following expression for the extensional viscosity: where and The asymptotic formula for extensional viscosity is given as The constant A can be expressed from Eq. ( 16) considering uniaxial extensional flow, that is g E ðm ¼ À0:5Þ, as follows: In this model, the expression for the uniaxial extensional viscosity at the high strain rate limit g E;U;1 corresponds to the molecular expression for a fully extended Fraenkel chain 36 and is given by the following formula: 23 where M e is the molar mass between entanglements, k max is the maximum chain stretch ratio, M c is the critical molar mass at which entanglements begin to occur, M is the molar mass, x represents the value by which the zero-shear rate viscosity, g 0 , scales with M x above M c (typically 3:5 6 0:2 for all linear and flexible molecules 37 ), and f eq f aligned represents the ratio of monomeric friction coefficients for equilibrium and fully aligned chains.Alternatively, the parameter A can also be determined from the known molecular parameters of the polymer melt by combining Eqs. ( 17) and ( 18) as follows: 23

B. Modified Giesekus model
6][27] The model is given as follows: where s is the extra-stress tensor, s p and s s are the polymer and solvent contribution to the stress tensor, g 1 is the solvent viscosity, g p is the polymer viscosity, g 0 is the zero-shear-rate viscosity, D is the strain rate tensor, k is the relaxation time, the symbol r represents the upper-convected time derivative, G is the modulus, and a is the parameter characterizing anisotropic hydrodynamic drag.The minimum and maximum anisotropies correspond to a ¼ 0 and a ¼ 1, respectively. 38The steady-state shear viscosity is given by the following equations: 39 The steady-state extensional viscosity is given by the following equation: It has been showed [40][41][42][43] that increasing the molecular weight of the oligomeric solvent leads to co-aligned orientation of short solvent molecules with long polymer chains, which reduces the monomeric friction coefficient and the extensional viscosity.In order to include this effect in viscoelastic models, let us assume that the solvent contribution to the stress tensor is given by the molecularized generalized Newtonian fluid model introduced in Sec.III A (with gðII D Þ being replaced by g 1 and A replaced by A G ) as follows: In this case, the expression for shear viscosity remains unchanged (because f ðII L Þ¼1), but the expression for extensional viscosity changes from Eq. ( 27) to the following expression: which gives the following asymptotic formula: The constant A G can be expressed from Eq. ( 30) for the case of uniaxial extensional flow, that is, g E ðm ¼ À0:5Þ, as follows: Since A G is defined as a positive value, g E;U;1 > 2 a ðg 0 À g 1 Þ.

C. Modified explicit Yao model
Yao proposed a constitutive equation in which the extra stress tensor is an explicit function of the objective velocity gradient with finite stretch and rotational recovery. 28The model is based on a modified stress law for a general elastic material that includes finite stretch and disentanglement Here, s is the extra stress tensor, G is the elastic modulus, G 0 is the initial elastic modulus (G 0 ¼ g 0 k0 ), n Àw represents the elastic modulus shape function controlled by the strain hardening coefficient w (0 w 1) 35 taking into account the effect of finite chain stretch on the elastic modulus (i.e., that G increases with increasing chain stretch), S 0 is the ceiling stretch (typically in the range of 1-3), 35 k 0 is the relaxation time at small strain, D is the strain rate tensor, B Ã is the generalized Finger tensor representing the accumulated elastic strain in the polymer melt considering both the effect of rotational relaxation controlled by index of rotational recovery n 0 (1 n 0 2) 35 and the effect of the finite stretch on the relaxation time by using the shape function n a0 controlled by the parameter a 0 with a typical value around 1 28 (i.e., a decrease in k 0 with increased chain stretch), and L is the objective velocity gradient.The explicit nature of the model represents the main advantage over previous implicit Yao models, 34,35 from which the current model is based.The explicit Yao model can be further generalized by including the solvent contribution via a molecularized Newtonian model (i.e., analogous to the modification describe above for the modified Giesekus model) to include the role of the chemical environment reported in Refs.40-43, which leads to the following expression for the extra stress tensor: where n is given by Eq. ( 34).The steady-state shear viscosity is given by the following equation: where n ¼ S0

S0þ
k P ffi ffi The steady-state extensional viscosity is given by the following equation: where The asymptotic formula for extensional viscosity is given as The constant A Y can be expressed from Eq. ( 40) considering uniaxial extensional flow, that is, g E ðm ¼ À0:5Þ, as follows: As can be seen, the expression for the Yao model constant A Y is the same as for the mGNF model constant A [compare Eqs. ( 17) and ( 41)].Note that the original Yao's expressions for both viscosities are recovered if g 1 ¼ 0 (i.e., g P ¼ g 0 and k P ¼ k 0 ) in Eqs. ( 36)- (38).It should also be mentioned that the explicit Yao model is viscoelastic; that is, it can predict the first and second normal stress differences, 28 while the mGNF model cannot.

IV. RESULTS AND DISCUSSION
Measured data for all three iPP samples combined with veryhigh-strain-rate-dependent shear and uniaxial extensional viscosities taken from our previous work [19][20][21][22] are shown in Figs.1-3 together with the mGNF model (Fig. 1), the original and modified explicit Yao model (Fig. 2), and original and modified Giesekus model (Fig. 3).The corresponding model parameters are given in Tables III-VII.
The mGNF model has a high ability to describe all the measured rheological data shown in Fig. 1.Here, g 0 , g 1 , k 1 , a, and n ¼ 0, g E;U;1 parameters were taken from our previous work 19,[21][22][23] and k 2 and b were determined from the fit of uniaxial and planar extensional viscosities (see Table III).As expected, the reptation-mode relaxation time k 1 and Rouse mode relaxation time k 2 increase with the weightaveraged molecular weight, M w , while the extensional strain hardening parameter b decreases with increasing M w.

23,44
A comparison between the single-mode Giesekus model predictions and the measured data is shown in Fig. 2 (left), and the corresponding model parameters are provided in Table IV.Also in this case, the parameters g 0 and g 1 were taken from our previous work 19,[21][22][23] (i.e., the polymer viscosity g p is known in advance because g p ¼ g 0 À g 1 ).The remaining model parameters k and a were determined by fitting the measured strain-rate dependent shear, uniaxial, and planar extensional viscosities.As can be seen in Fig. 2 (left), the model captures the shape of shear viscosity vs shear rate, even when single-mode approach was used.6][47][48] Weak extensional thinning in uniaxial extensional viscosity at moderate strain rates can only be achieved if a is sufficiently high.At very high extensional strain rates, the model predicts a constant extensional viscosity given as The asymptotic formulas for the uniaxial and planar extensional viscosities for the case of the most extensional thinning (i.e., when a ¼ 1) are thus given as follows: As can be seen, none of the above extensional viscosities ever drop below 2g 0 þ g 1 .A closer analysis of Eq. (42) shows that the Giesekus model provides extensional thinning for uniaxial and planar extensional viscosity only when 2=3 < a 1 and 1=2 < a 1, respectively.The dependence of both extensional viscosities on the extensional strain rate for different values of a is shown in Fig. 4 using the parameters for the 76 k polypropylene melt.This explains why the Giesekus model is unable to describe the measured extensional strain thinning followed by extensional strain thickening in both extensional viscosities for the given iPP melts shown in Fig. 2 (left).On the other hand, the Giesekus's reptation-mode relaxation time k and the anisotropy parameter a (causing a decrease in g E;1 as a increases) correctly increase with M w as expected.
The modified Giesekus model shows improved behavior in both extensional viscosities at very high extensional strain rates, as seen in Fig. 2 (right).The model describes extensional thinning for both extensional viscosities at moderate strain rates in the same way as the original Giesekus model (by setting 2=3 < a 1 with a minimum achievable uniaxial and planar extensional viscosity of 2g 0 þ g 1 and 2g 0 þ 2g 1 , respectively), but it can also handle extensional strain hardening at very high strain rates (i.e., at values higher than the reciprocal value of the Rouse time) due to a modified expression for the solvent contribution to the stress tensor.The ability of the model to describe the extensional viscosities increases with the reduced M w of the tested iPPs, because in that case the minimum attainable extensional viscosities predicted by the original Giesekus model approach the experimentally determined minimum extensional viscosities.It is clear that the proposed modification of the "solvent" viscosity (to account for the effect of the chemical environment) can significantly improve the ability of the viscoelastic constitutive equations to describe the extensional rheology at very high strain rates.It should be mentioned that the behavior of the Giesekus model at moderate extensional strain rates can be further improved by modifying the expression for the polymer contribution to the stress tensor by incorporating the finite extensibility of the polymer molecules into dumbbell kinetic theory using the Peterlin approximation as shown in Refs.45-48.The trend between the modified Giesekus's reptation-mode relaxation time k and M w remains the same as for the original Giesekus and mGNF models.Since k 2 is of the same order as in the mGNF model, but its value increases (and does not decrease) with decreasing M w , the k 2 is also correlated with Rouse time in this case, but inversely.The effect of M w on the strain-hardening-related parameters a and b is not clear in this case, as the ability of the modified Giesekus model to fit the extensional viscosity data varies significantly with M w (the ability of the model to fit experimental data increases with reduced M w , see Fig. 2 (right) and Table V for the fitting errors evaluated using the root mean square error, RMSE), which does not allow a direct physical interpretation of both parameters.
A comparison between the experimental data and the original explicit Yao model predictions is shown in Fig. 3 (left).As can be seen, the model is unable to describe the secondary Newtonian plateau in the shear viscosity curve, fails to describe the uniaxial viscosity data, and unrealistically predicts that the uniaxial and planar extensional viscosities cross each other at moderate strain rates and become infinite at very high strain rates.The corresponding model parameters (six in total) are listed in Table VI, where g 0 was taken from Refs.19 and 21-23 and the ceiling stretch S 0 was adjusted to 1 to ensure that model predicts extensional thinning at moderate strain rates 35 (which is typical flow behavior for linear polymer melts such as the tested iPPs).The Yao relaxation time is found to increase correctly with M w , but its value is about one order of magnitude higher than the reptation-mode relaxation time, that is, on the order of milliseconds, which is interestingly the same order as the longest relaxation time for the given iPP.For example, the longest relaxation time for iPP with M w ¼ 76 kg/mol was found to be 5.36 ms at 230 C 20 and the Yao's relaxation time for the same material and temperature is 1.91 ms (see Table VI).The index of rotational recovery n 0 was found to be equal Here, g S is the shear viscosity, g E,U is the uniaxial extensional viscosity, g E,P is the planar extensional viscosity, _ c is the shear rate, and _ e is the extensional strain rate.to 1, independent of M w , which is a typical value when no second normal stress difference, N 2 , is available during the fitting procedure (in which case the model predicts N 2 ¼ 0 28 ).Yao pointed out that n 0 > 1 to correctly predict N 2 .There is no simple relationship between the strain hardening coefficient w and M w (see Table VI), which can be explained by the loss of physical meaning of w due to the poor ability of the model to represent extensional viscosities at a very wide strain rate range.
The modified Yao model removes the shortcomings of the original Yao model, as can be deduced from Fig. 3 (right) and Table VII; that is, it correctly predicts both the occurrence of the secondary Newtonian plateau in the shear viscosity curve and the uniaxial and planar extensional viscosities over the entire strain rate range.Here, g 0 , g 1 , and g E;U;1 parameters were taken from our previous work 19,[21][22][23] (i.e., polymer viscosity g p is known in advance because g p ¼ g 0 À g 1 ), S 0 and n 0 were kept the same, equal to 1 as in the case of the original Yao model and remaining five parameters (k P , a 0 , w, k 2 , and b) were determined from the fit of shear, uniaxial, and planar extensional viscosities (see Table VII).As can be seen, the conclusions formulated for the Yao relaxation time k P remain unchanged for the modified version of the model and both strain hardening parameters w and b have a direct physical meaning as they both increase with reduced M w as expected. 44k 2 is of the same order as in the case of the mGNF model, but its value increases with decreasing M w ; that is, k 2 is also inversely related to the Rouse time as in the case of the mGiesekus model.
It should be mentioned that the molecularized generalized Newtonian fluid (mGNF) constitutive equation as well as the explicit Yao model can only be used to model steady-state flows.On the other hand, the proposed modification (i.e., taking into account the role of the oligomeric solvent by replacing the commonly used Newton's law with a simplified mGNF equation in the fully viscoelastic constitutive equations) demonstrated here on the Giesekus model enables the modeling of transient flow types, since the extra-stress tensor in both the original and the modified Giesekus model is given as the sum of the time-dependent polymer stress and the time-independent stress contribution from the oligomeric solvent.
A second note is related to the high strain rate viscosity data (taken from Refs.19 and 49), which are visualized in Figs.1-3 together with the low strain rate data measured in this work.In our previous work, we demonstrated through finite element method (FEM) simulations that for the given iPP melts, die design used and experimental conditions, the effect of viscous dissipation and pressure canceled each other out (for more details, see the section "effect of viscous dissipation and pressure" in our previous paper 19 ).We have shown that, first, the secondary Newtonian viscosity, g 1 (occurring above shear rates of 2 Â 10 6 s À1 ), depends linearly on the weight-averaged molecular weight, M w , for both entangled linear and branched polypropylene melts, indicating that polymer chains became fully disentangled at the    19,49 Second, the free volume increases at the secondary Newtonian plateau due to complete chain disentanglement (most likely due to the coalescence of shaped free volume cavities), as indicated by a large change in the flow activation energies of PPs (E 0 > E 1 ). 19,49Moreover, it was experimentally demonstrated by Takahashi et al. 50that there is no wall slip at very high shear rates for PP at 230 C (i.e., for the same polymer and flow conditions utilized in our work), which can also be explained by the complete chain disentanglements and significant increase in the free volume (i.e., no cohesive slip appears to occur and the adsorption force of the iPP chains to the wall surface appears to be high enough to suppress adhesive slip at very high strain rates).

V. CONCLUSIONS
In this work, the Cogswell methodology was used to determine strain-rate-dependent planar and uniaxial extensional viscosities for various polypropylene melts using novel rectangular and circular orifice (zero-length) dies.The obtained experimental data were combined with shear and uniaxial extensional viscosity data determined at very high strain rates.The ability of the recently proposed mGNF constitutive equation to describe the measured data was shown to be very    high.The Giesekus and recently introduced explicit Yao viscoelastic constitutive equations were modified to include the effect of chemical environment (i.e., the role of the oligomeric solvent influencing the monomeric coefficient and extensional viscosities) using a simplified version of the mGNF constitutive equation for the solvent contribution to the stress tensor.The proposed generalization of Newton's law for an oligomeric solvent has been shown to significantly improve the ability of both constitutive equations to describe measured experimental data, especially at very high extensional strain rates using parameters with a clear physical meaning.

FIG. 1 .
FIG. 1.Comparison between measured strain rate-dependent shear (open symbols), uniaxial (red full triangles and squares) and planar extensional viscosities (blue full circles), and mGNF model fits (curves) at 230 C for three linear isotactic polypropylenes [76 k (a), 64 k (b), and 56 k (c)].Here, g S is the shear viscosity, g E,U is the uniaxial extensional viscosity, g E,P is the planar extensional viscosity, _ c is the shear rate, and _ e is the extensional strain rate.

FIG. 2 .
FIG. 2. Comparison between measured strain rate-dependent shear (open symbols), uniaxial (red full triangles and squares) and planar extensional viscosities (blue full circles), and model fits (curves) at 230 C for three linear isotactic polypropylenes [76 k (a) and (d), 64 k (b) and (e), and 56 k (c) and (f)].The original and modified singlemode Giesekus model predictions are shown in [(a)-(c)] and [(d)-(f)], respectively.Here, g S is the shear viscosity, g E,U is the uniaxial extensional viscosity, g E,P is the planar extensional viscosity, _ c is the shear rate, and _ e is the extensional strain rate.

FIG. 3 .
FIG. 3. Comparison between measured strain rate-dependent shear (open symbols), uniaxial (red full triangles and squares) and planar extensional viscosities (blue full circles), and model fits (curves) at 230 C for three linear isotactic polypropylenes [76 k (a) and (d), 64 k (b) and (e), and 56 k (c) and (f)].The original and modified explicit Yao model predictions are provided in [(a)-(c)] and [(d)-(f)], respectively.Here, g S is the shear viscosity, g E,U is the uniaxial extensional viscosity, g E,P is the planar extensional viscosity, _ c is the shear rate, and _ e is the extensional strain rate.

FIG. 4 .
FIG. 4. Effect of the anisotropic parameter a on the normalized uniaxial (a) and planar (b) extensional viscosities predicted by the Giesekus model for 76 k iPP melt at 230 C.

TABLE I .
Basic characteristics of the isotactic polypropylenes used.

TABLE III .
Parameters of the mGNF model.

TABLE IV .
Parameters of the single-mode original Giesekus model.

TABLE V .
Parameters of the single-mode modified Giesekus model.

TABLE VI .
Parameters of the original explicit Yao model.

TABLE VII .
Parameters of the modified explicit Yao model.