Next Article in Journal
The β Form in PVDF Nanocomposites with Carbon Nanotubes: Structural Features and Properties
Next Article in Special Issue
Evolutionary Multi-Objective Optimization of Extrusion Barrier Screws: Data Mining and Decision Making
Previous Article in Journal
Metal-Organic Frameworks Applications in Synergistic Cancer Photo-Immunotherapy
Previous Article in Special Issue
Modeling and Experimental Studies on Polymer Melting and Flow in Injection Molding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization for the Contrary-Rotating Double-Screw Extrusion of Plastics

by
Andrzej Nastaj
and
Krzysztof Wilczyński
*
Polymer Processing Department, Faculty of Mechanical and Industrial Engineering, Warsaw University of Technology, Narbutta 85, 02-524 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(6), 1489; https://doi.org/10.3390/polym15061489
Submission received: 28 January 2023 / Revised: 10 March 2023 / Accepted: 13 March 2023 / Published: 16 March 2023

Abstract

:
A novel computer optimization system for the contrary-rotating double-screw extrusion of plastics was developed in this study. The optimization was based on the process simulation performed with the use of the global contrary-rotating double-screw extrusion software TSEM. The process was optimized using the GASEOTWIN software developed for this purpose using genetic algorithms. Several examples of optimization of the contrary-rotating double screw extrusion process parameters, i.e., the extrusion throughput, and minimize the plastic melt temperature and the plastic melting length.

Graphical Abstract

1. Introduction

Nowadays, computer simulations facilitate the development of polymer processing techniques and enable the prediction of the process course based on process parameters, i.e., material data, geometry data, and operating data. Simulations, however, do not allow researchers to conduct this process in reverse, i.e., to compute the process parameters to obtain the optimal output parameters.
Extrusion is the most used mass technique in the polymer processing industry. It is used for manufacturing films, sheets, pipes, and profiles, as well as for specialized operations, e.g., for compounding, reinforcing, pelletizing, etc. Extruders are divided into single-screw and double (or twin)-screw machines. Double-screw machines can be co-rotating or contrary-rotating (or counter-rotating).
Underwood [1] and Verbraak and Meijer [2] were the first to perform experiments to optimize extrusion. The main drawback of this approach is the number of experiments needed. Optimization by experiments is time-consuming, expensive, and nothing warranted that global optimum has been found.
Therefore, process simulations with statistical support seem to be better. The first optimizations using simulations were performed by Tadmor and Klein [3], Maddock and Smith [4], as well as by Helmy and Parnaby [5]. Later, Potente and Krell [6] used the REX software [7,8,9], Thibodeau and Lafleur [10,11] used the software of Ecole Polytechnique de Montreal [12,13], and the authors used the SSEM (Single Screw Extrusion Model) software. The drawback of these statistical approaches was the number of simulations needed, and a danger of selecting the local optima.
Covas and Gaspar-Cunha proposed a new strategy for extrusion optimization based on genetic algorithms. They developed procedures for the optimization of single-screw extrusion [14,15,16] and later for co-rotating double-screw extrusion [17,18,19]. The authors also used genetic algorithms for the optimization of single-screw extrusion [20]. The concept of genetic algorithms was also applied for scaling up the extrusion [21,22], as well as for the optimization of injection molding [23,24]. Nastaj and Wilczyński developed the optimization procedures for starve-fed single-screw extrusion [25,26] using the original models of this process [27,28,29,30].
Genetic algorithms are characterized by the following features:
-
Optimization variables are coded;
-
Searching for a solution starts from some population, which means the probability of getting stuck at a local extreme is low;
-
The rules of selection are probabilistic;
-
Objective functions are used, and derivatives are not necessary.
So far, there is a lack of optimization procedures for contrary-rotating double-screw extrusion, although computer models of this process have been recently proposed [31,32,33].
In this paper, a novel computer optimization system for the contrary-rotating double-screw extrusion of plastics was developed. Optimization was based on the process simulations performed with the use of the global contrary-rotating double-screw extrusion software TSEM [34]. The process was optimized using the GASEOTWIN software developed for this purpose using genetic algorithms. Several examples of the parameter optimization of the contrary-rotating double-screw extrusion process were studied in order to maximize the flow rate, i.e., the extrusion throughput, and minimize the plastic melt temperature and the plastic melting length, which are presented in the following sections.

2. Contrary-Rotating Double-Screw Extrusion

Double (or twin)-screw extrusion is divided into co-rotating extrusion (the screws are rotating in the same direction), and contrary-rotating (or counter-rotating) extrusion (the screws are rotating in the opposite direction (Figure 1)).
Contrary-rotating double-screw extruders are mainly applied for processing thermally sensitive plastic, e.g., polyvinyl chloride (PVC). Co-rotating double-screw extruders have specialized applications, e.g., for compounding, mixing, filling, and reinforcing plastics. The scheme of contrary-rotating double-screw extrusion is depicted in Figure 1.
Contrary-rotating double-screw extruders, compared with single-screw extruders, provide better feeding capability, as they can feed the machine with materials in the form of powder or with materials exhibiting slip properties. In contrary-rotating extruders, the material essentially does not flow from one screw to the other, as in co-rotating extruders. In the case of contrary-rotating machines, there is a co-rotating movement in the inter-screw gap, so high shear stresses are not produced in this area, as in co-rotating machines.
Contrary-rotating extruders are usually fed with dosing (metered feeding or starve feeding). The flow of plastics in these machines is completely different from the flow in the single-screw and co-rotating extruders. This flow results from a positive displacement mechanism, which does not appear in other extruders. The degree of positive displacement is dependent on the degree of screw meshing. It takes place most fully in closely intermeshing counter-rotating machines.
The material in the contrary-rotating extruder was transported in a C-shaped chamber (Figure 1b), and some leakage flows were also observed (Figure 1a), i.e., the calendering flow Qc, the flight flow Qf, the pressure (tetrahedron) flow Qt, and the side flow Qs. The C-chamber was formed by six surfaces, namely the screw root surface, the inner barrel surface, the side surfaces of the screw flights (twice), and the front surfaces of the screw flights (twice). There were distinguished leakage flows between these surfaces, i.e., the calendering flow Qc between the screw root and the screw flight; the flight flow Qf over the screw flights, through the clearance between the screw flight and the barrel; the back pressure inter-screw flow Qt (the tetrahedron flow) through the tetrahedral clearance between the flight flanks of screws, in the radial direction; and the side flow Qs through the side gap between the flight flanks of screws, in the tangential direction.

3. Modeling of Contrary-Rotating Double-Screw Extrusion

Some fundamental books present the state of the art in the modeling of polymer extrusion, e.g., [35,36,37,38,39], as well as some papers, e.g., [40,41,42,43]. Wilczyński et al. summarized this body of research in a review paper [44].
The fundamentals of contrary-rotating double-screw extrusion were first discussed a hundred years ago [45,46,47,48], and many designs of these machines were later developed [49,50].
Contrary-rotating extruders are entirely different from single-screw extruders, as well as from co-rotating machines. They were first presented by Kiesskalt [46] and Schenkel [49] as positive displacement pumps. Doboczky [51] and Janssen [52] proposed the flow-pumping characteristics of these extruders, and they considered the leakage flows. White and Adewale [53] performed the modeling of the flow in these machines while taking the level of screw intermeshing into consideration. Numerical FEM simulations were carried out by Li and Manas-Zloczower [54] and by Kajiwara et al. [55]. Hong and White [56,57] developed a FAN analysis for the non-Newtonian flow in these machines. They presented the concept of screw characteristics, which enabled them to model the flow for various screw designs and compute the pressure, fill factors, and temperature profiles.
Investigations of solid plastic conveying and plastic melting in contrary-rotating machines are very limited [51,52]. Wilczyński and White were the first who experimentally investigated and modeled the melting process [58,59]. Further studies were presented by Wang and Min [60,61] and by Wilczyński et al. [62].
White et al. [56,57,58,59] developed a theory that allowed for the prediction of both the pumping capacity and plastic melting profiles in these machines. The first model of this process was presented by Wilczynski et al. [31] to predict solid conveying, melting, and melt flow. Jiang et al. [63] developed a model for contrary-rotating double-screw extruders with flood feeding.
Lewandowski et al. [33] proposed an integrated model for the entire double-screw extruder, including the die. This model is based on combining melt-conveying models with melting and solid-conveying models. Three-dimensional, non-Newtonian FEM computations for melt conveying were performed to determine the pumping characteristics of the screws, which were included in the global model. To our knowledge, this was the first and, up to now, the only available global model of contrary-rotating double-screw extruders based on three-dimensional non-Newtonian FEM computations. This model was later extended, and it is now part of the software called Multi-Screw System [34].
The model we used in the studies presented here was based on our previous experimental investigations, which were discussed in the literature [58,62]. In these investigations, the “screw pulling-out technique” was applied to observe the plastic flow along the screws. The screws were removed from the barrel, the plastic contained in the screw channel was estimated, and the plastic samples were scrapped off from each screw. In the hopper section, there were solid granules (pellets) freely transported along the screws. The pellets were collected at the bottom part of the barrel adjacent to the pushing flights of the screws. The pellets were heated with the heat from the barrel and by being dragged into the gap between the screws. They formed some kind of solid bed, which decreased in length along the screws. The solid bed moved along the screw as a part of the C-chamber, and it was dragged into the calendering gap. High pressure was developed in the gap region, where friction forces initiated the melting process. Afterward, melting proceeded from the hot barrel, and a melt layer was formed between the barrel and the solid bed, which was scrapped off using the screw flights. The molten plastics flowed in a starved manner. The area close to the die was pressurized, and the screw became fully filled with plastic. Finally, the plastic flowed through the die under the pressure developed in the screw channel.
Based on these observations, models were developed for melting in both those regions, i.e., in the calendering gap and in the melt layer. Melting was initiated between the screws by the friction on the pellets, which was caused by the calendering stresses between the screws. The melting action at the barrel was induced by the barrel temperature being higher than the melting point and was propagated by the viscous dissipation heating of the melt film produced.
The model we used was experimentally validated for various screw configurations, various operating conditions, and various materials, namely polypropylene (PP), low-density polyethylene (LDPE), and high-density polyethylene (HDPE). This was presented in the literature, e.g., [31,33,34,59].
An example of computations with the use of TSEM software (and experimentation) is shown in Figure 2 for the screw system depicted in Figure 3. The experiment was performed for the extrusion of polypropylene (PP) with the screw speed of N = 100 rpm and the mass flow rate (throughput) of G = 10 kg/h. The dimensionless extrusion characteristics are presented with pressure (P—navy blue line) and temperature (T—green line) profiles, while the profile of plastic melting is represented by the solid-bed profile (SBP—blue line), and the profile of screw filling is represented by the fill factor (FF—red line). As can be seen, the pressure was only built in the region of fully filled screws. It is also interesting to note that the screws were fully filled with the material only at the final short section before the die. In the rest of the extruder, the screws were only 10–30% filled.
Another example of computations with the use of TSEM software (and experimentation) is shown in Figure 4 for the screw system depicted in Figure 5. This experiment was performed for the extrusion of low-density polyethylene (LDPE) at the screw speed of N = 80 rpm with the mass flow rate (throughput) of G = 8 kg/h. In this case, it is also seen that the pressure is built only in the region of fully filled screws. However, the screws were also fully filled with the material in the region of shearing elements (ZSS), and the pressure also developed in this region. In the rest of the extruder, the screws were only 20–30% filled.

4. Optimization Procedure

The GASEOTWIN software (Genetic Algorithm Screw Extrusion Optimization) has been developed to solve the optimization issue. Simulations with the use of the TSEM (Twin-Screw Extrusion Model) software [44] were the source of optimization data. The GASEOTWIN software is characterized by TSEM-specific data exchange between the simulation program and the optimization program.
The GASEOTWIN software, combined with the TSEM simulation software, enables process optimization with any number of optimized parameters and various criteria of optimization criteria. The accuracy of searching for the optimized response surface is determined by the number of divisions of the data range, which results from the length of the series of these numbers in a binary form. In GASEOTWIN software, the length of binary series is adjustable, and its maximum length is equal to 255 characters.
A “roulette wheel” concept was used for selection with an elitist strategy (Figure 6). In this strategy, the best dataset automatically continues to be used in the next generation, protecting the algorithm from losing the dataset with the highest value of objective functions. The algorithm stops if there is no improvement in the value of the objective function for 100 iterations. An action of the “roulette wheel” is presented in Figure 4. The areas of the “roulette wheel” defining the genotypes were proportional to the objective functions determined by these genotypes. As an example, the genotype Ge8 yielded the highest objective function F8 = 0.9996 and occupied a surface equal to 16.62% of the entire surface of this “roulette wheel”. The genotype Ge2 had the lowest objective function F2 = 0.1772, occupying a surface equal to 2.95% of the entire surface of the “roulette wheel”.
The optimization procedures were determined using GA parameters, i.e., the number of optimized variables, the size of the initial population, the length of chromosomes, the probability of crossover and mutation, and the points of crossover.

5. Optimization

5.1. Research Program

Optimization was performed for the contrary-rotating double-screw extrusion of polypropylene (PP). The calculations were carried out for a Lesitritz LSM 30/34 modular intermeshing contrary-rotating double-screw extruder. This is the machine intended for compounding. We considered the flow in the screw configurations depicted in Figure 3. The screws had a 34 mm diameter with a 30 mm distance between the center lines.
Polypropylene (PP) has a density of 0.904 g/cm3 (solid) and 0.739 g/cm3 (melt), and a melt flow rate (MFR) of 2.7 g/10 min (230 °C, 2.16 kg). The rheological flow properties of PP were determined on the basis of capillary rheometry and modeled using the Klein equation as follows:
l n η = A 0 + A 1 l n γ ˙ + A 11 l n 2 γ ˙ + A 12 T l n 2 γ ˙ + A 2 T + A 22 T 2
where η is the viscosity; γ ˙ is the shear rate; T is the temperature; and A0, A1, A11, A12, A2, and A22 are the model parameters (A0 = 14.0587, A1 = −0.4535, A11 = −0.0281, A12 = 2.71 × 10−4, A2 = −0.0316, and A22 = 4.45 × 10−5).
Optimization was carried out to maximize the extrusion output Qmax (kg/h), minimize the plastic temperature at the die outlet Tout (°C), and minimize the plastic melting length Lpl (mm). Contrary-rotating double-screw extruders are often used for the extrusion of thermally sensitive polymers; thus, minimizing the plastic temperature is important. Moreover, minimizing the temperature means minimizing power consumption. Minimizing the plastic melting length provides a sufficiently large flow space for a good mixing process of the plasticized material.
The optimized parameters were the screw rotation, the cylinder temperatures, and the extrusion throughput. These are presented in Table 1, where their range is also shown. Optimization was carried out for various weighted criteria, which were divided into two groups. In the first group (Table 2, items 1, 2, and 3), there was one criterion with a dominant weight (wi = 0.8), whereas in the second group (Table 2, items 4, 5, and 6), there were two criteria with both dominant and equal weight (wi = 0.4).
The global objective function is defined as
F i = w Q · Q i _ n o r m + w T o u t · T o u t   i _ n o r m + w L p l · L p l   i _ n o r m
where the output parameters (optimization criteria) are normalized as
Q i _ n o r m = Q i Q m i n Q m a x Q m i n
T o u t   i _ n o r m = T o u t   m a x T o u t   i T o u t   m a x T o u t   m i n
L p l   i _ n o r m = L p l   m a x L p l   i L p l   m a x L p l   m i n
where Fi is the global objective function, Qi_norm is the normalized flow rate (extrusion throughput), wQ is the weight of flow rate, Tout i_norm is the normalized plastic temperature at the die outlet, wTout is the weight of plastic temperature at the die outlet, Lpl i_norm is the normalized plastic melting length, wLpl is the weight of plastic melting length, and i is the number of the next value from the dataset.

5.2. Results

The optimization results are presented for various weights of the optimization criteria in Table 2.
Based on the obtained results, it can be concluded that, despite the different weights of the optimization criteria, the most common solutions were the sets of operating parameters that were similar in value (Table 2, items 1, 2, and 4). They had the following values: N = 227.5 rpm, TII = 180 °C, TIII = 180 °C, TIV = 180 °C and 187.5 °C, and Q = 65.63 kg/h and 62.50 kg/h. The optimization results showed that the operating parameters of the process both in the case when the highest weight was assumed for the mass flow rate (extrusion throughput; wQ = 0.8) and in the case when the highest weight was assumed for the lowest temperature of the material at the outlet of the die (wTout = 0.8) were also close together. It follows that the selection of either of these sets of operating parameters (Table 2, items 1 and 2) allows us to obtain the course of the extrusion process with relatively the highest throughput and the lowest plastic melt temperature.
Process simulations were performed for each set of optimal parameters (N, TII, TIII, TIV, and Q). The simulation results using the optimal parameters for various weights of optimization criteria are presented in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18. The overall process characteristics (extrusion characteristics) obtained with the optimal parameters are depicted in Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12, whereas the process variables, i.e., the pressure profiles (P—navy blue line), the temperature profiles (T—green line), the solid-bed profiles (SBP—blue lines), and the fill-factor profiles (FF—red line), using the optimal process parameters (N, TII, TIII, TIV, and Q) are depicted in Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18.
It is interesting to note that, as shown in Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18, the simulation results for optimal sets (3) and (6) differed from other simulations. These simulations were carried out with the weights of wQ = 0.1, wTout = 0.1, and wLpl = 0.8 as well as wQ = 0.2, wTout = 0.4, and wLpl = 0.4, both indicating the cases in which the weight of the polymer melting length wLpl had the dominating value.
An example of the optimization results (for criteria weights wQ = 0.8, wTout = 0.1, and wLpl = 0.1) as shown on the screen of the GASEOTWIN program is depicted in Figure 19. The parameters of optimization are also seen, as well as the values of optimal parameters.

6. Conclusions

A novel computer optimization system for contrary-rotating double-screw extrusion of plastics was developed. Optimization was based on the process simulation performed with the use of the global contrary-rotating double-screw extrusion software TSEM. The process was optimized using the GASEOTWIN software, developed for this purpose using genetic algorithms.
An example of using this system was presented to show the possibilities of the developed system. Optimization was carried out in order to maximize the flow rate, i.e., the extrusion throughput, minimize the plastic melting temperature at the die outlet, and minimize the plastic melting length. These parameters are important when using double-screw extruders. Contrary-rotating extruders are often used for the extrusion of thermally sensitive plastics; thus, minimizing the plastic temperature is important. Moreover, minimizing the plastic temperature means minimizing power consumption. Minimizing the plastic melting length provides a sufficiently large flow space for the good mixing of the plasticized material. Additionally, these optimization criteria were used with various weights, which depend on the manufacturer’s requirements and are decided by the manufacturer. In this example, the optimized parameters were the screw rotation, the cylinder temperatures, and the extrusion throughput, that is, the basic operating parameters of the process. Of course, optimization can be also performed for the optimization of the geometry of screws.
To our knowledge, this is the first optimization software for contrary-rotating double-screw extrusion that uses process simulations. It enables the optimization of extrusion process parameters using various optimization criteria with various weights.

Author Contributions

Conceptualization, A.N. and K.W.; methodology, A.N. and K.W.; software, A.N.; validation, A.N.; investigation, A.N.; writing—original draft preparation, A.N. and K.W.; writing—review and editing, K.W.; visualization, A.N.; supervision, K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Underwood, W.M. Experimental Method for Designing Extrusion Screws. Chem. Eng. Prog. 1962, 58, 59. [Google Scholar]
  2. Verbraak, C.P.J.M.; Meijer, H.E.H. Screw Design in Injection Molding. Polym. Eng. Sci. 1989, 29, 479–487. [Google Scholar] [CrossRef]
  3. Tadmor, Z.; Klein, I. Engineering Principles of Plasticating Extrusion; Van Nostrand Reinhold Co.: New York, NY, USA, 1970; ISBN 978-044-215-635-0. [Google Scholar]
  4. Maddock, B.H.; Smith, D.J. Extruder Design by Computer Printout. SPE J. 1972, 28, 12–17. [Google Scholar]
  5. Helmy, H.A.A.; Parnaby, J. Computer—Aided Optimal Melt Screw Design. Polym. Eng. Sci. 1976, 16, 437–449. [Google Scholar] [CrossRef]
  6. Potente, H.; Krell, B. Screw optimization by means of DOE and multiple regression. In Proceedings of the 55th Annual Technical Conferences of the Society of Plastics Engineers (ANTEC 1997), Toronto, ON, Canada, 27 April–2 May 1997; pp. 135–139. [Google Scholar]
  7. Potente, H.; Hanhart, W.; Schöppner, V. Potential Applications for Computer-aided Extruder Design. Int. Polym. Process. 1993, 8, 335–344. [Google Scholar] [CrossRef]
  8. Potente, H.; Hanhart, W.; Reski, T. Design and Processing Optimization of Extruder Screws. Polym. Eng. Sci. 1994, 34, 937–945. [Google Scholar] [CrossRef]
  9. Potente, H.; Schöppner, V.; Ujma, A. Successful Simulation of Wall—Slipping Plastics. J. Polym. Eng. 1997, 17, 153–170. [Google Scholar] [CrossRef]
  10. Thibodeau, C.A.; Lafleur, P.G. Computer Design and Screw Optimization. In Proceedings of the 58th Annual Technical Conference of the Society of Plastics Engineers (ANTEC 2000), Orlando, FL, USA, 7–11 May 2000; pp. 276–282. [Google Scholar]
  11. Thibodeau, C.A.; Lafleur, P.G. Computer Design and Screw Optimization. In Proceedings of the Polymer Processing Society 16th Annual Meeting (PPS-16), Shanghai, China, 18–23 June 2000; p. 15. [Google Scholar]
  12. Amellal, K.; Lafleur, P.G. Computer Simulation of Conventional and Barrier Screw Extruders. Plast. Rub. Compos. Pro. 1993, 19, 227–239. [Google Scholar]
  13. Vincelette, A.R.; Guerrero, C.S.; Carreau, P.J.; Lafleur, P.G. A Model for Single-Screw Plasticating Extruders. Int. Polym. Process. 1989, 4, 232–241. [Google Scholar] [CrossRef]
  14. Covas, J.A.; Cunha, A.G.; Oliveira, P. An Optimization Approach to Practical Problems in Plasticating Single Screw Extrusion. Polym. Eng. Sci. 1999, 39, 443–456. [Google Scholar] [CrossRef]
  15. Covas, J.A.; Gaspar-Cunha, A. The Use of an Optimization Approach to the Design of Extrusion Screw. In Proceedings of the Polymer Processing Society 16th Annual Meeting (PPS-16), Shanghai, China, 18–23 June 2000. [Google Scholar]
  16. Gaspar-Cunha, A.; Covas, J.A. The Design of Extrusion Screw: An Optimization Approach. Int. Polym. Process. 2001, 16, 229–240. [Google Scholar] [CrossRef] [Green Version]
  17. Gaspar-Cunha, A.; Covas, J.A.; Vergnes, B. An Optimization Methodology for Setting the Operating Conditions in Twin-Screw Extrusion. In Proceedings of the Polymer Processing Society 18th Annual Meeting (PPS-18), Guimaraes, Portugal, 18–21 June 2002. [Google Scholar]
  18. Gaspar-Cunha, A.; Poulesquen, A.; Vergnes, B.; Covas, J.A. Optimization of Processing Conditions for Polymer Twin-Screw Extrusion. Int. Polym. Process. 2002, 17, 201–213. [Google Scholar] [CrossRef]
  19. Gaspar-Cunha, A.; Covas, J.; Vergnes, B. Defining the Configuration of Co-Rotating Twin-Screw Extruders with Multiobjective Evolutionary Algorithms. Polym. Eng. Sci. 2005, 45, 1159–1173. [Google Scholar] [CrossRef]
  20. Nastaj, A.; Wilczyński, K. Process Optimization for Single Screw Extrusion of Polymeric Materials—Simulation Studies. Polimery 2018, 4, 297–304. [Google Scholar] [CrossRef]
  21. Covas, J.A.; Gaspar-Cunha, A. Extrusion Scale-up: An Optimization-based Methodology. Int. Polym. Process. 2009, 24, 67–82. [Google Scholar] [CrossRef]
  22. Covas, J.A.; Gaspar-Cunha, A. A Scaling-up Methodology for Co-rotating Twin-Extruders. In Proceedings of the 27th Annual Meeting of the Polymer Processing Society (PPS-27), Marrakesh, Morocco, 10–14 May 2011; pp. 1–6. [Google Scholar]
  23. Fernandes, C.; Pontes, A.J.; Viana, J.C.; Gaspar-Cunha, A. Using Multiobjective Evolutionary Algorithms in the Optimization of Operating Conditions of Polymer Injection Molding. Polym. Eng. Sci. 2010, 50, 1667–1678. [Google Scholar] [CrossRef]
  24. Fernandes, C.; Pontes, A.J.; Viana, J.C.; Gaspar-Cunha, A. Using Multi-objective Evolutionary Algorithms for Optimization of the Cooling System in Polymer Injection Molding. Int. Polym. Process. 2012, 27, 213–223. [Google Scholar] [CrossRef] [Green Version]
  25. Nastaj, A.; Wilczyński, K. Optimization for Starve Fed/Flood Fed Single Screw Extrusion of Polymeric Materials. Polymers 2020, 12, 149. [Google Scholar] [CrossRef] [Green Version]
  26. Nastaj, A.; Wilczyński, K. Optimization and Scale-Up for Polymer Extrusion. Polymers 2021, 13, 1547. [Google Scholar] [CrossRef] [PubMed]
  27. Wilczyński, K.J.; Nastaj, A.; Lewandowski, A.; Wilczyński, K. A Composite Model for Starve Fed Single Screw Extrusion of Thermoplastics. Polym. Eng. Sci. 2014, 54, 2362–2374. [Google Scholar] [CrossRef]
  28. Wilczyński, K.J.; Lewandowski, A.; Nastaj, A.; Wilczyński, K. Modeling for Starve Fed/Flood Fed Mixing Single-Screw Extruders. Int. Polym. Process. 2016, 31, 82–91. [Google Scholar] [CrossRef]
  29. Wilczyński, K.J.; Lewandowski, A.; Nastaj, A.; Wilczyński, K. A Global Model for Starve-Fed Nonconventional Single-Screw Extrusion of Thermoplastics. Adv. Polym. Technol. 2017, 36, 23–35. [Google Scholar] [CrossRef]
  30. Wilczyński, K.J.; Nastaj, A.; Wilczyński, K. A Computer Model for Starve-Fed Single-Screw Extrusion of Polymer Blends. Adv. Polym. Technol. 2018, 37, 2142–2151. [Google Scholar] [CrossRef]
  31. Wilczyński, K.; Jiang, Q.; White, J.L. A Composite Model for Melting, Pressure and Fill Factor Profiles in a Metered Fed Closely Intermeshing Counter-Rotating Twin Screw Extruder. Int. Polym. Process. 2007, 22, 198–203. [Google Scholar] [CrossRef]
  32. Wilczyński, K.; Lewandowski, A. Study on the Polymer Melt Flow in a Closely Intermeshing Counter-Rotating Twin Screw Extruder. Int. Polym. Process. 2014, 29, 649–659. [Google Scholar] [CrossRef]
  33. Lewandowski, A.; Wilczyński, K.J.; Nastaj, A.; Wilczyński, K. A Composite Model for an Intermeshing Counter-Rotating Twin-Screw Extruder and its Experimental Verification. Polym. Eng. Sci. 2015, 55, 2838–2848. [Google Scholar] [CrossRef]
  34. Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich, Germany, 2021; ISBN 978-1-56990-660-6. [Google Scholar]
  35. Rauwendaal, C. Polymer Extrusion, 5th ed.; Carl Hanser Verlag: Munich, Germany, 2014; ISBN 978-1-56990-516-6. [Google Scholar]
  36. White, J.L.; Potente, H. Screw Extrusion. Science and Technology; Hanser Publishers: Munich, Germany, 2003; ISBN 978-3-446-19624-7. [Google Scholar]
  37. Tadmor, Z.; Gogos, C.G. Principles of Polymer Processing, 2nd ed.; John Wiley & Sons Inc.: New York, NY, USA, 2006; ISBN 978-0-471-38770-1. [Google Scholar]
  38. Osswald, T.; Hernandez-Ortiz, J.P. Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich, Germany, 2006; ISBN 978-3-446-40381-9. [Google Scholar]
  39. Agassant, J.F.; Avenas, P.; Carreau, P.J.; Vergnes, B.; Vincent, M. Polymer Processing. Principles and Modelling, 2nd ed.; Carl Hanser Verlag: Munich, Germany, 2017; ISBN 978-1-56990-605-7. [Google Scholar]
  40. Altinkaynak, A.; Gupta, M.; Spalding, M.A.; Crabtree, S.L. Melting in a Single Screw Extruder: Experiments and 3D Finite Element Simulations. Int. Polym. Process. 2011, 26, 182–196. [Google Scholar] [CrossRef]
  41. Ariffin, A.; Ahmad, M.S.B. Review: Single Screw Extruder in Particulate Filler Composite. Polym.-Plast. Technol. 2011, 50, 395–403. [Google Scholar] [CrossRef]
  42. Teixeira, C.; Gaspar-Cunha, A.; Covas, J.A. Flow and Heat Transfer Along the Length of a Co-rotating Twin Screw Extruder. Polym.-Plast. Technol. 2012, 51, 1567–1577. [Google Scholar] [CrossRef]
  43. Malik, M.; Kalyon, D.M.; Golba, J.C., Jr. Simulation of Co-Rotating Twin Screw Extrusion Process Subject to Pressure-Dependent Wall Slip at Barrel and Screw Surfaces: 3D FEM Analysis for Combinations of Forward- and Reverse-Conveying Screw Elements. Int. Polym. Process. 2014, 29, 51–62. [Google Scholar] [CrossRef]
  44. Wilczyński, K.; Nastaj, A.; Lewandowski, A.; Wilczyński, K.J.; Buziak, K. Fundamentals of Global Modeling for Polymer Extrusion. Polymers 2019, 11, 2106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Leistritz, P.; Burghauser, F. German Patent. 1927, 453, 727.
  46. Kiesskalt, S. Untersuchungen an einer Kapsel Pumpe. VDI Z. 1927, 71, 453. [Google Scholar]
  47. Montelius, C.O.J. U.S. Patent. 1929, 1, 802.
  48. Montelius, C.O.J. Der Kogelige Danske Videnscabernes Selscabs Skrifter. Tek. Tidskr. 1933, 6, 61–63. [Google Scholar]
  49. Schenkel, G. Kunststoffe Extrudertechnik; Hanser: Munich, Germany, 1963. [Google Scholar]
  50. White, J.L. Twin Screw Extrusion. Technology and Principles; Hanser: Munich, Germany, 1990. [Google Scholar]
  51. Doboczky, Z. Theoretische and wirkliche Ausstoßleistung der Doppelschnecken Extruder. Plastverarbeiter 1965, 19, 395–400. [Google Scholar]
  52. Janssen, L.P.B.M. Twin Screw Extrusion; Elsevier: Amsterdam, The Netherlands, 1978. [Google Scholar]
  53. White, J.L.; Adewale, A. A Unified View of Modeling Flow in Counter-rotating Twin Screw Extruders. Int. Polym. Process. 1993, 8, 210–217. [Google Scholar] [CrossRef]
  54. Li, T.; Manas Zloczower, I. Flow Field Analysis of an Intermeshing Counter-Rotating Twin Screw Extruder. Polym. Eng. Sci. 1994, 34, 551–558. [Google Scholar] [CrossRef]
  55. Kajiwara, T.; Nagashima, Y.; Naakano, Y.; Funatsu, K. Numerical Study of Twin Screw Extruders by ThreeDimensional Flow Analysis—Development of Analysis Technique and Evaluation of Mixing Performance for Full Flight Screws. Polym. Eng. Sci. 1996, 36, 2142–2152. [Google Scholar] [CrossRef]
  56. Hong, M.H.; White, J.L. Fluid Mechanics of Intermeshing Counter-Rotating Twin Screw Extruders. Int. Polym. Process. 1998, 13, 342–346. [Google Scholar] [CrossRef]
  57. Hong, M.H.; White, J.L. Simulation of Flow in an Intermeshing Modular Counter-Rotating Twin Screw Extruder: Non-Newtonian and Non-Isothermal Behavior. Int. Polym. Process. 1999, 14, 136–143. [Google Scholar] [CrossRef]
  58. Wilczyński, K.; White, J.L. Experimental Study of Melting in an Intermeshing Counter-Rotating Twin Screw Extruder. Int. Polym. Process. 2001, 16, 257–262. [Google Scholar] [CrossRef]
  59. Wilczyński, K.; White, J.L. Melting Model for Intermeshing Counter-Rotating Twin-Screw Extruders. Polym. Eng. Sci. 2003, 43, 1715–1726. [Google Scholar] [CrossRef]
  60. Wang, D.; Min, K. In-Line Monitoring and Analysis of Polymer Melting Behavior in an Intermeshing Counter-Rotating Twin-Screw Extruder by Ultrasound Waves. Polym. Eng. Sci. 2005, 45, 998–1010. [Google Scholar] [CrossRef]
  61. Wang, D.; Min, K. Experiments and Analysis of Effect of Calender Gaps on Melting of PVC Powders in an Intermeshing Counter-Rotating Twin-Screw Extruder. Int. Polym. Process. 2006, 21, 17–23. [Google Scholar] [CrossRef]
  62. Wilczyński, K.; Lewandowski, A.; Wilczyński, K.J. Experimental Study of Melting of LDPE/PS Polyblend in an Intermeshing Counter-Rotating Twin Screw Extruder. Polym. Eng. Sci. 2012, 52, 449–458. [Google Scholar] [CrossRef]
  63. Jiang, Q.; White, J.L.; Yang, J. A Global Model for Closely Intermeshing Counter-Rotating Twin Screw Extruders with Flood Feeding. Int. Polym. Process. 2010, 25, 223–235. [Google Scholar] [CrossRef]
Figure 1. The plastic flow in a contrary-rotating double-screw extruder: (a) leakage flows, Qc—calendering flow, Qf—flight flow, Qt—pressure (tetrahedral) flow, Qs—side flow; (b) C-shaped chamber, 1—side surface of the screw flight, 2—barrel surface, 3—front surface of the screw flight, 4—surface of the screw root (adopted with permission from Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich 2021 [34]).
Figure 1. The plastic flow in a contrary-rotating double-screw extruder: (a) leakage flows, Qc—calendering flow, Qf—flight flow, Qt—pressure (tetrahedral) flow, Qs—side flow; (b) C-shaped chamber, 1—side surface of the screw flight, 2—barrel surface, 3—front surface of the screw flight, 4—surface of the screw root (adopted with permission from Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich 2021 [34]).
Polymers 15 01489 g001
Figure 2. General extrusion characteristics, simulation data using TSEM model—twin-screw extrusion, Q = 10 kg/h, screw speed N = 100 rpm: SBP—solid-bed profile.
Figure 2. General extrusion characteristics, simulation data using TSEM model—twin-screw extrusion, Q = 10 kg/h, screw speed N = 100 rpm: SBP—solid-bed profile.
Polymers 15 01489 g002
Figure 3. Screw configurations.
Figure 3. Screw configurations.
Polymers 15 01489 g003
Figure 4. General extrusion characteristics, simulation data using TSEM model—twin-screw extrusion, Q = 8 kg/h, screw speed N = 80 rpm: SBP—solid-bed profile (adopted with permission from Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich 2021 [34]).
Figure 4. General extrusion characteristics, simulation data using TSEM model—twin-screw extrusion, Q = 8 kg/h, screw speed N = 80 rpm: SBP—solid-bed profile (adopted with permission from Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich 2021 [34]).
Polymers 15 01489 g004
Figure 5. Screw configurations (adopted with permission from Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich 2021 [34]).
Figure 5. Screw configurations (adopted with permission from Wilczyński, K. Rheology in Polymer Processing. Modeling and Simulation; Carl Hanser Verlag: Munich 2021 [34]).
Polymers 15 01489 g005
Figure 6. A “roulette wheel”.
Figure 6. A “roulette wheel”.
Polymers 15 01489 g006
Figure 7. Extrusion characteristics (weights of optimization criteria: wQ = 0.8, wTout = 0.1, and wLpl = 0.1).
Figure 7. Extrusion characteristics (weights of optimization criteria: wQ = 0.8, wTout = 0.1, and wLpl = 0.1).
Polymers 15 01489 g007
Figure 8. Extrusion characteristics (weights of optimization criteria: wQ = 0.1, wTout = 0.8, and wLpl = 0.1).
Figure 8. Extrusion characteristics (weights of optimization criteria: wQ = 0.1, wTout = 0.8, and wLpl = 0.1).
Polymers 15 01489 g008
Figure 9. Extrusion characteristics (weights of optimization criteria: wQ = 0.1, wTout = 0.1, and wLpl = 0.8).
Figure 9. Extrusion characteristics (weights of optimization criteria: wQ = 0.1, wTout = 0.1, and wLpl = 0.8).
Polymers 15 01489 g009
Figure 10. Extrusion characteristics (weights of optimization criteria: wQ = 0.4, wTout = 0.4, and wLpl = 0.2).
Figure 10. Extrusion characteristics (weights of optimization criteria: wQ = 0.4, wTout = 0.4, and wLpl = 0.2).
Polymers 15 01489 g010
Figure 11. Extrusion characteristics (weights of optimization criteria: wQ = 0.4, wTout = 0.2, and wLpl = 0.4).
Figure 11. Extrusion characteristics (weights of optimization criteria: wQ = 0.4, wTout = 0.2, and wLpl = 0.4).
Polymers 15 01489 g011
Figure 12. Overall process characteristics for twin-screw extrusion (weights of optimization criteria: wQ = 0.2, wTout = 0.4, and wLpl = 0.4).
Figure 12. Overall process characteristics for twin-screw extrusion (weights of optimization criteria: wQ = 0.2, wTout = 0.4, and wLpl = 0.4).
Polymers 15 01489 g012
Figure 13. Pressure profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 1–3).
Figure 13. Pressure profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 1–3).
Polymers 15 01489 g013
Figure 14. Solid-bed profiles (SBPs) using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 1–3).
Figure 14. Solid-bed profiles (SBPs) using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 1–3).
Polymers 15 01489 g014
Figure 15. Fill-factor profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 1–3).
Figure 15. Fill-factor profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 1–3).
Polymers 15 01489 g015
Figure 16. Pressure profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 4–6).
Figure 16. Pressure profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 4–6).
Polymers 15 01489 g016
Figure 17. Solid-bed profiles (SBPs) using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 4–6).
Figure 17. Solid-bed profiles (SBPs) using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 4–6).
Polymers 15 01489 g017
Figure 18. Fill-factor profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 4–6).
Figure 18. Fill-factor profiles using optimal process parameters for various weights of optimization criteria: wQ, wTout, and wLpl (Table 2, results 4–6).
Polymers 15 01489 g018
Figure 19. Optimization results for twin-screw extrusion for weights of optimization criteria equal to wQ = 0.8, wTout = 0.1, and wLpl = 0.1.
Figure 19. Optimization results for twin-screw extrusion for weights of optimization criteria equal to wQ = 0.8, wTout = 0.1, and wLpl = 0.1.
Polymers 15 01489 g019
Table 1. Research program.
Table 1. Research program.
Screw Rotation
N, rpm
Cylinder Temperature TI, °CCylinder Temperature TII, °CCylinder Temperature TIII, °CCylinder Temperature TIV, °CExtrusion Output Q, kg/h
40 ÷ 240 180180 ÷ 240180 ÷ 240180 ÷ 2401 ÷ 100
Table 2. The objective functions (maxima) and the optimized parameters with various weights of optimization criteria.
Table 2. The objective functions (maxima) and the optimized parameters with various weights of optimization criteria.
Criterion Weights
wQ = 0.8
wTout = 0.1
wLpl = 0.1
wQ = 0.1
wTout = 0.8
wLpl = 0.1
wQ = 0.1
wTout = 0.1
wLpl = 0.8
wQ = 0.4
wTout = 0.4
wLpl = 0.2
wQ = 0.4
wTout = 0.2
wLpl = 0.4
wQ = 0.2
wTout = 0.4
wLpl = 0.4
No123456
Screw rotation
N, rpm
227.5227.5152.5227.524096.3
Cylinder temperature
TII, °C
180180219.4180180196.9
Cylinder temperature
TIII, °C
180180219.4180180180
Cylinder temperature
TIV, °C
187.5180230.6180187.5187.5
Extrusion throughput
Q, kg/h
65.6362.509.3865.6365.6312.50
Temperature of plastic at die outlet Tout, °C176.25181.49206.42175.75181.58197.24
Length of plastic melting
Lpl, mm
960920160960920240
Objective function Fi0.62430.85630.83730.66730.49340.5819
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nastaj, A.; Wilczyński, K. Optimization for the Contrary-Rotating Double-Screw Extrusion of Plastics. Polymers 2023, 15, 1489. https://doi.org/10.3390/polym15061489

AMA Style

Nastaj A, Wilczyński K. Optimization for the Contrary-Rotating Double-Screw Extrusion of Plastics. Polymers. 2023; 15(6):1489. https://doi.org/10.3390/polym15061489

Chicago/Turabian Style

Nastaj, Andrzej, and Krzysztof Wilczyński. 2023. "Optimization for the Contrary-Rotating Double-Screw Extrusion of Plastics" Polymers 15, no. 6: 1489. https://doi.org/10.3390/polym15061489

APA Style

Nastaj, A., & Wilczyński, K. (2023). Optimization for the Contrary-Rotating Double-Screw Extrusion of Plastics. Polymers, 15(6), 1489. https://doi.org/10.3390/polym15061489

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop