Homogeneity of multilayers produced with a Static Mixer

Size: px
Start display at page:

Download "Homogeneity of multilayers produced with a Static Mixer"

Transcription

1 Homogeneity of multilayers produced with a Static Mixer J.C.vanderHoeven,R. Wimberger-Friedl,H.E.H.Meijer Philips Research Laboratories, Prof. Holstlaan 4, 5656 AA Eindhoven, The Netherlands DPI, P.O. Box 513; 5600 MB Eindhoven, The Netherlands Abstract: A multiflux static mixer can be used to produce multilayered structures. The flow is repeatedly cut, stretched and stacked by mixing elements in the channel of such a device. In the standard design, however, the obtained layer thicknesses are inhomogeneous. The causes for the multiflux static mixer s deviation from ideal behaviour are identified by 3D numerical simulations as unequal pressure drops in the separating flows. Changes in the arrangement of the elements are proposed and their effects are verified by simulations and experiments. A significant improvement of the layer homogeneity is achieved by introducing additional elements with separating walls at the inlets and at the outlets of the mixing elements. Keywords Static mixer, multilayer, layer homogeneity, extrusion, finite element analysis, Fidap, SEPRAN, layer thickness distribution Introduction To achieve the most efficient mixing process when mixing highly viscous fluids, the application of the bakers transformation is essential. In this transformation a fluid is successively stretched, cut and stacked, in the same way that a baker rolls and folds his dough. The bakers transformation can be realized in a continuous process, for instance by using a multiflux static mixer. In the channel of this device, several mixing elements are present, in each of which the flow undergoes a bakers transformation. The geometry of the multiflux static mixer was first described by Sluijters [1]. It is the most compact design one can think of for this purpose. It has a low pressure drop and is easily produced at low cost. The operating principle of a mixing element is demonstrated in figure, showing several cross sections through this element, using a white and a coloured fluid. In cross section a the two input flows are kept separate by a horizontal wall, in cross section b the flows have been stacked and cut in two by a vertical wall. In cross sections c, d, e and f the flows are stretched and finally stacked (acceleration and deceleration), resulting in a four layered flow. In a second element, the number of layers can be doubled to eight and so on. The final number of layers equals m 2 n,wheremis the number of layers fed into the mixer and n is the number 1

2 of mixing elements used. Van der Sanden et.al. [2], Mueller et.al. [3] and Baer et.al. [4] used the multiflux static mixer to make multilayers for research purposes. Schilo and Ostertag [5] experimentally compared the mixing performance of the multiflux static mixer and the Kenics static mixer. Numerical investigations into the Kenics mixer performance have been done by among others: Bryde and Sawley [6], Hobbs et.al. [7], Hobbs and Muzzio [8, 9, 10], Avalosse and Crochet [11], Ling and Zhang [12] and Arimond and Erwin [13]. The Ross mixer, see Schrenk [14] for instance, is yet another static mixer design. In stead of cutting the flow in two flows, the Ross mixer cuts the flow in four before stretching and stacking. Schrenk [15] used a Ross mixer to make multilayers with sub-micrometer thicknesses which had interesting optical properties. Ebeling et.al. [16] used co-extruded multilayer sheets, supplied by The Dow Chemical Company, for their research into the effect of the peel rate and temperature on the delamination toughness of polycarbonate styrene-acrylonitrile micro-layers. It is important to have control over the layer thicknesses and to have straight interfaces for such applications. This is, however, not generally achieved. Both, the thicknesses of the individual layers, as well as the thicknesses of the layers relative to each other, are inhomogeneous. This investigation studies relevant processes and parameters which determine the layer distribution in the extrudate of the multiflux static mixer. It is not the mixing quality that is considered, but the ability of the mixer to produce multilayers with straight interfaces and layers with equal thicknesses. Minagawa and White [17] performed an experimental study on the influence of the viscosity ratio and the die cross-section on the interface shape in co-extrusion. Schrenk et.al. [18] used a critical interfacial shear stress criterion for the prediction of the onset of interfacial instabilities. Mavridis [19] added to this a contribution of the elasticity difference between the layers. Wilson and Khomami [20, 21] studied interfacial instabilities in co-extrusion of polypropylene and high-density polyethylene. The problems described above are less relevant for the present investigation, since the materials used here have equal rheological properties. This approach, in which interfacial tension effects will also be ignored, facilitates a study of the basic functioning of the multiflux mixer. The flow in the static mixer was simulated with 2D and 3D numerical models, indicating several causes of the layer inhomogeneities. The layers in the static mixer were visualised and compared with experimental results for validation of the models. Ways of improving the homogeneity of the layer thickness distribution, henceforth referred to as the layer distribution, by changing the geometry of the mixer channel, will be investigated, both, numerically and experimentally. Experimental The experimental setup is equipped with two velocity controlled pistons to feed the mixer. In this way, the fluxes can be controlled accurately, relative to feeding the mixer with extruders. The cylinders and the mixer are provided with heating elements and thermocouples for controlling the temperature of the polymer melt. The mixer is assembled from separate rectangular steel blocks that leave a quadrilateral channel of mm 3. This channel can be filled with separate mixing elements of 1 cm 3 volume 2

3 which are supported by the channel walls only. This allows a random variation of all kinds of elements within the length of the channel. Different extrusion dies can be used at the exit. The polymers used in this work were grey and white polycarbonate, PC, (Lexan 121R, General Electric Plastics) and blue and white ABS (Ronfalin FA-50, DSM). Figure 2 depicts the viscosities of the materials as determined using a plate-plate rheometer (Rheometric Scientific SR-5000) under the indicated conditions. Before insertion into the feeding cylinders of the static mixer, the polymers were compression molded into rods. PC was extruded at 280 C and ABS at 240 C, both with a total volume flux of 0.38 cm 3 s 1. ABS can be used for a frozen flow investigation by removing the material from the mixing channel in a non-destructive way. As ABS is rubbery at around 90 C, the multilayered flow can be isolated by disassembling the static mixer at the end of an experiment. The frozen flow can be sliced in order to study the thickness of the layers and the interface shapes. PC is used for its Newtonian behaviour in the low shear regions. The experiments in which PC is used will closely resemble the 3D numerical simulations, in which a constant viscosity is used. As the frozen PC flow could not be isolated from the mixer, only the layer distribution in the PC extrudate of the mixer is studied in the experiments. Numerical simulations For numerical models with two or more fluids with different rheological properties, techniques described by Galaktionov et.al. [22], Karagiannis et.al. [23], Mavridis et.al. [24] and Gifford [25], amongst others, may be useful. Calculation of moving boundaries using the pseudo concentration method is reported by, for instance, Haagh [26] and Thompson [27]. In the experiments performed in this work, however, the fluids are two polymer melts, which only have different colours. Having confirmed that the rheological properties of the melts are not influenced by the pigments, the flow in the static mixer is modelled with a single fluid only. The following assumptions are used: Both fluids have the same rheological properties The fluids do not react or diffuse There is no interfacial tension The flow is isothermal There is no slip along solid walls Gravity can be neglected Assuming a parabolical velocity profile, the velocity of the fluid in the smallest cross section of the mixer has a maximum value U max =21 mm s 1. The minimum height in the mixer channel is H min = 4.5 mm, the density of the melts is ρ 1200 kg m 3 and the viscosities of the polymer melts are η PC 2000 Pa s at 240 Candη ABS Pa s at 180 C (see also figure 2). This results in maximum values of the Reynolds numbers Re and for the PC and the ABS melt, respectively, and inertia is consequently negligible. The governing equations are given by: η u = p 3 (1)

4 u=0 (2) These equations were solved with the finite element analysis packages Fidap 1 and SEPRAN 2. In the Fidap models, 9 node quadrilateral elements are used for the 2D models and 27 node bricks for the 3D models of the static mixer. The 3D SEPRAN model uses spectral elements. The Powerlaw, Carreau and Bingham viscosity models used in the 2D Fidap simulations are given by equation 3, 4 and 5, respectively. η( γ)= η 0 γ n 1 (3) η( γ)= η 0 (1+λ 2 γ 2 ) 1 n 2 (4) τ = η γ + τ y τ τ y (5) Viscous heating is estimated using equation 6. When the pressure drop over the mixer with three mixing elements is estimated by using the model given by Sluiters [1] and the specific heat C p = J kg 1 K 1, the temperature increase of the melt for PC is 0.1 K. Viscous heating is therefore neglected. P v = ρ C p T v T = P ρ C p (6) The layers in the Fidap model are visualised with the help of particle tracking. First the velocity field is calculated. In a separate calculation, the trajectories of 640 mass and size-less particles, released at the inlets, are determined. The velocity field in the mixer is thus not influenced by the particles. The layer distribution in any cross section through the mixer channel can be visualised by calculating the intersections of the particle trajectories with the cross section plane. The particles released in the upper inlet can be given a different colour from the particles released in the lower inlet, in order to visualise the layers. Related methods were used by Bryde and Sawley [6], Hobbs and Muzzio [8, 9, 10], Avalosse and Crochet [11], Ling and Zhang [12], and many others. In the SEPRAN model an adaptive front tracking technique was implemented. Details of this technique are described by Galaktionov et.al. [22]. Results and discussion Introduction In the use of the multiflux device for the production of regular multilayered stacks, the layer distribution after leaving the extrusion die is important. The transition from a pressure driven flow to a free flow affects the layer distribution. Similarly, the transitions from separated to recombined flows at the ends of separating walls parallel to the layer-stack are also of importance. Figure 3 shows the simulation of a 2D multilayered flow. Obviously, the homogeneous layer distributions in the two inlets are disturbed by the recombination (see layer distribution at cross section 1) and again by 1 Supported by FLUENT, 2 Ingenieursburo SEPRA, Leidschendam, The Netherlands, 4

5 the extrusion process 3 (see cross section 2). From this simple calculation it can be concluded that the layer distribution in the extrudate of the static mixer can not be made homogeneous when the layer distribution in the mixer channel is homogeneous. In an extruded homogeneous layer-stack all layers have equal cross section surfaces. Also, the velocity is constant over the cross section. In other words, the volume fluxes of the individual layers are equal. The layer distribution in the extrudate of the multiflux static mixer should thus be made homogeneous by making the volume fluxes of the individual layers in the mixer equal. Typical mixer configuration Figure 4 shows, both, the FEM mesh of the static mixer as used in the Fidap simulations and the typical mixer configuration. The two inlet flows are separated by a horizontal wall. This element of the static mixer is called an H-element. Subsequently, there are two types of mixing elements. In the first element, the flow at the righthand side from top to bottom is led downwards (and vice versa at the lefthand side), while in the second element the flow at the lefthand side is led downwards. The element types are abbreviated by E R and E L, respectively. Before leaving the static mixer, the flow is led through an open element, referred to as an O-element. The configuration of the static mixer is denoted by giving the abbreviations of the adjacent element types. For example, the static mixer in figure 4 has an H-E R E L E R -O configuration. In the cross section between the inlet H-element and the first mixing element (cross section 1, see figure 5), four openings of equal size can be distinguished. It is also clear, however, that when the cross sections perpendicular to the flow velocity are considered, there are two larger openings (perpendicular to the mixer axis) and two smaller openings (not perpendicular to the mixer axis). The consequences of this are shown in figure 6: in cross section 1 the volume fluxes through the smaller openings are smaller, resulting in unequal layer thicknesses in the extrudate. The four openings perpendicular to the flow velocity related to cross section 2 are of equal size (see figure 5). This is the result of the L and R-type alternation of adjacent elements. In the model, the calculated volume fluxes through these openings are equal. Figure 6 demonstrates that the differences between the flow velocities through the four openings in cross section 1 are much larger than in cross sections 2 and 3. If all elements in the mixer were of either L or R-type only, the homogeneity of the final layer thickness distribution would be much worse. Figure 6 also demonstrates that the shapes of the four velocity profiles in cross section 2 (and cross section 3) are not equal. This could also result in layer thickness differences. The velocity profiles of the flows entering the O-element (cross section 4) are asymmetrical as well. This causes a kind of twisting behaviour of the extrudate in experiments. Figure 7 shows a cross section through the extrudate of the static mixer in the H-E R E L E R -O configuration, obtained from extrusion of white and grey PC without a die. Apart from the deformed cross-section due to extrudate swell (the flow does not solidify immediately after extrusion from the square channel), the photo demonstrates that the layer distribution is not only strongly inhomogeneous, but is twisted as well. 3 The extrusion process is modelled as a simple transition from no slip at the wall to full wall slip. 5

6 Figures 8a and 9a show cross sections through the frozen ABS flows, isolated from the mixer, halfway through the first and the second mixing elements, respectively, showing the layer distributions. In figures 8b and 9b the layers in the same cross sections visualised with particle tracking in the 3D Fidap model of the mixer are given. The layer distributions calculated with the SEPRAN model are given in figures 8c and 9c. The differences between the experimental and the numerical layer distributions are remarkably small. Clearly, the numerical simulations using a constant viscosity, give very good predictions of the layer distribution in the mixer channel. Improved mixer configurations Considering the results shown so far, an obvious way of improving the layer distribution is by regulating the asymmetry in the flow-fields of the separate arms of the elements. The size of the bigger openings in cross section 1 of the static mixer (figures 4 to 6) can be adapted, for instance with a so-called extra resistance. Figure 10 shows a 2D FEM model of the static mixer with an extra resistance, representing a side view of the inlet H-element and a part of the first mixing element. Equal volume fluxes through the upper and lower inlet, Q 1 /Q 2 = 1, are desired at equal inlet pressures p 1 =p 2. In figure 11, the values of the calculated flux ratios Q 1 /Q 2 are plotted for several heights of the resistance h r /h mix for a constant viscosity and for viscosities modelled with the powerlaw (see equation 3). It is clear from this result that the height of the extra resistance h r necessary to result in equal fluxes depends on the rheological properties of the fluid. Thus, such a design is not a robust solution. Another way of balancing the fluxes is by shifting the inclining, compressing wall further downstream, in this way increasing the effective cross section of the smaller opening. Practically, this can be achieved by inserting an extra element with a vertical midplane (called an I-element) between the inlet H-element and the first mixing element. Figure 12 shows a 2D FEM model of this arrangement. In figure 13, the calculated flux ratios Q 1 /Q 2 at equal inlet pressures p 1 =p 2 are plotted for several lengths l I /h mix of the I-element with a constant viscosity, two powerlaw viscosities, a Carreau viscosity model (see equation 4) and a Bingham viscosity model (see equation 5). Regardless of the rheological properties of the fluids, the fluxes become equal when the length of the I-element is at least one times the height of a mixing element l I /h mix 1. The same modification can be used to avoid the twisting of the extrudate at the exit. When an H-element is used between the last mixing element and the O-element, the velocity profiles can develop into symmetrical profiles before the two flows are stacked. The length of this H-element should also be at least equal to the height of the mixing element l H /h mix 1. Figure 14 shows a cross section through the extrudate of the static mixer in the H-I-E L E R E L -H-O configuration, obtained from extrusion of white and grey PC without a die. The use of the I and the H-element clearly improves the layer quality. The layer distribution is still not completely homogeneous, however, but is much better than in the extrudate of the mixer in the H-E R E L E R -O configuration (compare figures 7 and 14). In addition, the layers are less twisted. The remaining deviations from an ideal distribution may be explained by the differences in the shapes of the velocity profiles in the cross sections between adjacent mixing elements (see also figure 6). The 6

7 influence of this on the layer homogeneity is demonstrated by placing an H-element at the outlet of each mixing element. The outflows of the mixing elements can develop, just like in the H-element between the last mixing element and the O-element of the mixer in the H-I-E L E R E L -H-O configuration. When this is done, I-elements should be used in front of all mixing elements, just like between the inlet H-element and the first mixing element of the mixer in the H-I-E L E R E L -H-O configuration. Note that in this H-IEH-IEH- IEH-O configuration there is no need to use alternating L and R-type mixing elements. Figure 15 shows the FEM mesh of the static mixer in the H-IEH-IE configuration. Figure 16 demonstrates that the shapes of the velocity profiles of the flow through the four openings on cross section 1 and 2, where the actual cutting takes place, are equal, indicating equal volume-flux divisions. Figure 17 shows a photograph of a cross section through the frozen outflow of the mixer in the H-IEH-IEH-IEH-O configuration, obtained from extrusion of white and grey PC without a die. Obviously, the symmetry of the velocity profiles are important for a perfect layer homogeneity, as the homogeneity of the layer distribution has improved with respect to the mixer in the H-I-E L E R E L -H-O configuration. A disadvantage of the H-IEH- IEH-IEH-O configuration is that it makes the mixer three times longer, compactness being one of the great advantages of the static mixer. The H-I-E L E R E L -H-O configuration can be used as a good compromise, as the layer quality is considerably better and the mixer length only increases by two element lengths. Conclusions The layer distribution in the mixer is different from the layer distribution in the extrudate, as a result of the difference between the velocity profiles in the static mixer and in the extrudate. It should be stressed that the layer distribution in the extrudate is determined by the volume fluxes of the individual layers in the static mixer. When these volume fluxes are unequal, the layer distribution in the extrudate will also be unequal. The flux division can be corrected by adding extra elements with a vertical wall in front of the mixing elements, and adding extra elements with a horizontal wall behind the mixing elements. The length of these extra elements should at least be equal to the total height of a mixing element. A disadvantage is that the improved mixer becomes three times longer than the standard mixer. A more compact modification is obtained when an element with a vertical wall is added before the first mixing element, an element with a horizontal wall is added behind the last mixing element and adjacent mixing elements are alternatingly of the L and R-type. However, the improvement of the layer distribution is somewhat less than that obtained with the former modification of the mixer geometry. Acknowledgements The authors would like to thank Hans de Bruin, who built the experimental setup and performed all experiments, and Peter Kruijt and Patrick Anderson who did the SEPRAN simulations. The help and advice of Wim Zoetelief and Ramin Badie are also gratefully acknowledged. 7

8 References [1] R. Sluijters, Chemische Techniek, 3, 33 (1965) [2] M.C.M. van der Sanden, L.G.C. Buijs, F.O. de Bie and H.E.H. Meijer, Polymer, 35, 2783 (1994) [3] C.D. Mueller, S. Nazarenko, T. Ebeling, T.L. Schuman, A. Hiltner and E. Baer, Polym. Eng. Sci., 37, 355 (1997) [4] E. Baer, D. Jarus, and A. Hiltner, Proc. SPE ANTEC Tech. Papers, 3947 (1999) [5] D. Schilo and K. Ostertag, Verfahrenstechnik, 6, 45 (1972) [6] O. Bryde and M.L. Sawley, Computers & Fluids, 28, 1 (1999) [7] D.M. Hobbs, P.D. Swanson and F.J. Muzzio, Chem. Eng. Sci., 53, 1565 (1998) [8] D.M. Hobbs and F.J. Muzzio, Fluid Mechanics and Transport Phenomena, 43, 3121 (1997) [9] D.M. Hobbs and F.J. Muzzio, Chem. Eng. Sci., 53, 3199 (1998) [10] D.M. Hobbs and F.J. Muzzio, AIChE Journal, 43, 3121 (1997) [11] Th. Avalosse and M.J. Crochet, AIChE Journal, 43, 588 (1997) [12] F.H. Ling and X. Zhang, Chem. Eng. Comm., 136, 119 (1995) [13] J. Arimond and L. Erwin, Chem. Eng. Commun., 37, 105 (1985) [14] W.J. Schrenk, R.K. Shastri, R.F. Ayres and D.J. Gosen, U.S. Patent (1992) [15] W.J. Schrenk, U.S. Patent (1993) [16] T. Ebeling, A. Hiltner and E. Baer, Polymer, 40, 1525 (1999) [17] N. Minagawa and J.L. White, Polym. Eng. Sci., 15, 825 (1975) [18] W.J. Schrenk, N.L. Bradley, Jr. T. Alfrey and H. Maack, Polym. Eng. Sci., 18, 620 (1978) [19] H. Mavridis and R.N. Shroff, Polym. Eng. Sci., 34, 559 (1994) [20] G.M. Wilson and B. Khomami, J. Non-Newtonian Fluid Mech., 45, 355 (1992) [21] G.M. Wilson and B. Khomami, J. Rheology, 37, 315 (1993) [22] A.S. Galaktionov, P.D. Anderson and G.W.M., Peters, Vol of Lecture Notes in Computer Science, (1997) [23] A. Karagiannis, A.N. Hrymak and J. Vlachopoulos, Rheol. Acta, 29, 71 (1990) [24] H. Mavridis, A.N. Hrymak and J. Vlachopoulos, AIChE Journal, 33, 410 (1987) [25] W.A. Gifford, Polym. Eng. Sci., 37, 315 (1997) [26] G.A.A.V. Haagh and F.N. van de Vosse, Int. J. Numer. Methods Fluids, 28, 1355 (1998) [27] E. Thompson, Int. J. Numer. Methods Fluids, 6, 749 (1986) 8

9 a. b. c. d. e. f. Figure 1: Several cross sections through the first element of the static mixer demonstrating its working principle by using a stratified white and coloured fluid. 1.0E E+05 η* [Pa s] PC (Grey) PC (White) ABS (Blue) ABS (White) 1.0E E E E E E E E+03 ω [rad/s] Figure 2: The melt viscosities measured with a Rheometric Scientific SR-5000 in parallel plate configuration, at 240 C for grey and white PC and at 180 C for blue and white ABS. 9

10 Inlet 2 Inlet 1 Y Free outflow X section 1 section 2 Figure 3: 2D simulation demonstrating the influence of stacking and extruding on the homogeneous layer distributions at the inlets. 10

11 Inlets (H-element) section 1 st 1 element section 2 nd 2 element section 3 rd 3 element section 4 Outlet (O-element) Figure 4: The FEM mesh demonstrating the geometry of a static mixer in the H-E R E L E R -O configuration. H-element section 1 st 1 element section 2 nd 2 element section 3 rd 3 element section 4 O-element Figure 5: Several cross sections perpendicular to the flow velocity through the static mixer channel demonstrating the different sizes of the openings. H-element section 1 1 st element section 2 2 nd element section 3 3 rd element section 4 O-element Figure 6: Several cross sections through the static mixer with three elements showing the contours of the dimensionless velocity. Black contours are related to the highest velocities, white contours to zero velocity. 11

12 Figure 7: A photograph of a cross section through the extrudate of the static mixer in the H- E R E L E R -O configuration, obtained from extrusion without a die at 280 C using white and grey PC. 12

13 z[mm] y[mm] a. b. c. Figure 8: The layer distributions in cross sections through the middle of the first element of the static mixer: a: Isolated frozen ABS flow, b: Fidap particle tracking, c: SEPRAN adaptive front tracking z[mm] y[mm] a. b. c. Figure 9: The layer distributions in cross sections through the middle of the second element of the static mixer: a: Isolated frozen ABS flow, b: Fidap particle tracking, c: SEPRAN adaptive front tracking. 13

14 p Q 1 1 p 2 Q 2 h h r h 1 2 h mix Figure 10: A 2D model of the static mixer with an extra resistance adjusting the openings in the channel to obtain equal volume fluxes Q 1 and Q 2 through the upper and lower inlet at equal inlet pressures p 1 =p 2. 1HZWRQLDQ 3RZHUODZÃQÃ Ã 3RZHUODZÃQÃ Ã 3RZHUODZÃQÃ Ã 4 ÃÃ4 K U ÃÃK PL[ Figure 11: The calculated ratios of the volume fluxes through the upper inlet and the lower inlet Q 1 /Q 2 of the 2D static mixer model with an extra resistance for several heights of the resistance h r /h mix using different viscosity models. 14

15 `I p Q 1 1 p 2 Q 2 h H-element I-element Mixing-element h 1 2 h mix Figure 12: A 2D model of the static mixer with a vertical wall (I-element) of length l I between the two inlets and the mixing element in order to obtain equal volume fluxes Q 1 and Q 2 through the upper and lower inlet at equal inlet pressures p 1 =p 2. 4 ÃÃ4 1HZWRQLDQ 3RZHUODZÃQÃ Ã 3RZHUODZÃQÃ Ã &DUUHDX %LQJKDP O, ÃÃK PL[ Figure 13: The calculated ratios of the volume fluxes through the upper and lower inlet Q 1 /Q 2 of the 2D static mixer model for several lengths l I /h mix of the vertical wall (I-element) using different viscosity models (Carreau model with λ =0.25 and n=0.25, Bingham model with τ y = 2.7 [kpa]. η 0 = 2000 [Pa s]) 15

16 Figure 14: A photograph of a cross section through the extrudate of the static mixer in the H- I-E R E L E R -H-O configuration, obtained from extrusion without a die at 280 C using white and grey PC. 16

17 Inlets I-element st 1 element H-element I-element nd 2 element Figure 15: The FEM mesh of the static mixer with an I-element in front of each element and an H-element behind each element, the H-IEH-IE configuration. Inlets section 1 I-element 1 st element H-element section 2 I-element 2 nd element Figure 16: The velocity profiles on several cross sections through the static mixer in the H-IEH-IE configuration, demonstrating equal volume fluxes through the four openings on the cross sections X 1 and X 2 where the cutting of the flow takes place. The shapes of these four velocity profiles are also equal, indicating that the flows of the individual layers being cut are equally divided. 17

18 Figure 17: A photograph of a cross section through the extrudate of the static mixer in the H- IEH-IEH-IEH-O configuration, obtained from extrusion without a die at 280 C using white and grey PC. 18

Multilayer Rheology Effects in Coextruded Structure Design

Multilayer Rheology Effects in Coextruded Structure Design 2008 Best Paper Multilayer Rheology Effects in Coextruded Structure Design Print (10)» Best Papers» 2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion

More information

THE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES

THE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES THE 3D VISCOELASTIC SIMULATION OF MULTI-LAYER FLOW INSIDE FILM AND SHEET EXTRUSION DIES Kazuya Yokomizo 1, Makoto Iwamura 2 and Hideki Tomiyama 1 1 The Japan Steel Works, LTD., Hiroshima Research Laboratory,

More information

Modeling of Anisotropic Polymers during Extrusion

Modeling of Anisotropic Polymers during Extrusion Modeling of Anisotropic Polymers during Extrusion Modified on Friday, 01 May 2015 10:38 PM by mpieler Categorized as: Paper of the Month Modeling of Anisotropic Polymers during Extrusion Arash Ahmadzadegan,

More information

Determining the Processability of Multilayer Coextruded Structures

Determining the Processability of Multilayer Coextruded Structures Determining the Processability of Multilayer Coextruded Structures Joseph Dooley The Dow Chemical Company, Midland, MI ABSTRACT Multilayer coextrusion is a process in which two or more polymers are extruded

More information

The Effect of Material Selection in Feed-Block Coextrusion of a Three Layer Film

The Effect of Material Selection in Feed-Block Coextrusion of a Three Layer Film The Effect of Material Selection in Feed-Block Coextrusion of a Three Layer Film Jiri Vlcek 1, Walter Kopytko 2, Martin Zatloukal 2 and Jiri Svabik 1 1. Compuplast International Inc., Zlin, Czech Republic

More information

Analysis and Optimization of Low-Pressure Drop Static Mixers

Analysis and Optimization of Low-Pressure Drop Static Mixers FLUID MECHANICS AND TRANSPORT PHENOMENA Analysis and Optimization of Low-Pressure Drop Static Mixers Mrityunjay K. Singh, Tae G. Kang, Patrick D. Anderson, and Han E. H. Meijer Materials Technology, Eindhoven

More information

CONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION

CONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION Second International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 6-8 December 1999 CONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION

More information

2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion Process

2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion Process 2009 Best Paper Understanding and Quantification of Die Drool Phenomenon During Polypropylene Extrusion Process Print (10)» 2010 Best Paper An Engineering Approach to the Correction of Rotational Flow

More information

Simulations with Dispersive Mixtures Model

Simulations with Dispersive Mixtures Model Simulations with Dispersive Mixtures Model P.J.A. Janssen, P.D. Anderson & G.W.M. Peters MT03.05 May 12, 2003 1 Introduction The combination of two different materials in a blend can yield interesting

More information

EFFECT OF TYPICAL MELT TEMPERATURE NON-UNIFORMITY ON FLOW DISTRIBUTION IN FLAT DIES

EFFECT OF TYPICAL MELT TEMPERATURE NON-UNIFORMITY ON FLOW DISTRIBUTION IN FLAT DIES EFFEC OF YPICAL MEL EMPERAURE NON-UNIFORMIY ON FLOW DISRIBUION IN FLA DIES Olivier Catherine, Cloeren Incorporated, Orange, X Abstract In this study, the influence of non-uniform incoming melt temperature

More information

SIMULATION OF POLYMERIC FLOW IN A TWIN-SCREW EXTRUDER: AN ANALYSIS OF ELONGATIONAL VISCOSITY EFFECTS

SIMULATION OF POLYMERIC FLOW IN A TWIN-SCREW EXTRUDER: AN ANALYSIS OF ELONGATIONAL VISCOSITY EFFECTS SIMULATION OF POLYMERIC FLOW IN A TWIN-SCREW EXTRUDER: AN ANALYSIS OF ELONGATIONAL VISCOSITY EFFECTS A. Shah and M. Gupta Mechanical Engineering-Engineering Mechanics Department Michigan Technological

More information

ANALYSIS ON PLANAR ENTRY CONVERGING FLOW OF POLYMER MELTS

ANALYSIS ON PLANAR ENTRY CONVERGING FLOW OF POLYMER MELTS Journal of Materials Science and Engineering with Advanced Technology Volume 2, Number 2, 2010, Pages 217-233 ANALYSIS ON PLANAR ENTRY CONVERGING FLOW OF POLYMER MELTS College of Industrial Equipment and

More information

TWO-DIMENSIONAL SIMULATIONS OF THE EFFECT OF THE RESERVOIR REGION ON THE PRESSURE OSCILLATIONS OBSERVED IN THE STICK-SLIP INSTABILITY REGIME

TWO-DIMENSIONAL SIMULATIONS OF THE EFFECT OF THE RESERVOIR REGION ON THE PRESSURE OSCILLATIONS OBSERVED IN THE STICK-SLIP INSTABILITY REGIME 1 TWO-DIMENSIONAL SIMULATIONS OF THE EFFECT OF THE RESERVOIR REGION ON THE PRESSURE OSCILLATIONS OBSERVED IN THE STICK-SLIP INSTABILITY REGIME Eleni Taliadorou and Georgios Georgiou * Department of Mathematics

More information

Measuring Slip at a Polymer/Polymer Interface during Three-Layer Flow

Measuring Slip at a Polymer/Polymer Interface during Three-Layer Flow Rapid Communication Nihon Reoroji Gakkaishi Vol.41, No.4, 235~239 (Journal of the Society of Rheology, Japan) 2013 The Society of Rheology, Japan Measuring Slip at a Polymer/Polymer Interface during Three-Layer

More information

Excerpt from the Proceedings of the COMSOL Users Conference 2006 Boston

Excerpt from the Proceedings of the COMSOL Users Conference 2006 Boston Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu

More information

Morphology Evolution in PS/LDPE Blends in a Twin Screw Extruder: Effects of Compatibilizer

Morphology Evolution in PS/LDPE Blends in a Twin Screw Extruder: Effects of Compatibilizer Korean J. Chem. Eng., 18(1), 33-39 (2001) Morphology Evolution in PS/LDPE Blends in a Twin Screw Extruder: Effects of Compatibilizer Do Young Moon*, Moo Hyun Kwon and O Ok Park *Chemical Division R&D Center,

More information

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická

More information

Shell Balances in Fluid Mechanics

Shell Balances in Fluid Mechanics Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell

More information

A NEW DISPERSIVE AND DISTRIBUTIVE STATIC MIXER FOR THE COMPOUNDING OF HIGHLY VISCOUS MATERIALS

A NEW DISPERSIVE AND DISTRIBUTIVE STATIC MIXER FOR THE COMPOUNDING OF HIGHLY VISCOUS MATERIALS A NEW DISPERSIVE AND DISTRIBUTIVE STATIC MIXER FOR THE COMPOUNDING OF HIGHLY VISCOUS MATERIALS Paul Gramann and Bruce Davis, The Madison Group: PPRC. Tim Osswald, University of Wisconsin-Madison Chris

More information

Particle removal in linear shear flow: model prediction and experimental validation

Particle removal in linear shear flow: model prediction and experimental validation Particle removal in linear shear flow: model prediction and experimental validation M.L. Zoeteweij, J.C.J. van der Donck and R. Versluis TNO Science and Industry, P.O. Box 155, 600 AD Delft, The Netherlands

More information

Microfluidics 1 Basics, Laminar flow, shear and flow profiles

Microfluidics 1 Basics, Laminar flow, shear and flow profiles MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,

More information

Lecture 3. Properties of Fluids 11/01/2017. There are thermodynamic properties of fluids like:

Lecture 3. Properties of Fluids 11/01/2017. There are thermodynamic properties of fluids like: 11/01/2017 Lecture 3 Properties of Fluids There are thermodynamic properties of fluids like: Pressure, p (N/m 2 ) or [ML -1 T -2 ], Density, ρ (kg/m 3 ) or [ML -3 ], Specific weight, γ = ρg (N/m 3 ) or

More information

Rheological Properties of ABS at Low Shear Rates: Effects of Phase Heterogeneity

Rheological Properties of ABS at Low Shear Rates: Effects of Phase Heterogeneity Malaysian Polymer Journal, Vol 4, No, p9-36, 9 Available online at wwwfkkksautmmy/mpj Rheological Properties of ABS at Low Shear Rates: Effects of Phase Heterogeneity Asif Ali Qaiser *, Yasir Qayyum and

More information

AN ANALYSIS OF THE EFFECT OF ELONGATIONAL VISCOSITY ONTHEFLOWINAFLATDIE

AN ANALYSIS OF THE EFFECT OF ELONGATIONAL VISCOSITY ONTHEFLOWINAFLATDIE AN ANALYSIS OF THE EFFECT OF ELONGATIONAL VISCOSITY ONTHEFLOWINAFLATDIE Y. Sun and M. Gupta Mechanical Engineering-Engineering Mechanics Department Michigan Technological University Houghton, MI 49931

More information

Final Polymer Processing

Final Polymer Processing 030319 Final Polymer Processing I) Blow molding is used to produce plastic bottles and a blow molding machine was seen during the Equistar tour. In blow molding a tubular parison is produced by extrusion

More information

Validation 3. Laminar Flow Around a Circular Cylinder

Validation 3. Laminar Flow Around a Circular Cylinder Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous

More information

2 Navier-Stokes Equations

2 Navier-Stokes Equations 1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1

More information

OPTICAL PROPERTY INVESTIGATION IN SEQUENTIAL MULTI-COMPONENT MOLDING

OPTICAL PROPERTY INVESTIGATION IN SEQUENTIAL MULTI-COMPONENT MOLDING OPTICAL PROPERTY INVESTIGATION IN SEQUENTIAL MULTI-COMPONENT MOLDING Chao-Tsai Huang 1, Meng-Chih Chen 1, Yuan-Rong Chang 1, Wen-Li Yang 1, Shi-Chang Tseng 2,and Huang-Yi Chang 2, 1. CoreTech System (Moldex3D)

More information

Corrections to flow data in polymer melts

Corrections to flow data in polymer melts Corrections to flow data in polymer melts Narongrit Sombatsompop Polymer PROcessing and Flow (P-PROF) Materials Technology, School of Energy & Materials King Mongkut s University of Technology Thonburi

More information

Analysis of Heat Transfer in Pipe with Twisted Tape Inserts

Analysis of Heat Transfer in Pipe with Twisted Tape Inserts Proceedings of the 2 nd International Conference on Fluid Flow, Heat and Mass Transfer Ottawa, Ontario, Canada, April 30 May 1, 2015 Paper No. 143 Analysis of Heat Transfer in Pipe with Twisted Tape Inserts

More information

ANALYSIS OF EXTRUSION PROBLEMS WITH FLOW SIMULATION

ANALYSIS OF EXTRUSION PROBLEMS WITH FLOW SIMULATION Compuplast International Inc. 2010 ANALYSIS OF EXTRUSION PROBLEMS WITH FLOW SIMULATION Vijay Kudchadkar Compuplast International Inc. Objectives Demonstrate the application of CAE/CFD/Simulation towards

More information

Flow Induced Molecular Weight Fractionation during Capillary Flow of Linear Polymer Melt

Flow Induced Molecular Weight Fractionation during Capillary Flow of Linear Polymer Melt Flow Induced Molecular Weight Fractionation during Capillary Flow of Linear Polymer Melt JAN MUSIL a,b and MARTIN ZATLOUKAL a,b a Centre of Polymer Systems, University Institute Tomas Bata University in

More information

On Relationship between PVT and Rheological Measurements of Polymer Melts

On Relationship between PVT and Rheological Measurements of Polymer Melts ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 3, 2005 On Relationship between PVT and Rheological Measurements of Polymer Melts Tomas Sedlacek, Peter Filip 2, Peter Saha Polymer Centre, Faculty

More information

Visualization of flow pattern over or around immersed objects in open channel flow.

Visualization of flow pattern over or around immersed objects in open channel flow. EXPERIMENT SEVEN: FLOW VISUALIZATION AND ANALYSIS I OBJECTIVE OF THE EXPERIMENT: Visualization of flow pattern over or around immersed objects in open channel flow. II THEORY AND EQUATION: Open channel:

More information

Investigation of Polymer Long Chain Branching on Film Blowing Process Stability by using Variational Principle

Investigation of Polymer Long Chain Branching on Film Blowing Process Stability by using Variational Principle Investigation of Polymer Long Chain Branching on Film Blowing Process Stability by using Variational Principle ROMAN KOLARIK a,b and MARTIN ZATLOUKAL a,b a Centre of Polymer Systems, University Institute

More information

SANDEEP TYAGI and ANUP K. GHOSH*

SANDEEP TYAGI and ANUP K. GHOSH* Linear Viscoelastic and Transient Behavior of Polypropylene and Ethylene Vinyl Acetate Blends: An Evaluation of the Linear Palierne and a Nonlinear Viscoelastic Model for Dispersive Mixtures SANDEEP TYAGI

More information

Performance evaluation of different model mixers by numerical simulation

Performance evaluation of different model mixers by numerical simulation Journal of Food Engineering 71 (2005) 295 303 www.elsevier.com/locate/jfoodeng Performance evaluation of different model mixers by numerical simulation Chenxu Yu, Sundaram Gunasekaran * Food and Bioprocess

More information

Authors: Correspondence: ABSTRACT: Keywords:

Authors: Correspondence: ABSTRACT: Keywords: Implementation of a material model with shear rate and temperature dependent viscosity Authors: Mathias Vingaard, Benny Endelt, Jesper declaville Christiansen Department of Production Aalborg University

More information

10.52 Mechanics of Fluids Spring 2006 Problem Set 3

10.52 Mechanics of Fluids Spring 2006 Problem Set 3 10.52 Mechanics of Fluids Spring 2006 Problem Set 3 Problem 1 Mass transfer studies involving the transport of a solute from a gas to a liquid often involve the use of a laminar jet of liquid. The situation

More information

Answers to questions in each section should be tied together and handed in separately.

Answers to questions in each section should be tied together and handed in separately. EGT0 ENGINEERING TRIPOS PART IA Wednesday 4 June 014 9 to 1 Paper 1 MECHANICAL ENGINEERING Answer all questions. The approximate number of marks allocated to each part of a question is indicated in the

More information

3D CFD ANALYSIS OF HEAT TRANSFER IN A SCRAPED SURFACE HEAT EXCHANGER FOR BINGHAM FLUIDS

3D CFD ANALYSIS OF HEAT TRANSFER IN A SCRAPED SURFACE HEAT EXCHANGER FOR BINGHAM FLUIDS 3D CFD ANALYSIS OF HEAT TRANSFER IN A SCRAPED SURFACE HEAT EXCHANGER FOR BINGHAM FLUIDS Ali S.* and Baccar M. *Author for correspondence Department of Mechanical Engineering, National Engineering School

More information

MODELING THREE-DIMENSIONAL NON-NEWTONIAN FLOWS IN SINGLE- SCREW EXTRUDERS

MODELING THREE-DIMENSIONAL NON-NEWTONIAN FLOWS IN SINGLE- SCREW EXTRUDERS MODELING THREE-DIMENSIONAL NON-NEWTONIAN FLOWS IN SINGLE- SCREW EXTRUDERS Christian Marschik Wolfgang Roland Bernhard Löw-Baselli Jürgen Miethlinger Johannes Kepler University Institute of Polymer Extrusion

More information

Experimental Investigation of the Development of Interfacial Instabilities in Two Layer Coextrusion Dies

Experimental Investigation of the Development of Interfacial Instabilities in Two Layer Coextrusion Dies SCREW EXTRUSION R. Valette 1 *, P. Laure 2, Y. Demay 1, J.-F. Agassant 1 1 Centre de Mise Forme des MatØriaux, Ecole Nationale des Mines de Paris, Sophia Antipolis, France 2 Institut Non-LinØaire de Nice,

More information

FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study

FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study Nityanand Pawar Mechanical Engineering, Sardar Patel College of Engineering, Mumbai, Maharashtra, India nitya.pawar@gmail.com

More information

PHYSICS 221 Fall 2013 EXAM 2: November 6, :15pm 10:15pm. Name (printed): Recitation Instructor: Section #:

PHYSICS 221 Fall 2013 EXAM 2: November 6, :15pm 10:15pm. Name (printed): Recitation Instructor: Section #: PHYSICS 221 Fall 2013 EXAM 2: November 6, 2013 8:15pm 10:15pm Name (printed): Recitation Instructor: Section #: INSTRUCTIONS: This exam contains 25 multiple choice questions, plus two extra credit questions,

More information

RHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko

RHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko RHEOLOGY Principles, Measurements, and Applications I -56081-5'79~5 1994 VCH Publishers. Inc. New York Part I. CONSTITUTIVE RELATIONS 1 1 l Elastic Solid 5 1.1 Introduction 5 1.2 The Stress Tensor 8 1.2.1

More information

An Adjustable Gap In-Line Rheometer

An Adjustable Gap In-Line Rheometer An Adjustable Gap In-Line Rheometer By D. M. Kalyon, H. Gokturk and I. Boz Highly Filled Materials Institute Hoboken, NJ 07030 Introduction The rheological behavior of polymer melts, and structured fluids

More information

Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section

Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section Article Nihon Reoroji Gakkaishi Vol.34, No.2, 105~113 (Journal of the Society of Rheology, Jaan) 2006 The Society of Rheology, Jaan Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight

More information

CPGAN # 006. The Basics of Filament Stretching Rheometry

CPGAN # 006. The Basics of Filament Stretching Rheometry Introduction Measurement of the elongational behavior of fluids is important both for basic research purposes and in industrial applications, since many complex flows contain strong extensional components,

More information

On the effects of Non-Newtonian fluids above the ribbing instability

On the effects of Non-Newtonian fluids above the ribbing instability On the effects of Non-Newtonian fluids above the ribbing instability L. Pauchard, F. Varela LÓpez*, M. Rosen*, C. Allain, P. Perrot** and M. Rabaud Laboratoire FAST, Bât. 502, Campus Universitaire, 91405

More information

150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces

150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with

More information

Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid

Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid E. Tejaswini 1*, B. Sreenivasulu 2, B. Srinivas 3 1,2,3 Gayatri Vidya Parishad College of Engineering

More information

Color Mixing in Extrusion: Simulations and Experimental Validation

Color Mixing in Extrusion: Simulations and Experimental Validation NSF GRANT # DI-0140412 NSF PROGRA NAE: aterials Processing and anufacturing Color ixing in Extrusion: Simulations and Experimental Validation Ica anas-zloczower Department of acromolecular Science Case

More information

Q1. Which of the following is the correct combination of dimensions for energy?

Q1. Which of the following is the correct combination of dimensions for energy? Tuesday, June 15, 2010 Page: 1 Q1. Which of the following is the correct combination of dimensions for energy? A) ML 2 /T 2 B) LT 2 /M C) MLT D) M 2 L 3 T E) ML/T 2 Q2. Two cars are initially 150 kilometers

More information

DISPERSIVE AND DISTRIBUTIVE MIXING CHARACTERIZATION IN EXTRUSION EQUIPMENT

DISPERSIVE AND DISTRIBUTIVE MIXING CHARACTERIZATION IN EXTRUSION EQUIPMENT DISPERSIVE AND DISTRIBUTIVE MIXING CHARACTERIZATION IN EXTRUSION EQUIPMENT Winston Wang and Ica Manas-Zloczower Department of Macromolecular Science Case Western Reserve University Abstract Mixing is a

More information

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute

More information

Stopper rod tilt and effect on nozzle flow. Rajneesh Chaudhary (Graduate Student) & Brian G. Thomas. Department of Mechanical Science and Engineering

Stopper rod tilt and effect on nozzle flow. Rajneesh Chaudhary (Graduate Student) & Brian G. Thomas. Department of Mechanical Science and Engineering Stopper rod tilt and effect on nozzle flow Rajneesh Chaudhary (Graduate Student) & Brian G. Thomas Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign, IL Project

More information

EXPERIMENTAL AND NUMERICAL STUDY ON THE FLOW VISUALIZATION IN A TRI-HELICAL STATIC MIXER

EXPERIMENTAL AND NUMERICAL STUDY ON THE FLOW VISUALIZATION IN A TRI-HELICAL STATIC MIXER 392 Journal of Marine Science and Technology, Vol. 19, No. 4, pp. 392-397 (2011) EXPERIMENTAL AND NUMERICAL STUDY ON THE FLOW VISUALIZATION IN A TRI-HELICAL STATIC MIXER Kuo-Tung Chang*, Jer-Huan Jang*,

More information

Rheometry. II.1 Introduction

Rheometry. II.1 Introduction II Rheometry II.1 Introduction Structured materials are generally composed of microstructures dispersed in a homogeneous phase [30]. These materials usually have a yield stress, i.e. a threshold stress

More information

A calibration method for a new type of rheometer

A calibration method for a new type of rheometer ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL 14, 2006 A calibration method for a new type of rheometer C Salas-Bringas 1, WK Jeksrud 1, O-I Lekang 1 and RB Schüller 2 1 Dep of Mathematical Sciences

More information

Temperature dependence of critical stress for wall slip by debonding

Temperature dependence of critical stress for wall slip by debonding J. Non-Newtonian Fluid Mech. 94 (2000) 151 157 Temperature dependence of critical stress for wall slip by debonding Yogesh M. Joshi a, Prashant S. Tapadia a, Ashish K. Lele a, R.A. Mashelkar b, a Chemical

More information

CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP

CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP 8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,

More information

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral

More information

CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE

CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,

More information

Last update: 15 Oct 2007

Last update: 15 Oct 2007 Unsteady flow of an axisymmetric annular film under gravity HOUSIADAS K, TSAMOPOULOS J -- PHYSICS OF FLUIDS --10 (10)2500-2516 1998 Yang XT, Xu ZL, Wei YM Two-dimensional simulation of hollow fiber membrane

More information

Transport by convection. Coupling convection-diffusion

Transport by convection. Coupling convection-diffusion Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.

More information

Experimental and Theoretical Investigation of Hydrodynamics Characteristics and Heat Transfer for Newtonian and Non-newtonian Fluids

Experimental and Theoretical Investigation of Hydrodynamics Characteristics and Heat Transfer for Newtonian and Non-newtonian Fluids International Journal of Energy Science and Engineering Vol. 2, No. 3, 2016, pp. 13-22 http://www.aiscience.org/journal/ijese ISSN: 2381-7267 (Print); ISSN: 2381-7275 (Online) Experimental and Theoretical

More information

Influence of steady shear flow on dynamic viscoelastic properties of un-reinforced and Kevlar, glass fibre reinforced LLDPE

Influence of steady shear flow on dynamic viscoelastic properties of un-reinforced and Kevlar, glass fibre reinforced LLDPE Bull. Mater. Sci., Vol. 27, No. 5, October 2004, pp. 409 415. Indian Academy of Sciences. Influence of steady shear flow on dynamic viscoelastic properties of un-reinforced and Kevlar, glass fibre reinforced

More information

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These

More information

AGITATION AND AERATION

AGITATION AND AERATION AGITATION AND AERATION Although in many aerobic cultures, gas sparging provides the method for both mixing and aeration - it is important that these two aspects of fermenter design be considered separately.

More information

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

Eddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):

Eddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr): AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress

More information

Analysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications

Analysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications Analysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications P. Steinert* 1, I. Schaarschmidt 1, R. Paul 1, M. Zinecker 1, M. Hackert-Oschätzchen 1, Th. Muschalek

More information

NUMERICAL ANALYSIS OF THE EFFECT OF EXTRUSION CONDITIONS ON FLOW IN SLIT DIE RHEOMETER

NUMERICAL ANALYSIS OF THE EFFECT OF EXTRUSION CONDITIONS ON FLOW IN SLIT DIE RHEOMETER Third International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia - December 3 NUMERICAL ANALYSIS OF THE EFFECT OF EXTRUSION CONDITIONS ON FLOW IN SLIT DIE RHEOMETER

More information

Rheology/Viscometry. Viscometry

Rheology/Viscometry. Viscometry Viscometry Capillary viscometry is considered as the most accurate method for the determination of the viscosity of Newtonian liquids. By this technique the time is measured a certain volume needs to flow

More information

Exercise sheet 5 (Pipe flow)

Exercise sheet 5 (Pipe flow) Exercise sheet 5 (Pipe flow) last edited June 4, 2018 These lecture notes are based on textbooks by White [13], Çengel & al.[16], and Munson & al.[18]. Except otherwise indicated, we assume that fluids

More information

DYNAMIC RHEOLOGICAL MEASUREMENTS OF AQUEOUS POLYESTER DISPERSION IN BATCH REACTOR AND TWIN SCREW EXTRUDER

DYNAMIC RHEOLOGICAL MEASUREMENTS OF AQUEOUS POLYESTER DISPERSION IN BATCH REACTOR AND TWIN SCREW EXTRUDER DYNAMIC RHEOLOGICAL MEASUREMENTS OF AQUEOUS POLYESTER DISPERSION IN BATCH REACTOR AND TWIN SCREW EXTRUDER Ali Goger *, Michael R. Thompson *, John L. Pawlak **, David J.W. Lawton *** * Department of Chemical

More information

Solvent casting flow from slot die

Solvent casting flow from slot die Korea-Australia Rheology Journal, 27(4), 325-329 (November 2015) DOI: 10.1007/s13367-015-0032-x www.springer.com/13367 Short Communication Semi Lee and Jaewook Nam* School of Chemical Engineering, Sungkyunkwan

More information

Experiments at the University of Minnesota (draft 2)

Experiments at the University of Minnesota (draft 2) Experiments at the University of Minnesota (draft 2) September 17, 2001 Studies of migration and lift and of the orientation of particles in shear flows Experiments to determine positions of spherical

More information

Fluid Mechanics Abdusselam Altunkaynak

Fluid Mechanics Abdusselam Altunkaynak Fluid Mechanics Abdusselam Altunkaynak 1. Unit systems 1.1 Introduction Natural events are independent on units. The unit to be used in a certain variable is related to the advantage that we get from it.

More information

CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR

CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM

More information

Expansions-contractions Flows

Expansions-contractions Flows III Expansions-contractions Flows III.1 Introduction Flows of viscoplastic materials through internal passages of abruptly varying cross sections are found in a wide variety of industrial and natural processes.

More information

ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics

ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics For Friday, October 26 th Start reading the handout entitled Notes on finite-volume methods. Review Chapter 7 on Dimensional Analysis

More information

Investigation of Non-Newtonian Fluids Flow Behavior in a Double Step Expansion Channel: Part 1

Investigation of Non-Newtonian Fluids Flow Behavior in a Double Step Expansion Channel: Part 1 Investigation of Non-Newtonian Fluids Flow Behavior in a Double Step Expansion Channel: Part 1 Qais A. Rishack Mechanical Engineering Department Engineering College, Basra University Basra, Iraq Khudheyer

More information

H. W. Müllner (Sp), J. Eberhardsteiner, Technische Universität Wien (A); W. Fidi, Semperit Technische Produkte Ges.m.b.H. & Co. KG, Wimpassing (A)

H. W. Müllner (Sp), J. Eberhardsteiner, Technische Universität Wien (A); W. Fidi, Semperit Technische Produkte Ges.m.b.H. & Co. KG, Wimpassing (A) Dienstag, 4. Juli 2006 Tuesday, July 4, 2006, 9.30 10.00 h Section A Rheologische Charakterisierung der Strangaufweitung von Kautschukmischungen mittels numerischer Simulationen Rheological Characterisation

More information

MECHANICAL CHARACTERISTICS OF STARCH BASED ELECTRORHEOLOGICAL FLUIDS

MECHANICAL CHARACTERISTICS OF STARCH BASED ELECTRORHEOLOGICAL FLUIDS 8 th International Machine Design and Production Conference 427 September 9-11, 1998 Ankara TURKEY ABSTRACT MECHANICAL CHARACTERISTICS OF STARCH BASED ELECTRORHEOLOGICAL FLUIDS E. R. TOPCU * and S. KAPUCU

More information

Wall-Slip of Highly Filled Powder Injection Molding Compounds: Effect of Flow Channel Geometry and Roughness

Wall-Slip of Highly Filled Powder Injection Molding Compounds: Effect of Flow Channel Geometry and Roughness Wall-Slip of Highly Filled Powder Injection Molding Compounds: Effect of Flow Channel Geometry and Roughness Berenika Hausnerovaa,b, Daniel Sanetrnika,b, Gordana Paravanovab a Dept. of Production Engineering,

More information

ON THE ROLE OF EXTENSIONAL RHEOLOGY, ELASTICITY AND DEBORAH NUMBER ON NECK-IN PHENOMENON DURING FLAT FILM PRODUCTION

ON THE ROLE OF EXTENSIONAL RHEOLOGY, ELASTICITY AND DEBORAH NUMBER ON NECK-IN PHENOMENON DURING FLAT FILM PRODUCTION ON THE ROLE OF EXTENSIONAL RHEOLOGY, ELASTICITY AND DEBORAH NUMBER ON NECK-IN PHENOMENON DURING FLAT FILM PRODUCTION Martin Zatloukal 1, Tomas Barborik 1 and Costas Tzoganakis 2 1 Polymer Centre, Faculty

More information

NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE TEMPERATURE DISTRIBUTION INSIDE OIL-COOLED TRANSFORMER WINDINGS

NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE TEMPERATURE DISTRIBUTION INSIDE OIL-COOLED TRANSFORMER WINDINGS NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE TEMPERATURE DISTRIBUTION INSIDE OIL-COOLED TRANSFORMER WINDINGS N. Schmidt 1* and S. Tenbohlen 1 and S. Chen 2 and C. Breuer 3 1 University of Stuttgart,

More information

Interaction of Non-Newtonian Fluid Dynamics and Turbulence on the Behavior of Pulp Suspension Flows

Interaction of Non-Newtonian Fluid Dynamics and Turbulence on the Behavior of Pulp Suspension Flows ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 13, 2005 Interaction of Non-Newtonian Fluid Dynamics and Turbulence on the Behavior of Pulp Suspension Flows Juha-Pekka T. Huhtanen, and Reijo J.

More information

Effect of radius ratio on pressure drop across a 90 bend for high concentration coal ash slurries

Effect of radius ratio on pressure drop across a 90 bend for high concentration coal ash slurries This paper is part of the Proceedings of the 11 International Conference th on Engineering Sciences (AFM 2016) www.witconferences.com Effect of radius ratio on pressure drop across a 90 bend for high concentration

More information

Simulation of pressure drop for combined tapered and nontapered die for polypropylene using ansys Polyflow

Simulation of pressure drop for combined tapered and nontapered die for polypropylene using ansys Polyflow IOSR Journal of Polymer and Textile Engineering (IOSR-JPTE) e-issn: 2348-019X, p-issn: 2348-0181, Volume 1, Issue 3 (May-Jun. 2014), PP 22-29 Simulation of pressure drop for combined tapered and nontapered

More information

ESTIMATION OF BINGHAM RHEOLOGICAL PARAMETERS OF SCC FROM SLUMP FLOW MEASUREMENT

ESTIMATION OF BINGHAM RHEOLOGICAL PARAMETERS OF SCC FROM SLUMP FLOW MEASUREMENT ESTIMATION OF BINGHAM RHEOLOGICAL PARAMETERS OF SCC FROM SLUMP FLOW MEASUREMENT L. N. Thrane, C. Pade and T. Svensson Danish Technological Institute, Concrete Centre, Taastrup, Denmark Abstract Different

More information

This chapter focuses on the study of the numerical approximation of threedimensional

This chapter focuses on the study of the numerical approximation of threedimensional 6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation

More information

Research Letter Processability of Pultrusion Using Natural Fiber and Thermoplastic Matrix

Research Letter Processability of Pultrusion Using Natural Fiber and Thermoplastic Matrix Research Letters in Materials Science Volume 27, Article ID 37123, 5 pages doi:1.1155/27/37123 Research Letter Processability of Pultrusion Using Natural Fiber and Thermoplastic Matrix Tham Nguyen-Chung,

More information

Multiscale simulations of the injection moulding of a diffractive optical element

Multiscale simulations of the injection moulding of a diffractive optical element Multiscale simulations of the injection moulding of a diffractive optical element Terje Tofteberg * and Erik Andreassen SINTEF Materials and Chemistry, Oslo, Norway * e-mail: terje.tofteberg@sintef.no

More information

Evaluations of Flow and Mixing Efficiency in the Kneading Disks of a Novel Tri-Screw Extruder

Evaluations of Flow and Mixing Efficiency in the Kneading Disks of a Novel Tri-Screw Extruder Journal of Applied Fluid Mechanics, Vol. 9, No. 1, pp. 51-6, 216. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Evaluations of Flow and Mixing Efficiency in the Kneading Disks

More information

Measurement of the Transverse and Longitudinal Viscosities of Continuous Fibre Reinforced Composites

Measurement of the Transverse and Longitudinal Viscosities of Continuous Fibre Reinforced Composites Measurement of the ransverse and Longitudinal Viscosities of Continuous Fibre Reinforced Composites P. Harrison,. Haylock and A.C. Long University of Nottingham - School of Mechanical, Materials & Manufacturing

More information

Fluid Mechanics. du dy

Fluid Mechanics. du dy FLUID MECHANICS Technical English - I 1 th week Fluid Mechanics FLUID STATICS FLUID DYNAMICS Fluid Statics or Hydrostatics is the study of fluids at rest. The main equation required for this is Newton's

More information