Next Article in Journal
Equilibrium, Thermodynamic, Reuse, and Selectivity Studies for the Bioadsorption of Lanthanum onto Sericin/Alginate/Poly(vinyl alcohol) Particles
Next Article in Special Issue
The Viscoelastic Swirled Flow in the Confusor
Previous Article in Journal
Composite Aramid Membranes with High Strength and pH-Response
Previous Article in Special Issue
Rheological Properties Related to Extrusion of Polyolefins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rheological Basics for Modeling of Extrusion Process of Wood Polymer Composites

by
Krzysztof Wilczyński
1,*,
Kamila Buziak
1,
Adrian Lewandowski
1,
Andrzej Nastaj
1 and
Krzysztof J. Wilczyński
2
1
Polymer Processing Department, Faculty of Production Engineering, Warsaw University of Technology, 85 Narbutta, 02-524 Warsaw, Poland
2
Politech Ltd., 86-031 Osielsko, Poland
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(4), 622; https://doi.org/10.3390/polym13040622
Submission received: 3 February 2021 / Revised: 13 February 2021 / Accepted: 17 February 2021 / Published: 19 February 2021

Abstract

:
Wood polymer composites are materials with pseudoplastic and viscoelastic properties. They have yield stress and exhibit slip during flow. Studies on extrusion and rheology, as well as on process modeling of these highly filled materials are limited. Extensive rheological and extrusion modeling studies on the wood polymer composite based on the polypropylene matrix were performed. Viscous and slip flow properties were determined (with Rabinowitsch, Bagley, and Mooney corrections) at broad (extrusion) range of shear rate and temperature, using a high-pressure capillary rheometer. Rheological models of Klein and power-law were used for flow modeling, and Navier model was applied for slip modeling. A novel global computer model of WPC extrusion with slip effects has been developed, and process simulations were performed to compute the extrusion parameters (throughput, power consumption, pressure, temperature, etc.), and to study the effect of the material rheological characteristics on the process flow. Simulations were validated experimentally, and were discussed with respect to both rheological and process modeling aspects. It was concluded that the location of the operating point of extrusion process, which defines the thermo-mechanical process conditions, is fundamentally dependent on the rheological materials characteristics, including slip effects.

Graphical Abstract

1. Introduction

Wood polymer composites consist of thermoplastic matrix and wooden fillers, like fiber or flour. These composites are durable, stiff, and strong, and relatively cheap compared to other competitive materials. They can be applied instead of wood since they are more resistant to decay. The wood polymer composites are mainly based on the high-density polyethylene, polypropylene, and polyvinyl chloride.
The wood polymer composites market has grown significantly in the last years due to the increasing building and automotive market. Worldwide production has increased from about 2.5 million tonnes in 2012 to about 4 million tonnes in 2015. European production has reached 350,000 tonnes in 2015 [1]. In 2019, the global wood plastic composite market size was estimated at USD 5.3 billion and is expected to register a growth rate of 11.4% over the forecast period 2000–2017 [2].
The basic technology of processing of wood polymer composites is extrusion for profile production. Processing of these highly filled composites differs substantially from polymer processing which results from different thermo-rheological material characteristics, structure, etc. Limited studies were performed on extrusion and rheology of wood polymer composites. These issues were recently discussed by the authors [3,4]. A fundamental research in this area was carried out by Xiao and Tzoganakis [5,6,7,8], Li and Wolcott [9,10,11], Błędzki et al. [12], Mohanty et al. [13], Klyosov [14], Oksman Niska [15], as well as by Vlachopoulos [16,17,18], and Zolfaghari [19].
Wood polymer composites have pseudoplastic properties, and these are also viscoelastics. The viscosity decreases when shear rate and temperature are increased, and when filler content increases, the viscosity is higher. They have yield stress and exhibit slip at the walls in extrusion [5,9,16]. The slip velocity depends on the shear rate and the filler content. When the shear rate increases, the slip velocity increases, which may lead to the plug flow [17]. Increasing of filler content also promotes the plug flow [19]. Recently, extensive FEM simulation studies on extrusion with slip effects and extrusion of viscoplastics have been performed by Lewandowski and Wilczyński [20,21].
The design of the extrusion process of wood polymer composites requires a knowledge of the flow mechanism of these highly filled materials. Only a few melting observations have been reported for wood polymer composites, however. Xiao and Tzoganakis [5,6,8] intensively investigated an extrusion of a high density polyethyelene (HDPE) composite. They observed that partially filled channels were seen until the melting process was completed. However, this has not been explained by them. This phenomenon may probably result from the screw/die interactions and insufficient pressure generation in the extruder, which leads to a starvation in the screw channel. The authors also observed a slight segregation of the components that made up the HDPE/WPC blend [8]. Close to the trailing flight of the screw, mostly HDPE was seen whereas the WPC was mostly against pushing flight. Experimental data were compared to simulations of the process using a commercially available software program (unknown). It was found that predictions of the pressure profiles generated in the extruder, in general, did not agree well with those measured.
Wood polymer composites exhibit characteristic surface tearing upon exiting from an extruder die even at low shear rates, which was reported by Hristov et al. [16,17] and Santi et al. [18]. In industrial practice, die cooling is used to produce smooth profiles, and surface tearing may be also eliminated at higher shear rates and high loading of wood flour. Recently, Santi et al. [18] modeled the flow of WPC through the die trying to explain these phenomena.
Recently, Wilczynski et al. [3] have performed extensive experimental studies on single screw extrusion of a polypropylene (PP) composite. Based on these studies, the authors developed the only available global model of extrusion of wood polymer composites [4]. According to these studies, the composite solid conveying and melting were strongly dependent on the wood filler content. The well-known classical mechanism of polymer melting developed by Tadmor et al., e.g., [22], was absent in the case of composites with high filler content (over 50%). The characteristic polymer melt “pool” was not observed at a pushing flight of the screw. However, in the case of composites with low filler content (less than 50%), this mechanism was seen. According to these observations, two models of melting have been proposed, a one-dimensional model for lower content of wood filler, and a two-dimensional model for higher content, which are schematically depicted in Figure 1 and described in [3,4]. In a classical two-dimensional model, when melting progresses, the volume of the solid material decreases by decreasing the solid width, i.e., X < W, Y = H (Figure 1a), while in one-dimensional model, the volume of the solid decreases by decreasing the solid height Y < H, X = W (Figure 1b).
When modeling extrusion process, the no-slip condition in the extruders and dies is assumed, which means the flowing polymer adheres to the wall, e.g., by White [23], Tadmor [24], or Rauwendaal [25], although this is not always obvious.
Several extrusion models have been developed using this assumption, e.g., by Agur [26], Vincelette [27], Potente [28], Wilczyński [29], and recently, e.g., [30]. The review papers on modeling of the extrusion process were delivered by Ilinca [31], Altinkaynak [32], and Malik [33]. The review contributions on global modeling of the extrusion process were delivered by Teixeira [34], Wilczyński [35,36], and Hyvärinen [37].
So far, any global model of extrusion of wood polymer composites that would include slip effects is not known. In this paper, slip effects are considered and discussed both for WPC rheological studies as well as for WPC extrusion modeling. To produce global models of the extrusion process which would include slip effects, the models for both the screw and the die need to be developed, which was discussed by Malik [33], Potente [38,39,40,41], Gooneie [42], and Duretek [43].

2. Material and Rheological Studies

The wood polymer composite (PP copo inj 4, Beologic) composed of the polypropylene (PP) matrix and 50% wood filler was applied in this research. The polypropylene density was equal to ρ = 0.9 g/cm3, the polymer melt flow index was equal to MFI = 25 g/10 min (at M = 2.16 kg, and T = 230 °C), and the melting temperature was equal to Tm = 160 °C. The composite melt density was equal to ρm = 0.95 g/cm3, the composite bulk density ρb = 0.4–0.6 g/cm3, and the polymer melt flow index was equal to MFI = 4 g/10 min (at M = 5 kg, and T = 190 °C).
The WPC viscous properties have been measured at various temperatures with the use of a capillary rheometer (RG-25, Goettfert, Buchen, Germany). The capillaries of diameter D = 1, 2, 4 mm, and the length/diameter ratio L/D = 0/1, 0/2, 0/4, 10/1, 20/2, and 40/4 were used, and the Rabinowitsch correction, the Bagley correction, and the Mooney correction have been applied [44]. The measurements have been carried out at the temperatures T = 180 °C, 190 °C, and 200 °C, and at the shear rates from γ ˙ = 1 s−1 to γ ˙ = 1000 s−1. The examples of viscosity characteristics at the temperatures T = 180 °C, 190 °C, 200 °C are shown in Figure 2. An obvious pseudoplastic behavior is seen, which means the viscosity decreases with increasing the shear rate. The viscosity also decreases with increasing the temperature. However, the effect of temperature on the viscosity is not significant. The results of Mooney analysis are depicted in Figure 3, where the shear rate vs. 1/D (Figure 3a) and the slip velocity vs. shear stress (Figure 3b) are shown. The effects of the corrections of Rabinowitsch, Bagley, and Mooney on the viscosity curve are presented in Figure 4. Rabinowitsch correction means an increase of the shear rate at the constant shear stress (the constant capillary pressure drop), which results in shifting the viscosity curve shown in Figure 4a. Bagley correction means a decrease of the shear stress (the pressure decrease due to the pressure inlet loss) at the constant shear rate, which results in a shift in the viscosity curve depicted in Figure 4b. Mooney correction means a decrease of the shear rate at the constant shear stress (the constant capillary pressure drop), which results in shifting the viscosity curve shown in Figure 4c.
Two characteristic regions of slip behavior are seen in Figure 3. At the low shear stress, the slip can be characterized as weak, which is followed by a sharp increase in the slip velocity at the high shear stress, which was also observed, e.g., by Hristov et al. [16].
It is worth noting that the viscosity of composites may be evaluated based on the viscosity of the matrix and the filler content using the well-known equations of Einstein/Batchelor or Krieger/Dougherty, as it was done by Le Moigne et al. [45] and Polychronopoulos et al. [46].

3. FEM Flow Modeling

Studies on extrusion modeling with slip effects were mainly based on the works of Hatzikiriakos and Dealy [47], and Hatzikiriakos et al. [48,49,50].
For the slip analysis, a power-law model is usually applied which approximates the actual slip behavior of various fluids, including molten polymers [50].
u s = β τ w b
where us is the slip velocity, β is the slip coefficient, b is the slip power-law exponent, τw is the shear stress.
In this paper, the CFD software Polyflow [51] was used for modeling to develop screw pumping characteristics which may be implemented into the global model of the extrusion process.
The modeling procedures to build the correct geometrical model of the flow, to select the proper boundary conditions and evolution schemes, as well as the meshing concepts and procedures, were discussed in [20].
The scheme of modeling and boundary conditions are depicted in Figure 5, where;
BC1—Inflow (Qin) = (Q0), the flow rate Q0 is imposed at the inlet of the domain,
BC2—Outflow, which means that vanishing normal forces and tangential velocities are imposed at the outlet (fn & vs) = (0, 0),
BC3—Slip at the die/barrel wall, i.e., vanishing normal velocity is imposed, while the tangential force is a function of (vs − vwall), i.e., fs = f (vsmn − vwall), (vn & fs) = (0, f (vs − vwall) & vwall = 0),
BC4—Slip at the screw, i.e., vanishing normal velocity is imposed, while the tangential force is a function of (vs − vwall), i.e., fs = f (vs − vwall), (vn & fs) = (0, f (vs − vwall) & vwall = (cartesian velocities vx & vy & vz = N).
It is worth noting that the conditions BC1 and BC2 imply that the pressure may be developed in the screw/die system. However, since the pressure at the screw element end is unknown (at the die exit it is equal to zero), the pressure gradient may be calculated in this case only.
The Ostwald-de Waele model was used for process modeling with Polyflow:
η = m γ ˙ n 1
where η is the viscosity, Pa·s, γ ˙ is the shear rate, s−1, m is the coefficient of consistency, Pa·sn, n is the flow exponent.
The coefficient of consistency was equal to m = 30,100 Pa·sn, and the flow exponent was equal to n = 0.25.
The Navier law [51] was used for slip modeling with Polyflow:
τ w = F s l i p U w a l l U s U s U w a l l e s l i p 1
where τw is the shear stress, Us is the fluid velocity in the tangential direction, Uwall is the wall velocity in the tangential direction (Uwall = 0 by default), and Fslip and eslip are the model parameters (the full slip appears when Fslip = 0, the model is linear when eslip = 1).
The parameters of the Navier equation were obtained using rheometric data presented in Figure 3, which were approximated using a power-law model:
τ w = F s l i p U s e s l i p
where τw is the shear stress, N/mm2, Us is the slip velocity, mm/s, and Fslip = 0.064 and eslip = 0.43 are the model parameters.
An influence of slip parameters (Fslip = 0.064, eslip = 0.43) on the polymer flow in the die is presented in Figure 6 where die characteristics are shown. The cylindrical, three-zone die of diameter D = 45 mm at the onset and diameter D = 5 at the end was used in the study. The lengths of die sections were equal to L = 180 mm. 30 mm, 60 mm (Figure 5a). Various flow rates were used for simulations: Q = 0.0005 kg/s, 0.001 kg/s, 0.0015 kg/s, 0.0021 kg/s, 0.0048 kg/s, 0.0072 kg/s. A non-linear dependence of the flow rate on the die pressure drop is clearly seen. An example of pressure and velocity distributions is shown in Figure 7. Slipping at the die surface is very low.
An influence of slip parameters (Fslip = 0.064, eslip = 0.43) on the polymer flow in the screw is presented in Figure 8 where screw characteristics are shown. The conventional screw geometry was used in the study with the screw element of diameter D = 45 mm, the depth of the screw channel H = 3 mm, and the length L = 180 mm (Figure 5b). Various flow rates were used for simulations: Q = 0.001 kg/s, 0.0021 kg/s, 0.0048 kg/s, 0.0072 kg/s. A non-linear dependence of the flow rate on the screw pressure increase is clearly seen. When screw speed increases, the flow rate increases at the constant pressure, and the pressure increases at the constant flow rate. The drag flow can be evaluated as the flow rate at the pressure change Δp = 0. An example of pressure and velocity distributions is shown in Figure 9. Slipping is clearly seen at the screw and barrel surfaces.

4. Extrusion Modeling

Extrusion is the system in which the extruder co-operates with the die. The flow in the extruder affects the flow in the die, and the flow in the die influences the flow in the extruder.
Any changes in the flow conditions in the extruder cause changes in the flow conditions in the die, and any changes in the flow conditions in the die cause changes in the flow conditions in the extruder. The global extrusion process modeling requires taking into account this extruder/die co-operation.
The global model of the extrusion process is composed of the elementary models which describe the flow in the screw, that is in the solid conveying zone, in the delay zone, in the melting zone, and finally in the metering zone, as well as in the die zone. The calculations are carried out in the screw increments, where the parameters of the process, e.g., temperature or pressure, at the exit of the current increment are the inlet data for the next increment.
The extrusion global models use various procedures of computation which should be suitable for the modeled process, which is discussed, e.g., in [30,52,53]. In the extrusion with polymer flood feeding, the forward procedure of calculations is used. The extrusion throughput is not known here, and it is calculated by iteration searching for the extruder operating point which defines the extrusion throughput and the die pressure (the extrusion pressure). In the extrusion with polymer metered feeding, the inverse backward procedure of calculations is used, the flow rate is set here and equals the flow rate in the dosing device. The issue of global modeling was recently discussed in details by the authors [36].
As it was said, the global model of the extrusion process is composed of elementary models for the screw and the model for the die. The polymer mass flow rate is constant in the screw/die increments, and the process parameters may be assumed locally constant in these increments. This modeling approach is called the lumped parameter modeling which is accurate enough for engineering practice and applications.
The parameters of extrusion process, e.g., flow rate, temperature, power, pressure, screw filling, at the onset of the current screw increment are equal to the parameters at the exit of the former increment, i.e.,
fi−1_out (z) = fi_in (z)
fi−1_out (z + Δz) = fi_in (z + Δz)
where fi_in (z)—input data (temperature, pressure, or solid material content) at the onset of i—increment, fi−1_out (z)—output data at the exit of (i − 1)—increment, z—coordinate along the screw/die axis, Δz—screw/die increment.
The extrusion model that is discussed here is an extended version of the previous models [4], with novel concepts for slip modeling using three-dimensional, non-newtonian flow characteristics for the screw/die system.
Summarizing the model presented here is composed of the elementary models:
  • Solid conveying section;
  • delay section;
  • melting section, the model is dependent on the filler content, a one-dimensional model for less than 50% filler is used, a two-dimensional model for over 50% filler is used (Figure 1);
  • screw melt conveying section including slip effects using three-dimensional, non-Newtonian flow characteristics;
  • die flow section with slip effects using three-dimensional, non-Newtonian flow characteristics.
The calculation procedure is depicted in Figure 10. The conventional process with polymer flood feeding is considered here, and thus the forward scheme of calculation is applied. The computations begin with an initial mass flow rate. Then the process is modeled using different melting models (one-dimensional or two-dimensional). Afterwards, the flow characteristics with slip effects are used for flow modeling in the screw/die system. The exit screw pressure is checked for the convergence with the die pressure. Then, the polymer flow rate is changed (increased or decreased) according to the result of comparison, which can be positive or negative. The calculations are reiterated until the convergence of calculation has been reached.
The screw pumping characteristics were discussed in detail in [20]. In this paper, these characteristics have been approximated by the formula
Q * = a b n c Δ p *
where Q* is the flow rate in the dimensionless form, Δp* is the pressure in the dimensionless form, and a, b, c are the model parameters (a = 1.003, b = 3.195, c = 0.561).

5. Process Simulations

A classical screw with diameter D = 45 mm and length to diameter ratio (L/D) = 27 was applied. The lengths of the screw zones were equal to LF = 485 mm (the feeding zone), LC = 320 mm (the compression zone), LM = 410 mm (the metering zone), and the total screw length was equal to L = 1215 mm. The depths of the screw channel in the zones were equal to HF = 8 mm (the feeding zone), HM = 3 mm (the metering zone), and the compression ratio which is defined as the ratio of the screw channel depth in the feeding zone to the screw channel depth in the metering zone was equal to CR = HF/HM = 2.66. The pitch of the screw was equal to the screw diameter, i.e., t = D, so the helix angle φ was equal to 17.65°. The width of the screw flight was equal to e = 5 mm. The die for cylindrical rods of diameter D = 5 mm was applied in the study. Simulations were performed for various screw speeds N = 20 rpm, 50 rpm, and 80 rpm.
The rheological equation of Klein was used for the viscosity description.
ln η = a 0 + a 1 ln γ ˙ + a 11 ln 2 γ ˙ + a 12 ln γ ˙ T + a 2 T + a 22 T 2
where: η is the viscosity, Pa·s, γ ˙ is the shear rate, s−1, T is the temperature, °C, and a0, a1, a11, a12, a2, a22 are the model parameters (a0 = 40.12992893, a1 = −0.85791717, a11 = −0.01421715, a12 = 0.00160931, a2 = −0.27500818, a22 = −0.00060042).
Results of computation are shown in Figure 11, Figure 12, Figure 13 and Figure 14. The overall extrusion process characteristics are depicted in Figure 11. These dimensionless characteristics are composed of the important process parameters, i.e., the pressure, the temperature, the solid bed content, and the power.
An effect of slip on the process flow is clearly seen in Table 1 and Figure 12, Figure 13 and Figure 14. Table 1 shows that the flow rate (extrusion throughput) increases when slip appears. The pressure decreases when slip appears (Figure 12), and melting is slower in this case, which may be caused by an increase in the flow rate (Figure 13). Power consumption is a bit higher in the case when slip appears, which may also result from an increase in the flow rate (Figure 14).

6. Experimental

Validation of simulations has been performed in this paper at the process simulation data, which were presented above. Various screw speeds N = 20 rpm, 50 rpm, and 80 rpm were applied, and the barrel/die temperature profile was set as T = 185 °C, 190 °C, 190 °C, 170 °C. The pressure was measured along the screw, and in the die.
The results of validation are presented in Figure 15, Figure 16 and Figure 17. The predictions were quite good, although the simulations were generally overestimated comparing to experiment, and some discrepancies in a few experimental points were observed between the simulation and experiment. For higher screw speeds, N = 50 rpm and N = 80 rpm, the influence of including the slip effects into the model is clearly seen. The computed pressure decreases, and the pressure profile shifts towards the experimental data. For the lower screw speed, N = 20 rpm, the slip effect on modeling is not clear. Two reasons for this may be considered, that is the effect of yield stress and the elongational viscosity (important in the conical die), which may occur at the low shear/elongation rates, and which is not included in the model.

7. Conclusions

So far, any global model of extrusion of wood polymer composites that would include slip effects was not known. In this paper, slip effects were considered and discussed both for WPC rheological studies as well as for WPC extrusion modeling.
The concept of the global model of single screw extrusion process of wood polymer composites has been discussed with respect to slip effects, extensively. Two aspects of slip effects have been considered, rheological and extrusion flow modeling. In the first case, the effect of slipping has to be removed from viscosity measurements. This can be performed by using the Mooney correction. In the second case, slipping in the extrusion process has to be considered. Both these effects can be called the rheological slip effect and the extrusion process slip effect.
Removing the slip effects from viscosity measurements results in increasing the viscosity since the shear rate decreases with the shear stress unchanged. This may have substantial effects in process modeling. From the other side, considering the slip effects in process modeling results in increasing the flow rate and pressure decrease. Both these effects should be considered when modeling the global extrusion process, and the final result is dependent on the current material and process parameters.
Modeling of the single screw extrusion process of wood polymer composites with slip effects is not a trivial problem. It is difficult to determine the slip parameters using the capillary rheometers, and then to model and simulate the extrusion process.
The rheological characteristics of wood polymer composites are generally not available in the material data bases, like CAMPUS or Autodesk MOLDFLOW. The (WPC) material characteristics are strongly dependent on the material structure, and the size of the fiber or flour. The smaller the size of the filler, the easier and more reliable the measurement. Proper preparation of the material is very important. The results are generally difficult to reproduce. When the material is more homogeneous, the results are more reproducible. The temperature flow range is narrow, and it is difficult to determine the melting point. Therefore, the concept of no flow temperature is used to determine the onset of the flow.
The concept of modeling of extrusion of wood polymer composites presented here requires validation in a broader range of material and process parameters, both in relation to the rheological characteristics of the material and the processing conditions. The research will be continued in this direction.
Finally, it is important to note that the concepts of flow modeling with slip effects presented in this paper, i.e., rheological/flow modeling issues, may be applied or extended to other material processing systems, e.g.,
-
Processing of the materials characterized by wall slip,
-
processing of filled/reinforced polymeric materials (by extrusion and injection molding),
-
food and cosmetics processing,
-
pharmaceutical industry [54],
-
3D printing [55,56,57].

Author Contributions

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

Funding

The authors would like to acknowledge support from National Science Center, Poland (DEC-2015/19/B/ST8/00948).

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. Nova-Institut GmbH. Available online: http://news.bio-based.eu/new-market-and-trend-report-published (accessed on 21 November 2017).
  2. Grand View Research 2019, Wood Plastic Composite Market Size, Share & Trends Analysis Report by Type (PE, PP, PVC), Application, Region, Segment Forecasts, 2020–2027. Available online: https://www.grandviewresearch.com/industry-analysis/wood-plastic-composites-market (accessed on 1 February 2020).
  3. Wilczyński, K.; Nastaj, A.; Lewandowski, A.; Wilczyński, K.J.; Buziak, K. Experimental Study for Extrusion of Polypropylene/Wood Flour Composites. Int. Polym. Process. 2015, 30, 113–120. [Google Scholar] [CrossRef]
  4. Wilczyński, K.; Buziak, K.; Wilczyński, K.J.; Lewandowski, A.; Nastaj, A. Computer Modeling for Single-Screw Extrusion of Wood-Plastic Composites. Polymers 2018, 10, 295. [Google Scholar] [CrossRef] [Green Version]
  5. Xiao, K.; Tzoganakis, C. Rheological Properties and Their Influence on Extrusion Characteristics of HDPE-Wood Composite Resins. SPE ANTEC 2002, 39, 252–257. [Google Scholar]
  6. Xiao, K.; Tzoganakis, C. Extrusion Characteristic of HDPE-Wood Composites. SPE ANTEC 2003, 39, 197–201. [Google Scholar]
  7. Xiao, K.; Tzoganakis, C. Rheological Properties of HDPE-Wood Composites. SPE ANTEC 2003, 39, 975–979. [Google Scholar]
  8. Xiao, K.; Tzoganakis, C. Extrusiom of HDPE-Wood Blends. SPE ANTEC 2004, 40, 448–453. [Google Scholar]
  9. Li, T.Q.; Wolcott, M.P. Rheology of Wood Plastics Melt. Part 1. Capillary Rheometry of HDPE Filled with Maple. Polym. Eng. Sci. 2005, 45, 549–559. [Google Scholar] [CrossRef]
  10. Li, T.Q.; Wolcott, M.P. Rheology of Wood Plastics Melt, Part 2: Effects of Lubricating Systems in HDPE/Maple Composites. Polym. Eng. Sci. 2006, 46, 464–473. [Google Scholar] [CrossRef]
  11. Li, T.Q.; Wolcott, M.P. Rheology of Wood Plastics Melt, Part 3: NonLinear Nature of the Flow. Polym. Eng. Sci. 2006, 46, 114–121. [Google Scholar] [CrossRef]
  12. Błędzki, A.K.; Reihmane, S.; Gassan, J. Thermoplastics Reinforced with Wood Fillers: A Literature Review. Polym.-Plast. Technol. Eng. 1998, 37, 451–468. [Google Scholar] [CrossRef]
  13. Mohanty, A.K.; Misra, M.; Drzal, L.T. Natural Fibers, Biopolymers, and Biocomposites; Taylor & Francis: Boca Raton, FL, USA, 2005; ISBN 978-084-931-741-5. [Google Scholar]
  14. Klyosov, A.A. Wood-Plastic Composites; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2007. [Google Scholar] [CrossRef]
  15. Oksman Niska, K.; Sain, M. Wood-Polymer Composites; Woodhead Publishing: Cambridge, UK, 2008. [Google Scholar]
  16. Hristov, V.; Takács, E.; Vlachopoulos, J. Surface Tearing and Wall Slip Phenomena in Extrusion of Highly Filled HDPE/Wood Flour Composites. Polym. Eng. Sci. 2006, 46, 1204–1214. [Google Scholar] [CrossRef]
  17. Vlachopoulos, J.; Hristov, V. Rheology of Wood Polymer Composites. In Proceedings of the Regional PPS Meeting, Pretoria, South Africa, 9–13 October 2006. [Google Scholar]
  18. Santi, C.R.; Hage, E.; Vlachopoulos, J.; Correa, C.A. Rheology and Processing of HDPE/Wood Flour Composites. Int. Polym. Process. 2009, 24, 346–353. [Google Scholar] [CrossRef]
  19. Zolfaghari, A.; Behravesh, A.H.; Shakouri, E.; Soury, E. Flow Balancing in Die Design of Wood Flour/HDPE Composite Extrusion Profiles with Consideration of Rheological Effect. Polym. Eng. Sci. 2010, 50, 543–549. [Google Scholar] [CrossRef]
  20. Lewandowski, A.; Wilczyński, K. Global Modeling of Single Screw Extrusion with Slip Effects. Int. Polym. Process. 2019, 34, 81–90. [Google Scholar] [CrossRef]
  21. Lewandowski, A.; Wilczyński, K. Global Modeling for Single Screw Extrusion of Viscoplastics. Int. Polym. Process. 2020, 35, 26–36. [Google Scholar] [CrossRef]
  22. Tadmor, Z.; Klein, I. Engineering Principles of Plasticating Extrusion; Van Nostrand Reinhold Co.: New York, NY, USA, 1970. [Google Scholar]
  23. White, J.L.; Potente, H. Screw Extrusion; Hanser: Munich, Germany, 2003. [Google Scholar]
  24. Tadmor, Z.; Gogos, C. Principles of Polymer Processing; Wiley & Sons: New York, NY, USA, 2006. [Google Scholar]
  25. Rauwendaal, C. Polymer Extrusion; Hanser: Munich, Germany, 2013. [Google Scholar] [CrossRef]
  26. Agur, E.E.; Vlachopoulos, J. Numerical Simulation of a Single-Screw Plasticating Extruder. Polym. Eng. Sci. 1982, 22, 1084–1094. [Google Scholar] [CrossRef]
  27. 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]
  28. Potente, H.; Hanhart, W.; Schöppner, V. Potential Applications for Computer-Aided Extruder Design. Int. Polym. Process. 1993, 8, 335–344. [Google Scholar] [CrossRef]
  29. Wilczyński, K. Single Screw Extrusion Model for Plasticating Extruders. Polym.-Plast. Technol. Eng. 1999, 38, 581–608. [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. Ilinca, F.; Hétu, J.-F. Three-Dimensional Finite Element Solution of the Flow in Single and Twin-Screw Extruders. Int. Polym. Process. 2010, 25, 275–286. [Google Scholar] [CrossRef] [Green Version]
  32. Altınkaynak, 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]
  33. Malik, M.; Kalyon, D.M.; Golba, J.C. 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]
  34. 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. Eng. 2012, 51, 1567–1577. [Google Scholar] [CrossRef]
  35. Wilczyński, K.; Nastaj, A.; Lewandowski, A.; Wilczyński, K.J. Multipurpose Computer Model for Screw Processing of Plastics. Polym.-Plast. Technol. Eng. 2012, 51, 626–633. [Google Scholar] [CrossRef]
  36. 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] [Green Version]
  37. Hyvärinen, M.; Jabeen, R.; Kärki, T. The Modelling of Extrusion Processes for Polymers—A Review. Polymers 2020, 12, 1306. [Google Scholar] [CrossRef]
  38. Potente, H.; Ridder, H.; Cunha, R.V. Global Concept for Describing and Investigation of Wall Slip Effects in the Extrusion Process. Macromol. Mater. Eng. 2002, 287, 836–842. [Google Scholar] [CrossRef]
  39. Potente, H.; Kurte-Jardin, M.; Klus, S.; Timmermann, K. Two Dimensional Description of Pressure-Throughput Behaviour of Newtonian Materials Considering Wall Slippage Effects. Int. Polym. Process. 2005, 20, 312–321. [Google Scholar] [CrossRef]
  40. Potente, H.; Timmermann, K.; Kurte-Jardin, M. Description of the Pressure/Throughput Behavior of a Single-Screw Plasticating Unit in Consideration of Wall Slippage Effects for Non-Newtonian Material and 1-D Flow. Int. Polym. Process. 2006, 21, 272–282. [Google Scholar] [CrossRef]
  41. Potente, H.; Bornemann, M.; Kurte-Jardin, M. Analytical Model for the Throughput and Drive Power Calculation in the Melting Section of Single Screw Plasticizing Units Considering Wall-Slippage. Int. Polym. Process. 2009, 24, 23–30. [Google Scholar] [CrossRef]
  42. Gooneie, A.; Schuschnigg, S.; Duretek, I.; Holzer, C. Numerical Simulations of the Flow of Wood Polypropylene Composites with Wall Slipping in a Profile Die: The Significance of Material Data. In AIP Conference Proceedings; AIP Publishing LLC: New York, NY, USA, 2015; p. 050014. [Google Scholar] [CrossRef]
  43. Duretek, I.; Schuschnigg, S.; Gooneie, A.; Langecker, G.R.; Holzer, C. Rheological Properties of Wood Polymer Composites and Their Role in Extrusion. J. Phys. Conf. Ser. 2015, 602, 012014. [Google Scholar] [CrossRef] [Green Version]
  44. Macosko, C.W. Rheology: Principles, Measurements and Applications; VCH Publishers: New York, NY, USA, 1994. [Google Scholar]
  45. Moigne, N.; Van den Oever, M.; Budtova, T. Dynamic and Capillary Shear Rheology of Natural Fiber-Reinforced Composites. Polym. Eng. Sci. 2013, 53, 2582–2593. [Google Scholar] [CrossRef]
  46. Polychronopoulos, N.D.; Charlton, Z.; Suwanda, D.; Vlachopoulos, J. Measurements and Comparison to Predictions of Viscosity of Heavily Filled HDPE with Natural Fibers. Adv. Polym. Technol. 2016, 37, 1161–1167. [Google Scholar] [CrossRef]
  47. Hatzikiriakos, S.G.; Dealy, J.M. Wall Slip of Molten High Density Polyethylenes. II. Capillary Rheometer Studies. J. Rheol. 1992, 36, 703–741. [Google Scholar] [CrossRef]
  48. Hatzikiriakos, S.G.; Dealy, J.M. Effects of Interfacial Conditions on Wall Slip and Sharkskin Melt Fracture of HDPE. Int. Polym. Process. 1993, 8, 36–43. [Google Scholar] [CrossRef]
  49. Hatzikiriakos, S.G. The Onset of Wall Slip and Sharkskin Melt Fracture in Capillary Flow. Polym. Eng. Sci. 1994, 34, 1441–1449. [Google Scholar] [CrossRef]
  50. Hatzikiriakos, S.G.; Mitsoulis, E. Slip Effects in Tapered Dies. Polym. Eng. Sci. 2009, 49, 1960–1969. [Google Scholar] [CrossRef]
  51. ANSYS Polyflow v.17. 2017. Available online: https://cfd.hu/doc/ansys_capabilities_170.pdf (accessed on 1 November 2020).
  52. 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]
  53. Wilczyński, K.J.; Buziak, K. Study on the Flow of Wood Plastic Composites in the Starve Fed Single Screw Extrusion. Polimery 2017, 62, 680–685. [Google Scholar] [CrossRef]
  54. Dhaval, M.; Sharma, S.; Dudhat, K.; Chavda, J. Twin-Screw Extruder in Pharmaceutical Industry: History, Working Principle, Applications, and Marketed Products: An In-depth Review. J. Pharm. Innov. 2020. [Google Scholar] [CrossRef]
  55. Rohani Rad, E.; Vahabi, H.; Formela, K.; Saeb, M.R.; Thomas, S. Injectable Poloxamer/Graphene Oxide Hydrogels with Well-Controlled Mechanical and Rheological Properties. Polym. Adv. Technol. 2019, 30, 2250–2260. [Google Scholar] [CrossRef]
  56. Schwab, A.; Levato, R.; D’Este, M.; Piluso, S.; Eglin, D.; Malda, J. Printability and Shape Fidelity of Bioinks in 3D Bioprinting. Chem. Rev. 2020, 120, 11028–11055. [Google Scholar] [CrossRef] [PubMed]
  57. La Galla, A.; Fiorio, R.; Erkoç, M.; Cardon, L.; D’hooge, D.R. Theoretical Evaluation of the Melting Efficiency for the Single-Screw Micro-Extrusion Process: The Case of 3D Printing of ABS. Processes 2020, 8, 1522. [Google Scholar] [CrossRef]
Figure 1. Melting mechanism: (a) Two-dimensional model (Tadmor model), (b) one-dimensional model: W—screw channel width, X—solid bed width, H—screw channel height, Y—solid bed height.
Figure 1. Melting mechanism: (a) Two-dimensional model (Tadmor model), (b) one-dimensional model: W—screw channel width, X—solid bed width, H—screw channel height, Y—solid bed height.
Polymers 13 00622 g001
Figure 2. Material characteristics (PP copo inj 4, Beologic): Viscosity curves.
Figure 2. Material characteristics (PP copo inj 4, Beologic): Viscosity curves.
Polymers 13 00622 g002
Figure 3. Material characteristics (PP copo inj 4, Beologic), Mooney analysis: (a) Shear rate vs. 1/D, (b) slip velocity vs. shear stress, D—capillary diameter.
Figure 3. Material characteristics (PP copo inj 4, Beologic), Mooney analysis: (a) Shear rate vs. 1/D, (b) slip velocity vs. shear stress, D—capillary diameter.
Polymers 13 00622 g003
Figure 4. Effect of rheometric corrections on viscosity curves (PP copo inj 4, Beologic): (a) Rabinowitsch correction, (b) Bagley correction, (c) Mooney correction.
Figure 4. Effect of rheometric corrections on viscosity curves (PP copo inj 4, Beologic): (a) Rabinowitsch correction, (b) Bagley correction, (c) Mooney correction.
Polymers 13 00622 g004
Figure 5. The boundary conditions for flow simulations: (a) Die flow, (b) screw flow.
Figure 5. The boundary conditions for flow simulations: (a) Die flow, (b) screw flow.
Polymers 13 00622 g005
Figure 6. Die characteristics (Fslip = 0.064, eslip = 0.43).
Figure 6. Die characteristics (Fslip = 0.064, eslip = 0.43).
Polymers 13 00622 g006
Figure 7. Pressure and velocity distributions for die flow (Fslip = 0.064, eslip = 0.43): Q = 0.001 kg/s. (a) pressure profile, (b) velocity distribution (in cross-sections indicated by arrows), and pressure distribution, (c) velocity profile.
Figure 7. Pressure and velocity distributions for die flow (Fslip = 0.064, eslip = 0.43): Q = 0.001 kg/s. (a) pressure profile, (b) velocity distribution (in cross-sections indicated by arrows), and pressure distribution, (c) velocity profile.
Polymers 13 00622 g007
Figure 8. Screw characteristics at varoius screw speed (Fslip = 0.064, eslip = 0.43).
Figure 8. Screw characteristics at varoius screw speed (Fslip = 0.064, eslip = 0.43).
Polymers 13 00622 g008
Figure 9. Pressure and velocity distributions for screw flow (Fslip = 0.064, eslip = 0.43): Q = 0.001 kg/s. (a) pressure profile, (b) pressure distribution, (c) velocity distribution, (d) velocity profile.
Figure 9. Pressure and velocity distributions for screw flow (Fslip = 0.064, eslip = 0.43): Q = 0.001 kg/s. (a) pressure profile, (b) pressure distribution, (c) velocity distribution, (d) velocity profile.
Polymers 13 00622 g009
Figure 10. Algorithm of calculation: Q—polymer flow rate (extrusion throughput), Qin—initial polymer flow rate, Δpscrew—pressure build-up in the screw, Δpdie—pressure drop in the die, p—pressure, T—temperature, SBP—solid bed profile, P—power.
Figure 10. Algorithm of calculation: Q—polymer flow rate (extrusion throughput), Qin—initial polymer flow rate, Δpscrew—pressure build-up in the screw, Δpdie—pressure drop in the die, p—pressure, T—temperature, SBP—solid bed profile, P—power.
Polymers 13 00622 g010
Figure 11. Overall extrusion process characteristics.
Figure 11. Overall extrusion process characteristics.
Polymers 13 00622 g011
Figure 12. Process simulation: Pressure profile.
Figure 12. Process simulation: Pressure profile.
Polymers 13 00622 g012
Figure 13. Process simulation: Melting profile.
Figure 13. Process simulation: Melting profile.
Polymers 13 00622 g013
Figure 14. Process simulation: Power consumption profile.
Figure 14. Process simulation: Power consumption profile.
Polymers 13 00622 g014
Figure 15. Validation of simulation: Pressure profile, N = 20 rpm.
Figure 15. Validation of simulation: Pressure profile, N = 20 rpm.
Polymers 13 00622 g015
Figure 16. Validation of simulation: Pressure profile, N = 50 rpm.
Figure 16. Validation of simulation: Pressure profile, N = 50 rpm.
Polymers 13 00622 g016
Figure 17. Validation of simulation: Pressure profile, N = 80 rpm.
Figure 17. Validation of simulation: Pressure profile, N = 80 rpm.
Polymers 13 00622 g017
Table 1. Research program.
Table 1. Research program.
Screw Speed, rpmThroughput (Predicted), kg/h
No SlipSlip
2007.9408.41
5016.0518.26
8027.8731.69
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wilczyński, K.; Buziak, K.; Lewandowski, A.; Nastaj, A.; Wilczyński, K.J. Rheological Basics for Modeling of Extrusion Process of Wood Polymer Composites. Polymers 2021, 13, 622. https://doi.org/10.3390/polym13040622

AMA Style

Wilczyński K, Buziak K, Lewandowski A, Nastaj A, Wilczyński KJ. Rheological Basics for Modeling of Extrusion Process of Wood Polymer Composites. Polymers. 2021; 13(4):622. https://doi.org/10.3390/polym13040622

Chicago/Turabian Style

Wilczyński, Krzysztof, Kamila Buziak, Adrian Lewandowski, Andrzej Nastaj, and Krzysztof J. Wilczyński. 2021. "Rheological Basics for Modeling of Extrusion Process of Wood Polymer Composites" Polymers 13, no. 4: 622. https://doi.org/10.3390/polym13040622

APA Style

Wilczyński, K., Buziak, K., Lewandowski, A., Nastaj, A., & Wilczyński, K. J. (2021). Rheological Basics for Modeling of Extrusion Process of Wood Polymer Composites. Polymers, 13(4), 622. https://doi.org/10.3390/polym13040622

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