Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section
|
|
- Lesley Lloyd
- 5 years ago
- Views:
Transcription
1 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 Channel with a Square Cross Section Shuichi TANOUE, Tomohiro NAGANAWA, and Yoshiyuki IEMOTO Deartment of Materials Science and Engineering, University of Fukui Bunkyo, Fukui-shi, Fukui , Jaan (Received : January 5, 2006) When a viscoelastic fluid flows in a straight channel with a noncircular cross section, a secondary flow is develoed. We discussed the effect of flow characteristics on the secondary flow of the viscoelastic fluid through a straight channel with a square cross section by the viscoelastic flow simulation using the finite element method. The Phan-Thien Tanner (PTT) model was emloyed as the constitutive equation. The second normal stress difference N 2 contributes to the secondary flow motion of a viscoelastic fluid through a straight channel with a square cross section. The degree of the secondary flow exressed by the maximum streamline function ψ max deends on the recoverable shear for the second normal stress difference N 2 /2σ, where σ is shear stress. According to the distribution of N 2 on the cross section, the flow direction of the secondary flow deends on the ressure dro generated by N 2. The elongational flow characteristics also contribute to the secondary flow state. On the whole, the degree of the secondary flow increases with the strain-hardening. However, the contribution of elongational flow characteristics to the secondary flow is smaller than that of N 2 /2σ. Key Words: Square channel flow / Viscoelastic flow simulation / Finite element method / Secondary flow 1. INTRODUCTION The secondary flow is the flow in the direction erendicular to the main flow. When a Newtonian fluid or a urely viscous fluid flows through a straight channel, no secondary flow is develoed. However, when a viscoelastic fluid flows in the straight channel with a noncircular cross section, a secondary flow is develoed. It is well known that the second normal stress difference N 2 of a viscoelastic fluid contributes to the develoment of secondary flow of the viscoelastic fluid through the straight channel with a noncircular cross section. Townsend et al. 1) discussed exerimentally the secondary flow in ie with non-circular cross section. They concluded that the major factor for the develoment of secondary flow is N 2. After that, theoretical studies using the numerical simulation were carried out by several researchers. Gervang and Larsen 2) discussed the secondary flow in straight ducts with rectangular cross section by numerical simulation using the Criminale-Ericksen-Filbey (CEF) model. 3) Xue et al. 4) investigated numerically the attern and strength of secondary flow in rectangular ies by the finite volume To whom corresondence should be addressed. tanoue@matse.fukui-u.ac.j, Tel: , Fax: method using the Phan-Thien Tanner (PTT) model. 5) Debbaut et al. 6, 7) studied the secondary flow exerimentally by observing the two layer flow of the same viscoelastic fluids and numerically by the finite element method using the Giesekus model. 8) They obtained the good agreement of their numerical results with exerimental observations. On the other hand, the secondary flow was focused by the study of co-extrusion rocess because the secondary flow contributes to the encasulation that the less viscous layer tends to encasulate the more viscous layer. Even if two layers are non-elastic fluids, the encasulation will occur. 9, 10) When two layers are viscoelastic fluids, this henomenon will occur more owerfully because a viscoelastic fluid is easy to develo the secondary flow even if mono-layer flows. Dooley et al. 11) studied the effect of the fluid viscoelasticity on the layer thickness of multilayer coextruded structure. Takase et al. 12) investigated numerically the effect of N 2 of test fluid on the encasulation in coextrusion rocess by viscoelastic flow simulation using the finite element method. They indicated that the degree of encasulation increases with N 2 of test fluid. Sunwoo et al. 13) discussed the viscoelastic coextrusion rocess by three-dimensional numerical simulation. They showed that the test fluid with large magnitude of N 2 rotruded into the test 105
2 Nihon Reoroji Gakkaishi Vol fluid with small magnitude of N 2 around the symmetric lane. By the develoment of comuter technology and calculation technique, we can do more detail investigations for the secondary flow on the viscoelastic flow through a straight channel by the viscoelastic flow simulation. In this study, we tried to discuss the effect of shear and elongational flow characteristics on the secondary flow state of the viscoelastic fluid through a straight channel with a square cross section by the viscoelastic flow simulation. The aims of this study are 1) the investigations of the effect of shear and elongational flow characteristics on the secondary flow state, and 2) the consideration of the mechanism of a secondary flow motion of the viscoelastic fluid through a straight channel with a square cross section. 2. CALCULATION THEORY 2.1 Governing Equations In this study, an isothermal, steady creeing viscoelastic flow of incomressible fluids with high viscosity is assumed. It is also assumed that the inertia is negligible because a viscoelastic fluid having high viscosity such as a olymer melt are assumed as a focused fluid and we want to consider urely the effect of shear and elongational flow characteristics on the secondary flow. The equation of continuity and the momentum equation are described as follows; (1) (2) where v is a velocity vector, a ressure, an extra stress tensor. The Phan-Thien Tanner (PTT) model 5) was used as a constitutive equation of the viscoelastic fluid. The basic form of the extra stress is shown as the following equations. (6) The notations and are resectively the uer and the lower convected time derivative defined as the following equations. (7) (8) In this study, the 1-mode tye model having one relaxation time was used because it is easy to understand the effect of flow characteristics on the secondary flow. We selected the PTT model as the constitutive equation in this study. One of the reasons is that we can consider simly the effect of N 2 and elongational flow characteristics on the secondary flow by changing one material constant (ξ or ε) of the PTT model. The detail story is described later. 2.2 Flow Modeling and Calculation Method In this study, we consider the viscoelastic flow including the secondary flow through a straight channel with a square cross section. Figure 1 shows the schematic diagram of a flow model. We used the rectangular coordinate system (x, y, z axes) for the calculations. The main flow direction was defined by the x axis, and then the secondary flow occurs on the y-z lane. In this study, we assumed that the flow geometry of the cross section, i.e., the shae of cross section on the y-z lane is the same along the main flow direction. Therefore we can assume that the ressure dro ( / x) and other hysical values such as velocity and stress excet for ressure are constant in the main flow direction. In short, we solved the fully develoed (3) (4) where η 0 is a zero shear viscosity, s is a material constant, D is a rate of deformation tensor and v is a velocity gradient tensor. The suerscrit t denotes the transose of a tensor. For the PTT model, E is given by the following equation. (5) where λ is a relaxation time, ε and ξ are material constants. G is a relaxation modulus defined as the following equation. Fig. 1. Schematic diagram of the flow model. 106
3 TANOUE NAGANAWA IEMOTO : Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section flow to the main flow direction (to x axis) in this study. The calculation region was the one-quarter region of the whole cross section of the straight channel shown in Fig. 1. The equation of continuity, the momentum equation and the constitutive equation were simultaneously solved for velocity comonents v, ressure and stress comonents E, using the Galerkin finite element method. In order to calculate the flow fields at the high shear rate, we used the under-relaxation method for the rate of deformation tensor roosed by Tanoue et al. 14) In addition, we need to calculate the ressure dro ( / x) as one of the variables. In order to calculate the ressure dro, we added the following equation of the volume flow rate Q to the governing equations: (9) where A is the cross-sectional area of a channel. Figure 2 shows the boundary conditions. We imosed the non-sli conditions on the wall, the symmetric conditions on the symmetric lanes, the x-z lane at y = 0 and the x-y lane at z = 0, and the ressure at the corner of the channel is equal to zero. The nine-node Lagrangian quadrilateral element was used in this study. The velocity and stress comonents were interolated in terms of bi-quadratic Lagrangian olynomials, while the ressure was interolated in terms of bi-linear olynomials. A set of nonlinear algebraic equations was solved by the Newton-Rahson method. the material constants, ξ = 0.2 and ε = 0.08 was defined as the base fluid. This base fluid is assumed to be the LDPE melt. The material constants for test fluids were set out as follows. 1) In Case 1, we set the material constants of the PTT model on ξ = 0.05, 0.2 and 0.5 at ε = 0.08 and s = 1/9, for consideration of the effect of the second normal stress difference on a secondary flow state. The relationshi between the normal stress differences and ξ for the PTT model is exressed as the following equation; (10) where N 1 is the first normal stress difference. Therefore the contribution of N 2 increases with ξ. Figure 3 shows the steadystate shear and lanar elongational flow characteristics of the PTT model in Case 1. Other arameters are s = 1/9 and ε = Test Fluids In this study, we focused the two flow characteristics affecting the secondary flow state of viscoelastic fluid through a straight channel: 1) the second normal stress difference, and 2) the elongational flow characteristics. The fluid having Fig. 2. Boundary conditions. Fig. 3. Steady-state flow characteristics curves of the PTT model in Case 1 (ε = 0.08 and s = 1/9). (a) Dimensionless shear viscosity η/η 0 and dimensionless rimary normal stress difference N 1 /G versus dimensionless shear rate λγ, (b) Dimensionless secondary normal stress difference N 2 /G versus dimensionless shear rate λγ, (c) Dimensionless lanar elongational viscosity η EP /η 0 versus dimensionless lanar elongation rate λε. 107
4 Nihon Reoroji Gakkaishi Vol The absolute values of the second normal stress difference N 2 increase with ξ. We can simly discuss the effect of N 2 on the flow state by changing ξ at constant ε. The first normal stress difference N 1 increases with a decrease of ξ. 2) In Case 2, we set the material constants of the PTT model on ε = 0.013, 0.08 and 0.4 at ξ = 0.2 and s = 1/9, for the consideration of the effect of the elongational flow characteristics on the secondary flow state. Figure 4 shows the steady-state shear and lanar elongational flow characteristics of the PTT model in Case 2. The strain-hardening increases with a decrease of ε. However, the shear flow characteristics are almost indeendent of ε. Therefore, we can simly discuss the effect of elongational flow characteristics on the flow state by changing ε at constant ξ. In this study, we assume the viscoelastic flow through a die on extrusion molding. Therefore, the focused region of dimensionless shear rate λγ will be less than at most The ψ R is the streamline function with dimensional unit. We used the half length of a channel width as the reference length H. Figure 6 shows ψ max as a function of the Weissenberg number We in the case using three kinds of finite focused region of dimensionless lanar elongation rate λε may be less than 1 because the velocity of secondary flow would be much smaller than that of main flow. The actual region of λε that we should focus in this study will be shown later. 3. RESULTS AND DISCUSSION 3.1 Mesh Refinement Before the detail discussion of the secondary flow, we should determine the mesh density for the calculation. In this consideration, we used the test fluid having the material constants of ξ = 0.2, ε = 0.08 and s = 1/9. Figure 5 shows the finite element meshes used in this consideration. The 3 x 3, 5 x 5 and 10 x 10 meshes were used. We estimated the maximum absolute value of dimensionless streamline function for the secondary flow ψ max. This dimensionless streamline function ψ was defined as the equation, ψ = ψ R /(<u>h), where <u> is the mean velocity of a main flow, H is the reference length and Fig. 4. Steady-state flow characteristics curve of the PTT model in Case 2 (ξ = 0.2 and s = 1/9). (a) Dimensionless shear viscosity η/η 0 and dimensionless rimary normal stress difference N 1 /G versus dimensionless shear rate λγ, (b) Dimensionless secondary normal stress difference N 2 /G versus dimensionless shear rate λγ, (c) Dimensionless lanar elongational viscosity η EP /η 0 versus dimensionless lanar elongation rate λε. Fig. 5. Finite element meshes. 108
5 TANOUE NAGANAWA IEMOTO : Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section element meshes shown in Fig.5. The Weissenberg number We is defined as the following equation. The curve in the case using the 5 x 5 mesh is almost the same as that in the case using the 10 x 10 mesh. While, the curve in case using the 3 x 3 mesh is a little different from other two cases. As a result, we selected the 5 x 5 mesh in this study. The relationshi between the absolute value of the maximum dimensionless streamline function ψ max and the Weissenberg number We will be discussed in the next section. 3.2 The Effect of the Second Normal Stress Difference on the Secondary Flow (Case 1) First, we discuss the effect of the second normal stress difference N 2, i.e., Case 1 on the secondary flow. Figure 7 Fig. 6. (11) Comarison of the absolute value of maximum dimensionless streamline function ψ max as a function of We between the finite element meshes. The material constants of the test fluid are ε = 0.08, ξ = 0.2 and s = 1/9. shows the rojection of the streamlines of secondary flow at We = 1. The numerals in each figure are the dimensionless values of streamline function estimated by dividing the streamline function by <u>h. There are two secondary flows on the y-z lane in each case. Both secondary flows are counter-rotating each other. In articular, the fluid flows in a siral in the straight channel. The absolute value of the streamline function indicates the degree of secondary flow. It is shown that the secondary flow occurred in each case. In addition, the degree of secondary flow increases with ξ, i.e., the ratio between the second normal stress difference N 2 and the first one N 1, N 2 /N 1. It imlies that the degree of secondary flow deends on the absolute value of N 2 /N 1. For comarison, there is no secondary flow for the fluid without N 2, i.e., ξ = 0. Obviously, N 2 contributes to the secondary flow in the channel flow of the viscoelastic fluid. Let s discuss why the secondary flow direction was determined as Fig.7. Figure 8(a) shows the contour of dimensionless second normal stress difference (τ yy )/ (η 0 <u>/h). The absolute values of (τ yy <u>/h) in the region near the wall (the region of y/h > ca. 0.6 or z/h > ca. 0.6) are larger than those near the center and corner of the channel. This is because the shear rate in the region near the wall is larger than that near the center and corner of the channel. In addition, the (τ yy <u>/h) distribution is shown symmetrical with resect to the line from the center (y/h, z/h ) = (0, 0) to the corner (y/h, z/h) = (1, 1). Obviously, the two secondary flows would occur symmetrically to the line from (y/h, z/h) = (0, 0) to (y/h, z/h) = (1, 1). The negative value of (τ yy <u>/h) exists along the wall at y/h = 1. In this region near the wall at y/h = 1, the comressive stress in y direction is imosed on the fluid element. This fluid element would try to swell in y direction. This leads to the ressure near the center of the Fig. 7. Contours of the steamline for secondary flow at We = 1 in Case 1. Other materials constants are ε = 0.08 and s = 1/9. 109
6 Nihon Reoroji Gakkaishi Vol channel. We can also discuss the similar story about the influence of the normal stress difference (τ yy <u>/h) along the wall at z/h = 1 on the increase of ressure near the center of the channel. As the results, the ressure on the center of the channel becomes high (See Fig.8(b).). While, the ressure gradient along the y or z direction on the symmetric line (y/h = 0 or z/h = 0) is larger than that along the line from (y/h, z/h) = (0, 0) to (y/h, z/h) = (1, 1). Then the fluid flows from the center of the channel toward the wall at z/h = 1 or y/h = 1. For those reasons, the secondary flow occurs on the cross section of the channel as shown in Fig.7. Next we consider the effect of We on the secondary flow state in Case 1. Figure 9 shows the rojection of the streamline of secondary flow at various We. The material constants of focused fluid in these figures are ξ = 0.2, ε = 0.08 and s = 1/9, i.e., these are for the base fluid. The degree of secondary flow at We = 1 is larger than that at We = 0.2. However, the degree of secondary flow at We = 10 is smaller than that at We = 1 although Fig. 8. Contours of the dimensionless stress difference (τ yy -τ zz <u>/h) and ressure /(η 0 <u>/h) at We = 1 in Case 1. The material constants of the fluid are ξ = 0.2, ε = 0.08 and s = 1/9. the absolute value of the second normal stress difference N 2 increases with shear rate, i.e., We. In order to discuss the degree of secondary flow quantitatively, Fig.10 shows the relationshi between ψ max and We. Because the large value of ψ max indicates the high ercentage of secondary flow velocity to main flow velocity, ψ max corresonds to the degree of secondary flow. The value of ψ max increases with ξ in the low We region (We < about 10). This tendency is the same as the calculation results by Xue et al. 4) In the high We region (We > about 10), it can be judged that ψ max is almost indeendent of ξ although ψ max at ξ = 0.05 is the largest at the same We. However, the value of ψ max decreases with increasing We at each material constant ξ in this region (We > about 10) although the contribution of N 2 increases with ξ. Therefore, other roerties may contribute to the secondary flow. We focused the recoverable shear for the second normal stress difference defined as N 2 /2σ, where σ is a shear stress. Fig. 9. Contours of the steamline of secondary flow at various We. The material constants of the fluid are ξ = 0.2, ε = 0.08 and s = 1/9. 110
7 TANOUE NAGANAWA IEMOTO : Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section The reason is that the shear viscosity would make some influences on the secondary flow in addition to N 2. Figure 11 shows the relationshi between N 2 /2σ and dimensionless shear rate λγ for the PTT model in Case 1. This can be calculated by using the steady-state shear flow characteristics shown in Fig.3. In the low shear rate region ( λγ < about 3 10), N 2 /2σ increases with λγ. However, N 2 /2σ shows the maximum at λγ = about 3 10 and decreases with λγ in the high shear rate region ( λγ > about 3 10). The shaes of the curves in Fig.11 are almost similar to those in Fig.10 for the same fluid although We at ψ max is different from λγ at the maximum N 2 /2σ. In addition, We at ψ max is roortional to λγ at the maximum N 2 /2σ. These imly that the degree of secondary flow would be related to not only N 2 but also N 2 /2σ. We calculated the average value of dimensionless shear rare on y-z lane at We that is the case of ψ max. The shear rate γ at each node on y-z lane is calculated using the second R invariant of the rate of deformation tensor. Here, We and at ψ max and λγ at the maximum N 2 /2σ are called We max, max and λγ max, resectively. In case of ξ = 0.05, We max = 5, max = 12 and λγ max = 13. In case of ξ = 0.2, We max = 1.5, max = 3.6 and λγ max = 5.6. In case of ξ = 0.5, We max = 0.8, = 2.0 and max λγ max = 3.6. In case of same ξ, is not max equal to λγ max. However, roortionally increases with max λγ max. This imlies that N 2 /2σ is one of the main factors for secondary flow although there are other factors for it. 3.3 The Effect of the Elongational Flow Characteristics on the Secondary Flow (Case 2) Next, we discuss the effect of the elongational flow characteristics, i.e., Case 2 on the secondary flow. Figure 12 shows the relationshi between ψ max and We. In the low We region (We < ca. 0.4), ψ max for the fluid with ε = is the smallest of all three test fluids. In the region of We > 0.4, ψ max for the fluid with ε = is almost the same as that for the fluid with ε = While, ψ max for the fluid with ε = 0.4 is the smallest of all three test fluids in the region of 0.4 < We < 2, and ψ max for the fluid with ε = 0.4 is the largest of all three test fluids in the region of We > 2. Each curve has the maximum oint, i.e., the degree of secondary flow decreases with an increase of We in the high We region. This tendency is similar to that in Case 1. Figure 13 shows the relationshi between N 2 /2σ and λγ for the PTT model in Case 2. This can be calculated by using the steady-state shear flow characteristics Fig. 10. The absolute value of maximum dimensionless streamline function ψ max as a function of We in Case 1. Other material constants are ε = 0.08 and s = 1/9. Fig. 12. The absolute value of maximum dimensionless streamline function ψ max as a function of We in Case 2. Other material constants are ξ = 0.2 and s = 1/9. Fig. 11. The recoverable shear for the second normal stress difference N 2 /2σ as a function of dimensionless shear rate λγ for the PTT model in Case 1. Fig. 13. The recoverable shear for the second normal stress difference N 2 /2σ as a function of dimensionless shear rate λγ for the PTT model in Case
8 Nihon Reoroji Gakkaishi Vol shown in Fig.4. The shaes of these curves are almost similar to that in Case 1 shown in Fig.11. This imlies that N 2 /2σ contributes to the degree of secondary flow. However the relationshi between N 2 /2σ and ψ max may not be totally corresonding. For examle, ψ max of the test fluid for the fluid with ε = 0.4 is the largest of all three fluids in We > ca. 2 although N 2 /2σ for the test fluid with ε = 0.4 is almost the same as other test fluids in the high shear rate region. This indicates that the elongational flow characteristics would make some influences of the secondary flow state. In order to consider quantitatively the effect of elongational flow characteristics on the secondary flow state, we consider the contour of the absolute value of lanar elongation rate on y-z lane. The absolute value of lanar elongation rate ε is defined as the follow equation. (12) We can consider the total value of lanar elongation rate on the secondary flow by use of ε. Figure 14 shows the contour of dimensionless lanar elongation rate λε on y-z lane for the test fluids at We = 1 and 100 in Case 2. Some elongational flow occurs along the stremaline on the secondary flow. The values of λε on y-z lane at We = 1 are lower than about 0.05 for each test fluid. The maximum value of λε on y-z lane for the test fluid with ε = and 0.08 is larger than that for the test fluid with ε = 0.4. On the other hand, the maximum value of λε on y-z lane for the test fluid with ε = and 0.08 is smaller than that for the test fluid with ε = 0.4 although the values of λε at We = 100 are lower than about 0.16 for each test fluid. According to the lanar elongational flow characteristics shown in Fig.15, the test fluids with ε = and 0.08 show the strong strain-hardening character in λε > ca It would not be easy to occur the accerelation flow in the region of λε > ca for the test fluid with ε = and This is one reason why the degree of secondary flow for the fluid with ε = 0.4 is the largest of all three test fluids at the same We in the high We region. However, the contribution of elongational flow characteristics to the secondary flow is smaller than that of second normal stress difference, i.e., the recoverable shear for second normal stress difference N 2 /2σ. Fig. 14. Contours of dimensionless lanar elongation rate λε on y-z lane in Case 2. The other material constants are ξ = 0.2 and s = 1/9. 112
9 TANOUE NAGANAWA IEMOTO : Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section to the secondary flow state. On the whole, the degree of secondary flow would increase with the strain-hardening. However, the contribution of elongational flow characteristics to the secondary flow is smaller than that of the recoverable shear N 2 /2σ. Fig. 15. Dimensionless lanar elongational viscosity η EP /η 0 versus dimensionless lanar elongation rate λε in low region in Case 2. λε 4. CONCLUSIONS We discussed the effect of flow characteristics on secondary flow states of the viscoelastic fluid through a straight channel with a square cross section by the viscoelastic flow simulation using the finite element method. The Phan-Thien Tanner (PTT) model was emloyed as the constitutive equation. The second normal stress difference N 2 contributes to the secondary flow of the viscoelastic fluid through a straight channel with a square cross section. The degree of secondary flow exressed by the absolute value of maximum streamline function ψ max deends on the recoverable shear for the second normal stress difference N 2 /2σ. In addition, according to the consideration of the distribution of N 2 on the cross section, the flow direction of secondary flow deends on the ressure dro generated by N 2. The elongational flow characteristics would also contribute REFERENCES 1) Townsend P, Walters K, Waterhouse WM, J Non-Newton Fluid Mech, 1, 176 (1976). 2) Gervang B, Larsen PS, J Non-Newton Fluid Mech, 39, 217 (1991). 3) Criminale WO, Ericksen JL, Filbey GL, Arch Rational Mech Anal, 1, 410 (1958). 4) Xue SC, Phan-Thien N, Tanner RI, J Non-Newton Fluid Mech, 59, 191 (1995). 5) Phan-Thien N, Tanner RI, J Non-Newton Fluid Mech, 2, 353 (1977). 6) Debbaut B, Avalosse T, Dooley J, Hughes K, J Non-Newton Fluid Mech, 69, 255 (1997). 7) Debbaut B, Dooley J, J Rheol, 43, 1525 (1998). 8) Giesekus H, J Non-Newton Fluid Mech, 11, 69 (1982). 9) Torres A, Hrymak AN, Vlachooulus J, Rheol Acta, 32, 513 (1993). 10) Gifford WA, Polym Engng Sci, 37, 315 (1997). 11) Dooley J, Hyun KS, Hughes K, Polym Engng Sci, 38, 1060 (1998). 12) Takase M, Kihara S, Funatsu K, Rheol Acta, 37, 634 (1998). 13) Sunwoo KB, Park SJ, Lee SJ, Ahn KH, Lee SJ, Rheol Acta, 41, 144 (2002). 14) Tanoue S, Kajiwara T, Iemoto Y, Funatsu K, Polym Engng Sci, 38, 409 (1998). 113
Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 67-73 Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section Mingkan Zhang, Xinrong Shen, Jianfeng
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 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 informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
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 informationAuthor's personal copy
J. Non-Newtonian Fluid Mech. 139 (2006) 1 20 Numerical investigation of the velocity overshoots in the flow of viscoelastic fluids inside a smooth contraction A. Afonso a, F.T. Pinho a,b, a Centro de Estudos
More informationPhan-Thien-Tanner Modeling of a Viscoelastic Fluid in the Stick-Slip Scenario. Andrew Griffith
Phan-Thien-Tanner Modeling of a Viscoelastic Fluid in the Stick-Slip Scenario Andrew Griffith Supervisor: Bruce A. Finlayson Department of Chemical Engineering University of Washington March 1, 7 Introduction
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 informationPressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids
American Journal of Alied Sciences (9): 19-195, 5 ISSN 1546-939 5 Science Publications Pressure-sensitivity Effects on Toughness Measurements of Comact Tension Secimens for Strain-hardening Solids Abdulhamid
More informationA new numerical framework to simulate viscoelastic free-surface flows with the finitevolume
Journal of Physics: Conference Series PAPER OPEN ACCESS A new numerical framework to simulate viscoelastic free-surface flows with the finitevolume method Related content - Gravitational collapse and topology
More informationMultilayer 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 informationA numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 6 (2016, 152 161 Research Article A numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows Lei Hou a, Yunqing Feng a,, Lin
More informationSimple constitutive models for linear and branched polymers
J. Non-Newtonian Fluid Mech. 116 (2003) 1 17 Simple constitutive models for linear and branched polymers Roger I. Tanner, Simin Nasseri School of Aerospace, Mechanical and Mechatronic Engineering, University
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 informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationDetermining 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 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 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 informationExact Solutions in Finite Compressible Elasticity via the Complementary Energy Function
Exact Solutions in Finite Comressible Elasticity via the Comlementary Energy Function Francis Rooney Deartment of Mathematics University of Wisconsin Madison, USA Sean Eberhard Mathematical Institute,
More informationNumerical Simulation of Particle Concentration in a Gas Cyclone Separator *
2007 Petroleum Science Vol.4 No.3 Numerical Simulation of Particle Concentration in a Gas Cyclone Searator * Xue Xiaohu, Sun Guogang **, Wan Gujun and Shi Mingxian (School of Chemical Science and Engineering,
More informationOutline. Motivation Governing equations and numerical methods Results: Discussion:
Bifurcation phenomena in strong extensional flows (in a cross-slot geometry) F. A. Cruz 1,*, R. J. Poole 2, F. T. Pinho 3, P.J. Oliveira 4, M. A. Alves 1 1 Departamento de Engenharia Química, CEFT, Faculdade
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
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 informationNumerical description of start-up viscoelastic plane Poiseuille flow
Korea-Australia Rheology Journal Vol. 1, No. 1, March 009 pp. 47-58 Numerical description of start-up viscoelastic plane Poiseuille flow Kwang Sun Park and Young don Kwon* Department of Chemical Engineering,
More informationStability analysis of constitutive equations for polymer melts in viscometric flows
J. Non-Newtonian Fluid Mech. 103 (2002) 221 250 Stability analysis of constitutive equations for polymer melts in viscometric flows Anne M. Grillet 1, Arjen C.B. Bogaerds, Gerrit W.M. Peters, Frank P.T.
More informationFinite Element Analysis of V-Bending of Polypropylene Using Hydrostatic-Pressure-Dependent Plastic Constitutive Equation*
Materials Transactions, Vol. 48, No. 1 (7). 6 to 664 #7 The Jaan Society for Technology of Plasticity Finite Element Analysis of V-Bending of Polyroylene Using Hydrostatic-Pressure-Deendent Plastic onstitutive
More informationRotating-surface-driven non-newtonian flow in a cylindrical enclosure
Korea-Australia Rheology Journal Vol. 22, No. 4, December 2010 pp. 265-272 Rotating-surface-driven non-newtonian flow in a cylindrical enclosure Yalin Kaptan 1, *, Ali Ecder 2 and Kunt Atalik 2 1 Hansung
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
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 information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationDirect Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids
Direct Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids by P. Y. Huang 1, J. Feng 2, H. H. Hu 3 and D. D. Joseph 1 1 Department of Aerospace Engineering
More informationSimulation of Die-swell Flow for Oldroyd-B Model with Feedback Semi-implicit Taylor Galerkin Finite Element Method
KMUTNB Int J Appl Sci Technol, Vol.8, No.1, pp. 55-63, (2015) Simulation of Die-swell Flow for Oldroyd-B Model with Feedback Semi-implicit Taylor Galerkin Finite Element Method Nawalax Thongjub and Vimolrat
More informationCESSATION 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 informationWeek 8 lectures. ρ t +u ρ+ρ u = 0. where µ and λ are viscosity and second viscosity coefficients, respectively and S is the strain tensor:
Week 8 lectures. Equations for motion of fluid without incomressible assumtions Recall from week notes, the equations for conservation of mass and momentum, derived generally without any incomressibility
More informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
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 informationδq T = nr ln(v B/V A )
hysical Chemistry 007 Homework assignment, solutions roblem 1: An ideal gas undergoes the following reversible, cyclic rocess It first exands isothermally from state A to state B It is then comressed adiabatically
More informationSimplifications to Conservation Equations
Chater 5 Simlifications to Conservation Equations 5.1 Steady Flow If fluid roerties at a oint in a field do not change with time, then they are a function of sace only. They are reresented by: ϕ = ϕq 1,
More informationChapter 8 Internal Forced Convection
Chater 8 Internal Forced Convection 8.1 Hydrodynamic Considerations 8.1.1 Flow Conditions may be determined exerimentally, as shown in Figs. 7.1-7.2. Re D ρumd μ where u m is the mean fluid velocity over
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 informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
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 informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationStudy of steady pipe and channel flows of a single-mode Phan-Thien Tanner fluid
J. Non-Newtonian Fluid Mech. 101 (2001) 55 76 Study of steady pipe and channel flows of a single-mode Phan-Thien Tanner fluid Manuel A. Alves a, Fernando T. Pinho b,, Paulo J. Oliveira c a Departamento
More informationInternational Journal of Mathematics Trends and Technology- Volume3 Issue4-2012
Effect of Hall current on Unsteady Flow of a Dusty Conducting Fluid through orous medium between Parallel Porous Plates with Temerature Deendent Viscosity and Thermal Radiation Harshbardhan Singh and Dr.
More informationA simple method for simulating general viscoelastic fluid flows with an alternate log-conformation formulation
J. Non-Newtonian Fluid Mech. 147 (2007) 189 199 A simple method for simulating general viscoelastic fluid flows with an alternate log-conformation formulation Oscar M. Coronado a, Dhruv Arora a, Marek
More informationAdiabatic Shear Bands in Simple and Dipolar Plastic Materials
Adiabatic Shear Bands in Simle and Diolar Plastic Materials T W \-1RIGHT us Army Ballistic Research Laboratory Aberdeen Proving Ground, MD 215 R C BATRA University of Missouri-Rolla Rolla, Missouri 6541
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 informationMODELING OF MECHANICAL RESPONSE IN CFRP ANGLE-PLY LAMINATES
HE 9 H INERNAIONAL ONFERENE ON OMPOIE MAERIAL Introduction It is known that the failure rocess in angle-ly laminate involves matrix cracking and delamination and that they exhibit nonlinear stress-strain
More informationIsogeometric analysis based on scaled boundary finite element method
IOP Conference Series: Materials Science and Engineering Isogeometric analysis based on scaled boundary finite element method To cite this article: Y Zhang et al IOP Conf. Ser.: Mater. Sci. Eng. 37 View
More informationHEAT TRANSFER ENHANCEMENT IN LAMINAR FLOW OF VISCOELASTIC FLUIDS THROUGH A RECTANGULAR DUCT
Congreso de Métodos Numéricos en Ingeniería 2009 Barcelona, 29 junio al 2 de julio 2009 SEMNI, España 2009 HEAT TRANSFER ENHANCEMENT IN LAMINAR FLOW OF VISCOELASTIC FLUIDS THROUGH A RECTANGULAR DUCT N.
More informationApplication of Automated Ball Indentation for Property Measurement of Degraded Zr2.5Nb
Journal of Minerals & Materials Characterization & Engineering, Vol. 10, No.7,.661-669, 011 jmmce.org Printed in the USA. All rights reserved Alication of Automated Ball Indentation for Proerty Measurement
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 ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE BARCELONA PLASTIC DAMAGE MODEL
XII International Conerence on Comutational Plasticity. Fundamentals and Alications COMPLAS XII E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) ANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE
More informationViscosity and Polymer Melt Flow. Rheology-Processing / Chapter 2 1
Viscosity and Polymer Melt Flow Rheology-Processing / Chapter 2 1 Viscosity: a fluid property resistance to flow (a more technical definition resistance to shearing) Remember that: τ μ du dy shear stress
More informationAn inverse problem by boundary element
1 An inverse problem by boundary element method T. Tran-Cong, T. Nguyen-Thien Faculty of Engineering and Surveying, University of Southern Queensland, Toowoomba, QLD 4350, Australia A.L. Graham Engineering
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More information5.5 The concepts of effective lengths
5.5 The concets of effective lengths So far, the discussion in this chater has been centred around in-ended columns. The boundary conditions of a column may, however, be idealized in one the following
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 informationOldroyd Viscoelastic Model Lecture Notes
Oldroyd Viscoelastic Model Lecture Notes Drew Wollman Portland State University Maseeh College of Engineering and Computer Science Department of Mechanical and Materials Engineering ME 510: Non-Newtonian
More informationUniversity Graz / Austria Institut für Chemie Volker Ribitsch
University Graz / Austria Institut für Chemie Volker Ribitsch 1 Rheology Oscillatory experiments Dynamic experiments Deformation of materials under non-steady conditions in the linear viscoelastic range
More informationCharacterization of Material Parameters
Proceedings of the World Congress on Engineering 29 Vol II WCE 29, July 1-3, 29, London, U.K. Characterization of Material Parameters S. M. Humayun Kabir, Tae-In Yeo, Sang-Ho Kim Abstract The resent work
More informationUnsteady Flow of a Dusty Conducting Fluid through porous medium between Parallel Porous Plates with Temperature Dependent Viscosity and Heat Source
Volume Issue3 3- June www.ijsret.org ISSN 78-88 Unsteady Flow of a Dusty Conducting Fluid through orous medium between Parallel Porous Plates with Temerature Deendent Viscosity and Heat Source Shalini
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 informationINTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL
Symosium on Ice (26) INTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL Mohamed O. ElSeify and Thomas G. Brown University of Calgary, Calgary, Canada ABSTRACT Current ice rubble load
More informationNumerical and experimental investigation on shot-peening induced deformation. Application to sheet metal forming.
Coyright JCPDS-International Centre for Diffraction Data 29 ISSN 197-2 511 Numerical and exerimental investigation on shot-eening induced deformation. Alication to sheet metal forming. Florent Cochennec
More informationParameters Optimization and Numerical Simulation for Soft Abrasive Flow Machining
Sensors & Transducers 014 by IFSA Publishing, S. L. htt://www.sensorsortal.com Parameters Otimization and Numerical Simulation for Soft Abrasive Flow Machining Qiaoling YUAN, Shiming JI, Donghui WEN and
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 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 informationNumerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries
Korea-Australia Rheology Journal Vol. 19, No. 4, December 2007 pp. 211-219 Numerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries Youngdon Kwon*
More informationThe effect of dynamic bending moments on the ratchetting behavior of stainless steel pressurized piping elbows
International Journal of echanical Engineering and Alications 2014; 2(2): 31-37 Published online ay 30, 2014 (htt://www.scienceublishinggrou.com/j/ijmea) doi: 10.11648/j.ijmea.20140202.12 The effect of
More informationEffect of Activation Energy and Crystallization Kinetics of Polyethylenes on the Stability of Film Casting Processes
Korea-Australia Rheology Journal Vol. 21, No. 2, June 2009 pp. 135-141 Effect of Activation Energy and Crystallization Kinetics of Polyethylenes on the Stability of Film Casting Processes Joo Sung Lee*
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationUNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume, Numer 4/0, pp. 9 95 UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
More informationEffect of accumulators on cavitation surge in hydraulic systems
Effect of accumulators on cavitation surge in hydraulic systems Donghyuk Kang, Satoshi Yamazaki, Shusaku Kagawa, Byungjin An, Motohiko Nohmi, Kazuhiko Yokota ISROMAC 6 International Symosium on Transort
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
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 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 informationElastic Model of Deformable Fingertip for Soft-fingered Manipulation
IEEE TRANSACTION ON ROBOTICS, VOL. 1, NO. 11, NOVEMBER 25 1 Elastic Model of Deformable Fingerti for Soft-fingered Maniulation Takahiro Inoue, Student Member, IEEE, and Shinichi Hirai, Member, IEEE Abstract
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
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 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 informationVIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES
Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationPhase transition. Asaf Pe er Background
Phase transition Asaf Pe er 1 November 18, 2013 1. Background A hase is a region of sace, throughout which all hysical roerties (density, magnetization, etc.) of a material (or thermodynamic system) are
More informationOn the numerical algorithm for isotropic kinematic hardening with the Armstrong Frederick evolution of the back stress
Comut. Methods Al. Mech. Engrg. 191 (2002) 3583 3596 www.elsevier.com/locate/cma On the numerical algorithm for isotroic kinematic hardening with the Armstrong Frederick evolution of the back stress Vlado
More informationNonlinear Static Analysis of Cable Net Structures by Using Newton-Raphson Method
Nonlinear Static Analysis of Cable Net Structures by Using Newton-Rahson Method Sayed Mahdi Hazheer Deartment of Civil Engineering University Selangor (UNISEL) Selangor, Malaysia hazheer.ma@gmail.com Abstract
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 informationISSN Article
Entropy 23, 5, 28-299; doi:.339/e5628 OPEN ACCESS entropy ISSN 99-43 www.mdpi.com/journal/entropy Article Analysis of Entropy Generation Rate in an Unsteady Porous Channel Flow with Navier Slip and Convective
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 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 informationModeling Strong Extensional Flows of Polymer Solutions and Melts
Modeling Strong Extensional Flows of Polymer Solutions and Melts Antony N. Beris CSCAMM Workshop Multiscale Modeling and Simulation of Complex Fluids April 19, 27 Acknowledgements: Dr. Joydeep Mukherjee,
More informationA J estimation scheme for surface cracks in piping welds
2th International Conference on Structural Mechanics in Reactor Technology (SMiRT 2) Esoo, Finland, August 9-14, 29 SMiRT 2-Division II, Paer 1672 A J estimation scheme for surface cracks in iing welds
More informationHypersonic flow: introduction
Hyersonic flow: introduction Van Dyke: Hyersonic flow is flow ast a body at high ach number, where nonlinearity is an essential feature of the flow. Also understood, for thin bodies, that if is the thickness-to-chord
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 informationSimilarity Approach to the Problem of Second Grade Fluid Flows over a Stretching Sheet
Applied Mathematical Sciences, Vol. 1, 2007, no. 7, 327-338 Similarity Approach to the Problem of Second Grade Fluid Flows over a Stretching Sheet Ch. Mamaloukas Athens University of Economics and Business
More informationApplication Of Optimal Homotopy Asymptotic Method For Non- Newtonian Fluid Flow In A Vertical Annulus
Application Of Optimal Homotopy Asymptotic Method For Non- Newtonian Fluid Flow In A Vertical Annulus T.S.L Radhika, Aditya Vikram Singh Abstract In this paper, the flow of an incompressible non Newtonian
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 information