Capillary Flow of Low-Density Polyethylene
|
|
- Muriel Jefferson
- 6 years ago
- Views:
Transcription
1 Capillary Flow of Low-Density Polyethylene Mahmoud Ansari, 1 Thanasis Zisis, 2 Savvas G. Hatzikiriakos, 1 Evan Mitsoulis 2 1 Department of Chemical and Biological Engineering, The University of British Columbia, Vancouver, BC, Canada 2 School of Mining Engineering and Metallurgy, National Technical University of Athens, Zografou, Athens, Greece The capillary flow of a commercial low-density polyethylene (LDPE) melt was studied both experimentally and numerically. The excess pressure drop due to entry (Bagley correction), the compressibility, the effect of pressure on viscosity, and the possible slip effects on the capillary data analysis have been examined. Using a series of capillary dies having different diameters, D, and length-to-diameter L/D ratios, a full rheological characterization has been carried out, and the experimental data have been fitted both with a viscous model (Carreau-Yasuda) and a viscoelastic one (the Kaye Bernstein, Kearsley, Zapas/Papanastasiou, Scriven, Macosko, or K-BKZ/PSM model). Particular emphasis has been given on the pressure-dependence of viscosity, with a pressure-dependent coefficient b p. For the viscous model, the viscosity is a function of both temperature and pressure. For the viscoelastic K- BKZ model, the time-temperature shifting concept has been used for the non-isothermal calculations, while the time pressure shifting concept has been used to shift the relaxation moduli for the pressure-dependence effect. It was found that only the viscoelastic simulations were capable of reproducing the experimental data well, while any viscous modeling always underestimates the pressures, especially at the higher apparent shear rates and L/D ratios. POLYM. ENG. SCI., 52: , ª 2012 Society of Plastics Engineers Correspondence to: Evan Mitsoulis; mitsouli@metal.ntua.gr Contract grant sponsor: Natural Sciences and Engineering Research Council (NSERC) of Canada; contract grant sponsor: NTUA; contract grant number: PEBE DOI /pen Published online in Wiley Online Library (wileyonlinelibrary.com). VC 2012 Society of Plastics Engineers INTRODUCTION Capillary rheometry is extensively used in both industry and academia to assess the rheological behavior of polymer melts at high shear rates before testing their processability in full industrial scale [1]. When such a flow is used and the raw data are collected, a number of important corrections should be applied before the rheological data can be compared with corresponding data from a rotational rheometer [2, 3]. First, capillary flow involves flow through a contraction of a certain angle, where there is a large pressure drop associated with such flow, known as entrance pressure or Bagley correction [1, 4]. This pressure is required in order to calculate the true shear stress, and also, frequently, the apparent extensional rheology of molten polymers, a method well practiced in industry [5 9]. Therefore, it is important to understand the origin of this excess pressure and consequently to be able to predict it. Many studies have previously attempted to examine the origin of entrance pressure and its prediction for low-density polyethylenes (LDPEs). Feigl and Öttinger [10] have simulated the axisymmetric contraction flow of an LDPE melt using a multimode Rivlin Sawyers integral constitutive model, which is a special case of the Kaye Bernstein, Kearsley, Zapas (K-BKZ) model. However, their numerical results well under-predicted the available experimental pressure data for entry flows. Barakos and Mitsoulis [11, 12] reported similar findings with the same K-BKZ model for the Bagley correction in capillary flow of an LDPE melt tested by a working group of the International Union of Pure and Applied Chemistry (IUPAC-LDPE). Béraudo et al. [13] used a multimode Phan-Thien/Tanner (PTT) constitutive relation and found that numerical predictions significantly under-estimated the experimental findings. Guillet et al. [14] also studied the entrance pressure losses for an LDPE melt both experimentally and numerically using the K-BKZ model. Again significant under-estimation was reported. Similar results were obtained by Hatzikiriakos and Mitsoulis [15], Mitsoulis et al. [16], and Mitsoulis and Hatzikiriakos [17] for different polymer melts with the K-BKZ model. All of these works fitted rheological data for the melts and performed isothermal, viscoelastic simulations with no pressure-dependence of the viscosity. Thus, it had become apparent that pressure underestimation was a problem, especially for the more elastic LDPE melts, which needed some extra attention and more careful work to be resolved. POLYMER ENGINEERING AND SCIENCE -2012
2 Along these lines, a recent study by the principal authors [18] examined three different polyethylenes, namely LDPE, a linear low-density polyethylene (LLDPE), and a high-density polyethylene (HDPE) in flow through tapered dies. The novelty was the pressuredependence of the viscosity, which was assumed a constant. In isothermal simulations, the problem of predicting the Bagley correction for LDPE was solved satisfactorily at extremely high shear rates (up to 1000 s 1 ) for the first time [18] by taking into account the effect of pressure on the viscosity on top of viscoelasticity. Thus, it became apparent that both viscoelasticity and pressure-dependence of the viscosity were significant in obtaining high entry pressures in agreement with experiments. It should be noted that such high entry pressures are important for other materials as well, such as pastes, as have been reported by many in the literature [19 22]. Since the effect of pressure on the viscosity was found to be so important, it was decided to re-examine more carefully experimental data on LDPE. It is well known that using capillary data from dies of various length-to-diameter (L/D) ratios, the dependence of viscosity on pressure can be assessed [2, 23 28]. As a first approximation, the following expression (Barus equation) can be used to determine the parameter, b p, known as the pressure coefficient of viscosity: Z ¼ Z p0 exp b p p ; (1) where Z is the viscosity at pressure p, and Z p0 is the viscosity at ambient pressure. Using a pressurized slidingplate rheometer, Dealy and coworkers [26, 29, 30] have determined the coefficient b p for various systems, including LDPE, poly-a-methystyrene-co-acrylonitrile (PaMSAN), and LLDPE. For the LDPE (of main interest in the present work), a value between Pa and Pa has been reported by various authors [24, 28, 31, 32]. This coefficient has also been reported to be a function of the shear rate, with b p of LDPE decreasing significantly with shear rate [24, 33]. It is, therefore, a partial objective of the present work to carefully re-examine the pressure-dependence of the viscosity for the LDPE melt at hand. Another issue is the effect of slip of the material along the die walls. Unlike Newtonian fluids, polymer melts slip over solid surfaces when the wall shear stress exceeds a critical value [34 38]. In particular, slip effects have been reported in the capillary flow of molten polyethylenes [18, 34, 35, 37]. Thus, in a comprehensive study of any melt, such as the one undertaken here, possible slip effects should be studied. This can be done by using at least three capillary dies having different diameters and the same L/D in order to keep the effect of pressure constant. If the flow curve shows a diameter dependency similar to that implied by the explicit relations derived by Mooney [39], slip is present and the Mooney method can be used to determine the slip velocity as a function of the wall shear stress. It is, therefore, another partial objective of the present work to examine whether the LDPE melt at hand slips at the wall. Another issue is the effect of the non-isothermality of flow. Any analysis for the determination of slip and the dependence of viscosity on pressure are complicated by viscous dissipation (viscous heating) effects. Viscous heating in the capillary extrusion of polymer melts has been the subject of several reviews and studies over the past decades [2, 3, 40 46]. Rosenbaum and Hatzikiriakos [47] have considered the analysis of the combined effects of viscous heating and wall slip using numerical analysis. They have reported that in cases where viscous heating and slip effects are significant, the traditional methods used to estimate the slip velocity (for example Mooney s technique) often fail to give accurate or even physically meaningful results. They have proposed, a data analysis procedure based on a mathematical model for the non-isothermal capillary flow of polymer melts coupled with heat transfer to analyze the capillary data instead. Pressure-dependence of viscosity and compressibility effects have been neglected in the analysis by Rosenbaum and Hatzikiriakos [47]. More recently, Laun [2, 3] considered the superimposed effect of dissipative heating and a pressuredependent viscosity in circular die flows at high shear rates. The problem was treated theoretically, and pertinent relations for capillary rheometry were derived. The applicability of the data treatment was demonstrated for the die flow of LDPE. However, compressibility and slip effects were not considered. It is, therefore, the main objective of this work to study in a comprehensive way the capillary flow of an LDPE melt both experimentally and numerically by considering all possible capillary effects, combined and separately. Namely, the effects of viscoelasticity, compressibility, pressure-dependence of viscosity, and viscous heating (non-isothermal flow) will be considered in order to assess the significance of each effect. It will be shown that viscoelasticity and the pressure-dependence of viscosity are primarily the reasons for high values of the entry pressure losses, especially at elevated shear rates, for LDPE melts. EXPERIMENTAL Materials A typical LDPE melt was used in this work in order to address the effects of pressure, temperature, and compressibility on its capillary flow. The viscosity of LDPE melts is known to be more sensitive to temperature and pressure compared to that of other polyethylenes, such as HDPE and/or LLDPE, and it is the main reason for being selected in this study. It is also known that branched polymers exhibit extensional strain-hardening effects, and therefore increased entrance pressure values are obtained 650 POLYMER ENGINEERING AND SCIENCE DOI /pen
3 TABLE 1. Properties of the LDPE resin used in this study. Sample ID Resin producer Melt index (g/10 min) (1908C) Density (g/cm 3 ) (258C) Zero-shear-rate viscosity (Pa s) 1608C LDPE DOW Chemicals ,230 compared to those obtained for their linear counterparts [15 17]. This is expected as the entry pressure has been theoretically associated with extensional viscosity [5 9]. The properties of the particular LDPE used in this study are summarized in Table 1. Rheological Testing Several rheological experiments were carried out in order to characterize the polymer in the molten state. First, an Anton Paar MCR-501 device was used in the parallel-plate geometry to determine its linear viscoelastic moduli over a wide range of temperatures from 130 to 2308C. Master curves were obtained using the time temperature superposition (TTS) principle, and the results are presented at the reference temperature of 1608C. These values are reported in Table 1 and are found to compare well with those obtained from linear viscoelastic measurements. The polymer was also rheologically characterized in simple extension using the SER-2 Universal Testing Platform from Xpansion Instruments [48, 49]. Details of sample preparation and operating principles of this rheometer to obtain reliable results are given in Delgadillo et al. [50, 51]. Uniaxial extension tests were performed at the reference temperature of 1608C. Capillary Flow Testing An Instron capillary rheometer (constant piston speed) was used to determine the entrance pressure (known also as Bagley method) [1] and the viscosity as a function of the wall shear stress, r W, and apparent shear rate, _c A ¼ 32Q=pD 3 at 1608C, where Q is the volumetric flow rate and D is the capillary diameter. A series of orifice dies (D ¼ cm and L/D ¼ 0) were also used to determine directly the entrance pressure (Bagley correction) as a function of the apparent shear rate for various contraction angles (15, 30, 45, 60, and 908). In addition, several dies having the same diameter (D ¼ cm) and various L/D ratios (L/D ¼ 2, 5, 16, 33) have been used to determine the effect of pressure on viscosity. These data were also used to determine the Bagley correction and compared with the one determined from the orifice die. The agreement found was satisfactory. In addition, capillary dies having the same L/D ratio of 20 and different diameters (D ¼ , , and cm) were used to examine slip effects for the case of LDPE by means of the Mooney technique [39]. GOVERNING EQUATIONS AND RHEOLOGICAL MODELING We consider the conservation equations of mass, momentum, and energy for weakly compressible fluids under non-isothermal, creeping, steady flow conditions. These are written as [52, 53]: u rrþrðru Þ ¼ 0; (2) 0 ¼ rp þrt; (3) rc p u rt ¼ kr 2 T þ t : ru; (4) where r is the density, u is the velocity vector, p is the pressure, s is the extra stress tensor, T is the temperature, C p is the heat capacity, and k is the thermal conductivity. For a compressible fluid, pressure and density are connected as a first approximation through a simple linear thermodynamic equation of state [53]: r ¼ r 0 ð1 þ b c pþ; (5) where b c is the isothermal compressibility with the density to be r 0 at a reference pressure p 0 (¼ 0). The viscous stresses are given for inelastic non-newtonian compressible fluids by the relation [52, 53]: t ¼ Zðjj ġ Þ ġ 2 3 ðru Þ I ; (6) where g ðjj _c Þ is the apparent non-newtonian viscosity, which is a function of the magnitude jjof _c the rate-ofstrain tensor _c ¼ru þru T, which is given by: jj¼ ġ rffiffiffiffiffiffiffi 1 2 II g ¼ 1 1=2 2 ġ : ġ ; (7) where II _c is the second invariant of _c II g ¼ ġ : ġ ¼ X X ġ ij ġ ij ; (8) i j The tensor I in Eq. 6 is the unit tensor. To evaluate the role of viscoelasticity in the prediction of Bagley correction, it is instructive to consider first purely viscous models in the simulations. Namely, the Carreau-Yasuda model was used to fit the shear viscosity data of the LDPE melt. The Carreau-Yasuda model is written as [1]: DOI /pen POLYMER ENGINEERING AND SCIENCE
4 listed in Table 2. We observe that the LDPE melt is very shear-thinning for shear rates above 1 s 1 due to the presence of long-chain branching [54]. The Carreau-Yasuda model fits the data well over the range of experiment results. For easy checks and simple analytical formulas, the high shear-rate range of the viscosity data can be fitted to the power-law model for the viscosity [1]: Z ¼ Kjġj n 1 ; (10) FIG. 1. The shear viscosity of the LDPE melt at 1608C fitted with the Carreau-Yasuda (Eq. 9) and the power-law (Eq. 10) models using the parameters listed in Table 2. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] Z ¼ Z 0;C ½1 þðljġjþ a Š n C 1 a ; (9) where Z 0,C is the Carreau zero-shear-rate viscosity, l is a time constant, n C is the Carreau power-law index, and a is the Yasuda exponent (¼ 2 for the simple Carreau model). The fitted viscosity of the LDPE melt by Eq. 9 is plotted in Fig. 1, while the parameters of the model are where K is the consistency index and n is the power-law index. These values are also given in Table 2. For the capillary flow simulations, the effect of pressure on viscosity should be taken into account as this becomes evident below. This effect is quite significant for LDPE compared to LLDPE and less significant for HDPE [31, 32]. This effect can be taken into account by multiplying the constitutive relation with a pressureshift factor, a p, defined by the Barus equation, that is [27, 32]: a p Z ¼ expðb Z P pþ (11) p0 where b p is the pressure coefficient and p is the absolute pressure, as discussed above. Viscoelasticity is included in the present work via an appropriate rheological model for the stresses. This is a K-BKZ equation proposed by Papanastasiou et al. [55] and modified by Luo and Tanner [56]. This is written as: t ¼ 1 Z t X N a k t t0 a exp C 1 t ðt 0 ÞþyC t ðt 0 Þ dt 0 ; (12) 1 y 1 l k¼1 k l k ða 3ÞþbI C 1þð1 bþi C where t is the current time, k k and a k are the relaxation times and relaxation modulus coefficients, N is the number of relaxation modes, a and b are material constants, and I C, I C 1 are the first invariants of the Cauchy-Green tensor C t and its inverse C 1 t, the Finger strain tensor. The material constant y is given by LDPE melt to rheologically characterize it. Figure 2 plots the master dynamic moduli G 0 and G 00 for LDPE at the reference temperature of 1608C. The model predictions obtained by fitting the experimental data to Eq. 12 with a spectrum of relaxation times, k k, and coefficients, a k, determined by a non-linear regression package [57], are also plotted. The parameters found from the fitting proce- N 2 ¼ y N 1 1 y (13) TABLE 2. Parameters for the LDPE resin obeying the Carreau-Yasuda model (Eq. 9) and the power-law model (Eq. 10) at 1608C. where N 1 and N 2 are the first and second normal stress differences, respectively. It is noted that y is not zero for polymer melts, which possess a non-zero second normal stress difference. Its usual range is between 0.1 and 0.2 in accordance with experimental findings [1]. As discussed above, experiments were performed in the parallel-plate and extensional rheometers for the Parameter Value g 0,C 617,230 Pa s k s n C a K 13,780 Pa s n n POLYMER ENGINEERING AND SCIENCE DOI /pen
5 FIG. 2. Experimental data (symbols) and model predictions of storage (G 0 ) and loss (G 00 ) moduli for the LDPE melt at 1608C using the relaxation times listed in Table 3. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] dure are listed in Table 3. The relaxation spectrum is used to find the average relaxation time, k, and zero-shear-rate viscosity, g 0, according to the formulas: PN a k l 2 k l ¼ k¼1 P N k¼1 Z 0 ¼ XN k¼1 a k l k (14) a k l k ; (15) FIG. 3. Experimental data (solid symbols) and model predictions of shear viscosity, g S, first normal stress difference, N 1, and elongational viscosity, g E, for the LDPE melt at 1608C using the K-BKZ model (Eq. 12) with the parameters listed in Table 3. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] The values of these parameters are k ¼ 130 s, g 0 ¼ 150,607 Pa s, indicating a highly elastic melt with a long average relaxation time. Figure 3 plots a number of calculated and experimental material functions for the LDPE melt at the reference temperature of 1608C. Namely, data for the shear viscosity, g S, the elongational viscosity, g E, and the first normal stress difference, N 1, are plotted as functions of corresponding rates (shear or extensional). The parameter b that controls the calculated elongational viscosity was fitted by using the extensional behavior of the melt. Figure 4 shows the extensional behavior of the LDPE at several Hencky strain rates at 1608C and the model predictions of Eq. 12 using multiple b- values listed in Table 3. It can be seen that the overall rheological representation of all material functions is excellent. Non-Isothermal Modeling For the non-isothermal calculations, it is necessary to derive a non-isothermal constitutive equation from the isothermal one. This is done by applying the time-temperature shifting concept as explained by Luo and Tanner [58]. This concept is based on the relative difference between the observer s time scale and the material s internal time scale. Some intuition can be obtained considering that as temperature raises so does the amount of molecular motion occur- TABLE 3. Relaxation spectra and material constants for the LDPE resin obeying the K-BKZ model (Eq. 12) at 1608C. a ¼ 7.336, y ¼ 0 k k k (s) a k (Pa) b k , , , FIG. 4. Extensional data for LDPE melt at 1608C and their best fit using the K-BKZ model (Eq. 12) with the parameters listed in Table 3. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] DOI /pen POLYMER ENGINEERING AND SCIENCE
6 ring in one unit of observer s time. Using n for the time measured by the material s own internal clock, the following relation holds between n and the observer s time t: dx ¼ dt a T ðtþ ; (16) where the denominator is the time-shifting factor [52]. The above relationship expresses the Morland Lee hypothesis [52, 58] in a differential form and can be used to obtain the following integral relation between the particle s elapsed time and the observing period: x ¼ Z t 0 a 1 T ½Tðt0 ÞŠdt 0 : (17) Thus as a fluid particle is tracked along a streamline segment Dl i, the particle s time corresponding to the residence time Dt 0 is given by Dx ¼ Dl i V i a T ðt i Þ ; (18) where V i is the particle velocity. The above Eq. 16 has been used by Luo and Tanner [58] to obtain the non-isothermal form of the K-BKZ equation (Eq. 12). This was done by replacing the observer s time t in Eq. 12 with the particle s time n given by Eq. 17. The resulting version of the constitutive equation becomes t ¼ 1 1 y Z xðtþ X N 1 k¼1 a k x x0 exp l k l k a h ða 3Þþb k I C 1þð1 b k ÞI C i tðxþ ðt0 ðx 0 ÞÞþyC tðxþ ðt 0 ðx 0 ÞÞ dx 0 : (19) C 1 For LDPE melts, the following temperature-shifting function has been found to be adequate [52, 59]: a T ðtþ ¼ Z ¼ exp E 1 Z 0 R g T 1 T 0 : (20) In the above, g 0 is a reference viscosity at T 0, E is the activation energy constant, R g is the ideal gas constant, and T 0 is a reference temperature (in K). The activation energy constant can be determined from zero-shear-rate viscosity data [59]. Another method can be applied to derive the non-isothermal constitutive equation from the isothermal one. This method is again based on the time temperature superposition temperature and simply consists of shifting the relaxation times k k from the temperature history within the material s time scale [60]. The equation used to shift the relaxation times in the material s history is given by [60, 61]: l k ðt 0 ðt 0 ÞÞ ¼ l k ðt 0 Þa T ðt 0 ðt 0 ÞÞ: (21) Both methods have been applied before for the process of film blowing [61] with virtually the same results. Test runs here confirmed this finding and the current results were derived with the n-approach. The viscoelastic stresses calculated by the above constitutive equation (Eq. 19) enter in the energy equation (Eq. 4) as a contribution to the viscous dissipation term. This assumes that all the elastic stored energy in the material is converted into heat and serves to raise the temperature of the melt. This treatment has been the standard practice in non-isothermal simulations so far (see, e.g., Barakos and Mitsoulis [62]) and it is followed in the present work as well. Peters and Baaijens [63] have devised a different approach to handle the splitting between the viscous contribution to energy dissipation and elastic stored energy. However, such a development is beyond the scope of this study. The thermal properties of the melt have been gathered from various sources. The thermal conductivity was given by Laun [3] to be 0.17 W/m K, and the heat capacity was equal to J/(kg K). The density is given by the equation [44]: rðtþ ¼0:8772 0:00059T; for T > 150 C; (22) where r is in g/cm 3 and T is in 8C. Its value at T 0 ¼ 1608C isr ¼ g/cm 3. The activation energy for this particular melt was found by applying the time temperature superposition to be equal to E ¼ 64,100 J/mol [50, 51]. Others have reported similar values [59]. All these properties and their provenance are gathered together in Table 4. In this table, we have also added the value for the coefficient of compressibility, b c [64] and the value of a constant pressurecoefficient of viscosity, b p [27]. The various thermal and flow parameters are combined to give appropriate dimensionless numbers [40, 44]. The relevant ones here are the Peclet number, Pe, and the Nahme-Griffith number, Na. These are defined as: Pe ¼ rc pur ; (23) k 654 POLYMER ENGINEERING AND SCIENCE DOI /pen
7 TABLE C. Values of the various parameters for the LDPE resin at Parameter Value Reference b c MPa 1 [64] b p MPa 1 [27] m 0.65 GPa np 1 (0.11 MPa np 1 ) This work n p 0.75 This work q g/cm 3 [44] C p 2.25 J/(g K) [3] k J/(s cm K) [3] E 64,100 J/mol [50,51] R g J/(mol K) [1] T C (433 K) [50,51] Na ¼ ZEU2 kr g T0 2 ; (24) where g ¼ f ðu=rþ is a nominal viscosity given by the Carreau-Yasuda model (Eq. 9) at a nominal shear rate of U/R, R g is the ideal gas constant, and U (¼ _c A R/4) is the average velocity in the capillary die. The Pe number represents the ratio of heat convection to conduction, and the Na number represents the ratio of viscous dissipation to conduction and indicates the extent of coupling between the momentum and energy equations. A thorough discussion of these effects in non-isothermal polymer melt flow is given by Winter [40]. With the above properties and a die radius R ¼ cm, the dimensionless thermal numbers are in the ranges: 3.3 \ Pe \ 658, \ Na \ 3, showing a relatively strong convection (Pe1), and a moderate coupling between momentum and energy equations (Na1). A value of Na [ 1 indicates temperature non-uniformities generated by viscous dissipation, and a strong coupling between momentum and energy equations. For a die radius R ¼ 0.04 cm, the corresponding ranges are: 5.2 \ Pe \ 1036, \ Na \ More details are given in Table 5. The fully developed values for temperature occur for L/D ¼! and give a maximum temperature rise DT max at the centreline of the tube according to the formula for a power-law fluid [65]: DT max ¼ K k nr 2 n þ 1 nþ1 Vmax : (25) 3n þ 1 n R According to Eq. 25 and for R ¼ cm, we get a DT max ¼ 29.28C. For R ¼ 0.04 cm, we get a DT max ¼ 38.38C. Thus, a significant temperature rise is expected for very long dies, but for shorter dies as the ones used here, the non-isothermal simulations are necessary to find out the extent to which temperatures rise. Pressure-Dependent Modeling Similarly with the time temperature superposition principle where the stresses are calculated at a different temperature using the shift factor a T, the time pressure superposition principle can be used to account for the pressure effect on the stresses. In both cases of viscous or viscoelastic models, the new stresses are calculated using the pressure-shift factor a p. For viscous models, Eq. 11 is used to modify the viscosity. For viscoelastic models, such as the K-BKZ model (Eq. 12), the pressure-shift factor modifies the relaxation moduli, a k, according to: a k ðpt ð 0 ÞÞ ¼ a k ðp 0 Þa p ðpt ð 0 ÞÞ: (26) This is equivalent to multiplying the stresses by a p, according to Eq. 12. It should be noted that a p is an exponential function of b p, which itself may depend on pressure p. Since in a flow field negative pressures may appear, especially around singularities as is the case in contraction flows, special care must be taken numerically to handle these functions for negative numbers. Failure to do so leads to nonsensical results and/or to divergence. It should be noted that in our recent work [18], the pressure dependence of the viscosity was not handled in this way, but simply the overall pressure results were multiplied by the a p -factor. Thus, this intrinsic handling of the pressure, which assumes different local values, is a major contribution of the present work. METHOD OF SOLUTION The solution to the above conservation and constitutive equations is carried out with two codes, one for viscous flows (u-v-p-t-h formulation) [66] and one for viscoelastic flows [62, 67]. The boundary conditions (BC) for the problem at hand are well known and can be found in our earlier publication [62]. Briefly, we assume no-slip and a constant temperature T 0 at the solid walls; at entry, a fully-developed velocity profile is imposed, corresponding to the flow rate at hand, and a constant temperature T 0 is assumed; at the outlet, zero surface traction and zero heat flux are assumed; at the centerline, symmetry is assumed. In one case, adiabatic BC are assumed on the solid walls to study the differences with the constant-wall temperature BC. The novelty in the present work is the addition of the Barus equation (Eq. 11) for the pressure-shift factor a p with either a constant b p or an experimentally found function b p ¼ b p (p), and the calculation of the viscoelastic TABLE 5. Range of the dimensionless parameters in the flow of LDPE resin at 1608C (die radius R cm). Apparent Shear Rate, _c A ðs 1 Þ Peclet Number, Pe Nahme Number, Na DOI /pen POLYMER ENGINEERING AND SCIENCE
8 stresses based on Eq. 26. Furthermore, the non-isothermal viscous and viscoelastic calculations were checked again against analytical solutions for the power-law fluid with viscous heating having either isothermal or adiabatic walls, and against previous numerical solutions for extrudate swell [62]. It is surprising that since the early work by Luo and Tanner [58], and Barakos and Mitsoulis [62], no other non-isothermal viscoelastic works seem to have appeared in the literature for integral constitutive equations, except for simplified 1-D models for some polymer processes dominated by extension [60, 61]. This despite the fact that polymer processing is almost always occurring under non-isothermal conditions. The viscous simulations are extremely fast and are used as a first step to study the whole range of parameter values and die designs. The viscoelastic simulations admittedly are much harder to do and they need good initial flow fields to get solutions at such elevated apparent shear rates as 1000 s 1. In our recent work [18], we explained how it was possible for the first time to do viscoelastic computations up to the highest apparent shear rates with good results. However, those calculations were carried out under isothermal conditions and with no implicit internal calculations for pressure-dependence. The final result for the pressures was multiplied by a p with a constant value of b p. In the current simulations, these limitations are lifted, and full viscoelastic non-isothermal simulations are carried out with a pressure dependence of the relaxation moduli for the first time. The iteration method is direct substitution (Picard scheme). The solution strategy starts at a given apparent shear rate from the viscous non-isothermal solution with pressure dependence, and then using this as an initial solution it continues for the non-isothermal viscoelastic solution without pressure dependence. Once a new solution is achieved, more iterations are performed for the pressure-dependence of the relaxation moduli. In our earlier efforts, only a constant b p value could be handled, but in the final efforts an experimentally found b p -function gave well behaved results for all cases, provided care was taken to handle numerically correctly the negative pressures in the field. The results for L/D ¼ 0 were obtained as in Ref. 18 by subtracting from the longest dies results obtained for the same dies without the reservoir. Thus, it was found again necessary to use a long die to fully unravel the high viscoelastic stresses of the LDPE melt. FIG. 5. Experimental data for pressure for LDPE melt at 1608C asa function of L/D for different apparent shear rates. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] EXPERIMENTAL RESULTS Entrance (End) Pressure Figure 5 plots the entrance pressure (or end pressure due to L/D ¼ 0) of LDPE at 1608C as a function of the die length L/D for an extended range of values of the apparent shear rate from 5 s 1 to 1000 s 1. This plot is also known as the Bagley plot. The data from an orifice die are also plotted. The data show a curvature with concavity pointing upwards consistent with Eq. 1 (effect of pressure on viscosity). These results are for a contraction angle 2a ¼ 908. Our recent work [18] has shown that the entry pressure, although significant, is nearly independent of contraction angle for angles 2a 458. The Flow Curves Figure 6 plots the flow curves of LDPE at 1608C obtained with capillary dies having lengths of L/D ¼ 2, 5, 16, and 33. The Bagley and Rabinowitsch corrections have been applied to the data. Also plotted are the linear FIG. 6. The flow curves of LDPE at 1608C determined in capillary rheometry with dies having L/D ratios from 2 to 33. The data at the small L/D s agree well with the LVE data obtained from a rotational rheometer at ambient pressure. The capillary data have been corrected by applying both the Bagley and Rabinowitsch corrections. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] 656 POLYMER ENGINEERING AND SCIENCE DOI /pen
9 FIG. 7. The pressure coefficient, b p, as a function of pressure. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] viscoelastic data at 1608C, which agree well with the data at L/D ¼ 2 and 5. The data clearly shows that there is an effect of pressure on viscosity. The flow curves have been corrected for these effects according to Barus Equation (Eq. 1) and the procedure discussed by others [23, 27, 28]. Based on the superposition procedure, the data for the higher L/D s are shifted horizontally to superpose with the LVE data. In this way the shifting parameter b super p (which is a shear-rate-dependent pressure coefficient) can be obtained. According to Couch and Binding [23] this parameter could be related to Barus pressure coefficient b p as: where n is the power-law index. b p ¼ nb super p ; (27) FIG. 9. Die diameter effect to check the possibility of slip at the wall. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] The calculated coefficients b p depend on the shear rate and the L/D (level of pressure). Duvdevani and Klein [33] and Liang [24] have also reported significant shear-rate dependence of the coefficient b p. Although the pressure dependence of b p has been reported in some studies [68, 69], it has been accepted for practical reasons that this parameter is independent of pressure [70 73]. In our case, b p shows a power-law dependency on pressure (see Fig. 7), according to the following equation: b p ¼ mp n p (28) where m ¼ 0.62 GPa n p 1 and n p ¼ 0.75, with the pressure p given in GPa. The pressure-corrected data for the LDPE are plotted in Fig. 8. The data superposes well and the capillary data agrees well with the LVE data. The values of b p plotted in Fig. 8 agree well with the reported values in the literature [23, 27, 28, 31]. The diameter dependence of the flow curve has also been studied in order to assess the slip behavior of LDPE. The results are shown in Fig. 9. No diameter dependence is evident as the flow curves for three significantly different diameters practically superpose. However, if viscous heating effects are significant, this may mask the diameter dependence as discussed by Rosenbaum and Hatzikiriakos [47]. FIG. 8. The pressure corrected flow curves of LDPE at 1608C determined in capillary rheometry with dies having L/D ratios from 2 to 33. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] NUMERICAL RESULTS Viscous Modeling It is instructive to perform first calculations with a purely viscous model, so that the effect of viscoelasticity will become evident later. The numerical simulations have been undertaken using the viscous Carreau-Yasuda model (Eq. 9). This constitutive relation is solved together with DOI /pen POLYMER ENGINEERING AND SCIENCE
10 FIG. 10. (a) A typical finite element grid for the simulations in an 18.75:1 tapered circular contraction with L/R ¼ 40 and 2a ¼ 908. The upper grid (M1) consists of 1700 elements and 7161 nodes, while the lower grid is created by subdivision of each M1 element into four subelements to form a denser grid for checking the results for grid-independence; (b) detailed grids near the die entry. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] the conservation equations of mass and momentum either for an incompressible or compressible fluid under isothermal or non-isothermal conditions (conservation of energy equation) with or without the effect of pressure-dependence of the viscosity. FIG. 11. The effect of compressibility on the pressure distribution for the LDPE melt at 1608C at various values of apparent shear rate. The effect is negligible. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] FIG. 12. The effect of pressure-dependent viscosity on the pressure distribution for the LDPE melt at 1608C at various values of apparent shear rate. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] For the finite element mesh arrangement, we have used our long experience with viscous and viscoelastic flows and chosen a grid that progressively adds more elements as one moves towards the singularity at the entrance to the die, while the elements become bigger as one moves away from the singularity. A typical finite element grid is shown in Fig. 10 for L/R ¼ 40. The domain represents an 18.75:1 tapered circular contraction with an entrance angle 2a ¼ 908. The grid consists of 1700 elements, 7161 nodes, and 22,261 unknown degrees of freedom (d.o.f.), while a four-times denser grid is also used, having been created by subdivision of each element into four sub-elements for checking purposes of grid-independent results. This checking consists of reporting the overall pressures in the system from the two meshes and making sure that the differences are less than 1% between the two results. Having fixed the Carreau-Yasuda model parameters and the problem geometry, the only parameter left to vary was the apparent shear rate in the die (_c A ¼ 4Q=pR 3 ). Simulations were performed for the whole range of experimental apparent shear rates, namely from 5 s 1 to as high as 1000 s 1. First, runs were carried out to study each effect separately, namely the effect of compressibility alone, the effect of a pressure-dependent viscosity alone, and the effect of a temperature-dependent viscosity alone, and these are compared with the base case of no effects at all (b c ¼ b p ¼ a T ¼ 0). These are shown in Figs , respectively. We note in Fig. 11 that compressibility plays a negligible role, as the results are virtually the same with or without compressibility. Admittedly, this is expected because the value for LDPE of b c ¼ MPa 1 is small [64]. On the other hand, a pressure-dependence of 658 POLYMER ENGINEERING AND SCIENCE DOI /pen
11 FIG. 13. The effect of temperature-dependent viscosity on the pressure distribution for the LDPE melt at 1608C at various values of apparent shear rate. The effect is small. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] FIG. 14. The effect of compressibility, pressure-, and temperature-dependent viscosity on the pressure distribution for the LDPE melt at 1608C at various values of apparent shear rate. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] FIG. 15. The effect of compressibility, pressure-, and temperature-dependent viscosity on the pressure distribution for the LDPE melt at 1608C at various values of apparent shear rate but for a longer die (L/D ¼ 40). [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] the viscosity, with a constant value of b p ¼ MPa 1 [27], plays a significant role (see Fig. 12), augmenting the pressure drop substantially, and this effect is stronger as the apparent shear rate increases. Finally, a temperaturedependent viscosity serves to lower somewhat the pressure in the system (see Fig. 13), but again this effect is negligible for small to intermediate apparent shear rates, and becomes more evident at the high end of shear rate values. The non-isothermal runs have been made with isothermal walls at T 0 ¼ 1608C. At the highest apparent shear rate of 1000 s 1 we have also assumed adiabatic walls, and this brings down the pressure appreciably (from 19 MPa to 15 MPa). In reality, the walls are neither isothermal nor adiabatic but somewhere in between; however they are much closer to the isothermal case [40, 44]. We note that the effect of a temperature-dependent viscosity is opposite to the effect of a pressure-dependent viscosity, and these two effects are in competition. The overall effect of the three contributions is shown in Fig. 14, where it can be seen that the pressure is higher than the base case for all apparent shear rates, driven mainly by the pressure-dependent viscosity effect. When the die is longer (L/D ¼ 40), the overall effects are more pronounced as evidenced in Fig. 15. Again, the main effect is the pressure-dependence of viscosity, which given a longer die has more length to show its influence. The pressure drop in the die is not linear but concave, although this is more evident at the highest apparent shear rate of 1000 s 1. The effect of temperature is not very strong, as evidenced in Fig. 16. For L/D ¼ 20, the centreline temperature distribution is shown for different apparent shear rates for isothermal walls and for the two cases of b p ¼ 0 and b p ¼ MPa 1. The maximum temperature rise at the centreline does not exceed 38C, although near the wall is 1 2 degrees higher due to a higher viscous dissipation contribution there. Collecting all the pressure results together for the five dies with L/D ¼ 0.2, 2, 5, 16, 33, and 5 apparent shear rates _c A gives in Fig. 17 the well known Bagley plot. The experimental results are shown as symbols while the numerical DOI /pen POLYMER ENGINEERING AND SCIENCE
12 FIG. 16. The effect of compressibility, pressure-, and temperature-dependent viscosity on the centreline temperature distribution for the LDPE melt at 1608C at various values of apparent shear rate (L/D ¼ 20). [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] results here and in the subsequent graphs are shown as lines. These simulation results have been obtained with a variable b p given by Eq. 28. The results are below the experimental data, and this is mainly due to the inability of a viscous model to capture the end correction (for L/D ¼ 0). A single value of b p ¼ MPa 1 gives somewhat higher values, especially for the highest L/D and _c A values, but it still cannot account for the L/D ¼ 0 end correction. As it was found out in Ref. 18, a purely viscous model for LDPE FIG. 18. Viscous vs. viscoelastic results for the axial pressure distribution with all effects present for the LDPE melt at 1608C at various values of apparent shear rate (L/D ¼ 33). [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] gives end corrections one order of magnitude lower than the experimental ones. It should be noted that underestimates, such as these, have been found in all previous viscous works, even though the effect of pressure was not taken into account in all of them. So the effect of pressure on the viscosity alone is not enough to correctly predict the pressure drop. It is then at this point that we turn our attention to the viscoelastic results. FIG. 17. The effect of L/D on the pressure for the LDPE melt at 1608C at various values of apparent shear rate (all effects accounted for). Simulations with the Carreau-Yasuda model (Eq. 9). [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] FIG. 19. Viscous vs. viscoelastic results for the axial temperature distribution along the centreline with all effects present for the LDPE melt at 1608C at various values of apparent shear rate (L/D ¼ 33). [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] 660 POLYMER ENGINEERING AND SCIENCE DOI /pen
13 by Eq. 28 was essential in capturing this behavior. A constant b p -value not only gave wrong results but also led to divergence in the most severe cases; it seems to lead to inadmissibly high pressures because of the exponential function in Eq. 1. FIG. 20. The effect of L/D on the pressure for the LDPE melt at 1608C at various values of apparent shear rate. Non-isothermal simulations with the K-BKZ model with pressure-dependence of the viscosity (variable b p ). [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com.] Viscoelastic Modeling It is again instructive to compare individual distributions of pressure and temperature between the viscous and the viscoelastic models. This is done in Fig. 18, where we show the pressure results from the Carreau-Yasuda and the K-BKZ models for the case of L/D ¼ 33 with D ¼ cm. The viscoelastic pressures are markedly higher than the viscous ones, and this becomes more apparent for elevated shear rates. The small kinks at z/r ¼ 0 correspond to the singularity at the entry of the die, which is more severe for viscoelastic fluids than for viscous ones. The straight lines in the die region are a consequence of the counter-effects of temperature and pressure, as is well known in rheometry of polymer melts. The corresponding results for the centreline temperature distributions are shown in Fig. 19 for the same conditions. Again the temperatures from the viscoelastic runs are higher and they reach in the highest conditions a centreline value of 1768C, an increase of 168C (cf. the ultimate increase of 38.38C for L/D ¼!). The higher temperatures are a consequence of the higher viscoelastic stresses, which give rise to more viscous dissipation, especially near the walls. Collecting again all the pressure results together for the five dies with L/D ¼ 0, 2, 5,16, 33 and different apparent shear rates _c A gives in Fig. 20 the corresponding viscoelastic Bagley plot. We now see an excellent agreement between experiments and simulations. To our knowledge this is the first time that such a good agreement is obtained for LDPE. It is a combination of temperature and pressure dependence that is responsible for that. The pressure dependence with the variable b p -value given CONCLUSIONS A commercial LDPE has been studied in entry flows through tapered dies with different L/D ratios with the purpose of predicting the Bagley correction. The experiments have shown a distinct pressure-dependence of viscosity with a pressure coefficient to be a power-law function of pressure. Full rheological characterization was carried out both with a viscous (Carreau-Yasuda) and a viscoelastic (K-BKZ) model. All necessary material properties data were collected for the simulations. The viscous model was found to underestimate the extrusion pressures. The viscoelastic model showed a very good agreement with the experimental results, which appears to be the first in the literature for the elastic LDPE melt. The simulations showed that: (a) compressibility is not important in these steady flows; (b) viscous dissipation is important, especially for the more severe conditions (high L/D and apparent shear rates); and (c) the pressure-dependence of viscosity is very important and its correct function has to be found experimentally. This is the first time that all these effects are taken into account in a viscoelastic simulation. As in our recent work [18], the present viscoelastic results are unique in reaching levels of apparent shear rates as high as 1000 s 1, and showed the importance of including long lengths to fully relax the viscoelastic stresses for highly viscoelastic polymer melts, such as this LDPE, in order to capture the end correction. REFERENCES 1. J.M. Dealy and K.F. Wissbrun, Melt Rheology and its Role in Plastics Processing Theory and Applications, Van Nostrand Reinhold, New York (1990). 2. H.M. Laun, Rheol. Acta, 42, 295 (2003). 3. H.M. Laun, Rheol. Acta, 43, 509 (2004). 4. E.B. Bagley, J. Appl. Phys., 28, 193 (1957). 5. F.N. Cogswell, Trans. Soc. Rheol., 16, 383 (1972). 6. F.N. Cogswell, Polymer Melt Rheology A Guide to Industrial Practice, Wiley, New York (1981). 7. D.M. Binding, J. Non-Newtonian Fluid Mech., 27, 173 (1988). 8. D.M. Binding, J. Non-Newtonian Fluid Mech., 41, 27 (1991). 9. M. Padmanabhan and C.W. Macosco, Rheol. Acta, 36, 144 (1997). 10. K. Feigl and H.C. Öttinger, J. Rheol., 38, 847 (1994). 11. G. Barakos and E. Mitsoulis, J. Rheol., 39, 193 (1995). 12. G. Barakos and E. Mitsoulis, J. Non-Newtonian Fluid Mech., 58, 315 (1995). DOI /pen POLYMER ENGINEERING AND SCIENCE
14 13. C. Béraudo, T. Coupez, A. Fortin, Y. Demay, B. Vergnes, and J.-F. Agassant, in Proceedings of XIIth International Congress on Rheology, A. Ait-Kadi, J.M. Dealy, D.F. James, and M.C. Williams, Eds., Laval University, Quebec City, Canada, 417 (1996). 14. J. Guillet, P. Revenue, Y. Béreaux, and J.-R. Clermont, Rheol. Acta, 35, 494 (1996). 15. S.G. Hatzikiriakos and E. Mitsoulis, Rheol. Acta, 35, 545 (1996). 16. E. Mitsoulis, S.G. Hatzikiriakos, K. Christodoulou, and D. Vlassopoulos, Rheol. Acta, 37, 438 (1998). 17. E. Mitsoulis and S.G. Hatzikiriakos, Rheol. Acta, 42, 309 (2003). 18. M. Ansari, A. Alabbas, E. Mitsoulis, and S.G. Hatzikiriakos, Int. Polym. Proc., 25, 287 (2010). 19. D.J. Horrobin and R.M. Nedderman, Chem. Eng. Sci., 53, 3215 (1998). 20. A.G. Gibson, Converging Dies, in Rheological Measurements, A.A. Collyer and D.W. Gregg, Eds., Chapman and Hall, London, 455 (1998). 21. A.B. Ariawan, S. Ebnesajjad, and S.G. Hatzikiriakos, Can. J. Chem. Eng., 80, 1153 (2002). 22. I. Ochoa and S.G. Hatzikiriakos, Powder Technol, 153, 108 (2005). 23. M.A. Couch and D.M. Binding, Polymer, 41, 6323 (2000). 24. J.-Z. Liang, Polymer, 42, 3709 (2001). 25. R. Cardinaels, P. Van Puyvelde, and P. Moldenaers, Rheol. Acta, 46, 495 (2007). 26. H.E. Park, S.T. Lim, H.M. Laun, and J.M. Dealy, Macromolecules, 47, 1023 (2008). 27. Y. Son, J. Polym. Res., 16, 667 (2009). 28. J. Aho and S. Syrjälä, J. App. Polym. Sci., 117, 1076 (2010). 29. F. Koran and J.M. Dealy, J. Rheol., 43, 1279 (1999). 30. H.E. Park and J.M. Dealy, Macromolecules, 36, 5438 (2006). 31. T. Sedlacek, M. Zatloukal, P. Filip, A. Boltizar, and P. Saha, Polym. Eng. Sci., 44, 1328 (2004). 32. E.S. Carreras, N. El Kissi, J.-M. Piau, F. Toussaint, and S. Nigen, Rheol Acta, 45, 209 (2006). 33. I.J. Duvdevani and I. Klein, Soc. Plast. Eng. Trans. J., 23, 41 (1967). 34. A.V. Ramamurthy, J. Rheol., 30, 337 (1986). 35. D. Kalika and M.M. Denn, J. Rheol., 31, 815 (1987). 36. S.G. Hatzikiriakos and J.M. Dealy, J. Rheol., 35, 497 (1991). 37. S.G. Hatzikiriakos and J.M. Dealy, J. Rheol., 36, 703 (1992). 38. L.A. Archer, Wall Slip: Measurement and Modelling Issues, in Polymer Processing Instabilities, S.G. Hatzikiriakos and K.B. Migler, Eds., Marcel Dekker, New York,73 (2005). 39. M. Mooney, J. Rheol., 2, 210 (1931). 40. H.H. Winter, Adv. Heat Transfer, 13, 205 (1977). 41. R.M. Ybarra and R.E. Eckert, AIChE J., 26, 751 (1980). 42. S.M. Dinh and R.C. Armstrong, AIChE J., 28, 294 (1982). 43. J.F. Milthorpe and R.I. Tanner, Int. J. Num. Methods. Eng., 24, 263 (1987). 44. E. Mitsoulis, R. Wagner, and F.L. Heng, Polym. Eng. Sci., 28, 291 (1988). 45. R.C. Warren, Viscous Heating. Rheological Measurement, Elsevier, London, 119 (1988). 46. Y.S. Ko and A.S. Lodge, Rheol. Acta, 30, 357 (1991). 47. E.E. Rosenbaum and S.G. Hatzikiriakos, AIChE J., 43, 598 (1997). 48. M.L. Sentmanat and Goodyear Co., U.S. Patent 6,578,413 (2003). 49. M.L. Sentmanat, Rheol. Acta, 43, 657 (2004). 50. O. Delgadillo-Velázquez, S.G. Hatzikiriakos, and M. Sentmanat, Rheol. Acta, 47, 19 (2008). 51. O. Delgadillo-Velázquez, S.G. Hatzikiriakos, and M. Sentmanat, J. Polym. Sci. Part B: Polym Phys., 46, 1669 (2008). 52. R.I. Tanner, Engineering Rheology, 2nd ed., Oxford University Press, Oxford (2000). 53. E. Mitsoulis and S.G. Hatzikiriakos, J. Non-Newtonian Fluid Mech., 157, 26 (2009). 54. J.M. Dealy and R.G. Larson, Structure and Rheology of Molten Polymers-from Structure to Flow Behavior and Back Again, Hanser, Munich (2006). 55. A.C. Papanastasiou, L.E. Scriven, and C.W. Macosco, J. Rheol., 27, 387 (1983). 56. X.-L. Luo and R.I. Tanner, Int. J. Num. Methods. Eng., 25, 9 (1988). 57. T. Kajiwara, G. Barakos, and E. Mitsoulis, Int. J. Polym. Anal. Charact., 1, 201 (1995). 58. X.-L. Luo and R.I. Tanner, Rheol. Acta, 26, 499 (1987). 59. J. Meissner, Pure Appl. Chem., 42, 551 (1975). 60. S.M. Alaie and T.C. Papanastasiou, Int. Polym. Proc., 8, 51 (1993). 61. M. Beaulne and E. Mitsoulis, J. Appl. Polym. Sci., 105, 2098 (2007). 62. G. Barakos and E. Mitsoulis, J. Non-Newtonian Fluid Mech., 62, 55 (1996). 63. G.W.M. Peters and F.P.T. Baaijens, J. Non-Newtonian Fluid Mech., 68, 205 (1997). 64. S.G. Hatzikiriakos and J.M. Dealy, Polym. Eng. Sci., 34, 493 (1994). 65. R.B. Bird, R.C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids, Vol. 1: Fluid Mechanics, 2nd ed., Wiley,New York (1987). 66. A. Hannachi and E. Mitsoulis, Adv. Polym. Technol., 12, 217 (1993). 67. X.-L. Luo and E. Mitsoulis, Int. J. Num. Methods. Fluids, 11, 1015 (1990). 68. V.-H. Karl, Angew. Makromol. Chem., 79, 11 (1979). 69. S.E. Kadijk and B.H.A.A. van den Brule, Polym. Eng. Sci., 34, 1535 (1994). 70. R.C. Penwell, R.S. Porter, and S. Middleman, J. Polym. Sci., A-2, 9, 731 (1971). 71. H.M. Laun, Rheol. Acta, 22, 171 (1983). 72. G. Hay, M.E. Mackay, K.M. Awati, and Y. Park, J. Rheol., 43, 1099 (1999). 73. A. Goubert, J. Vermant, P. Moldenaers, A. Göttfert, and B. Ernst, Appl. Rheol., 11, 26 (2001). 662 POLYMER ENGINEERING AND SCIENCE DOI /pen
Entry Flow of Polyethylene Melts in Tapered Dies
REGULAR CONTRIBUTED ARTICLES M. Ansari 1, A. Alabbas 1, S. G. Hatzikiriakos 1, E. Mitsoulis 2 * 1 Department of Chemical and Biological Engineering, University of British Columbia, Vancouver, BC, Canada
More informationVISCOELASTIC SIMULATIONS WITH INTEGRAL MODELS AT EXTREMELY HIGH SHEAR RATES
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 VISCOELASTIC SIMULATIONS WITH INTEGRAL MODELS AT EXTREMELY HIGH SHEAR RATES Evan Mitsoulis School of Mining Engineering
More informationCapillary Extrusion and Swell of a HDPE Melt Exhibiting Slip
Capillary Extrusion and Swell of a HDPE Melt Exhibiting Slip MAHMOUD ANSARI Department of Chemical and Biological Engineering, The University of British Columbia, Vancouver, BC, V6T 1Z3, Canada EVAN MITSOULIS
More informationTubing Extrusion of a Fluoropolymer Melt
REGULAR CONTRIBUTED ARTICLES E. Mitsoulis 1 *, S. G. Hatzikiriakos 2 1 School of Mining Engineering and Metallurgy, National Technical University of Athens, Athens, Greece 2 Department of Chemical and
More informationTWO-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 informationAnalysis of Melt Spinning Master-Curves of Low Density Polyethylene
Analysis of Melt Spinning Master-Curves of Low Density Polyethylene Ji-Zhao Liang, 1 Lei Zhong, 1 Kejian Wang 2 1 Research Division of Green Function Materials and Equipment, School of Mechanical and Automotive
More informationCONTRIBUTION 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 informationThe Effect of Rheology in Polymer Processing: A Simulation Point of View
ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 10, 2002 The Effect of Rheology in Polymer Processing: A Simulation Point of View Evan Mitsoulis School of Mining Engineering and Metallurgy, National
More informationJournal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 157 (2009) 26 33 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm Steady flow simulations
More informationSimulation 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 informationEvaluation and comparison of routes to obtain pressure coefficients from high-pressure capillary rheometry data
DOI 10.1007/s00397-006-0148-5 ORIGINAL CONTRIBUTION Evaluation and comparison of routes to obtain pressure coefficients from high-pressure capillary rheometry data Ruth Cardinaels Peter Van Puyvelde Paula
More informationModelling the Rheology of Semi-Concentrated Polymeric Composites
THALES Project No 1188 Modelling the Rheology of Semi-Concentrated Polymeric Composites Research Team Evan Mitsoulis (PI), Professor, NTUA, Greece Costas Papoulias (Research Student), NTUA, Greece Souzanna
More informationAN 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 informationMemory Phenomena in Extrudate Swell Simulations for Annular Dies
Memory Phenomena in Extrudate Swell Simulations for Annular Dies X.-L. LUO and E. MITSOULIS, Department of Chemical Engineering, University of Ottawa, Ottawa, Ontario, Canada, KIN 9B4 Synopsis Streamline
More informationANALYSIS 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 informationJournal of Non-Newtonian Fluid Mechanics
Journal of Non-Newtonian Fluid Mechanics 177-178 (2012) 97 108 Contents lists available at SciVerse ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: http://www.elsevier.com/locate/jnnfm
More informationModeling of Non Isothermal Film Blowing Process for Non Newtonian Fluids by Using Variational Principles
Modeling of Non Isothermal Film Blowing Process for Non Newtonian Fluids by Using Variational Principles Modified on Friday, 01 May 2015 10:21 PM by mpieler Categorized as: Paper of the Month Modeling
More informationInvestigation 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 informationRheological evaluation of melt blown polymer melt
Rheological evaluation of melt blown polymer melt Jiri rabek and Martin Zatloukal Citation: AIP Conf. Proc. 1526, 237 (2013); doi: 10.1063/1.4802618 View online: http://dx.doi.org/10.1063/1.4802618 View
More informationH. 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 informationMeasuring 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 informationViscoelastic Flows in Abrupt Contraction-Expansions
Viscoelastic Flows in Abrupt Contraction-Expansions I. Fluid Rheology extension. In this note (I of IV) we summarize the rheological properties of the test fluid in shear and The viscoelastic fluid consists
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationvs. Chapter 4: Standard Flows Chapter 4: Standard Flows for Rheology shear elongation 2/1/2016 CM4650 Lectures 1-3: Intro, Mathematical Review
CM465 Lectures -3: Intro, Mathematical //6 Chapter 4: Standard Flows CM465 Polymer Rheology Michigan Tech Newtonian fluids: vs. non-newtonian fluids: How can we investigate non-newtonian behavior? CONSTANT
More informationOn 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 informationEFFECT 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 informationStability analysis of a three-layer film casting process
Korea-Australia Rheology Journal Vol. 19, No. 1, March 2007 pp. 27-33 Stability analysis of a three-layer film casting process Joo Sung ee 1, Dong Myeong Shin, Hyun Wook Jung* and Jae Chun Hyun Department
More informationShear rheology of polymer melts
Shear rheology of polymer melts Dino Ferri dino.ferri@versalis.eni.com Politecnico Alessandria di Milano, 14/06/2002 22 nd October 2014 Outline - Review of some basic rheological concepts (simple shear,
More informationEVALUATION OF NONLINEAR DIFFERENTIAL MODELS FOR THE SIMULATION OF POLYMER MELTS
1 th Fall Rubber Colloquium EVALUATION OF NONLINEAR DIFFERENTIAL MODELS FOR THE SIMULATION OF POLYMER MELTS Jochen Kroll, Stefan Turek, Patrick Westervoß Institute of Applied Mathematics (LS III), TU Dortmund
More informationHEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS
HEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS Paulo J. Oliveira Departamento de Engenharia Electromecânica, Universidade da Beira Interior Rua Marquês D'Ávila e Bolama, 600
More informationDIFFERENCE IN THERMOFORMING PROCESSABILITY OBSERVED FOR THREE HIGH IMPACT POLYSTYRENES
Page 1 of 5 DIFFERENCE IN THERMOFORMING PROCESSABILITY OBSERVED FOR THREE HIGH IMPACT POLYSTYRENES Caroline Woelfle, Kurt Koppi, Stephane Costeux, Todd Hogan, Joe Dooley, Ronald Van Daele, Alexander De
More informationAuthor's personal copy
Computers & Fluids 57 (2012) 195 207 Contents lists available at SciVerse ScienceDirect Computers & Fluids journal homepage: www.elsevier.com/locate/compfluid A study of various factors affecting Newtonian
More informationAuthors: 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 informationTemperature 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 informationTHE SUBORDINATION OF THE THREE- DIMENSIONAL FLOW INSTALLATION IN THE CONVERGING CHANNEL ON RHEOLOGICAL CHARACTERISTICS OF POLYMER STREAM
International Journal of Civil Engineering and Technology (IJCIET) Volume 9, Issue 13, December 2018, pp. 949-956, Article ID: IJCIET_09_13_095 Available online at http://www.iaeme.com/ijciet/issues.asp?jtype=ijciet&vtype=9&itype=13
More informationMeasurement and Prediction of Fluid Viscosities at High Shear Rates
Chapter 5 Measurement and Prediction of Fluid Viscosities at High Shear Rates Jeshwanth K. Rameshwaram and Tien T. Dao Additional information is available at the end of the chapter http://dx.doi.org/10.5772/54282
More informationCPGAN # 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 informationExcerpt 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 informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [HEAL-Link Consortium] On: 19 November 2008 Access details: Access Details: [subscription number 772725613] Publisher Taylor & Francis Informa Ltd Registered in England
More informationInfluence 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 informationModeling the Rheology and Orientation Distribution of Short Glass Fibers Suspended in Polymeric Fluids: Simple Shear Flow
Modeling the Rheology and Orientation Distribution of Short Glass Fibers Suspended in Polymeric Fluids: Simple Shear Flow Aaron P.R. berle, Donald G. Baird, and Peter Wapperom* Departments of Chemical
More informationMadrid, 8-9 julio 2013
VI CURSO DE INTRODUCCION A LA REOLOGÍA Madrid, 8-9 julio 2013 NON-LINEAR VISCOELASTICITY Prof. Dr. Críspulo Gallegos Dpto. Ingeniería Química. Universidad de Huelva & Institute of Non-Newtonian Fluid Mechanics
More informationON 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 informationOn the performance of enhanced constitutive models for polymer melts in a cross-slot flow
J. Non-Newtonian Fluid Mech. 82 (1999) 387±427 On the performance of enhanced constitutive models for polymer melts in a cross-slot flow Gerrit W.M. Peters *, Jeroen F.M. Schoonen, Frank P.T. Baaijens,
More informationQuasi-1D Modeling of Polymer Melt Die Swell in Short Dies. JAE-HYEUK JEONG and ARKADY I. LEONOV
Quasi-D Modeling of Polymer Melt Die Swell in Short Dies JAE-HYEUK JEONG and ARKADY I. LEONOV Department of Polymer Engineering The University of Aron Aron, Ohio 4435-030 Abstract This paper describes
More informationMechanism of gross melt fracture elimination in the extrusion of polyethylenes in the presence of boron nitride
Rheol Acta (2004) 43: 624 633 DOI 10.1007/s00397-004-0359-6 ORIGINAL CONTRIBUTION Martin Sentmanat Savvas G. Hatzikiriakos Mechanism of gross melt fracture elimination in the extrusion of polyethylenes
More informationStress Overshoot of Polymer Solutions at High Rates of Shear
Stress Overshoot of Polymer Solutions at High Rates of Shear K. OSAKI, T. INOUE, T. ISOMURA Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan Received 3 April 2000; revised
More informationTHE 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 informationFlow 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 informationA two-dimensional numerical study of the stick slip extrusion instability
J. Non-Newtonian Fluid Mech. 146 (2007) 30 44 A two-dimensional numerical study of the stick slip extrusion instability Eleni Taliadorou a,1, Georgios C. Georgiou a,, Andreas N. Alexandrou b,2 a Department
More informationModeling 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 informationViscosity overshoot in the start-up of uniaxial elongation of low density polyethylene melts
Downloaded from orbit.dtu.dk on: Mar 11, 2019 Viscosity overshoot in the start-up of uniaxial elongation of low density polyethylene melts Rasmussen, Henrik K.; Nielsen, Jens Kromann; Bach, Anders; Hassager,
More informationPerformance 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 information2009 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 informationOn the congruence of some network and pom-pom models
Korea-Australia Rheology Journal Vol 8, No, March 2006 pp 9-4 On the congruence of some network and pom-pom models Roger I Tanner* School of Aerospace, Mechanical and Mechatronic Engineering, University
More informationLecture 2: Constitutive Relations
Lecture 2: Constitutive Relations E. J. Hinch 1 Introduction This lecture discusses equations of motion for non-newtonian fluids. Any fluid must satisfy conservation of momentum ρ Du = p + σ + ρg (1) Dt
More informationNon-linear Viscoelasticity FINITE STRAIN EFFECTS IN SOLIDS
FINITE STRAIN EFFECTS IN SOLIDS Consider an elastic solid in shear: Shear Stress σ(γ) = Gγ If we apply a shear in the opposite direction: Shear Stress σ( γ) = Gγ = σ(γ) This means that the shear stress
More informationRheology of cellulose solutions. Puu Cellulose Chemistry Michael Hummel
Rheology of cellulose solutions Puu-23.6080 - Cellulose Chemistry Michael Hummel Contents Steady shear tests Viscous flow behavior and viscosity Newton s law Shear thinning (and critical concentration)
More informationPolymer rheology at high shear rate for microinjection moulding
Polymer rheology at high shear rate for microinjection moulding Cheima Mnekbi, Michel Vincent, Jean-François Agassant To cite this version: Cheima Mnekbi, Michel Vincent, Jean-François Agassant. Polymer
More informationRHEOLOGY 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 informationRheological 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 informationJournal of Non-Newtonian Fluid Mechanics
Journal of Non-Newtonian Fluid Mechanics 66 (20) 262 27 Contents lists available at SciVerse ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: http://wwwelseviercom/locate/jnnfm
More informationMeasuring the rheology of thermoplastic polymer melts
Measuring the rheology of thermoplastic polymer melts Using rotational and capillary rheometry to characterize polymer melts RHEOLOGY AND VISCOSITY Introduction Rheology is the science of studying the
More informationARTICLE IN PRESS. J. Non-Newtonian Fluid Mech. 154 (2008) Contents lists available at ScienceDirect. Journal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 154 (2008) 77 88 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm Numerical simulation of
More informationFAILURE AND RECOVERY OF ENTANGLED POLYMER MELTS IN ELONGATIONAL FLOW
FAILURE AND RECOVERY OF ENTANGLED POLYMER MELTS IN ELONGATIONAL FLOW Yogesh M. Joshi and Morton M. Denn Benjamin Levich Institute for Physico-Chemical Hydrodynamics and Department of Chemical Engineering
More informationContraction flow measurements of extensional properties
Contraction flow measurements of extensional properties Mats Stading ) and Leif Bohlin ) Chalmers University of Technology and SIK, The Swedish Institute for Food and Biotechnology, PO Box 54, S-4 9 Göteborg,
More informationRolling of bread dough: Experiments and simulations
food and bioproducts processing 8 7 (2 0 0 9) 124 138 Contents lists available at ScienceDirect Food and Bioproducts Processing journal homepage: www.elsevier.com/locate/fbp Rolling of bread dough: Experiments
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationDYNAMIC AND TRANSIENT TESTING OF ASPHALT BINDER AND PAVING MIX
66 6th RILEM Symposium PTEBM'03, Zurich, 2003 DYNAMIC AND TRANSIENT TESTING OF ASPHALT BINDER AND PAVING MIX L. Zanzotto, O.J. Vacin and J. Stastna University of Calgary, Canada Abstract: A commercially
More informationWall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes
Excerpt from the Proceedings of the COMSOL Conference 9 Boston Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes Daoyun Song *1, Rakesh K. Gupta 1 and Rajendra P. Chhabra
More informationMODELING 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 informationEffect of Molecular Weight on Viscosity of Polypropylene Melt by Capillary Rheometry
RESEARCH ARTICLE Effect of Molecular Weight on Viscosity of Polypropylene Melt by Capillary Rheometry Daniel N. Njoroge Department of Mechanical and Manufacturing Engineering Aalborg University Fibigerstraede
More information3D 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 informationCessation of Couette and Poiseuille flows of a Bingham plastic and finite stopping times
J. Non-Newtonian Fluid Mech. 29 2005) 7 27 Cessation of Couette and Poiseuille flows of a Bingham plastic and finite stopping times Maria Chatzimina a, Georgios C. Georgiou a,, Ioannis Argyropaidas b,
More informationIndex. Boundary integral method, 27 Boundary location method, 255 Brinkman number, 14, 158
Index ABFIND,275 Adaptive control, 304, 312 Adaptive process model, 315 Air drag, 221 Alternating Direction Implicit (ADI) technique, 30 Axisymmetric flow, 48 die entry, 58 Boundary integral method, 27
More informationQuantitative prediction of transient and steady-state elongational viscosity of nearly monodisperse polystyrene melts
Downloaded from orbit.dtu.dk on: Sep 27, 2018 Quantitative prediction of transient and steady-state elongational viscosity of nearly monodisperse polystyrene melts Wagner, Manfred H.; Kheirandish, Saeid;
More informationLast 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 informationCorrections 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 informationChapter 6 Molten State
Chapter 6 Molten State Rheology ( 流變學 ) study of flow and deformation of (liquid) fluids constitutive (stress-strain) relation of fluids shear flow shear rate ~ dγ/dt ~ velocity gradient dv 1 = dx 1 /dt
More informationNumerical study of 2D heat transfer in a scraped surface heat exchanger
Computers & Fluids 33 (2004) 869 880 www.elsevier.com/locate/compfluid Numerical study of 2D heat transfer in a scraped surface heat exchanger K.-H. Sun a, *, D.L. Pyle a, A.D. Fitt b, C.P. Please b, M.J.
More informationAn 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 informationCompressible viscous flow in slits with slip at the wall
Compressible viscous flow in slits with slip at the wall Georgios C. Georgiou Department of Mathematics and Statistics, University of Cyprus, Kallipoleos 75, F!O. Box 537, Nicosia, Cyprus Marcel J. Crochet
More informationCHAPTER 3. CONVENTIONAL RHEOMETRY: STATE-OF-THE-ART. briefly introduces conventional rheometers. In sections 3.2 and 3.
30 CHAPTER 3. CONVENTIONAL RHEOMETRY: STATE-OF-THE-ART This chapter reviews literature on conventional rheometries. Section 3.1 briefly introduces conventional rheometers. In sections 3.2 and 3.3, viscometers
More informationA Phenomenological Model for Linear Viscoelasticity of Monodisperse Linear Polymers
Macromolecular Research, Vol. 10, No. 5, pp 266-272 (2002) A Phenomenological Model for Linear Viscoelasticity of Monodisperse Linear Polymers Kwang Soo Cho*, Woo Sik Kim, Dong-ho Lee, Lee Soon Park, Kyung
More informationRheology and Constitutive Equations. Rheology = Greek verb to flow. Rheology is the study of the flow and deformation of materials.
Rheology and Constitutive Equations Rheology = Greek verb to flow Rheology is the study of the flow and deformation of materials. The focus of rheology is primarily on the study of fundamental, or constitutive,
More informationLecture 2. Simple shear devices. Simple shear devices 2. Simple shear devices 3. Moving plate. Velocity V. Force F. Area A. height h.
Lecture 2 Rheometry Simple shear devices Steady shear viscosity Normal stresses Oscillating shear Extensional viscosity Scalings Nondimensional parameter Simple shear devices Conceptual device for simple
More informationMorphology 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 informationChoking of liquid flows
J. Fluid Mech. (989), vol. 99, pp. 563-568 Printed in Great Britain 563 Choking of liquid flows By S. M. RICHARDSON Department of Chemical Engineering & Chemical Technology, Imperial College, London SW7.
More informationCessation of annular Poiseuille flows of Bingham plastics
J. Non-Newtonian Fluid Mech. 142 2007) 135 142 Cessation of annular Poiseuille flows of Bingham plastics Maria Chatzimina a, Christos Xenophontos a, Georgios C. Georgiou a,, Ioannis Argyropaidas b, Evan
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationFree (open) boundary condition: some experiences with viscous flow simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2012; 68:1299 1323 Published online 6 June 2011 in Wiley Online Library (wileyonlinelibrary.com)..2608 Free (open) boundary
More informationSolvent 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 informationRheometry. 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 informationCONVERGING FLOW ON-LINE RHEOMETRY FOR AN ENGINEERING EXTENSIONAL VISCOSITY OF UPVC.
CONVERGING FLOW ON-LINE RHEOMETRY FOR AN ENGINEERING EXTENSIONAL VISCOSITY OF UPVC. H. J. Ettinger, J. F. T. Pittman*, J. Sienz Centre for Polymer Processing Simulation and Design, C2EC, School of Engineering,
More informationThe effect of temperature gradients on the sharkskin surface instability in polymer extrusion through a slit die
Rheol Acta (2006) 45: 943 950 DOI 10.1007/s00397-006-0086-2 ORIGINAL CONTRIBUTION Erik Miller Sung Jin Lee Jonathan P. Rothstein The effect of temperature gradients on the sharkskin surface instability
More informationNumerical simulation of large amplitude oscillatory shear of a high-density polyethylene melt using the MSF model
J. Non-Newtonian Fluid Mech. 130 (2005) 63 76 Numerical simulation of large amplitude oscillatory shear of a high-density polyethylene melt using the MSF model P. Wapperom a,, A. Leygue b, R. Keunings
More informationEvan Mitsoulis. School of Mining Engineering and Metallurgy. National Technical University of Athens. R.R. Huilgol
Cessation of Couette and Poiseuille ows of a Bingham plastic and nite stopping times Maria Chatzimina, Georgios C. Georgiou Department of Mathematics and Statistics, University of Cyprus P.O. Box 2537,
More informationOscillatory shear rheology of polymers
Oscillatory shear rheology of polymers Dino Ferri dino.ferri@versalis.eni.com Politecnico Alessandria di Milano, 14/6/22 9 th June 215 Outline - elastic solids and viscous liquids (revision) - viscoelastic
More informationFiber spinning and draw resonance: theoretical background
Fiber spinning and draw resonance: theoretical background dr. ir. Martien Hulsen March 17, 2005 1 Introduction The fiber spinning process is nicely described in the book by Dantig & Tucker ([1], Chapter
More informationLaminar 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