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 D
b, 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 D
b to nearly the end of the screws. Similarly, the diagrams compare the solutions of our semi-numerical modeling approach and the experimental results (
to
) 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:
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.