Next Article in Journal
Total Life Cycle of Polypropylene Products: Reducing Environmental Impacts in the Manufacturing Phase
Previous Article in Journal
Synthesis and Nonisothermal Crystallization Kinetics of Poly(Butylene Terephthalate-co-Tetramethylene Ether Glycol) Copolyesters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Network Analysis to Flow Systems with Alternating Wave Channels: Part B. (Superimposed Drag-Pressure Flows in Extrusion)

1
Institute of Polymer Extrusion and Compounding, Johannes Kepler University Linz, 4040 Linz, Austria
2
Kunststofftechnik Paderborn, Paderborn University, 33098 Paderborn, Germany
*
Author to whom correspondence should be addressed.
Polymers 2020, 12(9), 1900; https://doi.org/10.3390/polym12091900
Submission received: 24 July 2020 / Revised: 12 August 2020 / Accepted: 20 August 2020 / Published: 24 August 2020
(This article belongs to the Section Polymer Processing and Engineering)

Abstract

:
Due to progress in the development of screw designs over recent decades, numerous high-performance screws have become commercially available in single-screw extrusion. While some of these advanced designs have been studied intensively, others have received comparatively less attention. We developed and validated a semi-numerical network-theory-based modeling approach to predicting flows of shear-thinning polymer melts in wave-dispersion screws. In the first part (Part A), we systematically reduced the complexity of the flow analysis by omitting the influence of the screw rotation on the conveying behavior of the wave zone. In this part (Part B), we extended the original theory by considering the drag flow imposed by the screw. Two- and three-dimensional melt-conveying models were combined to predict locally the conveying characteristics of the wave channels in a discretized flow network. Extensive experiments were performed on a laboratory single-screw extruder, using various barrel designs and wave-dispersion screws. The predictions of our semi-numerical modeling approach for the axial pressure profile along the wave-dispersion zone accurately reproduce the experimental data. Removing the need for time-consuming numerical simulations, this modeling approach enables fast analyses of the conveying behavior of wave-dispersion zones, thereby offering a useful tool for design and optimization studies and process troubleshooting.

Graphical Abstract

1. Introduction

From a quantitative point of view, extrusion is the most important process in the polymer industry. Every year, this technique converts more than 114 million tons of polymeric materials to products ranging from sheets to profiles and from pipes to films [1]. To meet the ever-increasing demands of the plastics industry, a trend towards use of high-output extruders has emerged over recent decades. These developments have also resulted in a variety of high-performance screws being designed to provide excellent melt homogeneity at high output rates and proper discharge temperatures. Despite their importance, only a few studies have examined the performance of wave-dispersion screws [2,3,4,5,6].
The concept of wave-dispersion screws was initially proposed by Kruder [7]: Replacing the metering zone of conventional extruder screws, a so-called wave section is implemented in which the channel depth oscillates over a plurality of cycles. Assuming the resin to be mostly in a molten state, the deep valley portions minimize heat generation in the screw channel, while the shallow peak portions ensure repeated intensive mixing.
The efficiency of the originally proposed concept depends strongly on the solids content in the extrudate, as solid resin particles may restrict the flow path at the wave peaks, causing the extrusion process to become unstable. To improve the performance, Kruder [8] developed the double-wave screw. This term refers to wave sections consisting of two parallel screw channels separated by a barrier flight that is undercut relative to the main flight. The wave cycles in the adjacent channels are out of phase, which means that, as one sub-channel decreases in depth, the other sub-channel increases. The helical offset is typically arranged such that a wave peak of one channel lies opposite a valley in the other and vice versa. In contrast to the original single-flighted design, this concept allows the polymer melt approaching a wave peak to pass across the undercut barrier, thereby minimizing the risk of blockage and pressure fluctuations. With the material being repeatedly transferred between the channels, the screw design promotes cross-channel mixing and improves melt homogeneity.
A modified concept of the double-wave screw, trademarked as energy-transfer screw, was developed by Chung and Barr [9]. In their design, the clearances of the screw flights are selectively interrupted such that the functions of the main and the secondary flights alternate. In contrast to the double-wave screw, in which the material is conveyed within the same pair of subchannels, the energy-transfer screw allows the material to flow only in the direction of the natural flow caused by the rotation of the screw. While this concept once again improves the mixing performance of the screw, it can limit the pressure build-up capacity. Several wave-dispersion screw patents followed these developments [10,11,12,13].
Due to their great mixing capability, wave-dispersion zones were originally located at the discharge end of the extruder. In recent years, the geometrical design was extended to the melting zone to break up solid agglomerates and mix it with molten material in the screw channel.

1.1. Analysis of Melt Conveying and Pressurization

The analysis of polymer-melt flow in single-screw extruders has been the subject of numerous theoretical studies. Most of these examined the flow in conventional metering sections, whereas wave-dispersion zones received less attention in comparison. The first model of melt conveying in single-screw extruders was published anonymously [14] and later extended by several studies [15,16,17,18]. These early analyses investigated the flow of a temperature-independent Newtonian fluid and provided exact analytical solutions for the drag- and pressure flows in the cross- and down-channel directions.
Taking the shear-thinning flow behavior of polymer melt into account increases the complexity substantially. Due to the dependency of viscosity on shear rate, the flow components are coupled, and even for a one-dimensional flow of a power-law fluid, exact closed-form analytical solutions remain elusive. To extend the flow analyses to shear-thinning fluids, numerical solutions have been presented for one- and two-dimensional flows [19,20,21,22]. The analysis gains an additional degree of complexity when transverse flow is included, which contributes to the magnitude of shear rate and to the viscosity of the polymer melt, thus indirectly affecting the flow rate. This effect is particularly pronounced in the case of wave-dispersion screws, where cross-channel mixing is systematically increased by the geometry of the wave channels.
With the advent of more advanced computers, numerous works [23,24,25,26,27] addressed full three-dimensional flows in conventional metering channels using methods based on computational fluid dynamics (CFD). By numerically solving the governing flow equations, these analyses provided detailed insights into the velocity and pressure characteristics of the recirculatory flow. Numerical simulations have also proven useful in predicting the conveying behavior of wave-dispersion screws [2,3,4]. Due to the increased complexity of the flow situation and the limited computational power at the time, only a few channel geometries and physical conditions were considered.
To estimate the effect of shear-thinning flow behavior on the pumping capability of conventional metering zones, several studies proposed analytical regression models [28,29,30,31,32,33]. Approximating the numerical results of numerous flow simulations, these melt-conveying models were designed to predict the throughput as a function of characteristic process parameters without resorting to numerical methods.

1.2. Research Approach

Due to the lack of analytical melt-conveying models available for wave-dispersion screws and to remove the need for computationally expensive and time-consuming CFD simulations, we propose a semi-numerical modeling approach to predicting the flow of shear-thinning polymer melts in wave-dispersion zones. Rather than numerically solving a full set of conservation and constitutive equations, our semi-numerical method iteratively solves a linearized set of network equations and is therefore considerably faster. Figure 1 illustrates the main characteristics of the research approach, which was split into two parts (Part A and Part B).
In Part A [34], we systematically reduced the complexity of the flow analysis by ignoring the influence of screw rotation on the conveying behavior of the wave zone. Focusing on three-dimensional pressure flows allowed us to develop and validate a semi-numerical modeling approach to predicting (i) the pressure distribution, (ii) the local down- and cross-channel flows, and (iii) the transverse mixing capability in double-wave-channel systems.
In this part, we incorporated the influence of screw rotation into the modeling procedure. To this end, two- and three-dimensional melt-conveying models, presented in [30,31,32], were used in the network calculation, replacing the governing pressure-flow equation in the original theory. Moreover, the structure of the flow network was modified to additionally include leakage flow over the main flights. Experimental studies were carried out using 35 mm single-screw extruders equipped with various barrels and double-wave screws. The results of the experiments were compared to those of our extended modeling approach.

2. Experimental Procedure

2.1. Materials

As in the first part of this research, two materials were analyzed: (i) a high-density polyethylene (melt-flow rate of 0.24 g/10 min measured at 190 °C/5 kg) and (ii) a polypropylene random copolymer (melt-flow rate of 8.0 g/10 min measured at 230 °C/2.16 kg). The first material (HDPE) is used in the manufacturing of pipes, whereas the second (PP-R) is applied in industrial film extrusion. The rheological flow behavior of the polymer melts was approximated by a Carreau-Yasuda model:
η c = a t η + a t ( η 0 η ) ( 1 + ( a t   λ   γ ˙ ) a ) n c 1 a
where η 0 and η are the zero- and infinite-shear viscosity, respectively, λ is the characteristic relaxation time, and n c is the Carreau-Yasuda power-law index. The temperature-shift factor was evaluated by an approximated Arrhenius model:
a t = e x p ( α ( T T 0 ) )
where α is the temperature coefficient and T 0 the reference temperature. Table 1 lists the Carreau-Yasuda parameters and the melt densities of the materials. For detailed information on data measurement, the reader is referred to Part A [34].

2.2. Equipment

Experimental studies were performed on a 35 mm SML single-screw extruder (SML, Lenzing, Austria). A schematic of the test setup is shown in Figure 2. The electrical drive unit provided a maximum mechanical power of 38.4 kW and a maximum screw speed of 434 rpm. At the extruder head, a throttle block equipped with an adjustable screw was installed to set the discharge pressure.
Along the extruder barrel, thermal energy was supplied by electrical heaters grouped into three heating zones: T1 to T3. Further, the temperature of the downstream threaded ring and the throttle block was controlled by two heating zones: T4 to T5. A water-cooled feed housing was used to avoid heat generation in the feeding zone, while the remaining barrel sections were cooled by forced air. Table 2 shows the temperature profiles used in the experimental part.
To allow various specific rates, we used three extruder barrels of diameter Db = 35 mm with different designs in the feeding zone of the screw. One of these was produced with an entirely smooth inner surface (barrel 1), whereas the others were constructed with short helical and axial grooves in the feeding zone (barrels 2 and 3). Table 3 lists characteristic geometrical properties of the barrels.
For each barrel, we investigated two double-wave screws of outer diameter DS = 34.85 mm and axial length L = 1090 mm. In total, these involved two screws (screws 1 and 2) for the smooth-bore and two screws (screws 3 and 4) for the grooved-barrel configurations. To adapt the screw designs to the diverse conveying characteristics of conventional and feed-controlled extrusion processes, different screw geometries were used in each case. For each barrel system, however, the screws were identical except for the geometry of the wave zone. Figure 3 illustrates the channel-depth profiles of the screws.
For barrel 1, we used two double-wave screws with a feeding and a metering zone of channel depth hf = 8.1 mm and hm = 3.6 mm, respectively. The beginning of the wave zone was shifted from an axial position of 6.4 Db for screw 1 to 14.2 Db for screw 2, while the end of each wave zone was located roughly at the same axial position. The oscillating channel-depth profiles were constructed such that the numbers of wave peaks and crests were equal. The distance between the channel depth maxima, however, was decreased in the second case.
For barrels 2 and 3, the feeding zone of the double-wave screws was constructed with a compression from hf,1 = 4.9 mm to hf,2 = 6.2 mm, and the wave zone was extended to the end of the screw. In contrast to screws 1 and 2, the wave profiles showed a fully alternating behavior over the entire length of the zone. By changing the frequencies of the channel depth maxima, the number of wave peaks and valleys was almost doubled for screw 4.
For all double-wave screws, a Maddock shearing element with axial length L = 165 mm was employed. Table 4 compares geometrical parameters of the double-wave zones analyzed.
We positioned six piezoresistive pressure transducers along the barrel to measure the pressure profile along the screw. Due to the individual designs of the screws, the wave zones included four pressure sensors ( p 1 to p 4 ) in the cases of screws 1 and 2, and five pressure sensors ( p 1 to p 5 ) in the cases of screws 3 and 4 (Figure 3). The last pressure transducer p 6 was positioned at the screw tip to measure the discharge pressure. Axial positions, measuring ranges, and response times of the pressure sensors are summarized in Table 5. Further, the temperature of the polymer melt at the discharge end of the extruder was evaluated by a melt-temperature sensor.

2.3. Procedure

For each screw-barrel configuration and material, experiments were performed by increasing the screw speed from 25 rpm up to 400 rpm. Due to mechanical limitations of the barrels, we were not able to reach the maximum screw speed in all experimental setups. With the single-screw extruder operating at steady state, the throughput, the axial pressure profile, and the melt temperature were measured using a data acquisition unit. For each operating point, the processing parameters were recorded for one minute. Mean and standard deviation were evaluated to determine parameter fluctuations and the stability of the process. Table A1 in the Appendix A summarizes the experimental data for each operating point.
At the beginning of each test run, the discharge pressure was set to 250 bar for the HDPE and to 150 bar for the PP-R by changing the position of the screw mounted within the throttle block. These values represent characteristic settings found in pipe and film extrusion. The discharge pressure was adjusted once before each test at predefined screw speeds. During the experimental procedure, the position of the throttle screw was kept constant.
In addition, we carried out solidification experiments for selected processing conditions using the HDPE to analyze the melting behavior of the screws. To this end, carbon black masterbatch was added to the white polymeric material, the extrusion process was run in stationary mode and then stopped abruptly, cooling down the extruder and thus solidifying the polymer melt. After pulling the screw from the processing unit, thin sections of the solidified material were removed from the helical ribbon. The solids and melt content along the wave zones were evaluated by image analysis. The white and black sections indicated solid material and polymer melt mixed with the masterbatch, respectively (Figure 4). Taking several samples from various axial positions, we assessed the melting progress by relating the area of the black spots to the total area of the cross section.

3. Modeling

3.1. Melt-Conveying Models

A novel feature of the theory presented here is the integration of two- and three-dimensional melt-conveying models (proposed in [30,31,32]) into the network calculation to describe the local conveying characteristics of the wave zone.
Replacing the governing pressure-flow equation in the original theory, these models predict fully developed isothermal flows of wall-adhering power-law fluids in two- and three-dimensional screw channels, as shown in Figure 5. The former screw channel is infinitely wide and the corresponding flow field shows velocity components in the cross- and down-channel directions. The latter screw channel includes the influence of the screw flights, thereby producing velocity components in all directions of the coordinate system.
The melt-conveying models were developed by means of a hybrid modeling procedure that incorporates analytical, numerical, and data-based modeling into one approach. Detailed information is given in [33]. For each flow situation, the governing flow equations were rewritten in dimensionless form to determine the independent dimensionless input parameters. These quantities were then varied to create a large set of physically independent design points, whose target variables were evaluated numerically. The numerical results of the parametric design studies were then approximated using symbolic regression based on genetic programming. The hybrid modeling approach yielded two analytical regression models that predict a dimensionless volume flow rate Π V as a function of the corresponding sets of independent dimensionless input parameters:
2 D :           Π V = f ( t D b ,                 n ,   Π p , z )
3 D :           Π V = f ( h w ,                 t D b ,   n ,   Π p , z )
The two-dimensional melt-conveying model has three independent dimensionless input parameters: (i) the screw pitch ratio t / D b , (ii) the power-law index n , and (iii) a dimensionless pressure gradient Π p , z . Taking the influence of the screw flights into account, the three-dimensional model additionally requires (iv) the aspect ratio of the screw channel h/w to be predefined. A comparison of the models is shown in Figure 6, which illustrates two- and three-dimensional screw characteristic curves for the wave zone of screw 4 (with t/Db = 1.86) and various power-law indices. To illustrate the effect of the screw flights on the flow, the aspect ratio in the three-dimensional model was set to h/w = 0.07 and h/w = 0.35, representing the screw channels at a peak and a valley, respectively.
Analyzing the conveying behavior at the wave peak (Figure 6a), the two- and three-dimensional screw characteristic curves show a similar behavior due to the limited influence of the screw flights. In contrast, a significant difference is evident at the wave valley (Figure 6b), where the screw flights reduce the flow rate for a variety of pressure gradients. The impact of the screw flights is more pronounced for negative dimensionless pressure gradients as found in overridden melt-conveying zones. A drawback of the three-dimensional model is that the screw flights are represented as being equally large on both sides of the channel. In wave-dispersion screws, the barrier flight is undercut relative to the main flight. This effect can be of significant size in the case of energy-transfer screws, in which the functions of the main and the secondary flight alternate along the functional zone. To consider the different clearances of the active and passive flights and their diverse impacts on the flow, the following correction factor is introduced:
f * = 2 h δ f δ b 2 h
where δ f and δ b are the clearances of the main and the barrier flight, respectively. The correction factor relates the channel depth covered by the flights to the total possible channel depth of 2 h. In the following screw calculation routine, it is used to weight the impacts of the two- and three-dimensional melt-conveying models in the output prediction:
Π V = f * Π V , 3 d + ( 1 f * ) Π V , 2 d

3.2. Network Analysis

In Part A [34], we developed a calculation routine based on network theory to analyze pressure-driven flows in double-wave zones. Omitting the influence of the drag flow, the theory predicts the pressure distribution, the local down- and cross-channel flows, and the transverse mixing capability. This section revisits the fundamentals of the modeling approach and highlights the novel features incorporated to include the drag force of the screw. The usefulness of network theory in the flow analysis of extrusion equipment has been demonstrated by several studies [5,35,36,37,38,39,40].
The main idea of our semi-numerical approach is as follows: To simplify the complex flow in wave-dispersion zones, the screw channel is subdivided into very small segments with constant geometrical and physical parameters. These sections are indicated by network elements, each of which consists of a source and a resistance connected in parallel that represent the local drag and pressure flows, respectively. With these properties, the network elements are designed to describe the local output-pressure gradient relationship. Figure 7 illustrates the basic structure of a network element, where m ˙ d is the drag flow, k the conductance, and p i n   and p o u t are the pressures at the surrounding nodes.
The mass flow rate of an element is given by a linear combination of the drag and pressure flow. The latter is defined by the conductance and the pressure difference between the nodes:
m ˙ = m ˙ d + m ˙ p = m ˙ d + k ( p i n p o u t )
In Part A [34], we linearized a nonlinear die-flow equation for power-law fluids to determine the properties (theoretical drag flow and conductance) of each network element. The modeling approach presented here considers the nonlinear conveying behavior of the wave channel by linearizing our melt-conveying models Π V = f ( t / D b ,   n ,   Π p , z ) and Π V = f ( h / w ,   t / D b ,   n ,   Π p , z ) with respect to the current operating point. For both models, (theoretical) drag flow and conductance are obtained from the initial value and the slope of the linearization, respectively (Figure 8).
The final element properties are then evaluated by weighting the results according to our correction factor in Equation (5). Note that the theoretical drag flow used in the network calculation deviates from the physical drag flow.
Analogously to electrical circuits, the network elements are connected via nodal points. The axial and down-channel positions of the nodal points result from the screw geometry and a predefined discretization distance. Figure 9 illustrates the flow network of an unwound double-wave zone consisting of two subchannels (channels 1 and 2), two main flights, and a barrier flight. We used down-channel elements to describe the flow along the wave channels and cross-channel elements to include transverse flow. Compared to the theory presented in Part A [34], the flow network was extended to consider both transverse mixing over the barrier flight and leakage flow over the main flight.
As in the previous part [34], the cross-channel network elements were initialized with three elements connected in series: (i) one element from the channel center to the pushing flight, (ii) one element over the flight, and (iii) one element from the trailing flight to the center of the channel one turn behind. This procedure allowed us to accurately model the stepwise changes in channel height between the subchannels. The elements were then replaced by an equivalent element. To identify the pairs of nodal points connected in the cross-channel direction, we introduced a parameter counting the number of down-channel elements between these nodes:
N o f f = r o u n d (   D b   π c o s ( φ b ) ( w 1 + w 2 + e a + e b )   t a n ( φ b ) D b   π N z   c o s ( φ b ) )
where w 1 , w 2 , e a , and e b are the widths of the first subchannel, the second subchannel, the main flight, and the barrier flight, respectively. The parameter additionally depends on the resolution of the network determined by the number of elements per revolution N z . In contrast to the down-channel flow, transverse flow is described by a Newtonian output-pressure-gradient relationship based on representative viscosities. The flow in the cross-channel direction is in fact two-dimensional since the rate-limiting influence of the screw flight can be ignored. Note that our two-dimensional melt-conveying model is only valid if 0.6 < t / D b < 2.4 or 10.8 ° < φ b < 37.4 ° [30]. In the transverse direction with 90 ° φ b , these conditions are not fulfilled as the cross-channel component dominates the flow. In order to consider the specific flow conditions over the screw flight, new melt-conveying models for shear-thinning polymer melts are currently under development.
The equivalent circuit diagram for an arbitrary nodal point with index i is shown in Figure 10. The indices “1” and “2” describe the first and the second subchannel, while the indices “a” and “b” refer to the main and the barrier flight.
For solving the flow network, nodal analysis is applied. For each non-reference node the sum of the incoming flows must equal the sum of outgoing flows (cf. Kirchhoff’s current law):
i m ˙ i = 0
The network equations for an arbitrary nodal point with index i are given by:
Channel 1:
k 1 ( i 1 )   p 1 ( i 1 ) + [ k 1 ( i 1 ) k 1 ( i ) k b ( i ) k a ( i + N o f f ) ]   p 1 ( i ) + k b ( i )   p 2 ( i ) + k 1 ( i )   p 1 ( i + 1 ) + k a ( i + N o f f )   p 2 ( i + N o f f ) = m ˙ d , 1 ( i 1 ) + m ˙ d , 1 ( i ) + m ˙ d , b ( i ) m ˙ d , a ( i + N o f f )
Channel 2:
k a ( i )   p 1 ( i N o f f ) + k 2 ( i 1 )   p 2 ( i 1 ) + k b ( i )   p 1 ( i ) + [ k 2 ( i 1 ) k b ( i ) k 2 ( i ) k a ( i ) ]   p 2 ( i ) + k 2 ( i )   p 2 ( i + 1 ) = m ˙ d , 2 ( i 1 ) m ˙ d , b ( i ) + m ˙ d , 2 ( i ) + m ˙ d , a ( i )
Taking all nodes of the flow network into account yields a system of linear equations that can be summarized in matrix form:
m ˙ d + K · p = m ˙
where m ˙ d is the drag flow vector, K the conductance matrix, and p the pressure vector, and m ˙ includes the boundary conditions. The solution of this equation system yields the nodal pressures:
p = K 1 · ( m ˙ m ˙ d )
Particular attention must be paid to the network equations for the first couple of nodes in each subchannel. For i = 1 the inlet mass flow rates in both channels m ˙ 1 , 0 and m ˙ 2 , 0 are known, while for i < N o f f there are no cross-channel connections over the main flight to channel 1. In addition, the equations must be modified for the last several nodes. For i = n + 1 the outlet pressures of both channels are given by p 1 ( n + 1 ) = p 1 , o u t and p 2 ( n + 1 ) = p 2 , o u t , while for i < n + 1 N o f f there are no cross-channel connections over the main flight to channel 2. The inlet mass-flow rates are given by:
m ˙ 1 , 0 =   m ˙ 0   w 1 ( 1 ) w 1 ( 1 ) + w 2 ( 1 )                 and                 m ˙ 2 , 0 =   m ˙ 0   w 2 ( 1 ) w 1 ( 1 ) + w 2 ( 1 )
where m ˙ 0 is the output of the single-screw extruder.

3.3. Screw Simulation

Our semi-numerical simulation routine is based on the procedure shown in Figure 11. A detailed discussion was given in Part A [34].
At the beginning, the main simulation settings were configured. These include: (i) input of screw geometry, material properties, and processing conditions (screw speed, mass flow rate, melt temperature, and reference pressure), (ii) selection of material models, and (iii) definition of parameters required for the solving process. Using the data obtained from the experimental tests, the inlet mass flow rates were defined according to Equation (14). For each wave section, the last experimental pressure value ( p 4 or p 5 ) was used to define a reference pressure along the wave zone. Our simulations considered an isothermal flow of an incompressible polymer melt. For this reason, the melt temperature was used to shift the viscosity data to the desired temperature level, while the density was constant. For all calculations, the number of network elements per revolution was set to N z = 200 .
In the next step, the screw channel was discretized into a network of small segments, and the calculation was initialized. The iteration scheme started by determining the properties of the network elements. At each nodal point, the geometrical parameters were calculated and assigned to the connected network elements. To evaluate the rheological element properties, the Carreau-Yasuda parameters were converted into equivalent power-law parameters for the local shear rate. For detailed information, the reader is referred to Part A [34].
In the subsequent step, linearized element properties were derived ( m ˙ d and k ), and the network equations for each element were built. The resulting system of linear equations was then solved, and the calculated pressure field was used to update the element flow rates for the next iteration. A simulation was considered converged if the pressure differences between the first and final nodes in each subchannel was smaller than p < 0.01 bar between two iteration loops.
With the simulation procedure having converged, the modeling procedure provided (i) axial pressure profiles along the subchannels and (ii) mass flow rates in each subchannel and over the screw flights.

4. Results and Discussion

4.1. Comparison of Experimental and Calculated Results

The semi-numerical modeling approach presented in this article was used to reproduce the pressure profiles measured in the experimental part. For each operating point, we predefined the screw speed, the output, and the melt temperature according to our experimental results and calculated the pressure profiles along the subchannels of the wave zone.
For convenience, these results were averaged to obtain mean pressure curves for the double wave-channel system. Pressures were evaluated over the entire length of the wave section, except for screw 1, where the wave zone started at an axial position of 6.4 Db. Since no experimental pressures were measured in the early sections of the screw, we started the calculation at 12.5 Db.
Figure 12 shows axial pressure profiles for barrel 1 (smooth-bore extruder) and screw 1 for various operating conditions. The calculated pressure profiles are in excellent agreement with the measured data, taking the averaging of the calculated results and the standard deviation of the experimental data into account. Especially in the early screw sections, this result seems surprising, as the melt-conveying models applied in our semi-numerical modeling approach assume the polymeric material to be entirely in a liquid state, and the influence of solid resin particles on the conveying characteristics is ignored. The melting profiles in Figure 13 clearly indicate that for the screw-barrel system under investigation (barrel 1 and screw 1) the melt content is in the range of 50% at the beginning of the wave zone and increases steadily over the length of the functional zone.
The high accuracy of the melt-conveying models in the early screw sections may be a result of the specific way in which wave-dispersion screws melt dispersed solids, which gives rise to particle clusters and individual particles trapped within the polymer melt (Figure 14). Since melt films are generated on each side of the screw channels, the viscous-drag conveying behavior replaces frictional-drag conveying behavior at an early stage. This change is an important prerequisite for the validity of our melt-conveying models. The melt films are subjected to high shear rates, thus greatly affecting the viscosity and the pressure behavior of the system. The solid bulk material in the channel centers, in contrast, shows minor viscosity differences and has therefore less influence on the pressure characteristics. This behavior is also represented by our melt-conveying models. In a fully melt-dominant description, however, the effect is explained not by the presence of solids, but by the reduced shear rates at the channel centers (especially in overridden functional zones), which cause the viscosities to reach the zero-viscosity plateau. The cross sections obtained from our solidification experiments are shown in Figure A1 in the Appendix A. For the smooth-bore system, melting completed after roughly 25 Db, whereas the helically grooved barrel showed a melt content of approximately 70% at the end of the wave zone. In the latter case, the shearing element converted the remaining solids.
Characteristic of smooth-bore extruders, the wave-dispersion zones analyzed in Figure 12 generate pressure for all setups, as indicated by a positive pressure gradient. With increasing screw speed, the slope of the pressure profiles, and hence the pressure build-up capacity, increases. The latter is affected by the geometry of the wave channels and the cross-channel flows induced by the drag force of the screw and the transverse pressure gradients.
Compared to the conventional extruder (barrel 1), the double-wave screws used in combination with the helically grooved and the axially grooved barrels (barrels 2 and 3) operate balanced or overridden with respect to the governing pressure characteristics. For these systems, the extrusion process is feed-controlled, and pressure development occurs early in the solids-conveying zone, while the downstream functional zones work pressure-neutral or consume pressure. Figure 15 and Figure 16 illustrate axial pressure profiles for barrels 2 and 3, respectively. For both configurations, pressure profiles are shown for screws 3 and 4 in which the wave sections extend from 11.6 Db to nearly the end of the screws. Similarly, the diagrams compare the solutions of our semi-numerical modeling approach and the experimental results ( p 1 to p 5 ) for various setups.
As expected, by constructing the barrels with grooves, the output was more than doubled from barrel 1 to 2, and again significantly raised from barrel 2 to 3. In addition, a substantial increase in pressure from 200 to 800 bar is evident. Our main intention in using different barrels was to allow the processing unit to operate at various specific output rates to increase the process windows for model validation. For all setups, the calculated pressures accurately reflect the experimental pressure characteristics. The less viscous PP-R shows nearly constant pressures over the entire length of the wave zone, whereas the more viscous HDPE is subjected to a negative pressure gradient in both barrel systems.
For all experimental setups, pressure fluctuations caused by the presence of solid resin particles were within an acceptable range, as shown by the standard deviations obtained for the experimental data.
Using our two- and three-dimensional melt-conveying models in combination with network theory enables fast prediction of the conveying characteristics of double-wave screws. The modeling approach being based on a fully melt-dominant description, the method applies several modeling assumptions that can further be optimized. To increase the accuracy of the approach, the presence of solid material and phase transition effects can be taken into account. The extrusion literature provides several solids-conveying and melting models that may be implemented into the network calculation. A full description of the transport phenomena, however, was not the objective of this research.
The presented simulation routine uses two- and three-dimensional melt-conveying models to analyze local conveying characteristics of the wave zone. These models consider flows of shear-thinning polymer melts in unwound screw channels and ignore the influence of channel curvature. While the flat-plate assumption is widely used in extrusion theory and has often proven useful in screw analyses, the applicability for wave-dispersion zones remains questionable. Especially in the deep valley portions of the wave zones, the channel curvature reaches a critical level with the ratio of channel depth to screw diameter significantly exceeding the values typically found in conventional metering zones. In a recent article [41], we showed that for deep screw channels the flat-plate model underestimates the flow rate for a variety of processing conditions.
Another simplification made in the development of the melt-conveying models is the assumption of the flow as isothermal. Due to the temperature-dependent viscosity of the polymer melts, the velocity and temperature fields in the screw channel are in fact coupled. Roland et al. [40] presented an extended modeling approach based on network theory to predicting the non-isothermal conveying characteristics of conventional metering zones.

4.2. Distribution of Mass-Flow Rates and Pressure

This section presents the results of applying our network-based simulation routine in order to analyze the mass flow rates and pressures in the subchannels of the double-wave zones. Figure 17 and Figure 18 illustrate the behavior for a conventional (barrel 1 and screw 2) and a feed-controlled extrusion process (barrel 2 and screw 4).
The distribution of the mass-flow rates is governed primarily by the oscillating channel-depth profiles of the wave channels. Material flowing down a channel toward a peak is forced to split flow due the diminishing cross-sectional area. Some material portions remain in the original channel, and others travel across the barrier flight to increase the mass flow rate in the adjacent channel. This process is repeated multiple times throughout the wave-dispersion zone.
Taking the channel-depth profiles in Figure 3 into account, the local mass flow rates reach a maximum when the channel under investigation decreases in cross-sectional area, while the adjacent channel increases in cross-sectional area. Depending on the ascending or descending nature of the channel depth, extreme values are found in the region where the cross-sectional areas in the subchannels are equally large. The exact position depends on the interaction of the flow channels, the transverse pressure gradient and the drag force of the screw. Once this position has been reached, the direction of cross-channel flow changes.
For both extrusion processes, the pressures in channel 2 are generally higher than in channel 1. This is caused by the different clearances of the main and the barrier flights. At the pushing side of channel 2, transverse flow is largely restricted by the small clearance of the main flight. In contrast, the increased clearance of the barrier flight allows the material to transfer between the channels. For a majority of operating conditions, the standard deviations of the pressures evaluated in the experimental part exceeded slightly the calculated pressure differences between the subchannels.
Due to compression of the wave zones in the axial direction, cross-channel mixing improves with increasing screw length, as indicated by the rising amplitudes in the mass-flow distributions. To evaluate the transverse mixing performance, the following mixing index is introduced:
κ = i | m ˙ f , i | m ˙ 0
where the numerator represents the sum of the cross-channel element flow rates (over the main and the barrier flights) and the denominator the total output. The higher the mixing parameter, the more pronounced the transverse flow and thus cross-channel mixing.
Figure 19 compares the influence of screw speed and screw design on the transverse mixing performance. For a given experimental setup (screw-barrel-material combination), cross-channel mixing improves only slightly with increasing screw speed. As the sum of cross-channel flow rates and the total throughput increase nearly proportionally if the screw speed is raised (setups 1 and 2), a similar mixing index is obtained in both cases, indicating that the total throughput passes across the screw flights more than twice along the wave zone. In contrast, a significant difference in mixing performance is evident if the screw design is modified. For a given operating point, the mixing index is almost doubled if the number of peaks and valleys along the wave zone is increased (setups 2 and 3). Further comparisons were carried out to investigate the influence of the viscosity behavior. Again, we observed no significant differences in transverse mixing between the lower and the higher viscous material.

5. Conclusions

We have presented and validated a semi-numerical modeling approach to analyzing flows of shear-thinning polymer melts in flow systems with alternating wave channels. In the first part (Part A), we developed the mathematical groundwork for the theory presented here by focusing exclusively on pressure-driven flows. Reducing the complexity of the mathematical analysis in Part A allowed us to propose a network-based modeling approach to predicting: (i) the pressure distribution, (ii) the local down- and cross-channel flows, and (iii) the transverse mixing capability of double wave-channel systems. Rather than solving the full set of conservation equations, the approach iteratively solves a system of linear equations. For a number of experimental setups, the semi-numerical modeling approach was shown to yield accurate results at low computational cost.
While the method proposed in Part A can be used to model flows in various types of extrusion dies with changing channel geometry, the research presented here incorporated the influence of the screw rotation into the network-based calculation. Two significant modifications were made to extend the original theory.
First, to describe the local conveying characteristics of the wave zone, two- and three-dimensional melt-conveying models for shear-thinning polymer melts were combined, using a novel correction factor. For each subchannel, the parameter relates the distance covered by the flight flanks on both sides of the channel to the total possible channel depth, thereby providing a measure for determining the impact of the screw flights in the output calculation. In wave-dispersion screws, the barrier flight is typically undercut with respect to the main flight to allow the material to transfer between the subchannels. Assuming the polymer melt to stick to the walls of the screw channel, the rate-limiting influence of the secondary flight is therefore reduced compared to that of the main flight.
Second, we adapted the structure of the flow network to include both transverse flow over the barrier flight and leakage flow over the main flight. A Newtonian output-pressure-gradient relationship based on representative viscosities was used to describe transverse flow. The cross-channel elements in the flow network being positioned perpendicular to the screw channel, new melt-conveying models for shear-thinning polymer melts are currently under development that consider the specific flow conditions over the screw flights.
The validity of the network-based modeling approach was confirmed experimentally for both conventional and feed-controlled extrusion processes by analyzing a number of double-wave screws and processing conditions. Our main intention in using extruder barrels of different designs was to increase the parameter space available for model validation. The results of our semi-numerical modeling approach for the axial pressure profiles along the wave zones accurately reproduced the experimental data. Even in the screw sections with increased quantities of solid particles, our melt-dominant analysis provided highly satisfactory results.
Our network-theory-based simulation routine enables fast prediction of the conveying behavior of double-wave screws. Taking the influence of the shear-thinning flow behavior of polymer melts into account, the modeling approach permits a simple description of the local flow phenomena by a linear superposition of a drag and a pressure flow. Due to its efficient solving process, the modeling approach is expected to be particularly useful in (i) design and optimization studies and (ii) process troubleshooting of wave-dispersion screws. In both cases, a fast and accurate analysis of the screw concept is crucial.

Author Contributions

Conceptualization, C.M., W.R., M.D., V.S. and G.S.; Methodology, C.M., W.R. and S.S.; Software, C.M. and W.R.; Validation, C.M., W.R. and S.S.; Formal Analysis, C.M. and W.R.; Investigation, C.M., W.R. and S.S.; Resources, C.M., W.R., M.D., V.S. and G.S.; Data Curation, C.M.; Writing-Original Draft Preparation, C.M.; Writing-Review & Editing, W.R., M.D., V.S. and G.S.; Visualization, C.M. and W.R.; Supervision, W.R. and G.S.; Project Administration, C.M. and G.S.; Funding Acquisition, C.M. and G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Austrian Science Fund (FWF): I 3327-N34. Further, this work was financially supported by the COMET Center CHASE and Pro2Future. The COMET Centre CHASE is funded within the framework of COMET—Competence Centers for Excellent Technologies by BMVIT, BMDW, and the Federal Provinces of Upper Austria and Vienna. The COMET Centre Pro2Future is funded within the framework of COMET—Competence Centers for Excellent Technologies by BMVIT, BMDW, and the Federal Provinces of Upper Austria and Styria. The COMET program is run by FFG.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Cross sections along the wave zone obtained from the solidification experiments for two screw-barrel configurations: barrel 1 and screw 1, and barrel 2 and screw 3. Both tests were carried out at a screw speed of N = 100 rpm, which resulted in a throughput of m ˙ = 20.5 kg/h and m ˙ = 47.0 kg/h for the first and the second screw-barrel configuration, respectively.
Figure A1. Cross sections along the wave zone obtained from the solidification experiments for two screw-barrel configurations: barrel 1 and screw 1, and barrel 2 and screw 3. Both tests were carried out at a screw speed of N = 100 rpm, which resulted in a throughput of m ˙ = 20.5 kg/h and m ˙ = 47.0 kg/h for the first and the second screw-barrel configuration, respectively.
Polymers 12 01900 g0a1aPolymers 12 01900 g0a1b
Table A1. Operating points. Mean values of process parameters measured in the experimental part.
Table A1. Operating points. Mean values of process parameters measured in the experimental part.
TestBarrelScrewMaterial N m ˙ T M   p 1   p 2   p 3   p 4   p 5   p 6  
----rpmkg/h°Cbarbarbarbarbarbar
111HDPE5010.0211.875.2102.6125.2137.5-184.2
211HDPE10020.5218.161.291.5125.3150.1-213.4
311HDPE15031.0222.650.976.1112.6149.7-232.2
411HDPE20041.0226.251.071.9110.7152.1-245.8
511HDPE25049.0230.656.073.8155.0158.1-270.1
611HDPE30058.0234.462.181.6122.1165-290.5
711HDPE35067.0239.365.582.7122.0170.6-292.5
811HDPE40076.5249.765.585.5128.1180.6-297.4
911PP-R5011.5237.045.155.464.867.8-85.7
1011PP-R10023.0237.035.453.064.074.9-103.9
1111PP-R15034.0238.437.958.973.483.1-114.5
1211PP-R20045.5238.535.459.473.187.9-123.5
1311PP-R25055.5239.032.958.771.490.2-129.7
1411PP-R30065.0238.424.250.567.089.5-136.6
1511PP-R35074.5239.428.254.496.992.9-140.8
1611PP-R40083.0243.129.155.771.495.8-148.1
1712HDPE5011.5208.0122.1158.5158.0173.9-185.7
1812HDPE10021.0217.177.0126.3130.2163.2-209.4
1912HDPE15031.0223.589.4124.5129.4158.5-227.9
2012HDPE20041.0227.2103.2129.4134.2156.0-243.9
2112HDPE25048.0233.3122.8140.4117.4145.7-253.7
2212HDPE30058.0238.3135.8153.0126.7153.5-265.8
2312HDPE35067.5244.1149.6167.3136.0161.7-277.8
2412HDPE40077.0250.1163.1180.5143.2172.7-288.9
2512PP-R5012.5238.873.097.381.398.6-104.9
2612PP-R10022.5239.152.090.181.2105.3-123.6
2712PP-R15033.0242.251.491.483.3108.2-135.0
2812PP-R20044.0240.736.689.380.0110.1-144.5
2912PP-R25054.5240.437.3105.489.8120.3-154.9
3012PP-R30064.5238.645.0107.199.7125.9-164.0
3112PP-R35076.0240.248.2123.1106.7133.0-170.7
3212PP-R40087.0240.968.4127.6111.7136.2-174.1
3323HDPE2511.0196.6109.6117.5114.5102.797.5182.8
3423HDPE5023.0199.4137.8157.7142.9125.3119.3192.4
3523HDPE7536.0192.1220.1248.8222.2205.0201.0227.6
3623HDPE10047.0181.7356.4392.5359.0335.7335.0294.5
3723HDPE12559.5175.6453.4495.0459.7429.8434.0328.7
3823HDPE15072.0171.5538.4585.9546.1508.2519.2340.6
3923PP-R5020.0229.764.890.281.585.180.3113.9
4023PP-R10042.0222.2138.3179.3154.5150.1140.2144.9
4123PP-R15062.0215.2415.9470.7436.7426.2414.4326.6
4223PP-R20082.5210.6685.8747.7717.7681.1674.5280.0
4324HDPE2510.5200.9370.8360.8371.6325.4327.77.8
4424HDPE5023.0202.0499.6480.2476.7414.6408.266.3
4524HDPE7533.0199.2672.8651.9641.5564.3556.1140.2
4624PP-R259.0231.5121.9128.0132.9119.0118.087.2
4724PP-R5020.0230.5159.2156.5163.0150.7146.2114.6
4824PP-R7531.0229.6207.5204.2200.3184.0175.7133.1
4924PP-R10041.0224.9261.9262.9258.0225.8218.1156.6
5024PP-R12551.0222.0340.5343.1342.7298.4284.0192.4
5124PP-R15062.0219.8423.9422.3434.4380.8364.3223.9
5224PP-R17571.5217.4547.1547.9564.8503.3486.2265.8
5324PP-R20083.5212.5685.7679.2709.8641.8621.8312.0
5433HDPE2515.0196.4197.5174.6181.9133.1139.649.4
5533HDPE5030.0193.6298.8274.9268.4195.3209.070.9
5633HDPE7546.0186.3459.9427.0422.3336.8361.3132.7
5733HDPE10059.5173.7617.6586.2578.1472.9511.5198.4
5833PP-R5023.0225.626.835.228.017.225.20.4
5933PP-R10046.0217.2115.8127.4116.382.796.99.4
6033PP-R15069.0211.7529.9543.2532.2470.9500.4168.3
6133PP-R20093.0207.9697.2714.9702.3622.7663.0213.2
6234HDPE2514.0197.7238.0208.5207.8149.4153.447.0
6334HDPE5029.0194.0371.6319.1314.5218.4225.369.5
6434HDPE7544.0188.8541.5475.8478.1362.0373.2109.5
6534HDPE10058.0186.0697.7613.8621.3490.1501.3157.3
6634PP-R5023.0227.033.534.834.726.034.718.4
6734PP-R10046.0218.4143.5122.7128.581.993.632.2
6834PP-R15068.0216.1305.6271.1287.9222.7234.073.7
6934PP-R20093.0213.9459.8427.8452.5391.1401.2148.9
7034PP-R250100.5214.6542.8517.8550.0484.0502.3185.3
7134PP-R300104.0216.4438.4434.5472.9421.7450.7163.4
7234PP-R350110.5219.3400.9405.1434.7402.8427.8167.3
Table A2. Nomenclature.
Table A2. Nomenclature.
a width of transition N o f f number of offset elements
a t temperature-shift factor N z number of down-channel elements
D b barrel diameter p pressure
D s outer screw diameter p o u t outlet pressure
e flight width p pressure vector
e a flight width (main flight) t screw pitch
e b flight width (barrier flight) T temperature
hchannel depth T 0 reference temperature
h f channel depth (feeding zone) w channel width
h m channel depth (metering zone) w 1 channel width (subchannel 1)
i number of parallel screw flights w 2 channel width (subchannel 2)
k linearized die conductance α temperature coefficient
K conductance matrix δ b clearance of barrier flight
L axial length δ f clearance of main flight
L1axial start of wave zone γ ˙ shear rate
L2axial end of wave zone Γ melt content
m ˙ mass-flow rate η c viscosity (Carreau-Yasuda model)
m ˙ 0 mass-flow rate (at inlet) η 0 zero-shear viscosity
m ˙ mass-flow vector η infinite-shear viscosity
m ˙ d theoretical drag flow κ mixing parameter
m ˙ d drag-flow vector λ characteristic relaxation time
m ˙ f cross-channel element flow rate Π p , z dimensionless pressure gradient
m ˙ p pressure flow Π V dimensionless volume-flow rate
n c power-law index (Carreau-Yasuda model) ρ m melt density
n power-law index φ b pitch angle
N screw speed

References

  1. Market Study: Plastic Extrusion. Available online: https://www.ceresana.com (accessed on 12 May 2020).
  2. Somers, S.A.; Spalding, M.A.; Dooley, J.; Hyun, K.S. Numerical Analysis of the Thermal Mixing Effects of an Energy Transfer (ET) Screw Section. In Proceedings of the Annual Technical Conference of the Society of Plastics Engineers, Boston, MA, USA, 7–11 May 1995; pp. 222–229. [Google Scholar]
  3. Fan, P.; Vlachopoulos, J.; Smith, N.; Sheth, H. Computer Simulation of Melt Flow in Wave Screws. In Proceedings of the Annual Technical Conference of the Society of Plastics Engineers, Atlanta, CA, USA, 26–30 April 1998; pp. 97–101. [Google Scholar]
  4. Perdiakoulias, J.; Svabik, J.; Womer, T. 3D FEM Analysis of a Wave Type Screw Channel. In Proceedings of the Annual Technical Conference of the Society of Plastics Engineers, Chicago, CA, USA, 22–24 June 2009; pp. 1180–1189. [Google Scholar]
  5. Pape, J. Grundlagen der Prozesssimulation von Einschneckenkonzepten zur Hochleistungsplastifizierung. Ph.D. Thesis, Universität Paderborn, Paderborn, Germany, 2006. [Google Scholar]
  6. Dörner, M.; Marschik, C.; Schöppner, V.; Steinbichler, G. Development of an Analytical Model to Describe the Disperse Melting in Wave-Dispersion Screws. Polymers 2020, 12, 946. [Google Scholar] [CrossRef] [Green Version]
  7. Kruder, G.A. Extruder Screws. U.S. Patent 3,870,284, 11 March 1975. [Google Scholar]
  8. Kruder, G.A. Extrusion Apparatus and Method. U.S. Patent 4,173,417, 6 November 1979. [Google Scholar]
  9. Chung, I.C.; Barr, R.A. Energy Efficient Extruder Screw. U.S. Patent 4,405,239, 20 September 1983. [Google Scholar]
  10. Kruder, G.A. Triple Channel Wave Screw. U.S. Patent 5,219,590, 15 June 1993. [Google Scholar]
  11. Medici, J.C.; Takashima, S.T.; Calland, W.N. Multichannel Wave Screw with a Slush Channel Barrier Section. U.S. Patent 6,132,075, 17 October 2000. [Google Scholar]
  12. Barr, R.A. Extruder Screw with Improved Energy Efficient Metling. U.S. Patent 6,599,004, 29 July 2003. [Google Scholar]
  13. Womer, T.W.; Buck, E.J.; Hudak, B.J. Apparatus for Plasticating Thermoplastics. U.S. Patent 6,672,753, 6 January 2004. [Google Scholar]
  14. Anonymous. Screw Viscosity Pumps. Engineering 1922, 114, 606–607. [Google Scholar]
  15. Rowell, H.S.; Finlayson, D. Screw Viscosity Pumps. Engineering 1928, 126, 249–387. [Google Scholar]
  16. Carley, J.F.; Mallouk, R.S.; McKelvey, J.M. Simplified Flow Theory for Screw Extruder. Ind. Eng. Chem. 1953, 45, 974–978. [Google Scholar] [CrossRef]
  17. Mohr, W.D.; Saxton, R.L.; Jepson, C.H. Theory of Mixing in the Single-Screw Extruder. Ind. Eng. Chem. 1957, 49, 1857–1862. [Google Scholar] [CrossRef]
  18. Mohr, W.D.; Mallouk, R.S. Flow, Power Requirement, and Pressure Distribution of Fluid in a Screw Extruder. Ind. Eng. Chem. 1959, 51, 765–770. [Google Scholar] [CrossRef]
  19. Rotem, Z.; Shinnar, R. Non-Newtonian Flow between Parallel Boundaries in Linear Movements. Chem. Eng. Sci. 1961, 15, 130–149. [Google Scholar] [CrossRef]
  20. Griffith, R.M. Fully Developed Flow in Screw Extruders. Ind. Eng. Chem. Fundam. 1962, 1, 180–187. [Google Scholar] [CrossRef]
  21. Kroesser, F.W.; Middleman, S. The Calculation of Screw Characteristics for the Extrusion of non-Newtonian Melts. Polym. Eng. Sci. 1965, 5, 230–234. [Google Scholar] [CrossRef]
  22. Zamodits, H.J.; Pearson, J.R.A. Flow of Polymer Melts in Extruders. Part I. The Effect of Transverse Flow and of a Superpose Steady Temperature Profile. Trans. Soc. Rheol. 1969, 13, 357–385. [Google Scholar] [CrossRef]
  23. Spalding, M.A.; Dooley, J.; Hyun, K.S.; Strand, S.R. Three-Dimensional Numerical Analysis of a Single-Screw Extruder. SPE ANTEC Tech. Pap. 1993, 39, 1533–1541. [Google Scholar]
  24. Ghoreishy, M.H.R.; Razavi-Nouri, M.; Naderi, G. Finite element analysis of a thermoplastic elastomer melt flow in the metering region of a single screw extruder. Comp. Mat. Sci. 2005, 34, 389–396. [Google Scholar] [CrossRef]
  25. Rauwendaal, C. Finite element studies of flow and temperature evolution in single screw extruder. Plast. Rubber Comp. 2013, 33, 390–396. [Google Scholar] [CrossRef]
  26. Kadyirov, A.; Gataullin, R.; Karaeva, J. Numerical Simulation of Polymer Solutions in a Single-Screw Extruder. Appl. Sci. 2019, 9, 5423. [Google Scholar] [CrossRef] [Green Version]
  27. Vachagina, E.K.; Kadyirov, A.; Karaeva, J. Simulation of Giesekus fluid flow in extruder using helical coordinate system. IOP Conf. Ser. Mater. Sci. Eng. 2020, 733, 012033. [Google Scholar] [CrossRef]
  28. White, J.L.; Potente, H. Screw Extrusion; Hanser Publishers: Munich, Germany, 2001; ISBN 978-3-446-19624-7. [Google Scholar]
  29. Rauwendaal, C. Polymer Extrusion, 5th ed.; Hanser Publishers: Munich, Germany, 2014; ISBN 978-3446217744. [Google Scholar]
  30. Roland, W.; Kommenda, M.; Marschik, C.; Miethlinger, J. Extended Regression Models for Predicting the Pumping Capability and Viscous Dissipation of Two-Dimensional Flows in Single-Screw Extrusion. Polymers 2019, 11, 334. [Google Scholar] [CrossRef] [Green Version]
  31. Marschik, C.; Roland, W.; Loew-Baselli, B.; Miethlinger, J. A heuristic method for modeling three-dimensional non-Newtonian flows of polymer melts in single-screw extruders. J. Non Newton. Fluid Mech. 2017, 248, 27–39. [Google Scholar] [CrossRef]
  32. Marschik, C.; Osswald, T.; Roland, W.; Loew-Baselli, B.; Miethlinger, J. A Heuristic Model for Predicting Three-Dimensional Non-Newtonian Flows in Metering Channels. In Proceedings of the ANTEC 2018—Technical Conference & Exhibition, Orlando, FL, USA, 7–10 May 2018. [Google Scholar]
  33. Marschik, C.; Roland, W.; Loew-Baselli, B.; Steinbichler, G. Application of Hybrid Modeling in Polymer Processing. In Proceedings of the ANTEC 2020—Technical Conference & Exhibition, San Aotonio, TX, USA, 30 March–2 April 2020. [Google Scholar]
  34. Marschik, C.; Dörner, M.; Roland, W.; Miethlinger, J.; Schöppner, V.; Steinbichler, G. Application of Network Analysis to Flow Systems with Alternating Wave Channels: Part A (Pressure Flows). Polymers 2019, 11, 1488. [Google Scholar] [CrossRef] [Green Version]
  35. Booy, M.L. A Network Flow Analysis of Extrusion Dies and Other Flow Systems. Polym. Eng. Sci. 1982, 22, 432–437. [Google Scholar] [CrossRef]
  36. Michaeli, W.; Pöhler, W. Die Verwendung der Netzwerktheorie bei der Werkzeugauslegung. Plaste Kautsch. 1993, 2, 61–65. [Google Scholar]
  37. Hopmann, C.; Michaeli, W. Extrusion Dies for Plastics and Rubber, 4th ed.; Hanser Publishers: Munich, Germany, 2016; ISBN 978-1-56990-623-1. [Google Scholar]
  38. Köpplmayr, T.; Miethlinger, J. Application of the Network Simulation Method to Flat Dies with Inverted Prelands. Int. Polym. Process. 2013, 28, 322–330. [Google Scholar] [CrossRef]
  39. Marschik, C.; Roland, W.; Miethlinger, J. A Network-Theory-Based Comparative Study of Melt-Conveying Models in Single-Screw Extrusion: A. Isothermal Flow. Polymers 2018, 10, 929. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Roland, W.; Marschik, C.; Hammer, A.; Steinbichler, G. Modeling The Non-Isothermal Conveying Characteristics in Single-Screw Extrusion by Application of Network Analysis. In Proceedings of the ANTEC 2020—Technical Conference & Exhibition, San Aotonio, TX, USA, 30 March–2 April 2020. [Google Scholar]
  41. Roland, W.; Marschik, C.; Loew-Baselli, B.; Miethlinger, J. The effect of channel curvature on the flow rate and viscous dissipation of power-law fluids. In Proceedings of the ANTEC 2018—Technical Conference & Exhibition, Orlando, FL, USA, 7–10 May 2018. [Google Scholar]
Figure 1. Research approach: In the first part (Part A), the pressure flow was decoupled from the drag flow. With the shear-thinning flow behavior being considered, a coupled network of down- and cross-channel flows was solved. In this part (Part B), the influence of the rotation is added.
Figure 1. Research approach: In the first part (Part A), the pressure flow was decoupled from the drag flow. With the shear-thinning flow behavior being considered, a coupled network of down- and cross-channel flows was solved. In this part (Part B), the influence of the rotation is added.
Polymers 12 01900 g001
Figure 2. Schematic diagram of the 35 mm plasticating single-screw extruder.
Figure 2. Schematic diagram of the 35 mm plasticating single-screw extruder.
Polymers 12 01900 g002
Figure 3. Comparison of channel-depth profiles: screw 1 (a), screw 2 (b), screw 3 (c), and screw 4 (d). The dotted vertical lines indicate the beginning and the end of the wave zone. The axial positions of the pressure transducers are shown at the upper x-axis.
Figure 3. Comparison of channel-depth profiles: screw 1 (a), screw 2 (b), screw 3 (c), and screw 4 (d). The dotted vertical lines indicate the beginning and the end of the wave zone. The axial positions of the pressure transducers are shown at the upper x-axis.
Polymers 12 01900 g003
Figure 4. Solidification experiment: The solids and melt content along the wave zone, indicated by the white and black sections, respectively, were determined by image analysis.
Figure 4. Solidification experiment: The solids and melt content along the wave zone, indicated by the white and black sections, respectively, were determined by image analysis.
Polymers 12 01900 g004
Figure 5. Representation of the screw channels in the two- and three-dimensional melt-conveying models. Ignoring the curvature of the screw, the flat-plate model with moving barrel was applied in both cases.
Figure 5. Representation of the screw channels in the two- and three-dimensional melt-conveying models. Ignoring the curvature of the screw, the flat-plate model with moving barrel was applied in both cases.
Polymers 12 01900 g005
Figure 6. Comparison between the two- and three-dimensional melt-conveying models: Screw characteristic curves for the wave zone of screw 4 with t/Db = 1.86, showing the pumping capability of a wave peak (a) and valley (b).
Figure 6. Comparison between the two- and three-dimensional melt-conveying models: Screw characteristic curves for the wave zone of screw 4 with t/Db = 1.86, showing the pumping capability of a wave peak (a) and valley (b).
Polymers 12 01900 g006
Figure 7. Basic network element for melt conveying.
Figure 7. Basic network element for melt conveying.
Polymers 12 01900 g007
Figure 8. Linearization of the nonlinear screw characteristic curve at operating point ( Π p , z | Π V ) .
Figure 8. Linearization of the nonlinear screw characteristic curve at operating point ( Π p , z | Π V ) .
Polymers 12 01900 g008
Figure 9. Network discretization of an unwound double-wave zone consisting of (i) two subchannels, (ii) two main flights, and (iii) a barrier flight. The network elements in the down- and cross-channel directions are connected via nodal points.
Figure 9. Network discretization of an unwound double-wave zone consisting of (i) two subchannels, (ii) two main flights, and (iii) a barrier flight. The network elements in the down- and cross-channel directions are connected via nodal points.
Polymers 12 01900 g009
Figure 10. Equivalent circuit diagram.
Figure 10. Equivalent circuit diagram.
Polymers 12 01900 g010
Figure 11. Flow chart of our semi-numerical modeling approach based on network theory.
Figure 11. Flow chart of our semi-numerical modeling approach based on network theory.
Polymers 12 01900 g011
Figure 12. Comparison between experimental data and calculated solutions for barrel 1 and screw 1. Axial pressure profiles along the wave zone for PP-R (ad) and HDPE (e,f). The continuous lines indicate the solutions of our semi-numerical modeling approach, while the circles illustrate experimental data.
Figure 12. Comparison between experimental data and calculated solutions for barrel 1 and screw 1. Axial pressure profiles along the wave zone for PP-R (ad) and HDPE (e,f). The continuous lines indicate the solutions of our semi-numerical modeling approach, while the circles illustrate experimental data.
Polymers 12 01900 g012
Figure 13. Melting profiles for HDPE for the smooth-bore extruder (barrel 1) (a) and the helically grooved-feed extruder (barrel 2) (b).
Figure 13. Melting profiles for HDPE for the smooth-bore extruder (barrel 1) (a) and the helically grooved-feed extruder (barrel 2) (b).
Polymers 12 01900 g013
Figure 14. Cross section of polymer ribbon of both channels extracted at an axial position of 15.5 L/Db during a solidification experiment with barrel 1 and screw 1. The white and black sections indicate solids and polymer melt, respectively.
Figure 14. Cross section of polymer ribbon of both channels extracted at an axial position of 15.5 L/Db during a solidification experiment with barrel 1 and screw 1. The white and black sections indicate solids and polymer melt, respectively.
Polymers 12 01900 g014
Figure 15. Comparison between experimental and calculated data for barrel 2. Axial pressure profiles along the wave zone for PP-R and screw 3 (a) and HDPE and screw 4 (b).
Figure 15. Comparison between experimental and calculated data for barrel 2. Axial pressure profiles along the wave zone for PP-R and screw 3 (a) and HDPE and screw 4 (b).
Polymers 12 01900 g015
Figure 16. Comparison between experimental and calculated data for barrel 3. Axial pressure profiles along the wave zone for PP-R and screw 3 (a) and HDPE and screw 4 (b).
Figure 16. Comparison between experimental and calculated data for barrel 3. Axial pressure profiles along the wave zone for PP-R and screw 3 (a) and HDPE and screw 4 (b).
Polymers 12 01900 g016
Figure 17. Distribution of mass-flow rates (a) and pressures (b) in subchannels 1 and 2 for barrel 1 (smooth-bore extruder) and screw 2. The dashed and dotted lines represent the solutions for subchannels 1 and 2, respectively.
Figure 17. Distribution of mass-flow rates (a) and pressures (b) in subchannels 1 and 2 for barrel 1 (smooth-bore extruder) and screw 2. The dashed and dotted lines represent the solutions for subchannels 1 and 2, respectively.
Polymers 12 01900 g017
Figure 18. Distribution of mass-flow rates (a) and pressures (b) in subchannels 1 and 2 for barrel 2 (helically grooved extruder) and screw 4. The dashed and dotted lines represent the solutions for subchannels 1 and 2, respectively.
Figure 18. Distribution of mass-flow rates (a) and pressures (b) in subchannels 1 and 2 for barrel 2 (helically grooved extruder) and screw 4. The dashed and dotted lines represent the solutions for subchannels 1 and 2, respectively.
Polymers 12 01900 g018
Figure 19. Comparison of mixing indices for various processing conditions: Influence of screw speed for a given screw-barrel-material configuration (a) and influence of screw design for a given operating point (b).
Figure 19. Comparison of mixing indices for various processing conditions: Influence of screw speed for a given screw-barrel-material configuration (a) and influence of screw design for a given operating point (b).
Polymers 12 01900 g019
Table 1. Material parameters ( ρ m is the melt density at 200 °C and 200 bar).
Table 1. Material parameters ( ρ m is the melt density at 200 °C and 200 bar).
ParameterUnitHDPEPP-R
η 0 Pas95,2302957
η Pas00
λ s32560.173
n c -0.2890.425
a -11
α 1/K0.01380.0162
T 0 K473.15473.15
ρ m kg/m3752730
Table 2. Temperature profiles in °C.
Table 2. Temperature profiles in °C.
Material T 1 T 2 T 3 T 4 T 5
HDPE40190200200200
PP-R40230240240240
Table 3. Comparison of barrel designs.
Table 3. Comparison of barrel designs.
NumberBarrelGroove Design
NumberWidthHeightLength *Pitch
Barrel 1Smooth----
Barrel 2Helically grooved65 mm2.8 mm140 mm105 mm
Barrel 3Axially grooved68 mm2.8 mm140 mm
* measured from the front end of the hopper.
Table 4. Comparison of double-wave zones.
Table 4. Comparison of double-wave zones.
Dimensions UnitScrew 1Screw 2Screw 3Screw 4
Outer diameterDSmm34.8534.8534.8534.85
Axial start of wave zone *L1mm225459.5407407
Axial end of wave zone *L2mm876.5880.01021.51056.5
Length of wave zone Δ L mm651.5420.5614.5649.5
Pitchtmm45.549.065.065.0
Pitch angle φ b °22.424.030.630.6
Flight widthemm3.23.53.53.5
Number of parallel flightsi-2222
Barrier undercut δ b mm0.60.61.451.45
Flight clearance δ f mm0.0750.0750.0750.075
* measured from the front end of the hopper.
Table 5. Positions, measuring ranges, and response times of pressure transducers.
Table 5. Positions, measuring ranges, and response times of pressure transducers.
Pressure SensorAxial PositionRangeResponse Time
L/Dbbarms
p 1 15.30–10008
p 2 18.30–10008
p 3 21.30–10008
p 4 25.30–10008
p 5 27.30–10008
p 6 34.70–5008

Share and Cite

MDPI and ACS Style

Marschik, C.; Roland, W.; Dörner, M.; Schaufler, S.; Schöppner, V.; Steinbichler, G. Application of Network Analysis to Flow Systems with Alternating Wave Channels: Part B. (Superimposed Drag-Pressure Flows in Extrusion). Polymers 2020, 12, 1900. https://doi.org/10.3390/polym12091900

AMA Style

Marschik C, Roland W, Dörner M, Schaufler S, Schöppner V, Steinbichler G. Application of Network Analysis to Flow Systems with Alternating Wave Channels: Part B. (Superimposed Drag-Pressure Flows in Extrusion). Polymers. 2020; 12(9):1900. https://doi.org/10.3390/polym12091900

Chicago/Turabian Style

Marschik, Christian, Wolfgang Roland, Marius Dörner, Sarah Schaufler, Volker Schöppner, and Georg Steinbichler. 2020. "Application of Network Analysis to Flow Systems with Alternating Wave Channels: Part B. (Superimposed Drag-Pressure Flows in Extrusion)" Polymers 12, no. 9: 1900. https://doi.org/10.3390/polym12091900

APA Style

Marschik, C., Roland, W., Dörner, M., Schaufler, S., Schöppner, V., & Steinbichler, G. (2020). Application of Network Analysis to Flow Systems with Alternating Wave Channels: Part B. (Superimposed Drag-Pressure Flows in Extrusion). Polymers, 12(9), 1900. https://doi.org/10.3390/polym12091900

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