BACKMIXING IN SCREW EXTRUDERS

Similar documents
The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

SECTION 2 Electrostatics

Fluid Physics 8.292J/12.330J % (1)

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Upper bound for ropelength of pretzel knots

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

SNAP Centre Workshop. Basic Algebraic Manipulation

True Nature of Potential Energy of a Hydrogen Atom

Principle Of Superposition

Parallel Vector Algorithms David A. Padua

CM3110 Transport Processes and Unit Operations I

MATH Exam 1 Solutions February 24, 2016

SOLUTIONS: ECE 606 Homework Week 7 Mark Lundstrom Purdue University (revised 3/27/13) e E i E T

U8L1: Sec Equations of Lines in R 2

Information-based Feature Selection

11TH INTERNATIONAL SYMPOSIUM ON PARTICLE IMAGE VELOCIMETRY - PIV15 Santa Barbara, California, Sept , 2015

Kinetics of Complex Reactions

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

Boundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary

The target reliability and design working life

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

Salmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations

The Generalized Newtonian Fluid - Isothermal Flows Constitutive Equations! Viscosity Models! Solution of Flow Problems!

CS / MCS 401 Homework 3 grader solutions

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation

EXPERIMENTAL INVESTIGATION ON LAMINAR HIGHLY CONCENTRATED FLOW MODELED BY A PLASTIC LAW Session 5

EVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD

Calculus 2 Test File Spring Test #1

OPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES


Math 176 Calculus Sec. 5.1: Areas and Distances (Using Finite Sums)

The Random Walk For Dummies

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

EDEXCEL NATIONAL CERTIFICATE UNIT 4 MATHEMATICS FOR TECHNICIANS OUTCOME 4 - CALCULUS

Chapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation

Random Variables, Sampling and Estimation

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

Math 116 Second Exam

The Method of Least Squares. To understand least squares fitting of data.

Calculus with Analytic Geometry 2

CENTRIFUGAL PUMP SPECIFIC SPEED PRIMER AND THE AFFINITY LAWS Jacques Chaurette p. eng., Fluide Design Inc. November 2004

7.) Consider the region bounded by y = x 2, y = x - 1, x = -1 and x = 1. Find the volume of the solid produced by revolving the region around x = 3.

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

Measures of Spread: Standard Deviation

Ray-triangle intersection

Multicomponent-Liquid-Fuel Vaporization with Complex Configuration

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim

Ma 530 Introduction to Power Series

R is a scalar defined as follows:

Optimization of the Brownie Pan

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

Overview. p 2. Chapter 9. Pooled Estimate of. q = 1 p. Notation for Two Proportions. Inferences about Two Proportions. Assumptions

Statistical Fundamentals and Control Charts

Session 5. (1) Principal component analysis and Karhunen-Loève transformation

Deterministic Model of Multipath Fading for Circular and Parabolic Reflector Patterns

CHAPTER 8 SYSTEMS OF PARTICLES

Streamfunction-Vorticity Formulation

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Calculus 2 Test File Fall 2013

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

Math 142, Final Exam. 5/2/11.

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

FIR Filter Design: Part II

6.) Find the y-coordinate of the centroid (use your calculator for any integrations) of the region bounded by y = cos x, y = 0, x = - /2 and x = /2.

Chapter 2 Feedback Control Theory Continued

Numerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

MATH 1080: Calculus of One Variable II Fall 2017 Textbook: Single Variable Calculus: Early Transcendentals, 7e, by James Stewart.

The "Last Riddle" of Pierre de Fermat, II

a b c d e f g h Supplementary Information

Math 113, Calculus II Winter 2007 Final Exam Solutions

Statistics 511 Additional Materials

Goodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen)

Infinite Sequences and Series

MATH 10550, EXAM 3 SOLUTIONS

Measuring Scales. Measuring Scales

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS

The Graetz Problem. R. Shankar Subramanian. , which can be larger or smaller than T 0

Mark Lundstrom Spring SOLUTIONS: ECE 305 Homework: Week 5. Mark Lundstrom Purdue University

C. C. Fu, Ph.D., P.E.

Monte Carlo Integration

Carleton College, Winter 2017 Math 121, Practice Final Prof. Jones. Note: the exam will have a section of true-false questions, like the one below.

Activity 3: Length Measurements with the Four-Sided Meter Stick

OPTIMAL PIECEWISE UNIFORM VECTOR QUANTIZATION OF THE MEMORYLESS LAPLACIAN SOURCE

SECOND-LAW ANALYSIS OF LAMINAR NON- NEWTONIAN GRAVITY-DRIVEN LIQUID FILM ALONG AN INCLINED HEATED PLATE WITH VISCOUS DISSIPATION EFFECT

Analysis of composites with multiple rigid-line reinforcements by the BEM

Statistical Intervals for a Single Sample

CALCULUS AB SECTION I, Part A Time 60 minutes Number of questions 30 A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAM.

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials

7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals

Probability, Expectation Value and Uncertainty

6.3 Testing Series With Positive Terms

1 Cabin. Professor: What is. Student: ln Cabin oh Log Cabin! Professor: No. Log Cabin + C = A Houseboat!

Orthogonal Gaussian Filters for Signal Processing

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

Transcription:

BACKMIXING IN SCREW EXTRUDERS Chris Rauwedaal, Rauwedaal Extrusio Egieerig, Ic. Paul Grama, The Madiso Group Abstract Mixig is a critical fuctio i most extrusio operatios. Oe of the most difficult mixig tasks is backmixig. A extrusio operatio where good backmixig is very importat is whe a low percetage color cocetrate, CC, is added to a virgi polymer. I this case, the iitial distace betwee the CC pellets may be mm or greater. If the fial striatio thickess eeds to be reduced to the micro level, the reductio of the striatio thickess eeds to be at least five orders of magitude this is quite a tough task! This paper will aalyze how the velocity profiles, axial mixig, ad residece time distributio are related. It will be show why simple coveyig screws have poor axial mixig capability. New mixer geometries that are specifically desiged to improve backmixig will be discussed. Results from extrusio experimets will be preseted. Cross Sectioal Mixig ad Axial Mixig Most aalyses of mixig focus o cross sectioal mixig, e.g. (). The cross sectioal mixig is determied mostly by the Couette shear rate betwee the rotatig screw ad statioary barrel. Usig the flat plate approximatio ad assumig pure drag flow, this shear rate ca be approximated by: πdn γ = H & () where D is the barrel diameter, N the rotatioal speed, ad H the chael depth. More accurate expressios for the shear rate have bee developed usig cylidrical coordiates (). Typical values of the Couette shear rate i sigle screw extruders rage from 5 to sec -. With a typical residece time i the melt coveyig zoe of about secods, the resultig total shear strai rages from about, to, uits. This meas that the striatio thickess i cross sectioal mixig is reduced by about three orders of magitude. I may cases, this is ot eough to achieve a level of mixig that appears uiform to the aked eye. plates, the dimesioless velocity φ=v/v max ca be writte as a fuctio of the dimesioless ormal coordiate ξ=y/h as follows: + φ ( ξ ) = ξ () Figure shows the velocity distributio for several values of the power law idex. The velocity profile for a Newtoia fluid (=.) is a parabola. The shear rate i the ceter of the chael is zero. As a result, o mixig will take place there. The ceter regio flattes as the power law idex reduces. I other words, the velocity profile becomes closer to a plug flow profile as the power law idex approaches zero. This meas that the low shear rate regio expads as the fluid becomes more shear thiig. Thus, the regio with poor mixig becomes larger whe the power law idex reduces. From this simple aalysis it becomes clear that the situatio for backmixig is substatially more difficult tha for cross sectioal mixig. For shear thiig fluids there is a cosiderable regio i the ceter of the chael where little or o axial mixig takes place. Residece Time Distributio The RTD ca be determied from the velocity profiles i the chael. The axial flow i a screw extruder is a pressure flow because there are o axial velocity compoets of the screw or barrel. As a first approximatio, the axial pressure flow ca be cosidered a flow betwee parallel plates. The velocity profile for pressure flow of a power law fluid betwee parallel plates is (4): WH H P V = ( + ) ml & (3) Where W is the width of the chael, H the height of the chael, P the pressure drop over axial legth L, the power law idex, ad m the cosistecy idex of the fluid. The velocity profile ca be writte as: Axial mixig or backmixig occurs by pressure flow (3). For pressure flow of a power law fluid betwee parallel

H H P v y) = ( + ) ml y H + ( (4) where ormal coordiate y rages from H/ to +H/. The velocity profile ca be writte as: v y) = v max y H + ( (5) The exteral RTD fuctio f(t)dt ca be determied from: dv& f ( t) dt = V& ( + ) y = ( + ) H H + dy Coordiate y ca be expressed as a fuctio of time by the followig substitutio: (6) L t = (7) v(y) With this substitutio, the exteral RTD fuctio ca be writte as: f ( t) dt ( + ) t ( + ) t + = (8) The miimum residece time t ca be determied from: L t = (9) v max The cumulative RTD fuctio F(t) ca be foud by itegratig the exteral RTD fuctio; it ca be expressed as: F t) + = t t t + + t ( () For a Newtoia fluid the RTD fuctio becomes: t = t + F( t) t t t t 3 dt () If we express the RTD as a fuctio of the dimesioless residece time θ, where θ is the actual residece time divided by the mea residece time, we get: F( θ ) = ( + + ) θ + + ( + ) θ () The expressio above for a power law fluid has ot bee published before. With this expressio the RTD ca plotted at several values of the power law idex, see figure. It is clear from figure that the RTD becomes arrower as the value of the power law idex reduces. This meas that backmixig reduces as the fluid becomes more shear thiig (lower power law idex). Figure cofirms what we have already see i the velocity profiles of figure. As the fluid becomes more shear thiig, the velocity profile becomes closer to plug flow ad, cosequetly, the RTD becomes arrower ad backmixig more problematic. RTD i Screw Extruders Pito ad Tadmor (5) developed expressios for the RTD i sigle screw extruders. The cumulative RTD fuctio ca be writte as: ( ξ + ( ξ ) + ξ 3 ) F( t) = F( ξ ) =.5 3 ξ (3) The dimesioless time θ (time divided by mea residece time) ca be expressed as a fuctio of the dimesioless ormal coordiate ξ: θ 3ξ + 3 6ξ ( ξ + + ξ 3ξ = (4) + ξ 3ξ ) The two expressios above were derived for a Newtoia fluid usig the flat plate approximatio cosiderig both dow- ad cross-chael velocity compoets. Figure 3 shows the RTD for a sigle screw extruder as well as for pressure flow of a Newtoia fluid betwee flat plates. The sigle screw RTD is arrower tha the flat plate RTD because of the recirculatio of the fluid i the screw chael. Fluid speds more time i the lower portio of the chael ξ=-/3 tha i the upper portio of the chael ξ=/3-. I order to properly determie the RTD of a extruder we have to cosider ot oly dow- ad cross-chael velocity compoets but also ormal velocity compoets that occur at the flight flaks. This will require a umerical aalysis, either FDA, FEA, or BEA. Eve

though the depth of the chael is usually quite small compared to the chael width, the residece time at the flight flaks is substatial because the ormal velocities are quite small. Joo ad Kwo (6) poited out limitatios of the Pito aalysis. As oe would expect, the Pito model uderpredicts the residece times relative to a full three-dimesioal aalysis, particularly with large values of the axial pressure gradiet. The Pito RTD for a sigle screw extruder is arrower tha the RTD for pressure flow betwee flat plate for a Newtoia fluid. The effect of shear thiig is to further arrow the RTD as discussed earlier. These two effects explai why backmixig is such a critical issue i screw extruders. Methods to Improve Backmixig A major cocer i backmixig is the fluid i the ceter regio of the screw chael where the axial shear strai is zero or close to zero. I a simple coveyig screw the fluid i the ier recirculatio regio will stay withi this regio util it reaches the ed of the screw. Whe this happes, the material flowig ito the die will be poorly mixed. Mixig pis ad slots i the screw flights will improve axial mixig because they achieve a short term splittig ad reorietatio of the fluid. The effect of mixig pis o backmixig is illustrated i figure 4. These results were obtaied usig a three-dimesioal BEM flow aalysis. It is clear that oe row of mixig pis has limited effect o axial mixig. Backmixig ca be improved by varyig the spacig betwee the pis, thus, itetioally creatig streams of differet axial velocities. The challege i improvig axial mixig is to efficietly trasfer fluid from the ier recirculatio regio to the outer regio ad vice versa. A simple but effective method of doig this is the iside-out mixer show i figure 5. The flight i this mixer is offset so that the material i the ceter regio is cut by the offset flight ad the pushed to the screw ad barrel surfaces by the ormal pressure gradiets that occur at the flight flak. Results of particle trackig usig BEA are show i figure 6. The redistributio of the material is show i figure 7. It shows how the fluid from the ceter regio is cut by the offset flight ad pushed to screw surface at the pushig side of the flight ad to the barrel surface at the trailig side of the flight. Figure 7 shows the redistributio of the fluid at several axial locatios from the flight offset. It is clear that a cosiderable axial distace is ecessary to brig about the redistributio of the material. Coclusios Backmixig is oe of the most difficult mixig tasks i screw extruders. This is due to the fact that the axial strai rates i the extrusio process are very low, particularly i the ceter regio of the chael. The axial velocity profile is close to plug flow, particularly for strogly shear-thiig fluids. Backmixig problems are especially severe whe addig a small percetage cocetrate to the extruder. Whe both the polymer ad cocetrate are i pellet form, the iitial striatio thickess ca be of the order of mm. If a fial striatio thickess of micro is required, the axial mixig has to achieve a reductio of striatio thickess of at least 5 orders of magitude. Cosiderig that the axial shear rate i melt coveyig is close to zero i the ceter regio of the chael, it is clear that axial mixig will be isufficiet uless efficiet mixig devices are used. The most efficiet way to improve axial mixig is to redistribute material from the ceter of the chael to the outer regio of the chael ad vice versa. Oe mixer that aims to achieve such redistributio is the iside-out mixer. Aother mixer that was desiged specifically to improve axial mixig is the CRD7 mixer, see figure 8. Aother method of reducig mixig problems with color cocetrates is to reduce the iitial striatio thickess. This ca be doe by reducig both the virgi plastic pellet size ad the color cocetrate pellet size. Graules are better tha pellets ad powder is better tha graules from a mixig poit of view. Smaller particle sizes may lead to other problems though, such as coveyig problems ad air etrapmet. Addig the colorat i liquid form ca also reduce the iitial striatio thickess. This is oe of the mai reaso liquid colorats are used. However, they ca create problems by formig a lubricatig layer o the barrel surface ad reducig the coveyig efficiecy of the extruder. Refereces. C. Rauwedaal, Polymer Mixig, A Self-Study Guide, Haser Publishers, Muich (998). C. Rauwedaal, T. A. Osswald, G. Tellez, ad P. J. Grama, Flow Aalysis i Screw Extruders Effect of Kiematic Coditios, Iteratioal Polymer Processig, XIII, 4, 37-333 (998) 3. Z. Tadmor ad C.G. Gogos, Priciples of Polymer Processig, Joh Wiley ad Sos, New York (979) 4. R.B. Bird, W.E. Stewart, ad E.N. Lightfoot, Trasport Pheomea, Wiley, New York (96) 5. G. Pito ad Z. Tadmor, Polym. Eg. Sci.,, 79 (97) 6. J.W. Joo ad T.H. Kwo, Polym. Eg. Sci., 33, 5, 959 (993)

Velocity Distributio Power Law Fluid. Dimesioless Velocity.8.6.4 =. =.5 =.5. -.5 - -.5.5.5 Dimesioless Normal Coordiate Figure, Velocity profiles for various power law idex values RTD of Power Law Fluid i Flat Plate Pressure Flow.9 Cumulative RTD Fuctio.8.7.6.5.4.3.. =.5 =. =.5.5.5.5 3 Dimesioless Time Figure, Residece time distributio curves for several power law idex values RTD Sigle Screw Extruder.9 Cumulative RTD Fuctio.8.7.6.5.4.3.. Sigle Screw Extruder Flat Plate Newtoia.5.5.5 3 Dimesioless Time Figure 3, RTD of sigle screw extruder ad flat plate

Figure 4, Particle trackig results i a mixig sectio with elogatioal mixig pis Figure 5, Solid model of Iside-Out mixer

Figure 6, Particle trackig i the Iside-Out mixer Cross sectio Cross sectio Figure 7, Redistributio i the Iside-Out mixer Trasport directio Figure 8, The CRD7 mixer