Boundary-Conforming Finite Element Methods for Twin-Screw Extruders: Unsteady - Temperature-Dependent - Non-Newtonian Simulations

We present a boundary-conforming space-time finite element method to compute the flow inside co-rotating, self-wiping twin-screw extruders. The mesh update is carried out using the newly developed Snapping Reference Mesh Update Method (SRMUM). It allows to compute time-dependent flow solutions inside twin-screw extruders equipped with conveying screw elements without any need for re-meshing and projections of solutions - making it a very efficient method. We provide cases for Newtonian and non-Newtonian fluids in 2D and 3D, that show mesh convergence of the solution as well as agreement to experimental results. Furthermore, a complex, unsteady and temperature-dependent 3D test case with multiple screw elements illustrates the potential of the method also for industrial applications.

[1]  Xianming Zhang,et al.  Numerical simulation and experimental validation of mixing performance of kneading discs in a twin screw extruder , 2009 .

[2]  Thierry Coupez,et al.  A Direct 3D Numerical Simulation Code for Extrusion and Mixing Processes , 2009 .

[3]  P. Anderson,et al.  Simulation of distributive mixing inside mixing elements of co-rotating twin-screw extruders , 2013 .

[4]  P. Cleary,et al.  Accuracy analysis of SPH for flow in a model extruder with a kneading element , 2018, Chemical Engineering Science.

[5]  C. F. Curtiss,et al.  Dynamics of Polymeric Liquids, Volume 1: Fluid Mechanics , 1987 .

[6]  J. Petrie,et al.  Application of the grid convergence index to a laminar axisymmetric sudden expansion flow , 2016 .

[7]  Philippe A. Tanguy,et al.  Adaptive finite element simulations of fluid flow in twin-screw extruders , 2003, Comput. Chem. Eng..

[8]  Patrick Patrick Anderson,et al.  Extended finite element method for viscous flow inside complex three‐dimensional geometries with moving internal boundaries , 2012 .

[9]  Kenneth E. Jansen,et al.  A better consistency for low-order stabilized finite element methods , 1999 .

[10]  M. Malik,et al.  An Integrated Approach for Numerical Analysis of Coupled Flow and Heat Transfer in Co-rotating Twin Screw Extruders , 2007 .

[12]  Johannes Khinast,et al.  Co-Rotating Twin-Screw Extruders: Detailed Analysis of Conveying Elements Based on Smoothed Particle Hydrodynamics. Part 2: Mixing , 2015 .

[13]  Jean-Francois Hetu,et al.  Immersed boundary finite elements for 3D flow simulations in twin-screw extruders , 2013 .

[14]  Volker John,et al.  On the Divergence Constraint in Mixed Finite Element Methods for Incompressible Flows , 2015, SIAM Rev..

[15]  P. Cleary,et al.  Effect of geometry and fill level on the transport and mixing behaviour of a co-rotating twin screw extruder , 2018, Computational Particle Mechanics.

[16]  T. Tezduyar,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. I: The concept and the preliminary numerical tests , 1992 .

[17]  M. Malik,et al.  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 , 2014 .

[18]  Marek Behr,et al.  On stabilized space‐time FEM for anisotropic meshes: Incompressible Navier–Stokes equations and applications to blood flow in medical devices , 2017 .

[19]  A. Murakami,et al.  Heat transfer behavior of melting polymers in laminar flow field , 2004 .

[20]  Marek Behr,et al.  The Shear-Slip Mesh Update Method , 1999 .

[21]  Yuri Bazilevs,et al.  An immersogeometric variational framework for fluid-structure interaction: application to bioprosthetic heart valves. , 2015, Computer methods in applied mechanics and engineering.

[22]  Thomas J. R. Hughes,et al.  NURBS-based isogeometric analysis for the computation of flows about rotating components , 2008 .

[23]  Curtiss,et al.  Dynamics of Polymeric Liquids , .

[24]  Thomas F. Irvine,et al.  Shear rate dependent thermal conductivity measurements of non-Newtonian fluids , 1997 .

[25]  Dong-Ryul Lee Shear rate dependence of thermal conductivity and its effect on heat transfer in a Non-Newtonian flow system , 1998 .

[26]  Pierre J. Carreau,et al.  Review of some useful rheological equations , 1979 .

[27]  Patrick Patrick Anderson,et al.  Adaptive non‐conformal mesh refinement and extended finite element method for viscous flow inside complex moving geometries , 2012 .

[28]  Andrew N. Hrymak,et al.  Numerical simulation of pressure and velocity profiles in kneading elements of a co‐rotating twin screw extruder , 2000 .

[29]  Tim A. Osswald,et al.  Polymer rheology : fundamentals and applications , 2014 .

[30]  L. Richardson The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[31]  F. Ilinca,et al.  Three-dimensional Finite Element Solution of the Flow in Single and Twin-Screw Extruders , 2010 .

[32]  Kazumori Funatsu,et al.  3-D non-isothermal flow field analysis and mixing performance evaluation of kneading blocks in a co-rotating twin srew extruder , 2001 .

[33]  G. Wang,et al.  Effect of Dynamic Center Region on the Flow and Mixing Efficiency in a New Tri-Screw Extruder Using 3D Finite Element Modeling , 2013 .

[34]  M. Kostic,et al.  Investigation of Thermal Conductivity of a Polymer Solution as Function of Shearing Rate , 1999, Heat Transfer: Volume 4.

[35]  S. Mittal,et al.  Computation of unsteady incompressible flows with the stabilized finite element methods: Space-time formulations, iterative strategies and massively parallel implementations , 1992 .

[36]  T. Tezduyar,et al.  Mesh Moving Techniques for Fluid-Structure Interactions With Large Displacements , 2003 .

[37]  Sadao Sato,et al.  Effects of filler on heat transmission behavior of flowing melt polymer composites , 2006 .

[38]  Marek Behr,et al.  Stabilized finite element method for flows with multiple reference frames , 2015 .

[39]  M. L. Booy,et al.  Geometry of fully wiped twin-screw equipment , 1978 .

[40]  Mukund V. Karwe,et al.  Velocity distributions and volume flow rates in the nip and translational regions of a co-rotating, self-wiping, twin-screw extruder , 2002 .