Estimation of the Shear Viscosity of Mixed-Polymer Materials for Screw Extrusion-Based Recycling Process Modeling

The scope of this work is the development of a method to estimate the temperature and shear rate-dependent viscosity of mixtures composed of two polymers. The viscosity curve of polymer mixtures is crucial for the modeling and optimization of extrusion-based recycling, which is the most efficient way to recycle polymeric materials. The modeling and simulation of screw extruders requires detailed knowledge of the properties of the processed material, such as the thermodynamic properties, the density, and the rheological behavior. These properties are widely known for pure materials; however, the incorporation of impurities, like other polymers in recycled materials, alters the properties. In this work, miscible, immiscible, and compatibilized immiscible polymer mixtures are considered. A new method based on shear stress is proposed and compared to the shear rate-based method. Several mixing rules are evaluated for their accuracy in predicting mixture viscosity. The developed methods allow the prediction of the viscosity of a compatibilized immiscible mixture with deviations below 5% and that of miscible polymer mixtures with deviations below 3.5%.


Introduction
The recycling of polymer films is not only ethically required but also enforced by the European Union in its Directive 94/62/EC, which was first formulated in 1994 and revised in 2018.It states that by end of December 2025, 50 wt% of packaging plastics must be recycled.This is raised to 55% by the end of 2030 [1].Plastic films are produced by either flat-film or blown-film extrusion processes.In these processes, single-screw extruders are commonly used.From a hopper polymer, bulk material is fed to the solid conveying zone of the single-screw extruder.Afterwards, the bulk material is compacted, compressed, plasticated, conveyed, mixed, and extruded through a die to receive the characteristic film shape [2][3][4][5][6].
The packaging of food and other perishable goods demands certain barrier properties from the packaging materials [7].Mostly, oxygen (O 2 ) and water vapor (H 2 O) are unwanted in the products, whereas carbon dioxide (CO 2 ) is often used as an inert gas inside the packaging.The diffusion of gases through porous solids can be described with Knudsen diffusion [8][9][10].The physical resistance that opposes the diffusivity is the barrier property [11].The permeation process can be explained by a combination of Henry's and Fick's laws [12].Since different polymers exhibit different barrier properties, a combination of different layers of polymers can suffice to reach the imposed limits of permeability.
These combinations are called multilayer films, which are usually created by using multiple extruders that are connected to a feedblock that arranges the layers and feeds them to a die.One of the major challenges that comes at the end of life of multilayer film packaging is the recycling process.Films can be dumped in landfills, energy recovery can be employed, or they can be recycled [13].The most common choice of use at the end of Polymers 2024, 16, 1339 2 of 23 life in Europe is still energy recovery or thermal recycling [14].Other routes are chemical recycling, feedstock recycling, or mechanical recycling.The enforcement of directives and regulations will tip the balance towards mechanical recycling, in which end-of-life products are collected, sorted, and reprocessed to new products.Contrary to chemical recycling [15], mechanical recycling does not allow the separation of different polymers of multilayer films.Those polymers are processed in the same material stream.The mixing of polymers has effects on most material properties, especially the rheological, mechanical, and thermal properties.Compatibilization of immiscible polymers can assist in hindering the deterioration of these properties [16,17].
Kallel et al. [18] discussed the effect of compatibilization of polyethylene-polystyrene (PE-PS) and polyethylene-polypropylene (PE-PP) blends.By adding maleic acid anhydridegrafted ethylene-propylene copolymer (EP-g-MAH) for PE-PP and styrene-ethylenebutadiene-styrene (SEBS) for PE-PS, respectively.They analyzed rheological behavior with a parallel-plate rheometer and concluded that at low shear rates, compatibilizers act as emulsifiers in the molten state.This effect is diminished at high shear rates, where curves tend to converge.This effect was also previously reported by Palierne [19] and Bousmina [20].Yoo et al. [21] showed the effect of maleic acid anhydride-grafted styrene-acrylonitrile (SAN-g-MAH) on the rheological properties of polypropylene carbonate-polymethyl methacrylate (PPC-PMMA) blends and reported a strong enhancement effect of the compatibilization, even at parallel-plate rheometry testing frequencies of 10 2 rad•s −1 , corresponding to a maximum shear rate of approximately 400 s −1 .Kozlowski and Mantia [22] reported a quite different effect when studying blends of PP and LCP with and without the addition of PPgMAH (maleic acid anhydride-grafted polypropylene).The shear viscosities derived from capillary rheometry showed a decrease with the addition of the compatibilizer.Holsti-Miettinen et al. [23] reviewed the effect of the chemical structure of compatibilizers on the properties of compatibilized PP-polyamide 6 (PA 6) blends.They used EBA-g-FA (fumaric acid-grafted ethylene-butylene acrylate), PPgMAH, SEBS-g-MAH and E-AE-GMA (ethylene-ethyl acrylate-glycidyl methacrylate terpolymer with glycidyl methacrylate) and tested the blends using a capillary rheometer.In a shear rate region between 10 2 and 10 4 s −1 , all the compatibilizer additions led to an increase in viscosity of the blend.SEBS-g-MAH had the strongest effect in this study, revealing even higher viscosities than the single components at low shear rates.
The calculation of extrusion characteristics demands rheological material data [3][4][5][6].By altering this behavior through the incorporation of impurities like other polymers, the process design must be adapted [6,[24][25][26].Since it would be far too time-consuming to test every material mixture before extrusion, material models for mixtures have been developed and are now more important than ever before.Although several mixing rules exist, there is no calculation method that considers the shear rate/shear stress-dependent behavior or polymers, until now.
Bingham [27] (1914) and Heitmiller et al. [28] (1964) formulated an approach with a reciprocal additive relation to calculate the shear viscosity of a mixture of two polymers considering Newtonian fluids with the same density.An expression for the mixture's viscosity η M B/H is derived by balancing the flow rates and formulated as: where v A is the volume fraction of component A, and η A and η B are the shear viscosities of component A and B, respectively.Heitmiller et al. [28] assume that the liquids are divided into many layers that are flowing adjacently without mixing.Lin [29] (1979) presented a comparable approach with an additional interaction parameter that considers the interfacial slip between the layers [30][31][32][33] or instabilities [34][35][36] that may influence the results.
Polymers 2024, 16, 1339 3 of 23 Tsenoglou [37] (1987) estimated the stress relaxation modulus of entangled polymer blends and formulated his mixing rule for a mixture's zero-shear viscosity η 0,M T,Orig : where v B is the volume fraction of component B. As we deal with mixtures of two materials, v B is 1 minus v A .η 0A and η 0B are the zero-shear viscosities of component A and B, respectively.The mixing rule of Tsenoglou was adapted as described in [38] to analyze if it allows to estimate the shear rate-dependent viscosity η M T of a mixture of polymers in a wide range of shear rates: Arrhenius [39] (1887) and Bersted et al. [40] (1981) proposed a mixing rule η M A/B for blends of branched and linear components, which is then assumed to be independent of the shear rate.The viscosity of the mixture is given as: where w A and w B are the weight fractions of components A and B. The model is validated with viscosity data of a mixture of high-and low-density polyethylene (HDPE and LDPE) from another study [41].Furthermore, the model is used to describe the viscosity of neat HDPE as a mixture as well.
Friedman and Porter [41] published a mixing rule that originates in the work of Kendall and Monroe [42].Their mixing rule for predicting the viscosity of a mixture η M KM/FP : uses a scaling exponent c between the viscosity and the molecular weight, which is usually around 3.4 [43,44].The simplest mixing model is a linear model η M lin : Walther [45] presented his double logarithmic mixing rule for predicting the viscosity of a mixture η M W in 1931.According to Gao and Li [46], the constant C is set to 0: Cragoe [47] presented his reciprocal logarithmic mixing rule for predicting the viscosity of a mixture η M Cr in 1933.According to Gao and Li [46], the constant L is set to 2000: When one of these mixing rules is used to predict the viscosity of polymer mixtures, the shear rate/shear stress-dependent viscosity of the pure materials must be considered.From the mentioned mixing rules, the rule of Bersted et al. [40] is the only one that was applied considering the shear rate.Contrary, Utracki and Wilkie [48] recommend estimating the viscosity of a mixture based on shear stress, but they do not show a method of how to do this in detail.In theory, a shear stress-based evaluation method should allow an estimation of the rheological behavior of blends with higher accuracy.
This paper presents a method to predict the shear rate/shear stress-dependent viscosity of a mixture of two polymers from the viscosity of the pure polymers and from the volume fraction.A new method to calculate the viscosity of a polymer mixture based on shear stress considering several mixing rules is proposed.The results are compared to measured values, which are also compared to a shear rate-dependent calculation method that is also considering the previously mentioned mixing rules.Accurately predicting the mixture viscosity can drastically speed up the determination of the material behavior for the simulation of polymer recycling processes.Moreover, the effect of compatibilization on the viscosity of PP-PA12 mixtures is analyzed.

Materials and Methods
The polymeric materials used are PP HB600TF and PP HD234CF (both Borealis AG, Vienna, Austria) and polyamide 12 (PA12) Grilamid L25 (EMS-CHEMIE HOLDING AG, Herrliberg, Switzerland).Both PP grades are typical film extrusion materials.Maleic acid anhydride-grafted PP (PPgMAH) Orevac CA 100 (Arkema Functional Polyolefines, Colombes Cedex, France) was used as a coupling agent between the PP HB600TF and the PA12 grade.The polymeric materials with the corresponding melt flow rate (MFR) or melt volume rate (MVR) and the recommended processing temperatures are given in Table 1.The polymer mixtures were produced using a co-rotating twin-screw extruder: Leistritz ZSE 27 Maxx (Leistritz Group, Nuremberg, Germany).The temperature profile used was held constant for all material blends (see Figure 1) and a standard screw configuration was used.
Polymers 2024, 16, x FOR PEER REVIEW 4 of 25 the volume fraction.A new method to calculate the viscosity of a polymer mixture based on shear stress considering several mixing rules is proposed.The results are compared to measured values, which are also compared to a shear rate-dependent calculation method that is also considering the previously mentioned mixing rules.Accurately predicting the mixture viscosity can drastically speed up the determination of the material behavior for the simulation of polymer recycling processes.Moreover, the effect of compatibilization on the viscosity of PP-PA12 mixtures is analyzed.

Materials and Methods
The polymeric materials used are PP HB600TF and PP HD234CF (both Borealis AG, Vienna, Austria) and polyamide 12 (PA12) Grilamid L25 (EMS-CHEMIE HOLDING AG, Herrliberg, Switzerland).Both PP grades are typical film extrusion materials.Maleic acid anhydride-grafted PP (PPgMAH) Orevac CA 100 (Arkema Functional Polyolefines, Colombes Cedex, France) was used as a coupling agent between the PP HB600TF and the PA12 grade.The polymeric materials with the corresponding melt flow rate (MFR) or melt volume rate (MVR) and the recommended processing temperatures are given in Table 1.
Table 1.Polymeric materials used, including their abbreviations for this work, the MFR/MVR, and the processing temperature, as given in the material data sheets.

Grade
Abbreviation MFR/MVR The polymer mixtures were produced using a co-rotating twin-screw extruder: Leistritz ZSE 27 Maxx (Leistritz Group, Nuremberg, Germany).The temperature profile used was held constant for all material blends (see Figure 1) and a standard screw configuration was used.The contaminating polymer was supplied by the side-feeder (Congrav OP5 dosing, Kubota Brabender Technologie GmbH, Duisburg, Germany) in each setup, whereas the base polymer was provided by the main dosing (Congrav M, Kubota Brabender Technologie GmbH, Duisburg, Germany).As a coupling agent, the PPgMAH masterbatch was added to the PP1 in a preceding compounding process.
The amount of coupling agent that was needed to obtain a satisfying adhesion between PP and PA12 was determined in a previous step.For this purpose, the PP base polymer PP1 was blended with 1, 2, 3, 4 and 5 wt% of the PPgMAH masterbatch.Afterwards, 20 wt% PA12 was added to those PP compounds with different amounts of coupling agent, resulting in the mixtures that are shown in Table 2.The contaminating polymer was supplied by the side-feeder (Congrav OP5 dosing, Kubota Brabender Technologie GmbH, Duisburg, Germany) in each setup, whereas the base polymer was provided by the main dosing (Congrav M, Kubota Brabender Technologie GmbH, Duisburg, Germany).As a coupling agent, the PPgMAH masterbatch was added to the PP 1 in a preceding compounding process.
The amount of coupling agent that was needed to obtain a satisfying adhesion between PP and PA12 was determined in a previous step.For this purpose, the PP base polymer PP 1 was blended with 1, 2, 3, 4 and 5 wt% of the PPgMAH masterbatch.Afterwards, 20 wt% PA12 was added to those PP compounds with different amounts of coupling agent, resulting in the mixtures that are shown in Table 2.The test specimens for the Charpy impact test in Table 2 were produced as prescribed by DIN EN ISO 20753:2019 [49] using an injection molding machine Engel Victory 330/80 Tech (Engel Austria GmbH, Schwertberg, Austria).The process settings were a melt temperature of 230 • C, a mold temperature of 40 • C, pressure holding time of 40 s and a cooling time of 13 s.The middle part of the multipurpose test specimen was separated with a custom-made cutter to be used for Charpy impact testing, following the DIN EN ISO 179 norm [50].The samples were conditioned at a temperature of 23 • C at 50% relative humidity for a duration of 88 h in a CLS climatic chamber 70/600 (Klima Temperatur Systeme GmbH, Jennersdorf, Austria) before conducting the impact tests using a Zwick 5113E pendulum impact tester (ZwickRoell GmbH & Co. KG, Ulm, Germany) equipped with a 50 J pendulum.This mechanical test was repeated five times for each blend.After deciding on which amount of coupling agent to use for the main set of rheological experiments, those blends were manufactured as well.Table 3 shows all the different blends used in the rheological experiments.To determine the morphology of the blends, field electron scanning electron microscopy (FE-SEM) was deployed.Pellets of the blends were submerged in liquid nitrogen, placed in a cotton cloth, and broken with a hammer.The resulting pieces were placed on a sample holder and sputtered with a thin layer of gold before analyzing them with a MIRA3 LHM FE-SEM (Tescan Orsay Holding, Brno, Czech Republic).
The rheological tests were performed using a high-pressure capillary rheometer Rheograph 6000 (Göttfert, Buchen, Germany) with the corresponding Labrheo software (version 4.3.3).The tests were conducted at a temperature of 230 and 250 • C deploying capillaries with a diameter of 1 mm and an L/D ratio of 5, 10, 20 and 30, respectively.
Depending on the L/D ratio of the capillary dies in use, one pressure transducer was used to detect the chamber pressure before the inlet of the capillary.The measured pressure loss was corrected using Bagley correction [51].The true shear viscosity curve was obtained employing Weissenberg-Rabinowitsch correction [52].The resulting true viscosity curves were fitted using a simplified Bird-CarreauQ -Yasuda (BCY) [53] model (Equation ( 9)) in combination with the least-square method and a numerical solver: where η 0 is the zero-shear viscosity, B is a time constant describing the reciprocal value of the shear rate where the viscosity changes from the Newtonian zero-shear rate plateau to shear thinning behavior, n is the power-law exponent and a is a material specific model parameter that describes the transition from the Newtonian zero-shear-rate plateau to shear thinning behavior in more detail.Bersted et al. [40] presented a shear rate-based method of applying their mixing rule, but they did not describe the method how to do this in detail.We propose performing the shear rate-based method in the following manner.The first step is to fit a viscosity model (e.g., that of BCY [53]) for the pure materials.Four additional evaluation steps are needed to predict the viscosity of a mixture as a function of the shear rate assuming that the viscosity models of the pure materials are known: 1.
Define several shear rate values for the calculation.

2.
Calculate the viscosity values of the pure polymers at the defined shear rate values.

4.
If needed, apply a viscosity model on the calculated viscosity of the mixture.

Shear Stress-Based Calculation of the Mixture Viscosity
Utracki and Wilkie [48] recommend the calculation of the viscosity of polymer mixtures based on shear stress, because higher accuracy is expected.The shear stress-based calculation method is more complicated than the shear rate-dependent calculation, as numerical calculations are needed.The following evaluation steps must be carried out after fitting a viscosity model (e.g., that of BCY [53]) for the pure materials.

1.
In the first step, several shear stress values for the calculation are defined.In this paper, the shear stress values of the mixtures from the rheological tests were used.

2.
The viscosity values of the pure polymers are calculated for those defined shear stresses.For this purpose, the basic relation between shear stress and shear rate is used.The shear rate value .
γ must be adapted until the calculated shear stress value equals the shear stress value that was defined in step 1.This can be done with a numerical solver, as presented in [6,54,55] 3.
After that, the mixing rules presented in the Introduction (linear mixing rule, Kendall and Monroe [42]/Friedman and Porter [41], Tsenoglou [37], Arrhenius [39]-Bersted et al. [40], Bingham [27]-Heitmiller et al. [28]) are employed to calculate the shear viscosity of the mixture η m as a function of the shear stress τ.The comparison of the results using the mixing models to the measured values is much more complicated for the shear stress-dependent calculation, as the shear rate must be calculated from the shear stress and the shear viscosity and deviates from the measured values.Thus, the next step is needed.γ M,Model can now be calculated from the viscosity of the mixture η M,Model and from the shear stress τ: The results of this calculation are several points that describe the viscosity of a mixture and the corresponding shear rate values for the different mixing rules.

5.
Then, a viscosity model is used to approximate the calculated viscosity values of the mixture η M,Model considering the shear stress values obtained in step one and the shear rate value .
γ M,Model calculated in step four.In this paper, the simplified BCY (see Equation ( 9)) was used.6.
To compare the results of the shear stress-dependent calculation, the viscosity models from step five must now be applied to calculate the viscosity of the mixture for the different models at the shear rates from the rheological tests.
The deviations of the mixture viscosity that are presented in this paper are calculated as the difference between the predicted model values and the values received from the rheological tests.A positive deviation shows an overestimation by the mixing rules and a negative deviation and underestimation of the predicted mixture viscosity.

Comparison of the Mixing Models
The different mixing rules from the literature are first compared against each other.A test evaluation for two materials with the viscosities 2239 and 9636 Pa s was performed to compare the predicted mixture viscosity values obtained with the different mixing models.As Figure 2 shows, the linear mixing rule predicts the highest mixture viscosity.The mixing rules of Kendall and Monroe [42] and Friedman and Porter [41] deliver viscosity values that are slightly below the values calculated with the linear model.The values calculated with Tsenoglous's [37] mixing rule are lower than those predicted according to Kendall and Monroe [42]/Friedman and Porter [41], and Arrhenius [39] and Bersted et al. [40].The viscosity values obtained employing the mixing rules of Cragoe [47] and Walther [45] show nearly the same behavior as those obtained according to Arrhenius [39]-Bersted et al. [40].Consequently, the mixing rule of Cragoe and Walther was not further considered.When applying the reciprocal additive relations of Bingham [27] and Heitmiller et al. [28], the lowest viscosity values were obtained.
The results of this calculation are several points that describe the viscosity of a mixture and the corresponding shear rate values for the different mixing rules.5.Then, a viscosity model is used to approximate the calculated viscosity values of the mixture , considering the shear stress values obtained in step one and the shear rate value , calculated in step four.In this paper, the simplified BCY (see Equation ( 9)) was used.6.To compare the results of the shear stress-dependent calculation, the viscosity models from step five must now be applied to calculate the viscosity of the mixture for the different models at the shear rates from the rheological tests.
The deviations of the mixture viscosity that are presented in this paper are calculated as the difference between the predicted model values and the values received from the rheological tests.A positive deviation shows an overestimation by the mixing rules and a negative deviation and underestimation of the predicted mixture viscosity.

Comparison of the Mixing Models
The different mixing rules from the literature are first compared against each other.A test evaluation for two materials with the viscosities 2239 and 9636 Pa s was performed to compare the predicted mixture viscosity values obtained with the different mixing models.As Figure 2 shows, the linear mixing rule predicts the highest mixture viscosity.The mixing rules of Kendall and Monroe [42] and Friedman and Porter [41] deliver viscosity values that are slightly below the values calculated with the linear model.The values calculated with Tsenoglous's [37] mixing rule are lower than those predicted according to Kendall and Monroe [42]/Friedman and Porter [41], and Arrhenius [39] and Bersted et al. [40].The viscosity values obtained employing the mixing rules of Cragoe [47] and Walther [45] show nearly the same behavior as those obtained according to Arrhenius [39]-Bersted et al. [40].Consequently, the mixing rule of Cragoe and Walther was not further considered.When applying the reciprocal additive relations of Bingham [27] and Heitmiller et al. [28], the lowest viscosity values were obtained.

Shear Viscosity of PP Contaminated with a Different PP Grade
An overview of the shear viscosity curves for the different mixtures of PP 1 and PP 2 is given in Figure 3.The lines represent the BCY fits of the measured values, and the measured values are shown as symbols.Nearly all measured values fall between the viscosities of the pure materials.The corresponding parameters of the fits can be seen in Table A1 in Appendix A. The measurement data can be found in Appendix A in Tables A2 and A3.All presented viscosity and shear rate values are so-called true values.

Shear Viscosity of PP Contaminated with a Different PP Grade
An overview of the shear viscosity curves for the different mixtures of PP1 and PP2 is given in Figure 3.The lines represent the BCY fits of the measured values, and the measured values are shown as symbols.Nearly all measured values fall between the viscosities of the pure materials.The corresponding parameters of the fits can be seen in Table A1 in Appendix A. The measurement data can be found in Appendix A in Tables A2  and A3.All presented viscosity and shear rate values are so-called true values.Figure 4 shows the relative deviation between the predicted and the measured shear viscosity of polypropylene PP1 mixed with different concentrations of PP2.The results of the shear rate-based calculation method are depicted on the left-hand side in charts a-e.The deviations obtained by employing the shear stress-based calculation method are presented on the right-hand side in charts f-j.For the prediction of mixture viscosity, the linear mixing rule, the rule of Bingham and Heitmiller, the rule of Arrhenius and Bersted, the rule of Kendall and Monroe/Friedman and Porter, and the rule of Tsenoglou are used.The mixtures with 10 wt% PP2 and with 90 wt% PP2 are described very well with all mixing models and both calculation methods.In the case of the other mixtures, the shear rate-based calculation method delivers negative deviations for nearly all measurement points.In this case, the linear model predicts the mixture viscosity best, which is demonstrated by the lowest deviations from the measured values.The shear stressdependent calculation is more accurate for the mixtures with 30, 50, and 70 wt% PP2, which can be seen in the lower deviations from the measured values.
It seems like low concentrations of another polymer (10%) lead to rather strong deviations that can be described well with the reciprocal additivity relationship of Bingham and Heitmiller et al. [28].Mixtures with a higher concentration of impurities lead to comparatively higher viscosities and can be described by the linear mixing rule in the best manner.This holds true for both the shear rate-based and the shear stress-based calculation methods.For the shear rate-based calculation, the measured values are even higher than the linear mixing rule predicts.Figure 4 shows the relative deviation between the predicted and the measured shear viscosity of polypropylene PP1 mixed with different concentrations of PP2.The results of the shear rate-based calculation method are depicted on the left-hand side in charts a-e.The deviations obtained by employing the shear stress-based calculation method are presented on the right-hand side in charts f-j.For the prediction of mixture viscosity, the linear mixing rule, the rule of Bingham and Heitmiller, the rule of Arrhenius and Bersted, the rule of Kendall and Monroe/Friedman and Porter, and the rule of Tsenoglou are used.The mixtures with 10 wt% PP 2 and with 90 wt% PP 2 are described very well with all mixing models and both calculation methods.In the case of the other mixtures, the shear rate-based calculation method delivers negative deviations for nearly all measurement points.In this case, the linear model predicts the mixture viscosity best, which is demonstrated by the lowest deviations from the measured values.The shear stress-dependent calculation is more accurate for the mixtures with 30, 50, and 70 wt% PP 2 , which can be seen in the lower deviations from the measured values.
It seems like low concentrations of another polymer (10%) lead to rather strong deviations that can be described well with the reciprocal additivity relationship of Bingham and Heitmiller et al. [28].Mixtures with a higher concentration of impurities lead to comparatively higher viscosities and can be described by the linear mixing rule in the best manner.This holds true for both the shear rate-based and the shear stress-based calculation methods.For the shear rate-based calculation, the measured values are even higher than the linear mixing rule predicts.Relative deviation in % Shear rate in s-1  Figure 5 shows the results of the impact tests of the PP compounds with different PPgMAH concentrations.PP 1 has approximately half the impact strength of PA12 (the values are presented as lines).Compounding PP 1 with 20 wt% PA12 without compatibilization leads to a decrease in the impact strength.Adding 1, 2 or 3 wt% PPgMAH to PP 1 before blending it with PA12, which results in 0.8, 1.6 or 2.4 wt% PPgMAH in the final mixture, has only minor effects on the impact strength.At a concentration of 3.2 wt% (4 wt% PPgMAH in PP 1 ), a sudden increase occurs, and the impact strength jumps to those of neat PA12.At a PPgMAH content of 4 wt% (5 wt% PPgMAH in PP 1 ), no further enhancement of the impact strength can be found, which is why a PPgMAH content of 4% was chosen for the rheological tests.Due to the enhanced impact strength, it is assumed that sufficient compatibilization between PP 1 and PA12 is obtained at this concentration of coupling agent.

Shear Viscosity of PP Contaminated with PA12
Figure 5 shows the results of the impact tests of the PP compounds with different PPgMAH concentrations.PP1 has approximately half the impact strength of PA12 (the values are presented as lines).Compounding PP1 with 20 wt% PA12 without compatibilization leads to a decrease in the impact strength.Adding 1, 2 or 3 wt% PPgMAH to PP1 before blending it with PA12, which results in 0.8, 1.6 or 2.4 wt% PPgMAH in the final mixture, has only minor effects on the impact strength.At a concentration of 3.2 wt% (4 wt% PPgMAH in PP1), a sudden increase occurs, and the impact strength jumps to those of neat PA12.At a PPgMAH content of 4 wt% (5 wt% PPgMAH in PP1), no further enhancement of the impact strength can be found, which is why a PPgMAH content of 4% was chosen for the rheological tests.Due to the enhanced impact strength, it is assumed that sufficient compatibilization between PP1 and PA12 is obtained at this concentration of coupling agent.Figure 6 shows SEM-images with a magnification of 5000 of cryogenically broken pellets of blends of PP1-PA12 with and without coupling agent and a different PA12 concentration.Especially at lower PA12 concentration, the effect of the compatibilizer well visible.Except for a concentration of 20 wt% PA12, the blends without coupling agent appear as one large phase with small, bubble-shaped inclusions.Those inclusions become larger and more numerous with higher PA12 concentrations.In the SEM images of the samples with compatibilization, the inclusions appear much smaller and better dispersed.At a PA12 concentration of 20 wt%, the phases are not easily distinguishable anymore.Figure 6 shows SEM-images with a magnification of 5000 of cryogenically broken pellets of blends of PP 1 -PA12 with and without coupling agent and a different PA12 concentration.Especially at lower PA12 concentration, the effect of the compatibilizer well visible.Except for a concentration of 20 wt% PA12, the blends without coupling agent appear as one large phase with small, bubble-shaped inclusions.Those inclusions become larger and more numerous with higher PA12 concentrations.In the SEM images of the samples with compatibilization, the inclusions appear much smaller and better dispersed.At a PA12 concentration of 20 wt%, the phases are not easily distinguishable anymore.Figures 7 and 8 show the measured shear viscosity curves of the uncompatibilized PP1-PA12 blends, the compatibilized PP1-PA12-PPgMAH blends and the pure materials at temperatures of 230 °C and 250 °C.Especially at higher shear rates, PA12 exhibits a far higher shear viscosity than PP1.The viscosity of the uncompatibilized blends does not systematically rise with increasing PA12 content.Especially at a temperature of 230 °C, the viscosity of the blend at lower shear rates is even lower than the viscosity of pure PP1.At higher shear rates, the blends have a viscosity that is slightly higher than or similar to the viscosity of pure PP1.The standard deviations of the measurements are very low and therefore are not shown in the charts.
The blends with PP1 with coupling agent show a different behavior.The measured viscosities of those blends are slightly higher than the viscosity of pure PP1 along the total shear rate range.The values of the mixtures are slightly above the values of pure PP1.This is not astonishing, as the highest concentration of PA12 is only 20 wt%.The parameters obtained by the mathematical fit of the curves, as well as their deviation from the measured data, are depicted in Table A4 in Appendix B. The measurement data can be found in Appendix B in Tables A5 and A6.All presented viscosity and shear rate values are so-called true values.
The previously described mixing rules were applied to the compatibilized immiscible mixtures in a fashion similar to the miscible ones.Figures 9 and 10 show the relative deviation of the mixing models.The mixing models were calculated using the measured values depicted in Figures 7 and 8.The graphs show that the Bingham-Heitmiller model yields the least relative deviation compared to the other mixing rules when using the shear rate-based calculation method.This holds true for all mixtures at both 230 °C and 250 °C.The Bingham-Heitmiller model also deviates the least at 250 °C when calculated through the shear stress-based method.At 230 °C, the Arrhenius-Bersted model is more accurate.In contrast to the miscible mixtures, the linear model results in the highest calculated relative deviation for the immiscible blends.Overall, the relative deviations are higher than those found in blends made from the two PP grades.Figures 7 and 8 show the measured shear viscosity curves of the uncompatibilized PP 1 -PA12 blends, the compatibilized PP 1 -PA12-PPgMAH blends and the pure materials at temperatures of 230 • C and 250 • C. Especially at higher shear rates, PA12 exhibits a far higher shear viscosity than PP 1 .The viscosity of the uncompatibilized blends does not systematically rise with increasing PA12 content.Especially at a temperature of 230 • C, the viscosity of the blend at lower shear rates is even lower than the viscosity of pure PP 1 .At higher shear rates, the blends have a viscosity that is slightly higher than or similar to the viscosity of pure PP 1.The standard deviations of the measurements are very low and therefore are not shown in the charts.
The blends with PP 1 with coupling agent show a different behavior.The measured viscosities of those blends are slightly higher than the viscosity of pure PP 1 along the total shear rate range.The values of the mixtures are slightly above the values of pure PP 1 .This is not astonishing, as the highest concentration of PA12 is only 20 wt%.The parameters obtained by the mathematical fit of the curves, as well as their deviation from the measured data, are depicted in Table A4 in Appendix B. The measurement data can be found in Appendix B in Tables A5 and A6.All presented viscosity and shear rate values are so-called true values.
The previously described mixing rules were applied to the compatibilized immiscible mixtures in a fashion similar to the miscible ones.Figures 9 and 10 show the relative deviation of the mixing models.The mixing models were calculated using the measured values depicted in Figures 7 and 8.The graphs show that the Bingham-Heitmiller model yields the least relative deviation compared to the other mixing rules when using the shear rate-based calculation method.This holds true for all mixtures at both 230 • C and 250 • C. The Bingham-Heitmiller model also deviates the least at 250 • C when calculated through the shear stress-based method.At 230 • C, the Arrhenius-Bersted model is more accurate.In contrast to the miscible mixtures, the linear model results in the highest calculated relative deviation for the immiscible blends.Overall, the relative deviations are higher than those found in blends made from the two PP grades.

Discussion
To obtain a better overview of the results, the mean absolute deviations of the individual mixing rules were determined for both the shear rate-based and the shear stress-based calculation method.The results are shown in Figures 11 and 12.In general, using the shear stress-based calculation method predicts more accurate mixture viscosities for nearly all employed mixing rules, which is demonstrated by the low deviations.Exceptions are the linear mixing rule and the Bingham-Heitmiller model in the case of the immiscible mixture at 230 • C.
individual mixing rules were determined for both the shear rate-based and the shear stress-based calculation method.The results are shown in Figures 11 and 12.In general, using the shear stress-based calculation method predicts more accurate mixture viscosities for nearly all employed mixing rules, which is demonstrated by the low deviations.Exceptions are the linear mixing rule and the Bingham-Heitmiller model in the case of the immiscible mixture at 230 °C.
Considering the mixture of PP1 and PP2 (see Figure 11), the shear rate-based calculation method with the linear model yields the best results.This deviates the least from the measured viscosity data.For all other models, the shear stress-based calculation method is more accurate compared to the shear rate-based calculation method due to the smaller averaged deviation.Nevertheless, the shear rate-based calculation of the linear model delivers the best results.
In summary, based on the results shown, it can be recommended to apply the shear rate-based calculated linear mixing rule for homogeneous mixtures.Looking at the PP1-PA12 blends, or rather the effects of the compatibilizer used, various points need to be addressed.First, in this study, the PPgMAH influences the viscosity of the pure PP slightly, as can be seen in Figures 7 and 8.Other authors [21][22][23] found a strong effect of the compatibilizer on the viscosity of immiscible polymer blends.Even though the viscosity was just improved slightly in this study, it was clearly affected by the compatibilizer.
The measured viscosity of the uncompatibilized mixture cannot be described by any of the considered mixing models, as the values are below the values of the lower viscous material.The model of Lin [29] could probably describe this behavior.
The measured viscosity of the compatibilized mixture is between the viscosities of the pure materials and can be modeled by a mixing model with low deviations, as can be seen in Figures 9,10    behavior that PP1 and PA12 may show because of the different polarity of these materials.Furthermore, the results show that the mixing rule of Arrhenius-Bersted also delivers good results when the shear stress-based calculation is applied.This also means that the mixing rules of Cragoe and Walther [45] would also lead to low deviations when they are applied in a shear stress-based calculation, as they behave similarly to Arrhenius's mixing rule.
In summary, based on the results shown, it can be recommended to apply the shear rate-based calculated mixing rule of Bingham-Heitmiller for heterogeneous mixtures.Kallel et al. [18] found a diminishing effect of the compatibilizer at high shear rates.As can be seen in Figures 9 and 10, all mixing models predict the viscosity of compatibilized immiscible blends clearly better at low shear rates.At high shear rates, the measured viscosity of an immiscible mixture is clearly lower than predicted by the mixing rules, so the effect observed by Kallel et al. [18]  Considering the mixture of PP 1 and PP 2 (see Figure 11), the shear rate-based calculation method with the linear model yields the best results.This deviates the least from the measured viscosity data.For all other models, the shear stress-based calculation method is more accurate compared to the shear rate-based calculation method due to the smaller averaged deviation.Nevertheless, the shear rate-based calculation of the linear model delivers the best results.
In summary, based on the results shown, it can be recommended to apply the shear rate-based calculated linear mixing rule for homogeneous mixtures.
Looking at the PP 1 -PA12 blends, or rather the effects of the compatibilizer used, various points need to be addressed.First, in this study, the PPgMAH influences the viscosity of the pure PP slightly, as can be seen in Figures 7 and 8.Other authors [21][22][23] found a strong effect of the compatibilizer on the viscosity of immiscible polymer blends.Even though the viscosity was just improved slightly in this study, it was clearly affected by the compatibilizer.
The measured viscosity of the uncompatibilized mixture cannot be described by any of the considered mixing models, as the values are below the values of the lower viscous material.The model of Lin [29] could probably describe this behavior.
The measured viscosity of the compatibilized mixture is between the viscosities of the pure materials and can be modeled by a mixing model with low deviations, as can be seen in Figures 9,10  It can be concluded that the Bingham-Heitmiller model allows good prediction of the viscosity of a compatibilized mixture through both methods.This does make sense, as the Bingham-Heitmiller model assumes that the single phases of a blend do not intermix, but rather flow in separate layers [30], behavior that PP 1 and PA12 may show because of the different polarity of these materials.Furthermore, the results show that the mixing rule of Arrhenius-Bersted also delivers good results when the shear stress-based calculation is applied.This also means that the mixing rules of Cragoe and Walther [45] would also lead to low deviations when they are applied in a shear stress-based calculation, as they behave similarly to Arrhenius's mixing rule.
In summary, based on the results shown, it can be recommended to apply the shear rate-based calculated mixing rule of Bingham-Heitmiller for heterogeneous mixtures.
Kallel et al. [18] found a diminishing effect of the compatibilizer at high shear rates.As can be seen in Figures 9 and 10, all mixing models predict the viscosity of compatibilized immiscible blends clearly better at low shear rates.At high shear rates, the measured viscosity of an immiscible mixture is clearly lower than predicted by the mixing rules, so the effect observed by Kallel et al. [18] is confirmed by the calculations and measurements.
As can be seen in Figures 7 and 8, the viscosity of a compatibilized immiscible mixture increases only slightly when 5% to 20% more viscous material is added.This effect can be described by the Bingham-Heitmiller model with low deviations.Theoretically, the effect should be much stronger if 5% to 20% less viscous material were added, as can be seen in Figure 2. A reduction in the viscosity of an uncompatibilized mixture of immiscible material, as shown in Figures 7 and 8, can occur when different materials slide against each other, as reported in several publications [30,31,[56][57][58][59][60].This effect can be counteracted by compatibilization with the PPgMAH, as seen in Figures 7 and 8b.
The mixing models discussed in this paper, as well as the calculation methods presented, may be used for the predictive modeling and simulation of the functional zones of a single-screw extruder in the scope of a recycling process.Especially for the melt delay zone, melting zone, and conveying (metering) zone, the possibility of predicting the shear viscosity of mixtures is of particular interest.

Conclusions
In this study, the viscosity of polymer blends was studied for recycling modeling purposes.Different mixing rules were employed to calculate the shear rate-dependent viscosity curve of a mixture from viscosity curves of the involved pure materials.For this purpose, a shear rate-based calculation method and a newly developed stress-based calculation method were compared.The predicted mixture viscosities were compared to experimentally determined values.The polymer blends were created from two PP grades and from PP and PA12 at varying proportions, with the latter mixtures also being blended with 4 wt% compatibilizer of the PP.To assess the accuracy of the mixing rules used, the calculated viscosity was compared to data of the same blends measured with a high-pressure capillary rheometer.
The main conclusions of this paper are that for the miscible PP1-PP2 blends, several mixing rules can be applied to achieve an estimation of the viscosity with a deviation of less than 5% compared to the measured viscosity.Using the new shear stress-based calculation method, even smaller deviations can be reached.Out of the subjected mixing rules, a simple linear model, surprisingly, deviated the least from the measurements.
For the immiscible blends made from PP and PA12, a much higher average absolute deviation was calculated for most models.Only the model proposed by Bingham and Heitmiller differs less than 5% from the measured viscosity at both applied temperatures.In contrast to the miscible blends, the linear model is by far the most inaccurate.
Based on the results shown, it is recommended to apply the shear rate-based linear mixing model for mixtures of different grades from the same polymer type and the shear rate-based Bingham-Heitmiller model for compatibilized mixtures of different polymers.The shear stress-based models partially deliver better results, but they are much more complicated to apply due to the more complex evaluation and higher number of evaluation steps.The viscosity of uncompatibilized immiscible mixtures cannot be estimated by the applied mixing models, as the viscosity of the blends is partially lower than the viscosity of the lower viscous incorporated material.The proposed mixing models can be used to estimate the viscosity of mixtures, which is crucial when designing an extrusion process for recycling or for recycled materials.zone, melting zone, and conveying (metering) zone, the possibility of predicting the shear viscosity of mixtures is of particular interest.

Conclusions
In this study, the viscosity of polymer blends was studied for recycling modeling purposes.Different mixing rules were employed to calculate the shear rate-dependent viscosity curve of a mixture from viscosity curves of the involved pure materials.For this purpose, a shear rate-based calculation method and a newly developed stress-based calculation method were compared.The predicted mixture viscosities were compared to experimentally determined values.The polymer blends were created from two PP grades and from PP and PA12 at varying proportions, with the latter mixtures also being blended with 4 wt% compatibilizer of the PP.To assess the accuracy of the mixing rules used, the calculated viscosity was compared to data of the same blends measured with a highpressure capillary rheometer.
The main conclusions of this paper are that for the miscible PP1-PP2 blends, several mixing rules can be applied to achieve an estimation of the viscosity with a deviation of less than 5% compared to the measured viscosity.Using the new shear stress-based calculation method, even smaller deviations can be reached.Out of the subjected mixing rules, a simple linear model, surprisingly, deviated the least from the measurements.
For the immiscible blends made from PP and PA12, a much higher average absolute deviation was calculated for most models.Only the model proposed by Bingham and Heitmiller differs less than 5% from the measured viscosity at both applied temperatures.In contrast to the miscible blends, the linear model is by far the most inaccurate.
Based on the results shown, it is recommended to apply the shear rate-based linear mixing model for mixtures of different grades from the same polymer type and the shear rate-based Bingham-Heitmiller model for compatibilized mixtures of different polymers.The shear stress-based models partially deliver better results, but they are much more complicated to apply due to the more complex evaluation and higher number of evaluation steps.The viscosity of uncompatibilized immiscible mixtures cannot be estimated by the applied mixing models, as the viscosity of the blends is partially lower than the viscosity of the lower viscous incorporated material.The proposed mixing models can be used to estimate the viscosity of mixtures, which is crucial when designing an extrusion process for recycling or for recycled materials.

Figure 1 .
Figure 1.Setup of the twin-screw extruder with the temperature profile used (all values in °C) along the cylinder and the die and the position of the dosing systems.

Figure 1 .
Figure 1.Setup of the twin-screw extruder with the temperature profile used (all values in • C) along the cylinder and the die and the position of the dosing systems.

Figure 2 .
Figure 2. Exemplary evaluation of the predicted mixture viscosity using different mixing models, assuming two polymeric materials without considering shear rate-dependent behavior.Calculations according to Equations (1) and (3)-(8).The viscosity of material A is 9636 Pa s and that of material B is 2239 Pa s.

Figure 2 .
Figure 2. Exemplary evaluation of the predicted mixture viscosity using different mixing models, assuming two polymeric materials without considering shear rate-dependent behavior.Calculations according to Equations (1) and (3)-(8).The viscosity of material A is 9636 Pa s and that of material B is 2239 Pa s.

Figure 3 .
Figure 3. Shear viscosity curves of mixtures composed of PP1 and PP2 with different concentrations and shear viscosity curves of the pure polymers including the BCY fits at a testing temperature of 200 °C.

Figure 3 .
Figure 3. Shear viscosity curves of mixtures composed of PP 1 and PP 2 with different concentrations and shear viscosity curves of the pure polymers including the BCY fits at a testing temperature of 200 • C.

Figure 5 .
Figure 5. Impact strength of the PP1-PA12 specimen with different concentration of the coupling agent PPgMAH and of the pure material PP1.The concentration of PA12 is always 20%.The signs show the median values of six measurements, the error bars show the maximum and the minimum value.The impact strength of the pure PP1 and the pure PA12 are shown as horizontal lines.

Figure 5 .
Figure 5. Impact strength of the PP 1 -PA12 specimen with different concentration of the coupling agent PPgMAH and of the pure material PP 1 .The concentration of PA12 is always 20%.The signs show the median values of six measurements, the error bars show the maximum and the minimum value.The impact strength of the pure PP 1 and the pure PA12 are shown as horizontal lines.

Figure 6 .
Figure 6.SEM images of cryogenically broken pellets of PP1-PA12 with and without coupling agent for different PA12 concentrations.

Figure 6 .
Figure 6.SEM images of cryogenically broken pellets of PP 1 -PA12 with and without coupling agent for different PA12 concentrations.

Figure 7 .
Figure 7. Shear viscosity curves of the PP1-PA12 blends and the pure materials including the BCY fits of uncompatibilized (a) and compatibilized (b) PP1 with different concentrations of PA12 at a testing temperature of 230 °C.

Figure 7 . 25 Figure 8 .
Figure 7. Shear viscosity curves of the PP 1 -PA12 blends and the pure materials including the BCY fits of uncompatibilized (a) and compatibilized (b) PP 1 with different concentrations of PA12 at a testing temperature of 230 • C.

Figure 11 .
Figure 11.Averaged absolute relative deviation over the total shear rate or shear stress range for the different mixing models applied to mixtures of PP1 with all concentrations of PP2.Averaged absolute values of the deviations presented in Figure 4.
and 12.At 230 °C, the shear stress-based Arrhenius-Bersted mixing model delivers the lowest averaged absolute deviations, but the shear rate-based calculation of the Bingham-Heitmiller model shows only slightly higher deviations.At 250 °C, the shear stress-based calculation of the Bingham-Heitmiller mixing model

Figure 11 .
Figure 11.Averaged absolute relative deviation over the total shear rate or shear stress range for the different mixing models applied to mixtures of PP 1 with all concentrations of PP 2 .Averaged absolute values of the deviations presented in Figure 4.

Figure 12 .
Figure 12.Averaged absolute relative deviations of the mixing models applied to compatibilized PP1-PA12 blends at 230 °C and at 250 °C.Averaged absolute values of the relative deviations presented in Figures 9 and 10.
Figure 12.Averaged absolute relative deviations of the mixing models applied to compatibilized PP1-PA12 blends at 230 °C and at 250 °C.Averaged absolute values of the relative deviations presented in Figures 9 and 10.
and 12.At 230 • C, the shear stress-based Arrhenius-Bersted mixing model delivers the lowest averaged absolute deviations, but the shear rate-based calculation of the Bingham-Heitmiller model shows only slightly higher deviations.At 250 • C, the shear stress-based calculation of the Bingham-Heitmiller mixing model delivers the lowest deviations, but the shear rate-based calculation of the Bingham-Heitmiller model also shows deviations below 5%.

Author
Contributions: Conceptualization, G.Z., C.K., M.L. and E.W.; methodology, G.Z. and C.K.; software, C.K. and E.W.; validation, G.Z.; formal analysis, C.K. and E.W.; investigation, M.L. and E.W.; resources, C.K.; data curation, C.K.; writing-original draft preparation, M.L. and E.W.; writingreview and editing, G.Z. and C.K.; visualization, C.K. and E.W.; supervision, G.Z.; project administration, G.Z. and C.K.; funding acquisition, G.Z.All authors have read and agreed to the published version of the manuscript.Funding: This research was carried out in the course of the project NaKuRe and was financed by research funds granted by the Government of Upper Austria over the course of the program FTI-Struktur, project Wi-2021-305611/21-Au/NaKuRe.Polymers 2024, 16, x FOR PEER REVIEW 18 of 25

Table 1 .
Polymeric materials used, including their abbreviations for this work, the MFR/MVR, and the processing temperature, as given in the material data sheets.
1Derived from the material data sheet.

Table 2 .
Blends of PP 1 , PA12 and PPgMAH for the preliminary mechanical trials.

Base Polymer Concentration Base Polymer (wt%) Contaminating Polymer Concentration Contaminating Polymer (wt%) Coupling Agent Concentration Coupling Agent (wt%)
Polymers 2024,16,x FOR PEER REVIEW 7 of 25 the shear stress and the shear viscosity and deviates from the measured values.Thus, the next step is needed.4. The shear rate , can now be calculated from the viscosity of the mixture , and from the shear stress : , = , is confirmed by the calculations and Averaged absolute relative deviations of the mixing models applied to compatibilized PP 1 -PA12 blends at 230 • C and at 250 • C. Averaged absolute values of the relative deviations presented in Figures 9 and 10.

Table A2 .
Measured viscosity data of the PP 1 and PP 2 blends.MTSRate, measured true shear rate; MTVisc, measured true viscosity.First part of the data.

Table A3 .
Measured viscosity data of the PP 1 and PP 2 blends.MTSRate, measured true shear rate; MTVisc, measured true viscosity.Second part of the data.

Table A4 .
Fit parameters for the BCY model of uncompatibilized and compatibilized PP 1 and PA12 blends.

Table A5 .
Measured viscosity data of the compatibilized PP 1 and PA12 blends at 230 • C. MTSRate, measured true shear rate; MTVisc, measured true viscosity.First part of the data.

Table A6 .
Measured viscosity data of the compatibilized PP 1 and PA12 blends at 230 • C. MTSRate, measured true shear rate; MTVisc, measured true viscosity.Second part of the data.

Table A7 .
Measured viscosity data of the compatibilized PP 1 and PA12 blends at 250 • C. MTSRate, measured true shear rate; MTVisc, measured true viscosity.First part of the data.