Non‐Isothermal Modeling of a Non‐Newtonian Fluid Flow in a Twin Screw Extruder Using the Fictitious Domain Method

Abstract

In this study, using a developed mathematical model, the non‐isothermal behavior of a non‐Newtonian fluid flow in the conveying elements of an intermeshing co‐rotating twin screw extruder (TSE) is simulated based on the combination of mixed finite elements and fictitious domain methods. The flow equations are solved employing the standard Galerkin method and a streamline‐upwind/Petrov–Galerkin technique is used in the solution scheme of the energy equation to reduce numerical oscillation. This model is combined with the Carreau rheological model to solve governing equations of continuity, momentum, and energy in a 3D Cartesian coordinate system. Using a developed mathematical model, the velocity and pressure fields are simulated for a non‐isothermal flow in a co‐rotating TSE. The shear rate distribution as a criterion for viscous dissipation and also the temperature distribution is calculated based on the simulated flow fields. The applicability of this model is verified by the comparison of experimentally measured pressures and mass flow rates with the simulation results for a high‐density polyethylene melt. This comparison shows a good correlation between experimental data and model predictions.