The following pages are available to you as for the Regression Model.

• Grundlagen der Strömungsuntersuchung in Schneckenkanälen
• Druck-Durchsatzmodell für Gewindeelemente
• Druck-Durchsatzmodell für dichtkämmende Gewindemischelemente
• Druck-Durchsatzmodell für tangierende Gewindemischelemente
• Druck-Durchsatzmodell für Zahnmischelemente

The conservation laws for mass, momentum and energy are the basis of all mathematical-physical descriptions of flow processes. To describe a fact at any place at any time the velocity profile and the thermodynamical properties pressure, density and temperature have to be known. In addition, the conservation equation has to be connected with an equation to describe the material behavior. For solving the system with the help of conservation laws, additionally simplified assumptions have to be made and boundary conditions have to be given. The quality of the simplifications and the boundary conditions determine here the accuracy of the solution. All steps described above together form the physical-mathematical model.

To describe the flow of polymer melts – as mentioned above – the conservation equations have to be solved. The conservation equations are here first generally given in differential form.

Continuity law

The continuity law is the mathematical formulation of a mass balance in a fixed control room. It states that the mass saved in a volume corresponds to the difference of inflowing and outflowing mass flows.

$$\frac{\partial \rho}{\partial t} + \nabla(\rho \cdot \vec{v}) = 0$$

Here $\rho$ is the density of the polymer melt, t the time and $\vec{v}$ the velocity vector.

Equation of motion (Principle of linear momentum)

The equation of motion is the difference between the momentum going in a volume element and that coming out of it plus the forces effecting the system (e.g. because of gravity).

$$\frac{\partial}{\partial t}(\rho \cdot \vec{v}) = -\nabla(\rho \cdot \vec{v} \cdot \vec{v}) - \nabla p - \nabla \tau + \rho \cdot \vec{a}$$

Law of conservation of energy

A formulation of the balancing of the heat quantity saved in the control room compared to inflowing and outflowing heat flows is the principle of conservation of energy.

$$\rho \cdot c \cdot \left(\frac{\partial T}{\partial t} + \vec{v} \cdot \nabla T\right) = -\nabla \vec{q} + \tau:\nabla \vec{v} + \phi$$

With polymer melts the viscosity depends at the same time on the temperature and the shear rate. The literature describes models which model the single dependencies separately. Finally, the basic approaches are joined to a comprehensive description approach.

Most frequently used approaches to describe the shear rate dependency of the viscosity:

Power law $$\eta(\dot{\gamma}) = K \cdot a_T \cdot \dot{\gamma}^{n-1}$$

Carreau-approach $$\eta(\dot{\gamma}) = A \cdot a_T \cdot (1 + (B \cdot a_T \cdot \dot{\gamma})^c)^{-1}$$

Approach according to Carreau-Michaeli $$\eta(\dot{\gamma}) = A \cdot a_T \cdot (1 + B \cdot a_T \cdot \dot{\gamma})^{-c}$$

Yasuda-approach $$\eta(\dot{\gamma}) = A \cdot a_T \cdot (1 + (B \cdot a_T \cdot \dot{\gamma})^c)^{-\frac{c}{a}}$$

Approaches to describe the temperature dependency of the viscosity:

simplified Arrhenius-approach $$\ln(a_T) = -\beta \cdot (T - T_{ref})$$

Arrhenius-approach $$\ln(a_T) = -\frac{E}{R} \left(\frac{1}{T} - \frac{1}{T_{ref}}\right)$$

WLF-approach $$\log(a_T) = \frac{8.86 \cdot (T_{ref} - T_S)}{101.6 + T_{ref} - T_S} - \frac{8.86 \cdot (T - T_S)}{101.6 + T - T_S}$$

Mathematically the power law is easy to manage. Yet, using this approach has the advantage that this model is able to describe the flow behavior only in certain sections of the flow curve sufficiently well. A description of the whole flow curve in general with the help of this law is not possible. The first three approaches offer better descriptions. Yet, these have the advantage that they are mathematically difficult to manage. In the following the power law is used because of these disadvantages. Its parameter is locally adapted to the Carreau-Michaeli-law. The basis is that the average shear rate is known in the viewed flow section.

The basis of the finite element simulations is the division of a continuum into infinitely many degrees of freedom due to finite elements with finitely many degrees of freedom [Chu89]. These elements are connected at nodes. The profiles of the functions searched for are described with trial functions which depict an approximation of the real profile. With the help of a suitable method the function values at the nodes are combined to a global equation system. Solving this system of equations is the main task of the finite element method. The accuracy can be improved by a smaller division of the continuum or by raising the degree of the trial function of the elements.

The analysis of existing models to specify the pressure-throughput behavior of rectangular and twin screw flow channels gives evidence that these models are only suitable to a limited extent to describe mixing elements. This shall be explained on the basis of the results of a three-dimensional finite element simulation of a threaded-mixing element.

Figure: Velocity profiles in a threaded mixing element.

The figure illustrates the variation of the axial velocities in a conveying threaded mixing element with reconveying grooves.

It can be seen that negative axial velocities predominantly emerge in the area of the grooves. Correspondingly, the volume throughput in the grooves orientates in the opposite direction as the transport direction of the element. Deriving from this as well as from the geometries of the mixing elements the following requirements for pressure-throughput models are to be specified:

• areas, in which $\dot{V} \leq 0$, must be able to be represented,
• geometries, whose pitch to screw diameter ratio is $\frac{t}{D_s} > 2$, must be able to be incorporated and
• geometry sections, whose ratio of channel width to channel height is $\frac{b}{h} < 2$, must be able to be described.

Existing throughput equations for rectangular channels (e.g. [Jun00]) or for twin screw channels (e.g. [Ans93], [Mel00]) cannot fulfil these requirements.

To specify the throughput behavior of the screw elements, first of all, the basic geometry (twin screw channel, rectangular channel and disk elements) is analyzed and modelled. Subsequently, the models are superposed and integrated into an overall model for the description of the throughput behavior of the screw elements.

Table: Dimensionless parameters for the characterization of the geometry.

The data pool for the modelling was established on the basis of finite element simulations for all elements. To limit the amount of finite element calculations the geometry of the elements was standardized beforehand. The dimensionless parameters, which were utilized, are listed in the table.

In addition to the calculations the Power Law was applied. To restrict the error deriving from the application of the power law the pressure-throughput behavior as well as the shear rate were approximated. This enables the specification of the power law parameters in the interested shear rate area of the flow curve.

The presented models constitute approximations of the finite elements results. To attain physically sensible results certain boundary conditions are imposed on the model:

• The characteristic pressure-throughput curves are linear for Newton liquids.
• The models for the rectangular and the twin screw channels shall be applicable for conveying, reconveying and neutral channels.

For the modeling of the pressure-throughput behavior of rectangular channels the following boundary condition also has to be considered:

• In case of a Newtonian fluid (n=1), the influence of the helix angle on the pressure-throughput behavior decreases with an increasing ratio of b/h.

$$\pi_V = \left(\frac{b}{h} \rightarrow \infty, n = 1\right) = f(\varphi_s) \tag{1}$$

Due to the closely intermeshing geometry this requirement cannot be transferred to twin screw channels, because channel width and channel height can-not be varied independently. They are rather dependent on the pitch, the axial distance, the screw diameter, and the number of flights (see figure). Equation (1), however, can only be valid, if the channel width is independent of the pitch.

Figure: Influence of the helix angle on the geometry of twin screw channels.

For the modelling of the pressure-throughput behavior of the channel geometries the requirement of the applicability for conveying, reconveying and neutral channels implies, that the parameters applied by Ansahl [Ans93] and Melisch [Mel00] to specify the pressure-throughput behavior:

• the dimensionless volume throughput

$$\pi_\dot{V} = \frac{\dot{V}}{\frac{1}{2} \cdot h \cdot b \cdot v_0 \cdot \cos(\varphi_s)} \tag{2}$$

• the dimensionless pressure gradient

$$\pi_p = \frac{\Delta p}{\Delta z} \cdot \frac{h^{1+n}}{6 \cdot K \cdot v_0^n \cdot \cos^n(\varphi_s)} \tag{3}$$

are incompatible with this requirement, as they are not defined in the neutral case. On account of this, slightly modified parameters are used in these investigations.

$$\pi_\dot{V} = \frac{\dot{V}}{h \cdot b \cdot v_0} \tag{4}$$

and

$$\pi_p = \frac{\Delta p}{\Delta z} \cdot \frac{h^{1+n}}{K \cdot v_0^n} \tag{5}$$

To describe the process behavior of disk elements the following parameters are applied:

$$\pi_\dot{V} = \frac{\dot{V}}{A_{fr} \cdot D_a \cdot n_0} \tag{6}$$

$$\pi_p = \frac{\Delta p}{L} \cdot \frac{D_a}{K \cdot n_0^n} \tag{7}$$

The subsequent modelling of the pressure-throughput behavior of rectangular and twin screw channels shall form the basis of the modelling of the most diverse screw elements. The models of both channel types are based on the same procedure. The characteristic curve family can be separated into 3 sections (see figure).

• Section I: Here we observe a negative volume throughput. This can emerge in single channel areas; e.g. for grooved threaded elements a back-flow can be detected in the grooves.
• Section II: In this area a positive volume throughput and a positive pressure gradient are effective.
• Section III: In this area a negative pressure gradient and a positive volume throughput can be observed.

Figure: Boundary conditions for the modelling of the pressure-throughput behavior.

For a proper modelling special attention has to be paid to supply a sufficiently good description of the axis intercept points as depicted in the figure, as the partial models for the separate sections are linked with each other on the basis of these points. In the following the intersection point of the characteristic curve with the y-axis will be denominated as (dimensionless) drag conveying capacity. Physically this point represents a condition, in which the channel is completely filled but pressure is just not built up yet. The intersection point of the characteristic curve with the x-axis represents an operating point, in which the drag volume flow rates are as large as the volume flow rates, which arise from a pres-sure back stream. The consequence is that on aggregate there is no flow and that in such a way the examined channel area works as a closed mixing chamber.

The modelling of the pressure-throughput behavior is based upon a multitude of finite element simulations, in which the dimensionless parameters were varied systematically and completely factorial. The variation boundaries were selected in such a way that in part they clearly go beyond the hitherto technically realized spheres, so that a model extrapolation can be avoided as far as possible, which can always be problematic in statistical models. In detail, the parameters were varied within the boundaries as specified in the tables.

Table: Variation of the parameters for the finite element simulations of rectangular channels.

Table: Variation of the parameters for the finite element simulations of twin screw channels.

Table: Variation of the parameters for the finite element simulations of disk elements.

Altogether, 61.727 finite element simulations resulted for rectangular channels, 94.392 simulations for twin screw channels and about 14.000 simulations for the disk elements.

To describe the pressure-throughput behavior of rectangular channels the characteristic curve family is separated into three sections (see figure). For this, all those models were applied as specified in Eqs. (1) to (3). The parameters are dependent on the geometry and on the exponent of the power law. A detailed specification of the model equation can be found in the appendix C1.1 in [Kre04].

Section I: $$\pi_\dot{V} = \left(1 - \frac{\pi_p}{\pi_p|_{\pi_\dot{V}=0}}\right) \cdot \left(C_{RE,I,1} + C_{RE,I,2} \cdot \pi_p + C_{RE,I,3} \cdot \pi_p^2\right) \tag{1}$$

Section II: $$\pi_\dot{V} = \left(1 - \frac{\pi_p}{\pi_p|_{\pi_\dot{V}=0}}\right) \cdot \left(\pi_\dot{V}|_{\pi_p=0} + C_{RE,II,1} \cdot \pi_p + C_{RE,II,2} \cdot \pi_p^2\right) \tag{2}$$

Section III: $$\pi_\dot{V} = \pi_\dot{V}|_{\pi_p=0} + C_{RE,III,1} \cdot \pi_p + C_{RE,III,2} \cdot \pi_p^3 + +C_{RE,III,3} \cdot \pi_p^5 \tag{3}$$

with $$\pi_p|_{\pi_\dot{V}=0} = C_{Re,S,1} \cdot \cos(\varphi_s) + C_{Re,S,2} \cdot \cos^3(\varphi_s) + C_{Re,S,3} \cdot \cos^5(\varphi_s) \tag{4}$$ $$\pi_\dot{V}|_{\pi_p=0} = C_{Re,S,1} \cdot \cos(\varphi_s) + C_{Re,S,2} \cdot \cos^3(\varphi_s) + C_{Re,S,3} \cdot \cos^5(\varphi_s) \tag{5}$$

In the figure the calculated dimensionless pressure-throughput characteristics are contrasted with the Polyflow results of varying exponents of the power law. It can be seen, that the model offers a sufficiently good description for a wide area.

Figure: Comparison of the calculated and the predicted pressure-throughput characteristics for varying exponents of the power law.

The characteristic curves of conveying and reconveying channels are related to each other in a centrosymmetrical way

$$\pi_{\dot{V},FE}(\pi_{p,Fe}) = -\pi_{\dot{V},RFe}(-\pi_{p,RFe}) \tag{6}$$

i.e. if the characteristic curves of the conveying channels are known, those of the reconveying channels can be deduced from this basis.

The characteristic curve to specify the transportation neutral elements correspond with the specification in zone III, as the drag conveying capacity equals zero and as, therefore, the element is always overrun. The figure displays the pressure-throughput characteristics for the three types of transport and varying exponents of the power law. It can be recognized that the above outlined transformation of the characteristics allow a sufficiently good description of the Poly-flow results.

Figure: Comparison of the calculated and the predicted pressure-throughput characteristics of conveying, reconveying and neutral rectangular channels.

In an analogous way to the modelling procedure of rectangular channels in the modelling of twin screw channels also three areas can be separated (see figure). For the description of the particular zones all those models are used as specified in Eqs. (7) to (9). The parameters are dependent on the geometry and on the exponent of the power law. A detailed specification of the model equation can be found in the appendix C2.1 in [Kre04].

Section I: $$\pi_\dot{V} = \left(1 - \frac{\pi_p}{\pi_p|_{\pi_\dot{V}=0}}\right) \cdot \left(C_{DSE,I,1} + C_{DSE,I,2} \cdot \pi_p + C_{DSE,I,3} \cdot \pi_p^2\right) \tag{7}$$

Section II: $$\pi_\dot{V} = \left(1 - \frac{\pi_p}{\pi_p|_{\pi_\dot{V}=0}}\right) \cdot \left(C_{DSE,II,1} + C_{DSE,II,2} \cdot \pi_p + C_{DSE,II,3} \cdot \pi_p^2\right) \tag{8}$$

Section III: $$\pi_\dot{V} = \pi_\dot{V}|_{\pi_p=0} + C_{DSE,III,1} \cdot \pi_p + C_{DSE,III,2} \cdot \pi_p^3 + +C_{DSE,III,3} \cdot \pi_p^5 \tag{9}$$

with $$\pi_\dot{V}|_{\pi_p=0} = \frac{1}{C_{DSE,S,1} + C_{DSE,S,2} \cdot \overline{CL}} \tag{10}$$

$$\pi_p|_{\pi_\dot{V}=0} = \left(C_{DSE,P,1} + C_{DSE,P,2} \cdot \sqrt{CL}\right) \tag{11}$$

In the figure the dimensionless pressure-throughput characteristics are illustrated for varying exponents of the power law. It can be recognized that the model yields a sufficiently good description for a wide area.

Figure: Comparison of the calculated and the predicted pressure-throughput characteristics for varying exponents of the Power Law.

The description model of the throughput behavior of twin screw channels can be applied to conveying and reconveying channels. A comparison of the model predictions and the Polyflow results for both conveying and reconveying channels is illustrated in the figure. For all three transportation types a sufficiently good description can be identified.

Figure: Comparison of the calculated and the predicted pressure-throughput characteristics for conveying, reconveying and transportation neutral twin screw channels.

Contrary to the previous geometries for the description of the throughput behavior of disk elements only one zone is examined. The motive for doing this is the fact that disk elements are neutral elements and, therefore, are always pressure consumers. Accordingly, a complete modelling can simply be attained on the basis of Section III of this element.

$$\pi_\dot{V} = C_{SE,1} \cdot \pi_p + C_{SE,2} \cdot \pi_p^2 + C_{SE,3} \cdot \pi_p^3 \tag{12}$$

Figure: Comparison of the calculated and the predicted pressure-throughput characteristics for varying exponents of the Power Law.

In the figure the dimensionless pressure-throughput characteristics are delineated for varying exponents of the power law. It can be recognized that the model gives a sufficiently good description for a wide area. The parameters in Equation (12) depend on the dimensionless geometry parameters as well as on the exponent of the power law. These are explained in detail in the appendix C3.1 in [Kre04].

In modeling the pressure-throughput behavior of threaded elements two distinct zones have to be differentiated, in which disparate mechanisms affect the transportation of the melt:

• transportation as a result of drag conveying mechanisms in the channel area and
• forced conveying in the intermeshing zone.

In the screw channels there can be observed insignificant, partly even negative, axial velocities. In the intermeshing zone of the screws, on the contrary, relatively high axial velocities are calculated. From this follows that at least for minor throughputs the transportation in the threaded element is dominated by the flow in the intermeshing zone. With these operating points a back-flow takes place in the screw channel.

To be able to represent this effect, first of all, the flow in the intermeshing zone and the flow in the screw channels are modelled separately before they are ultimately linked with each other.

$$\dot{V}_{GE} = \dot{V}_{intermeshing\ area,GE} + \dot{V}_{channel,GE} \tag{1}$$

For the description of the flow in the channel area a channel model is used, in which the screw and the barrel are uncoiled along the barrel wall. Additionally, a kinematic reversal is carried out, i.e. the rotating screw is represented by stationary, parallel channels above which moves the unwounded barrel (see figure).

Figure: Conveyor channel model for double-flighted conveying threaded elements [Ans93]

The velocity $v_0$ of the uncoiled barrel (see figure) equals the peripheral velocity at the outside diameter of the screw, which rotates with a speed of $n_0$.

$$v_0 = \pi \cdot D_s \cdot n_0 \tag{2}$$

The modeling of the throughput behavior is now based on a volume flow balance at the control rooms ABC in Figure.

$$\dot{V}_{channel,GE} = k \cdot \dot{V}_{z,GE} + \dot{V}_{x,GE} \tag{3}$$

For the description of the flow in the conveyor channel $\dot{V}_z$ the approach specified in chapter is applied.

$$\dot{V}_{z,GE} = \pi_{V,z,GE} \cdot (\overline{h} \cdot b_{max} \cdot v_0) \tag{4}$$

It has to be taken into consideration that the calculation of $\pi_{V,z}$ is performed in sections:

$$\pi_{V,z,GE} = \begin{cases} \left(1 - \frac{\pi_{p,z}}{\pi_{p,z}|_{\pi_V=0}}\right) \cdot \left(C_{DSE,I,1} + C_{DSE,I,2} \cdot \pi_{p,z} + C_{DSE,I,3} \cdot \pi_{p,z}^2\right) & \pi_{V,z} < 0 \\ \left(1 - \frac{\pi_{p,z}}{\pi_{p,z}|_{\pi_V=0}}\right) \cdot \left(\pi_V|_{\pi_{p,z}=0} + C_{DSE,II,1} \cdot \pi_{p,z} + C_{DSE,II,2} \cdot \pi_{p,z}^2\right) & \text{für } 0 \leq \pi_{V,z} \leq \pi_V|_{\pi_{p,z}=0} \\ \pi_V|_{\pi_{p,z}=0} + C_{DSE,III,1} \cdot \pi_{p,z} + C_{DSE,III,2} \cdot \pi_{p,z}^3 + +C_{DSE,III,3} \cdot \pi_{p,z}^5 & \pi_V|_{\pi_{p,z}=0} \leq \pi_{V,z} \end{cases} \tag{5}$$

with

$$\pi_{p,z} = \frac{\Delta p}{Z_{Tv}} \cdot \frac{\overline{h}^{1+n}}{K \cdot v_0^n} \tag{6}$$

Accordingly, for the following measures it appears appropriate to formulate Equation (5) in a general form:

$$\pi_{\dot{V},z,GE} = A_{DSE,0} + A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2 + A_{DSE,3} \cdot \pi_{p,z}^3 + A_{DSE,4} \cdot \pi_{p,z}^5 \tag{7}$$

For the description of the leakage flow $\dot{V}_x$ the approach for rectangular channels is used as specified in Chapter Flow in the Channel.

$$\dot{V}_{x,GE} = \pi_{V,x} \cdot (s_R \cdot b_{flight} \cdot v_0) \tag{8}$$

It has to be taken into consideration that the calculation of $\pi_{V,x}$ is performed in sections:

$$\pi_{V,x,GE} = \begin{cases} \left(1 - \frac{\pi_{p,x}}{\pi_{p,x}|_{\pi_V=0}}\right) \cdot \left(C_{RE,I,1} + C_{RE,I,2} \cdot \pi_{p,x} + C_{RE,I,3} \cdot \pi_{p,x}^2\right) & \pi_{V,x} < 0 \\ \\ \left(1 - \frac{\pi_{p,x}}{\pi_{p,x}|_{\pi_V=0}}\right) \cdot \left(\pi_V|_{\pi_{p,x}=0} + C_{RE,II,1} \cdot \pi_{p,x} + C_{RE,II,2} \cdot \pi_{p,x}^2\right) & \text{für } 0 \leq \pi_{V,x} \leq \pi_V|_{\pi_{p,x}=0} \\ \\ \left(\pi_V|_{\pi_{p,x}=0} + C_{RE,III,1} \cdot \pi_{p,x} + C_{RE,III,2} \cdot \pi_{p,x}^3 + +C_{RE,III,3} \cdot \pi_{p,x}^5\right) & \pi_V|_{\pi_{p,x}=0} \leq \pi_{V,x} \end{cases}$$

with

$$\pi_{p,x} = \frac{\Delta p}{e_{max}} \cdot \frac{s_R^{1+n}}{K \cdot v_0^n} \tag{10}$$

By virtue of clarity a general notation is also chosen for Equation (9):

$$\pi_{\dot{V},x} = A_{RE,0} + A_{RE,1} \cdot \pi_{p,x} + A_{RE,2} \cdot \pi_{p,x}^2 + A_{RE,3} \cdot \pi_{p,x}^3 + A_{RE,4} \cdot \pi_{p,x}^5 \tag{11}$$

Are the equations (7) and (11) solved in terms of the dimension affected volume flow rates and are these then applied to Equation (3), one will receive:

$$\dot{V}_{channel,GE} = k \cdot (\overline{h} \cdot b_{max} \cdot v_0) \cdot (A_{DSE,0} + A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2$$ $$+ A_{DSE,3} \cdot \pi_{p,z}^3 + A_{DSE,4} \cdot \pi_{p,z}^5) + (s_R \cdot b_{steg} \cdot v_0) \cdot (A_{RE,0}$$ $$+ A_{RE,1} \cdot \pi_{p,x} + A_{RE,2} \cdot \pi_{p,x}^2 + A_{RE,3} \cdot \pi_{p,x}^3 + A_{RE,4} \cdot \pi_{p,x}^5) \tag{12}$$

Is this equation standardized to $k \cdot (\overline{h} \cdot b_{max} \cdot v_0)$, then follows:

$$\pi_{V,GE} = \frac{\dot{V}_{channel,GE}}{k \cdot (\overline{h} \cdot b_{max} \cdot v_0)} = (A_{DSE,0} + A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2$$ $$+ A_{DSE,3} \cdot \pi_{p,z}^3 + A_{DSE,4} \cdot \pi_{p,z}^5) + \frac{(s_R \cdot b_{steg} \cdot v_0)}{k \cdot (\overline{h} \cdot b_{max} \cdot v_0)} \cdot (A_{RE,0} \tag{13}$$ $$+ A_{RE,1} \cdot \pi_{p,x} + A_{RE,2} \cdot \pi_{p,x}^2 + A_{RE,3} \cdot \pi_{p,x}^3 + A_{RE,4} \cdot \pi_{p,x}^5)$$

As a next measure it proves necessary to build a linkage of the pressure gradients in channel direction and above the screw flight. For this purpose the approach proposed by Ansahl [Ans93] is reverted to. According to Ansahl the following equation applies:

$$\Delta p_x = \frac{\Delta p \cdot b_{max} + e_{max}}{\tan(\varphi_s)} \tag{14}$$

Is Equation (14) inserted into Equation (10) then follows:

$$\pi_{p,x} = \frac{\Delta p \cdot b_{max} + e_{max} \cdot s_R^{1+n}}{\Delta z \cdot e_{max} \cdot \tan(\varphi_s) \cdot K \cdot v_0^n} = \pi_{geo,GE} \cdot \pi_{p,z} \tag{15}$$

with

$$\pi_{geo,GE} = \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \left(\frac{s_R}{\overline{h}}\right)^{1+n} \tag{16}$$

From Equation (15) applied to (13) derives the equation to calculate the pressure-throughput behavior of the channel area of threaded elements for twin screw extruders: $$\pi_{V,GE} = \left(A_{DSE,0} + \pi_1 \cdot A_{RE,0}\right)$$ $$+ \left(A_{DSE,1} + \pi_1 \cdot \pi_{geo} \cdot A_{RE,1}\right) \cdot \pi_{p,z}$$ $$+ \left(A_{DSE,2} + \pi_1 \cdot \pi_{geo}^2 \cdot A_{RE,2}\right) \cdot \pi_{p,z}^2$$ $$+ \left(A_{DSE,3} + \pi_1 \cdot \pi_{geo}^3 \cdot A_{RE,3}\right) \cdot \pi_{p,z}^3$$ $$+ \left(A_{DSE,4} + \pi_1 \cdot \pi_{geo}^5 \cdot A_{RE,4}\right) \cdot \pi_{p,z}^5 \tag{17}$$

with

$$\pi_1 = \frac{(s_R \cdot b_{steg} \cdot v_0)}{k \cdot (\overline{h} \cdot b_{max} \cdot v_0)} \tag{18}$$

In the literature diverse approaches to describe the flow in the intermeshing zone have been developed. In the following the approach of Booy [14] is used.

$$\dot{V}_{intermeshing\ area,GE} = A_{intermeshing\ area} \cdot t \cdot n_0 \tag{19}$$

On the basis of measurements and flow simulations Bakalis [15] demonstrated, that the volume flow rates calculated by means of Equation (5.90) are well concordant with the measurements and the simulation results.

Pressure - Throughput Model

As a final step the volume throughput in the conveyor channel model (Eq. 20) is to be superposed with the volume throughput in the intermeshing zone (Eq. 19):

$$\dot{V}_{GE} = \dot{V}_{intermeshing,GE} + \dot{V}_{channel,GE} = \dot{V}_{intermeshing\ area,GE} + \pi_{V,channel,GE} \cdot (\overline{h} \cdot b_{max} \cdot v_0) \tag{20}$$

In the figure the theoretical pressure-throughput characteristics are illustrated for conveying and reconveying threaded elements.

Figure: Pressure-throughput behavior of conveying and reconveying threaded elements.

In the figure the model predictions are compared with the experimental results. The experiments were performed with a model extruder using silicon oil as melt.

Figure: Comparison of calculated and measured pressure gradients in threaded elements.

It can be recognized that the model can offer a sufficiently good description of the measured pressure gradients.

In the figure the variation of the axial velocities in a conveying, closely intermeshing thread mixing element are illustrated on two different levels. Out of the variation of the velocities it can be noticed that the flow in the grooves and in the intermeshing zone of both screws significantly differs from the flow in the channel area. In the intermeshing zone axial velocities can be detected, which are considerably higher than those in the channel area. In the grooves, however, another effect can be observed. As the grooves are of a back-flow design (in this case), in this area a back-flow is to be recognized accordingly. In the modelling, at first, a distinction is made between the intermeshing zone and the channel zone.

$$\dot{V}_{dGME} = \dot{V}_{channel,dGME} + \dot{V}_{intermeshing\ area,dGME} \tag{1}$$

The influence of the grooves on the behavior of the elements is then taken account of in the description of the separate areas.

Figure: Calculated velocity fields in a tightly intermeshing threaded mixing element.

The channel model for tightly intermeshing threaded mixing elements can be compared to that of threaded elements. The flow in the grooves was accounted for by an additional volume throughput in the balance (Eq. 2) of the marked control room ABC illustrated in see figure.

Figure: Channel model for a reconveying, intermeshing threaded mixing element.

A volume flow balance in the control room as depicted in the figure results in the following equation:

$$\dot{V}_{channel,dGME} = k \cdot \dot{V}_{z,dGME} + \dot{V}_{x,dGME} + l \cdot \dot{V}_{groove,dGME} \tag{2}$$

For the description of the flow in channel direction and in the grooves all models developed in chapter Rectangular channels and chapter Twin Screw Channels respectively are used in a general form:

$$\pi_{\dot{V},z,dGME} = A_{DSE,0} + A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2 + A_{DSE,3} \cdot \pi_{p,z}^3 + A_{DSE,4} \cdot \pi_{p,z}^5 \tag{3}$$

$$\pi_{\dot{V},x,dGME} = A_{RE,0} + A_{RE,1} \cdot \pi_{p,x} + A_{RE,2} \cdot \pi_{p,x}^2 + A_{RE,3} \cdot \pi_{p,x}^3 + A_{RE,4} \cdot \pi_{p,x}^5$$

$$\pi_{\dot{V},groove,dGME} = A_{RE,N,0} + A_{RE,N,1} \cdot \pi_{p,N} + A_{RE,N,2} \cdot \pi_{p,N}^2 + A_{RE,N,3} \cdot \pi_{p,N}^3 + A_{RE,4} \cdot \pi_{p,N}^5$$

with

$$\dot{V}_{z,dGME} = \pi_{\dot{V},z} \cdot (\bar{h} \cdot b_{max} \cdot v_0) \tag{6}$$

$$\dot{V}_{x,dGME} = \pi_{\dot{V},x} \cdot (h_N \cdot b_N \cdot v_0)$$

$$\dot{V}_{groove,dGME} = \pi_{\dot{V},N} \cdot (h_N \cdot b_N \cdot v_0)$$

and

$$\pi_{p,z} = \frac{\Delta p}{\Delta z} \cdot \frac{\bar{h}^{1+n}}{K \cdot v_0^n} \tag{9}$$

$$\pi_{p,x} = \frac{\Delta p}{x} \cdot \frac{s_R^{1+n}}{K \cdot v_0^n} \tag{10}$$

$$\pi_{p,N} = \frac{\Delta p}{x} \cdot \frac{h_N}{K \cdot v_0^n} \tag{11}$$

and

$$\pi_{\dot{V},z,dGME} = \frac{\dot{V}_z}{\bar{h} \cdot b_{max} \cdot v_0} \tag{12}$$

$$\pi_{\dot{V},x,dGME} = \frac{\dot{V}_x}{s_R \cdot b_{Steg} \cdot v_0} \tag{13}$$

$$\pi_{\dot{V},x,dGME} = \frac{\dot{V}_{Nut}}{h_N \cdot b_N \cdot v_0} \tag{14}$$

If now the Equations (3) to (8) are inserted into Equation (2), one receives:

$$\dot{V}_{channel,dGME} = k \cdot (\bar{h} \cdot b_{max} \cdot v_0) \cdot A_{DSE,0} + (A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2$$

$$+A_{DSE,3} \cdot \pi_{p,z}^3 + A_{DSE,4} \cdot \pi_{p,z}^5) + (s_R \cdot b_{steg} \cdot v_0) \cdot (A_{RE,0}$$

$$+A_{RE,1} \cdot \pi_{p,x} + A_{RE,2} \cdot \pi_{p,x}^2 + A_{RE,3} \cdot \pi_{p,x}^3 + A_{RE,4} \cdot \pi_{p,x}^5) \tag{15}$$

$$+l \cdot (h_N \cdot b_N \cdot v_0) \cdot (A_{RE,N,0} + A_{RE,N,1} \cdot \pi_{p,N} + A_{RE,N,2} \cdot \pi_{p,N}^2$$

$$+A_{RE,N,3} \cdot \pi_{p,N}^3 + A_{RE,4} \cdot \pi_{p,N}^5)$$

Is this equation standardized to $k \cdot (\bar{h} \cdot b_{max} \cdot v_0)$, then follows:

$$\pi_{\dot{V},channel,dGME} = \frac{\dot{V}_{channel,dGME}}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} = A_{DSE,0} + (A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2$$

$$+A_{DSE,3} \cdot \pi_{p,z}^3 + A_{DSE,4} \cdot \pi_{p,z}^5) + \frac{(s_R \cdot b_{thread} \cdot v_0)}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} \cdot (A_{RE,0}$$

$$+A_{RE,1} \cdot \pi_{p,x} + A_{RE,2} \cdot \pi_{p,x}^2 + A_{RE,3} \cdot \pi_{p,x}^3 + A_{RE,4} \cdot \pi_{p,x}^5) \tag{16}$$

$$+ \frac{l \cdot (h_N \cdot b_N \cdot v_0)}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} \cdot (A_{RE,N,0} + A_{RE,N,1} \cdot \pi_{p,N} + A_{RE,N,2} \cdot \pi_{p,N}^2$$

$$+A_{RE,N,3} \cdot \pi_{p,N}^3 + A_{RE,4} \cdot \pi_{p,N}^5)$$

To obtain a comprehensive description of the process behavior the pressure gradient in conveyor channel direction is to be linked with the pressure gradients above the radial clearance and with the pressure gradient in the groove. For this purpose the approach proposed by Ansahl [6] is reverted to (see chap. Linkage of the Volume Flow Rates). The linkage of the pressure gradients in the conveyor channel and in the groove can be carried out analogously. It only has to be taken into account that the grooves are inserted in a specific angle in the screw flank, so that the length, in which the pressure gradient is effective, in the grooves differs from the length of the screw flight. Consequently, the pressure gradients can be interlinked as follows.

$$\pi_{p,x} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \frac{s_R^{1+n}}{K \cdot v_0^n} = \pi_{geo,GE} \cdot \pi_{p,z} \tag{17}$$

$$\pi_{p,x} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \frac{h_N^{1+n}}{K \cdot v_0^n} = \pi_{geo,N} \cdot \pi_{p,z} \tag{18}$$

with

$$\pi_{geo,GE} = \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \left(\frac{s_R}{\bar{h}}\right)^{1+n} \tag{19}$$

$$\pi_{geo,N} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \left(\frac{h_N}{\bar{h}}\right)^{1+n} \tag{20}$$

Out of Equation (19), (20) and (16) one receives:

$$\pi_{\dot{V},channel,dGME} = (A_{DSE,0} + A_{DSE,1} \cdot \pi_{p,z} + A_{DSE,2} \cdot \pi_{p,z}^2 + A_{DSE,3} \cdot \pi_{p,z}^3$$

$$+A_{DSE,4} \cdot \pi_{p,z}^5) + \pi_1 \cdot (A_{RE,0} + A_{RE,1} \cdot (\pi_{p,z} \cdot \pi_{geo})$$

$$+A_{RE,2} \cdot (\pi_{p,z} \cdot \pi_{geo})^2 + A_{RE,3} \cdot (\pi_{p,z} \cdot \pi_{geo})^3 + A_{RE,4} \cdot (\pi_{p,z} \cdot \pi_{geo})^5)(21)$$

$$+\pi_2 \cdot (A_{RE,N,0} + A_{RE,N,1} \cdot (\pi_{p,z} \cdot \pi_{geo,N}) + A_{RE,N,2} \cdot (\pi_{p,z} \cdot \pi_{geo,N})^2$$

$$+A_{RE,N,3} \cdot (\pi_{p,z} \cdot \pi_{geo,N})^3 + A_{RE,N,4} \cdot (\pi_{p,z} \cdot \pi_{geo,N})^5)$$

with

$$\pi_1 = \frac{(s_R \cdot b_{thread} \cdot v_0)}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} \tag{22}$$

$$\pi_2 = \frac{l \cdot (h_N \cdot b_N \cdot v_0)}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} \tag{23}$$

the equation to calculate the pressure-throughput behavior of the channel area of threaded elements:

$$\pi_{\dot{V},GE} = (A_{DSE,0} + \pi_1 \cdot A_{RE,0} + \pi_2 \cdot A_{RE,N,0})$$

$$+(A_{DSE,1} + \pi_1 \cdot \pi_{geo} \cdot A_{RE,1} + \pi_2 \cdot \pi_{geo,N} \cdot A_{RE,N,1}) \cdot \pi_{p,z}$$

$$+(A_{DSE,2} + \pi_1 \cdot \pi_{geo}^2 \cdot A_{RE,2} + \pi_2 \cdot \pi_{geo,N}^2 \cdot A_{RE,N,2}) \cdot \pi_{p,z}^2 \tag{24}$$

$$+(A_{DSE,3} + \pi_1 \cdot \pi_{geo}^3 \cdot A_{RE,3} + \pi_2 \cdot \pi_{geo,N}^3 \cdot A_{RE,N,3}) \cdot \pi_{p,z}^3$$

$$+(A_{DSE,4} + \pi_1 \cdot \pi_{geo}^5 \cdot A_{RE,4} + \pi_2 \cdot \pi_{geo,N}^5 \cdot A_{RE,N,4}) \cdot \pi_{p,z}^5$$

The flow in the intermeshing zone is made up of a proportion of chamber transport and a transport in the grooves (see figure).

For the description of the throughput behavior in the intermeshing zone it is assumed that the throughput behavior is substantially dominated by a chamber flow. The grooves in the intermeshing zone alone result in an influence of the pressure gradient on the flow.

$$\dot{V}_{intermeshing\,area,dGME} = \overline A_{intermeshing\,area} \cdot t \cdot n_0 + l_{zw} \cdot \pi_{V,intermeshing\,area} \cdot (h_N \cdot b_N \cdot v_{zw}) \tag{25}$$

To specify the groove flow the description of the rectangular channels is applied:

$$\pi_{\dot{V},intermeshing,dGME} = A_{RE,Z,0} + A_{RE,Z,1} \cdot \pi_{p,zw} + A_{RE,Z,2} \cdot \pi_{p,zw}^2 + A_{RE,Z,3} \cdot \pi_{p,zw}^3 + A_{RE,Z,4} \cdot \pi_{p,zw}^5 \tag{26}$$

Figure: Axial velocities in the intermeshing zone of a closely intermeshing threaded mixing element.

For the linkage of the pressure gradients in the intermeshing zone and in the conveyor channel area it is assumed that the back pressure in the intermeshing zone is negligible and that the same approach can thus be used as above. Consequently, the following equation applies:

$$\pi_{p,zw} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \frac{h_N^{1+n}}{K \cdot v_0^n} = \pi_{geo,zw} \cdot \pi_{p,z} \tag{27}$$

with

$$\pi_{geo,zw} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \left(\frac{h_N}{\bar{h}} \cdot \frac {v_0}{v_{zw}} \right)^{1+n}$$

By the agency of the presented equations one equation can be generated to specify the pressure-throughput behavior of tightly intermeshing threaded mixing elements:

$$\pi_{\dot{V},GE} = \left(A_{DSE,0} + \pi_1 \cdot A_{RE,0} + \pi_2 \cdot A_{RE,N,0} + \pi_3 \cdot A_{RE,Z,0} + \frac{\overline A_{intermeshing\,area} \cdot t \cdot n_0}{k \cdot \bar{h} \cdot b_{max} \cdot v_0}\right)$$

$$+(A_{DSE,1} + \pi_1 \cdot \pi_{geo} \cdot A_{RE,1} + \pi_2 \cdot \pi_{geo,N} \cdot A_{RE,N,1} + \pi_3 \cdot \pi_{geo,zw} \cdot A_{RE,Z,1}) \cdot \pi_{p,z}$$

$$+(A_{DSE,2} + \pi_1 \cdot \pi_{geo}^2 \cdot A_{RE,2} + \pi_2 \cdot \pi_{geo,N}^2 \cdot A_{RE,N,2} + \pi_3 \cdot \pi_{geo,zw}^2 \cdot A_{RE,Z,2}) \cdot \pi_{p,z}^2 \tag{29}$$

$$+(A_{DSE,3} + \pi_1 \cdot \pi_{geo}^3 \cdot A_{RE,3} + \pi_2 \cdot \pi_{geo,N}^3 \cdot A_{RE,N,3} + \pi_3 \cdot \pi_{geo,zw}^3 \cdot A_{RE,Z,3}) \cdot \pi_{p,z}^3$$

$$+(A_{DSE,4} + \pi_1 \cdot \pi_{geo}^5 \cdot A_{RE,4} + \pi_2 \cdot \pi_{geo,N}^5 \cdot A_{RE,N,4} + \pi_3 \cdot \pi_{geo,zw}^5 \cdot A_{RE,Z,4}) \cdot \pi_{p,z}^5$$

with

$$\pi_3 = \frac{l_{zw} \cdot (h_N \cdot b_N \cdot v_{zw})}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} \tag{30}$$

The figure displays a comparison of the measured dimensionless pressure gradient and those calculated by means of the model for a conveying and for a reconveying threaded mixing element. Despite the multiplicity of simplifications it can be observed that the model is in a position to describe the experimental results with sufficient accurateness.

Figure: Comparison of the experimentally investigated and the calculated dimensionless pressure gradients.

The description of non-intermeshing mixing elements is based upon the approach already presented in the previous chapter for tightly intermeshing mixing elements. As a result of the reduced outside diameter it proves possible to build a pair of elements with completely disparate geometries. The description here-after is limited to elements, which consist of a conveying and a reconveying element with a similar geometry. The only distinction between the elements shall be the direction of the pitch.

When examining the flow in the intermeshing zone of these elements it directly attracts the attention that the axial velocities strongly deviate along the element pair. In those areas in the figure marked with an „a“ a maximum axial velocity of about +40 mm/s can be identified and in those marked with a „b“ an axial velocity of about -40 mm/s. Out of this the simplification can be deduced that the sum of the volume flow rates in two sequenced intermeshing zones approximately equals zero.

Figure: Distribution of the axial velocities in non-intermeshing threaded mixing elements.

As the volume flow rates consequently equalize in sequenced intermeshing zones, it follows, that the intermeshing zone does not make any contribution, or at least to a negligible degree, to the pressure-throughput behavior. Therefore, in the following this area will be disregarded.

The unwinding of the non-intermeshing threaded mixing element into the section results in the conveyor channel model as depicted in the figure.

Figure: Conveyor channel model for non-intermeshing threaded mixing elements.

The balance of the volume flow rates in the control room ABCD yields:

$$\dot{V}_{tGME} = (k \cdot \dot{V}_{z,fe} + \dot{V}_{x,fe} + l \cdot \dot{V}_{groove,fe}) + (k \cdot \dot{V}_{z,rfe} + \dot{V}_{x,rfe} + l \cdot \dot{V}_{groove,rfe}) \tag{1}$$

with

$$\dot{V}_{z,fe} = \pi_{\dot V,z,fe} \cdot (\bar{h} \cdot b_{max} \cdot v_0) \tag{2}$$

$$\dot{V}_{x,fe} = \pi_{\dot V,x,fe} \cdot (s_R \cdot b_{thread} \cdot v_0) \tag{3}$$

$$\dot{V}_{groove,fe} = \pi_{\dot V,N,fe} \cdot (h_N \cdot b_N \cdot v_0) \tag{4}$$

$$\dot{V}_{z,rfe} = \pi_{\dot V,z,rfe} \cdot (\bar{h} \cdot b_{max} \cdot v_0) \tag{5}$$

$$\dot{V}_{x,rfe} = \pi_{\dot V,x,rfe} \cdot (s_R \cdot b_{thread} \cdot v_0) \tag{6}$$

$$\dot{V}_{groove,rfe} = \pi_{\dot V,N,rfe} \cdot (h_N \cdot b_N \cdot v_0) \tag{7}$$

The dimensionless volume flow rates are defined as follows:

$$\pi_{\dot{V},z,fe} = A_{DSE,fe,0} + A_{DSE,fe,1} \cdot \pi_{p,z,fe} + A_{DSE,fe,2} \cdot \pi_{p,z,fe}^2 \tag{8}$$

$$+A_{DSE,fe,3} \cdot \pi_{p,z,fe}^3 + A_{DSE,fe,4} \cdot \pi_{p,z,fe}^5$$

$$\pi_{\dot{V},x,fe} = A_{RE,fe,0} + A_{RE,fe,1} \cdot \pi_{p,x,fe} + A_{RE,fe,2} \cdot \pi_{p,x,fe}^2 \tag{9}$$

$$+A_{RE,fe,3} \cdot \pi_{p,x,fe}^3 + A_{RE,fe,4} \cdot \pi_{p,x,fe}^5$$

$$\pi_{\dot{V},N,fe} = A_{RE,N,fe,0} + A_{RE,N,fe,1} \cdot \pi_{p,N} + A_{RE,N,fe,2} \cdot \pi_{p,N}^2 \tag{10}$$

$$+A_{RE,N,fe,3} \cdot \pi_{p,N}^3 + A_{RE,N,fe,4} \cdot \pi_{p,N}^5$$

$$\pi_{\dot{V},z,rfe} = A_{DSE,rfe,0} + A_{DSE,rfe,1} \cdot \pi_{p,z,rfe} + A_{DSE,rfe,2} \cdot \pi_{p,z,rfe}^2 \tag{11}$$

$$+A_{DSE,rfe,3} \cdot \pi_{p,z,rfe}^3 + A_{DSE,rfe,4} \cdot \pi_{p,z,rfe}^5$$

$$\pi_{\dot{V},x,rfe} = A_{RE,rfe,0} + A_{RE,rfe,1} \cdot \pi_{p,x,rfe} + A_{RE,rfe,2} \cdot \pi_{p,x,rfe}^2 \tag{12}$$

$$+A_{RE,rfe,3} \cdot \pi_{p,x,rfe}^3 + A_{RE,rfe,4} \cdot \pi_{p,x,rfe}^5$$

$$\pi_{\dot{V},N,rfe} = A_{RE,N,rfe,0} + A_{RE,N,rfe,1} \cdot \pi_{p,N,rfe} + A_{RE,N,rfe,2} \cdot \pi_{p,N,rfe}^2 \tag{13}$$

$$+A_{RE,N,rfe,3} \cdot \pi_{p,N,rfe}^3 + A_{RE,N,rfe,4} \cdot \pi_{p,N,rfe}^5$$

with the parameters defined in the table.

Simplifying assumed that the isobars operate at right angle to the screw axis the following equation applies:

$$\frac{\Delta p_{fe}}{\Delta z} = \frac{\Delta p_{rfe}}{\Delta z} \Leftrightarrow \pi_{p,z,fe} = \pi_{p,z,rfe} \tag{14}$$

The pressure gradients in channel direction and at right angle to it are again linked in an analogous manner to the approaches introduced so far.

$$\Delta p_x = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{\tan(\varphi_s)} \tag{15}$$

** Dimensionless parameters for the description of the pressure-throughput behavior of non-intermeshing threaded mixing elements: **

Consequently, the dimensionless factors can also be defined for the linkage of the volume flow rates:

$$\pi_{p,x,fe} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \frac{s_R^{1+n}}{K \cdot v_0^n} = \pi_{geo,x} \cdot \pi_{p,z} \tag{17}$$

$$\pi_{p,N,fe} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \frac{h_N^{1+n}}{K \cdot v_0^n} = \pi_{geo,N} \cdot \pi_{p,z} \tag{18}$$

$$\pi_{p,x,rfe} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \frac{s_R^{1+n}}{K \cdot v_0^n} = \pi_{geo,x} \cdot \pi_{p,z} \tag{19}$$

$$\pi_{p,N,rfe} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \frac{h_N^{1+n}}{K \cdot v_0^n} = \pi_{geo,N} \cdot \pi_{p,z}$$

with:

$$\pi_{geo,x} = \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \left(\frac{s_R}{\bar{h}}\right)^{1+n} \tag{20}$$

$$\pi_{geo,N} = \frac{\Delta p}{\Delta z} \cdot \frac{b_{max} + e_{max}}{e_{max} \cdot \tan(\varphi_s)} \cdot \cos\left(\frac{\pi}{2} - \varphi_s - \varphi_N\right) \cdot \left(\frac{h_N}{\bar{h}}\right)^{1+n} \tag{21}$$

Are these equations inserted into the balance (Eq. 1), one receives the model to specify the pressure-throughput behavior in the conveyor channel area:

$$\pi_{\dot{V},tGME} = \frac{\dot{V}_{tdGME}}{k \cdot (\bar{h} \cdot b_{max} \cdot v_0)} = (A_{DSE,fe,0} + A_{DSE,rfe,0} + \pi_1 \cdot (A_{RE,fe,0} + A_{RE,rfe,0}))$$

$$+\pi_2 \cdot (A_{RE,N,fe,0} + A_{RE,N,rfe,0})) + (A_{DSE,fe,1} + A_{DSE,rfe,1}$$

$$+\pi_1 \cdot \pi_{geo,x} \cdot (A_{RE,fe,1} + A_{RE,rfe,1}) + \pi_3 \cdot \pi_{geo,N} \cdot (A_{RE,N,fe,1} + A_{RE,N,rfe,1})) \cdot \pi_{p,z}$$

$$+(A_{DSE,fe,2} + A_{DSE,rfe,2} + \pi_1 \cdot \pi_{geo,x}^2 \cdot (A_{RE,fe,2} + A_{RE,rfe,2})$$

$$+\pi_2 \cdot \pi_{geo,N}^2 \cdot (A_{RE,N,fe,2} + A_{RE,N,rfe,2})) \cdot \pi_{p,z}^2 + (A_{DSE,fe,3} + A_{DSE,rfe,3}$$

$$+\pi_1 \cdot \pi_{geo,x}^3 \cdot (A_{RE,fe,3} + A_{RE,rfe,3}) + \pi_2 \cdot \pi_{geo,N}^3 \cdot (A_{RE,N,fe,3} + A_{RE,N,rfe,3})) \cdot \pi_{p,z}^3$$

$$+(A_{DSE,fe,4} + A_{DSE,rfe,4} + \pi_1 \cdot \pi_{geo,x}^5 \cdot (A_{RE,fe,4} + A_{RE,rfe,4})$$

$$+\pi_2 \cdot \pi_{geo,N}^5 \cdot (A_{RE,N,fe,4} + A_{RE,N,rfe,4})) \cdot \pi_{p,z}^5 \tag{22}$$

The comparison of the model predictions with the experimental results for distinct polymers confirms, that the model yields sufficient precision despite its multitude of simplifications.

Figure: Comparison of the calculated and the measured dimensionless pressure gradients for a non-intermeshing threaded mixing element.

For the calculation of the pressure-throughput behavior in turbine mixing elements the following model assumptions are necessary:

• The clearances in the intermeshing zone of the turbine mixing element are constant ($s_F \approx s_W$)
• The flow in the grooves can be described sufficiently accurate with a flow in the rectangular channels.

Figure: Real turbine mixing elements geometry and the geometrical substitution model.

The description of the process behavior is based upon the superimposition of the volume flow rates in the grooves and in the disk area.

$$\dot{V}_{ZE} = \dot{V}_{disc} + l_N \cdot \dot{V}_{groove} \tag{1}$$

To specify the volume flow rate the approach is applied as introduced in chapter Pressure-Throughput Model for Threaded Elements:

$$\dot{V}_{disc} = A_{free,SE} \cdot D_a \cdot n_0 \cdot (A_{SE,1} \cdot \pi_{p,SE} + A_{SE,2} \cdot \pi_{p,SE}^2 + A_{SE,3} \cdot \pi_{p,SE}^5) \tag{2}$$

The modeling of the throughput behavior in the grooves is performed by virtue of the approaches presented in chapter Pressure-Throughput Model for Threaded Elements:

$$\dot{V}_{groove} = h_N \cdot b_N \cdot v_0 \cdot (A_{RE,0} + A_{RE,1} \cdot \pi_{p,N} + A_{RE,2} \cdot \pi_{p,N}^2 + A_{RE,3} \cdot \pi_{p,N}^3 + A_{RE,4} \cdot \pi_{p,N}^5) \tag{3}$$

Out of this results the following equation to describe the pressure-throughput behavior of toothed mixing elements:

$$\pi_{\dot V,ZE} = A_{ZE,0} + A_{ZE,1} \cdot \pi_p + A_{ZE,2} \cdot \pi_p^2 + A_{ZE,3} \cdot \pi_p^3 + A_{ZE,4} \cdot \pi_p^5 \tag{4}$$

with:

$$\pi_{\dot V,ZE} = \frac{\dot{V}_{ZE}}{A_{fr} \cdot D_a \cdot n_0} \tag{5}$$

$$\pi_p = \frac{\Delta p}{L} \cdot \frac{D_a}{K \cdot n_0^n} \tag{6}$$

and

$$A_{ZE,0} = A_{RE,0} \cdot i_N \cdot \pi_{geo,V} \tag{7}$$

$$A_{ZE,1} = A_{SE,1} + A_{RE,1} \cdot i_N \cdot \pi_{geo,V} \cdot \pi_{geo,p} \tag{8}$$

$$A_{ZE,2} = A_{RE,2} \cdot i_N \cdot \pi_{geo,V} \cdot \pi_{geo,p}^2 \tag{9}$$

$$A_{ZE,3} = A_{SE,2} + A_{RE,3} \cdot i_N \cdot \pi_{geo,V} \cdot \pi_{geo,p}^3 \tag{10}$$

$$A_{ZE,4} = A_{SE,3} + A_{RE,4} \cdot i_N \cdot \pi_{geo,V} \cdot \pi_{geo,p}^5 \tag{11}$$

The two factors $\pi_{geo,V}$ and $\pi_{geo,p}$ are both necessitated for the linkage of the two volume flow rates. They are defined as:

$$\pi_{geo,V} = \frac{h_N \cdot b_N \cdot v_0}{A_{free} \cdot D_a \cdot n_0} \tag{12}$$

and:

$$\pi_{geo,p} = \sin(\varphi_{groove}) \cdot \pi^n \cdot \frac{D_a^{1+n}}{h^{1+n}} \tag{13}$$

The figure shows a comparison of the model predictions with the experimental results. It can be recognized that the model is in a position to describe the experimental results with satisfactory precision.

Figure: Comparison of experimental and approximated dimensionless pressure gradients of conveying turbine mixing elements.