Evan Mitsoulis. School of Mining Engineering and Metallurgy. National Technical University of Athens. R.R. Huilgol
|
|
- Winfred Parker
- 6 years ago
- Views:
Transcription
1 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, 678 Nicosia, CYPRUS Tel.: Fax: georgios@ucy.ac.cy Evan Mitsoulis School of Mining Engineering and Metallurgy National Technical University of Athens Heroon Polytechniou 9, 57 8 Zografou, Athens, GREECE R.R. Huilgol School of Informatics and Engineering Flinders University of South Australia G.P.O. Box 2, Adelaide, SA 5, AUSTRALIA Short title: CESSATION OF BINGHAM FLOWS
2 Abstract We solve the one-dimensional cessation Couette and Poiseuille ows of a Bingham plastic using the regularized constitutive equation proposed by Papanastasiou and employing nite elements in space and a fully implicit scheme in time. The numerical calculations conrm previous theoretical ndings that the stopping times are nite when the yield stress is nonzero. The decay of the volumetric ow rate, which is exponential in the Newtonian case, is accelerated and eventually becomes linear as the yield stress is increased. In all ows studied, the calculated stopping times agree very well with the theoretical upper bound estimates. KEYWORDS : Couette Flow Poiseuille Flow Bingham Plastic Papanastasiou Model Cessation. 2
3 Introduction In viscometric ows, one can bring the uid to a halt by setting the moving boundary to rest, in the case of Couette ows, or by reducing the applied pressure gradient tozeroinpoiseuille ows. In a Newtonian uid, the corresponding velocity elds decay to zero in an innite amount of time []. In a Bingham plastic, the velocity elds go to zero in a nite time, which emphasizes the role of the yield stress [2]. Glowinski [3] and Huilgol and co-workers [2, 4] provided explicit theoretical nite upper bounds on the time for a Bingham material to rest in various ows, such as the plane and circular Couette ows and the plane and axisymmetric Poiseuille ows. In each case, the theoretical bound depends on the density, the viscosity, the yield stress, and the least eigenvalue of the Laplacian operator on the ow domain. More recently, Huilgol [5] has also derived upper bounds for the cessation of round Poiseuille ow of more general viscoplastic uids. The objective of the present work is to compute numerically the stopping times and make comparisons with the theoretical upper bounds provided in the literature for the cessation of three ows of a Bingham uid: (a) the plane Couette ow (b) the plane Poiseuille ow and (c) the axisymmetric Poiseuille ow. Instead of the ideal Bingham-plastic constitutive equation, we employ the regularized equation proposed by Papanastasiou [6], in order to avoid the need of determining a priori the yielded and unyielded regions in the ow domain. It should be noted that preliminary results for the case of the plane Poiseuille ow can also be found in Ref. [7]. The paper is organized as follows. In section 2, we discuss the regularized Papanastasiou equation for a Bingham plastic. In section 3, we present the dimensionless forms of the governing equations for the three ows of interest along with the corresponding theoretical 3
4 stopping times. In section 4, we present and discuss representative numerical results for all ows. The numerical stopping times agree very well with the theoretical upper bounds. Some discrepancies are observed only for low Bingham numbers when the growth parameter in the Papanastasiou model is not suciently high. Finally, section 5 contains the conclusions of this work. 2 Constitutive equation Let u and denote the velocity vector and the stress tensor, respectively, and _ denote the rate-of-strain tensor, _ ru +(ru) T () where ru is the velocity-gradient tensor, and the superscript T denotes its transpose. The magnitudes of _ and are respectively dened as follows: _ = r 2 II _ = r 2 _ : _ and = r 2 II = r 2 : (2) where II stands for the second invariant of a tensor. In tensorial form, the Bingham model is written as follows: 8 >< >: _ = = _ + _ (3) where is the yield stress, and is a constant viscosity. In anyow of a Bingham plastic, determination of the yielded ( ) and unyielded ( ) regions in the ow eld is necessary, which leads to considerable computational diculties in the use of the model. These are overcome by using the regularized constitutive equation 4
5 proposed by Papanastasiou [6]: = [ ; exp(;m _)] _ + _ (4) where m is a stress growth exponent. For suciently large values of the regularization parameter m, the Papanastasiou model provides a satisfactory approximation of the Bingham model, while at the same time the need of determining the yielded and the unyielded regions is eliminated. The model has been used with great success in solving various steady and time-dependent ows (see, for example, [8, 9] and references therein). 3 Flow problems and governing equations The governing equations along with the boundary and initial conditions of the three timedependent, one-dimensional Bingham-plastic ows of interest are discussed below. The theoretical upper bounds for the stopping times are also presented. 3. Cessation of plane Couette ow The geometry of the plane Couette ow is shown in Fig. a. The steady-state solution is given by u s x(y) = ; y V (5) H where V is the speed of the lower plate (the upper one is kept xed) and H is the distance between the two plates. We assume that at t=, the velocity u x (y t) isgiven by the above prole and that at t= + the lower plate stops moving. To non-dimensionalize the x-momentum equation, we scale the lengths by H, the velocity by V, the stress components by V=H, and the time by H 2 =, where is the constant density of the uid. With these 5
6 scalings, the x-momentum equation : (6) The dimensionless form of the Papanastasiou model is reduced to yx = Bn [ ; exp(;m _)] (7) where _=j@u x =@yj, is the Bingham number, and is the dimensionless growth parameter. Bn H V M mv H (8) (9) The dimensionless boundary and initial conditions are as follows: u x ( t)= t> u x ( t)= t u x (y ) = ; y y 9 >= > : () In the case of a Newtonian uid (Bn=), the analytical solution of the time-dependent ow, governed by Eqs. (6), (7) and (), is known []: u x (y t) = 2 X k= k sin (ky) e;k2 2 t : () Hence, the ow ceases theoretically at innite time. If the uid is a Bingham plastic (Bn > ), however, the ow comes to rest in a nite amount of time, as demonstrated by Huilgol et al. [2], who provide the following upper bound for the dimensionless stopping time: " # T f = jju x (y )jj 2 2 Bn (2) 6
7 where jju x (y )jj = Z u 2 x(y ) dy =2 : (3) 3.2 Cessation of plane Poiseuille ow The geometry of the plane Poiseuille ow is depicted in Fig. b. The steady-state solution is given by u s x (y) = 8 >< >: 2 s (H ; y ) 2 y s (H 2 ; y 2 ) ; where (;@p=@x) s is the pressure gradient, and (H ; y) y y H (4) y = <H (5) (;@p=@x) s denotes the point at which the material yields. Note that ow occurs only if (;@p=@x) s > H. The volumetric ow rate is given by Q = 2W s H 3 " ; 3 2 y H + y 3 # (6) 2 H where W is the width of the plates (in the z-direction). We assume that at t= the velocity u x (y t)isgiven by the steady-state solution (4) and that at t= + the pressure gradient isvanished, or reduced to (;@p=@x) < (;@p=@x) s, in which case the ow is expected to stop. The evolution of the velocity is again governed by the x-momentum equation. Using the same scales as in the plane Couette ow, with V denoting now the mean velocity in the slit, the dimensionless form of the x-momentum equation is = f (7) 7
8 where f denotes the dimensionless pressure gradient. The dimensionless form of the constitutive equation is given by Eq. (7). The dimensionless steady velocity prole becomes: u s x(y) = 8 >< 2 f s ( ; y ) 2 y y (8) >: 2 f s ( ; y 2 ) ; Bn ( ; y) y y where y = Bn f s : (9) It turns out that y is the real root of the cubic equation: y 3 ; y +2=: (2) Bn It is clear that a steady ow in the channel occurs only if f s >Bn. The dimensionless boundary and initial conditions for the time-dependent problem ( t)= t u x ( t)= t u x (y ) = u s x(y) y 9 >= > : (2) In the case of Newtonian ow (Bn=), the time-dependent solution is given by [] u x (y t) = 48 3 X k= (;) k+ (2k (2k ; ) cos ; ) y exp 3 2 " ; (2k ; )2 2 4 t # (22) which indicates that the ow stops only after an innite amount of time. In the case of a Bingham plastic (Bn > ), Huilgol et al. [2] provide the following estimate for the stopping time: T f = 4 2 " # jju x (y )jj f<bn: (23) Bn ; f The above estimate is valid when f < Bn(or, equivalently, when f < f s ) otherwise, the ow will not stop. 8
9 3.3 Cessation of axisymmetric Poiseuille ow The geometry of the axisymmetric Poiseuille ow is depicted in Fig. c. The steady-state solution is given by u s z(r) = 8 >< >: 4 s (R ; r ) 2 r r (R 2 ; r 2 ) ; (R ; r) r r R (24) where (;@p=@z) s is the pressure gradient, and the yield point isgiven by r = 2 <R: (25) (;@p=@z) s The volumetric ow rate is given by Q = s R 4 " ; 4 3 r R + r 4 # : (26) 3 R We assume that at t= the velocity u z (r t)isgiven by the steady-state solution and that at t= + the pressure gradient isvanished, or reduced to (;@p=@z) < (;@p=@z) s. Scaling the lengths by the tube radius R, the velocity by the mean velocity V, the pressure and the stress components by V=R, and the time by R 2 =, we obtain the dimensionless form of the z-momentum = f (r rz) (27) where f is the dimensionless pressure gradient. The dimensionless form of the constitutive equation is given by rz = Bn [ ; exp(;m _)] (28) where _=j@u z =@rj, Bn R V (29) 9
10 and M mv R : (3) The dimensionless steady velocity prole takes the form 8 u s z(r) = >< 4 f s ( ; r ) 2 r r (3) >: 4 f s ( ; r 2 ) ; Bn ( ; r) r r where r satises r = 2Bn f s (32) and r 4 ; r +3=: (33) Bn Note that a steady ow in the tube occurs only if f s > 2Bn. The growth of r with Bn is illustrated in Fig. 2, in which steady-state velocity proles calculated for various Bingham numbers are shown. The dimensionless boundary and initial conditions ( t)= t u z ( t)= t u z (r ) = u s z(r) r 9 >= > : (34) The time-dependent solution for Newtonian ow (Bn=) is given by [] u z (r t)=6 X k= J (a k r) a 3 k J (a k ) e;a2 k t (35) where J and J are respectively the zeroth- and rst-order Bessel functions of the rst kind, and a k, k=, 2, :::are the roots of J. In the case of a Bingham plastic (Bn > ), Glowinski
11 [3] provides the following estimate for the stopping time: T f = jju z (r )jj + 2Bn ; f f<2bn : (36) where jju z (r )jj = 2 Z u 2 z(r ) rdr =2 (37) and is the smallest (positive) eigenvalue of the problem: d r dw r dr dr + w = w () = w() = : (38) It is easily found that =a 2 ' 5:783, where a is the least root of J (x), with the corresponding eigenfunction being given by w (x)=j (a x). Therefore, T f = a 2 +a 2 jju z (r )jj f<2bn : (39) 2Bn ; f The estimate (39) holds only when f<2bn (or, equivalently, when f<f s ) otherwise, the ow will not stop. 4 Numerical results Since there are no analytical solutions to the ows under study, in the case of the Bingham plastic or the Papanastasiou model, we have used a numerical method, namely the nite element method with quadratic (P 2 -C ) elements for the velocity. For the spatial discretization, we used the Galerkin form of the momentum equation. For the time discretization, we used the standard fully-implicit (Euler backward-dierence) scheme. At each time step, the nonlinear system of discretized equations was solved using the Newton method. In the case of Couette ow, a 2-element mesh rened near the two plates has been used. In Poiseuille ows, the mesh consisted of 8 elements and was rened near the wall. The code has been
12 tested by solving rst the Newtonian ows and making comparisons with the analytical solutions. In all three problems, the agreement between the theory and the calculations was excellent. Cessation of plane Couette ow Figures 3-5, show the evolution of the velocity forbn= (Newtonian uid), 2 and 2, respectively. The growth parameter has been taken to be M =2. The numerical solution in Figure 3 compares very well with the analytical solution () for the Newtonian ow. The numerical solutions for Bingham ow (Figs. 4 and 5), show that a small unyielded region, where the velocity is at, appears near the moving plate. Note that for high Bingham numbers (i.e., Bn > 5), very small time steps (of the order of ;9 )were necessary in order to get convergence in the early stages of cessation. The size of the unyielded region increases as the time proceeds. Its left limit initially moves to the right but at higher times starts moving to the left, as the ow approaches complete cessation. Figure 6 shows the evolution of the volumetric ow rate, Q(t) = Z u x (y t) dy (4) for various Bingham numbers. These curves show the dramatic eect of the yield stress, which accelerates the cessation of the ow. In the Newtonian case (Bn=) and for small Bingham numbers the decay of the volumetric ow rate is exponential, at least initially. At higher Bingham numbers, the decay ofq becomes polynomial and eventually linear. The times at which Q= ;3 and ;5 are plotted as functions of the Bingham number in Fig. 7. The two times coincide for moderate or large Bingham numbers, which indicates that the ow indeed stops at a nite time. In order to make comparisons with the theoretical upper bound (2) we consider as numerical stopping time the time at which Q= ;5. The comparison between 2
13 calculations and theory, provided in Fig. 8, shows a verygoodagreement for moderate and higher Bingham numbers. The small discrepancies observed for low Bingham numbers are due to the fact that the value of M is not suciently high, as discussed below. For very small Bn, the eect of M is not crucial, since the material is practically Newtonian, which explains why the calculations fall again below the theoretical upper bound, as they should. Cessation of plane Poiseuille ow Figures 9-2 show the evolution of u x (y t) forbn= (Newtonian case),, 5 and 2. In Fig. 3, we see the evolution of the calculated volumetric ow rate for various Bingham numbers. As in plane Couette ow, the decay of the volumetric ow rate is exponential for small Bingham numbers and becomes polynomial at higher Bn values. Figure 4 shows plots of the times required for the volumetric ow ratetobecome ;3 and ;5 versus the Bingham number for both plane and round Poiseuille ows (with M =2). Before proceeding to the comparisons with the theoretical estimate (23), let us investigate the eect of the growth parameter M on the calculated stopping times. As demonstrated in Fig. 5, which shows results obtained with M=2 and M=5, the calculated stopping times are not so sensitive tom when the Bingham number is moderate or high, i.e. Bn. For smaller Bingham numbers, i.e. ;3 Bn, the time required for the volumetric ow rate to become ;5 is reduced as M is increased. For very small Bingham numbers, the uid is essentially Newtonian, and therefore the value of M has no eect on the calculations. Hence, in order to get converged results in the range ;3 Bn, the value of M has to be increased further. However, our studies showed that when M =, convergence diculties are observed of Bn > :. One way toresolve the problem is to reduce the time step. This might be good for extending the calculations to a slightly higher Bn. Beyond this critical Bn value, the required time step is very small and the accumulated round-o errors are so high so 3
14 that the error in the calculated stopping time is higher than that corresponding to a smaller value of M. As a conclusion, decreasing the time step is not the best way to obtain converged results for M>5. If results for small (but not vanishingly small) Bingham numbers and large M are necessary, then continuation in M must be used at each time step. According to our numerical experiments, using such a continuation will increase the computational time by at least times. Since we arenotinterested in so small values of Bn, such calculations have not been pursued. A comparison between theory and calculations is provided in Fig. 6 for the case f= (i.e., when the imposed pressure gradient is set to zero). Again, the agreement between theory and computations is excellent for moderate and high Bingham numbers. Again, the small discrepancies observed for small values of the Bingham number are due to the fact that M is not suciently high. We have also examined the case in which the imposed pressure gradient f is not zero. An ideal Bingham plastic stops after a nite time if f Bn and reaches a new steady-state if f>bn(with the volumetric ow rate corresponding to the new value of f). This is not the case with a regularized Bingham uid. Since M is nite, the ow will reach a new steady-state in which the volumetric ow ratemay be small but not zero. To illustrate this eect, we considered the case in which Bn= and M =5 and carried out simulations for dierent values of f. In Fig. 7a, we see the evolution of the volumetric ow rate for dierent values of f. Figure 7b is a zoom of the previous gure showing that indeed the volumetric ow rate reaches a nite value when f 6=. Thisvalue may be reduced further by increasing the value of M. The new volumetric ow rate is plotted against f in Fig. 8. Finally, in Fig. 9 we compare the times required to reach Q= ;3 with the theoretical estimate (23). For smaller values of Q, the numerical results move closer to the theoretical curve, but in a smaller range of f, as it is easily deduced from Fig. 23. The deviations between theory and 4
15 experiment become larger as the value of the imposed pressure gradient is increased. These, however, can be further reduced by increasing the value of M. Cessation of axisymmetric Poiseuille ow The results for the axisymmetric Poiseuille ow arevery similar to those obtained for the planar case. Figures 2-23 show the evolution of u x for M=2 and Bn= (Newtonian case),, 5 and 2. In Fig. 24, we zoom near the wall in order to see how the velocity prole changes when Bn=2. It is clear that a second unyielded region of a smaller size appears near the wall, in which the velocity is zero. In Fig. 25, we see the evolution of the calculated volumetric ow rate (scaled by 2), Q(t) = Z u z (r t) rdr (4) for M =2 and various Bingham numbers. As in plane Poiseuille ow, the cessation of the ow is accelerated as the Bingham number is increased. The calculated stopping times for Q= ;5 and f=, plotted versus the Bingham number in Fig. 26, agree well with the theoretical estimate (39), with the small discrepancies observed for low Bn excepted. 5 Conclusions The Papanastasiou modication of the Bingham model has been employed in order to solve numerically the cessation of plane Couette, plane Poiseuille, and axisymmetric Poiseuille ows of a Bingham plastic. The nite element calculations showed that the volumetric ow rate decreases exponentially for low, polynomially for moderate, and linearly for high Bingham numbers. Unlike their counterparts in a Newtonian uid, the corresponding times for complete cessation are nite, in agreement with theory. The numerical stopping times are 5
16 found to be in very good agreement with the theoretical upper bounds provided in Refs. [3, 2], for moderate and higher Bingham numbers. Some minor discrepancies observed for rather low Bingham numbers can be reduced by increasing the regularization parameter introduced by the Papanastasiou model. A noteworthy dierence between the predictions of the ideal and the regularized Bingham model is revealed when the imposed pressure gradient is nonzero and below the critical value at which a nonzero steady-state Poiseuille solution exists. In contrast to the ideal Bingham ow, which reaches complete cessation at a nite time, the regularized ow reaches a velocity prole corresponding to a small but nonzero volumetric ow rate. The value of the latter may be reduced by increasing the value of the regularization parameter M but will always be nonzero. Acknowledgement This research was partially supported by the Research Committee of the University of Cyprus. References [] T. Papanastasiou, G. Georgiou and A. Alexandrou, Viscous Fluid Flow, CRC Press, Boca Raton, 999. [2] R.R. Huilgol, B. Mena, and J.M. Piau, Finite stopping time problems and rheometry of Bingham uids, J. Non-Newtonian Fluid Mech. 2 (22) 97. [3] R. Glowinski, Numerical Methods for Nonlinear Variational Problems, Springer-Verlag, New York,
17 [4] R.R. Huilgol, and B. Mena, On kinematic conditions aecting the existence and nonexistence of a moving yield surface in unsteady unidirectional ows of Bingham uids, J. Non-Newtonian Fluid Mech. 23 (24) 25. [5] R.R. Huilgol, Variational inequalities in the ows of yield stress uids including inertia: Theory and applications, Phys. Fluids 4 (22) 269. [6] T.C. Papanastasiou, Flows of materials with yield, J. Rheology 3 (987) 385. [7] M. Chatzimina, G.C. Georgiou, E. Mitsoulis, and R.R. Huilgol, Finite stopping times in Couette and Poiseuille ows of viscoplastic uids, in: Proceedings of the XIVth Int. Cong. Rheol., Seoul, Korea, pp. NFF22--NF22-4. [8] Y. Dimakopoulos and J. Tsamopoulos, Transient displacement of a viscoplastic material by air in straight and suddenly constricted tubes, J. Non-Newtonian Fluid Mech. 2 (23) 43. [9] E. Mitsoulis and R.R. Huilgol, Entry ows of Bingham plastics in expansions, J. Non- Newtonian Fluid Mech. 22 (24) 45. 7
18 FIGURE CAPTIONS. Flow problems under study: (a) cessation of plane Couette ow (b) cessation of plane Poiseuille ow and (c) cessation of axisymmetric Poiseuille ow. 2. Steady velocity distributions for various Bingham numbers in round Poiseuille ow M=2. 3. Evolution of the velocity in cessation of plane Couette ow of a Newtonian uid. Comparison of the analytical (solid lines) with the numerical (dashed lines) solutions. 4. Evolution of the velocity in cessation of plane Couette ow of a Bingham uid with Bn=2 and M=2. 5. Evolution of the velocity in cessation of plane Couette ow of a Bingham uid with Bn=2 and M=2. 6. Evolution of the volumetric ow rate during the cessation of plane Couette ow of a Bingham uid with M=2 and various Bingham numbers. 7. Calculated times for Q= ;3 and ;5 in cessation of plane Couette ow of a Bingham uid with M=2. 8. Comparison of the computed stopping time (Q= ;5 ) in cessation of plane Couette ow of a Bingham uid with the theoretical upper bound M=2. 9. Evolution of the velocity in cessation of plane Poiseuille ow of a Newtonian uid.. Evolution of the velocity in cessation of plane Poiseuille ow of a Bingham uid with Bn= and M=2. 8
19 . Evolution of the velocity in cessation of plane Poiseuille ow of a Bingham uid with Bn=5 and M=2. 2. Evolution of the velocity in cessation of plane Poiseuille ow of a Bingham uid with Bn=2 and M=2. 3. Evolution of the volumetric ow rate during the cessation of plane Poiseuille ow of a Bingham uid with M=2 and various Bingham numbers. 4. Calculated times for Q= ;3 and ;5 in cessation of plane and round Poiseuille ows of Bingham uids with M=2. 5. Calculated times for Q= ;3 and ;5 in cessation of plane Poiseuille ow of Bingham uids with M=2 (dashed) and M=5 (solid). 6. Comparison of the computed stopping time (Q= ;5 ) in cessation of plane Poiseuille ow of a Bingham uid with the theoretical upper bound f= and M=5. 7. (a) Evolution of the volumetric ow rate for various values of the imposed pressure gradient f (b) Zoom of the same plot showing that a nite volumetric ow rate is reached when f> plane Poiseuille ow with Bn= and M=5. 8. Volumetric ow rates reached with the regularized Papanastasiou model versus the imposed pressure gradient f plane Poiseuille ow, Bn= and M=5. 9. Comparison of the times required to reach Q= ;3 in cessation of plane Poiseuille ow of a regularized Bingham uid with the theoretical estimate of Huilgol et al. (22) for an ideal Bingham uid Bn= and M=5. 2. Evolution of the velocity in cessation of round Poiseuille ow of a Newtonian uid. 9
20 2. Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn= and M= Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn=5 and M= Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn=2 and M= Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn=2 and M=2 (zoom near the wall). 25. Evolution of the volumetric ow rate during the cessation of round Poiseuille ow of a Bingham uid with M=2 and various Bingham numbers. 26. Comparison of the computed stopping time (Q= ;5 ) in cessation of round Poiseuille ow of a Bingham uid with the theoretical upper bound f= and M=5. 2
21 u x = u x = H u x (y ) u x (y t) y u x =V u x = x t= t> (a) H u x (y ) u x (y t) y = t= t> (b) R u z (r ) u z (r t) r = t= t> (c) Figure : Flow problems under study: (a) cessation of plane Couette ow (b) cessation of plane Poiseuille ow and (c) cessation of axisymmetric Poiseuille ow. 2
22 2.5 u z 2.5 Bn= (Newtonian uid) Bn= Bn=3 Bn= Bn= (Solid) r Figure 2: Steady velocity distributions for various Bingham numbers in round Poiseuille ow M =2..2 u x t= y Figure 3: Evolution of the velocity in cessation of plane Couette ow of a Newtonian uid. Comparison of the analytical (solid lines) with the numerical (dashed lines) solutions. 22
23 .2 u x.8.6 t= y Figure 4: Evolution of the velocity in cessation of plane Couette ow of a Bingham uid with Bn=2 and M =2..2 u x.8.6 t= y Figure 5: Evolution of the velocity in cessation of plane Couette ow of a Bingham uid with Bn=2 and M =2. 23
24 .6 Q Bn= t Figure 6: Evolution of the volumetric ow rate during the cessation of plane Couette ow of a Bingham uid with M =2 and various Bingham numbers. T f Q= ;5 Q= ;3.... Bn Figure 7: Calculated times for Q= ;3 and ;5 in cessation of plane Couette ow of a Bingham uid with M =2. 24
25 T f Theory. Numerical (Q= ;5 )... Bn Figure 8: Comparison of the computed stopping time (Q= ;5 ) in cessation of plane Couette ow of a Bingham uid with the theoretical upper bound M =2. 2 u x.5 t= y Figure 9: Evolution of the velocity in cessation of plane Poiseuille ow of a Newtonian uid. 25
26 2 u x.5.5 t= y Figure : Evolution of the velocity in cessation of plane Poiseuille ow of a Bingham uid with Bn= and M =2. 2 u x.5 t= y Figure : Evolution of the velocity in cessation of plane Poiseuille ow of a Bingham uid with Bn=5 and M =2. 26
27 2 u x.5.5 t= y Figure 2: Evolution of the velocity in cessation of plane Poiseuille ow of a Bingham uid with Bn=2 and M =2. Q Bn= t Figure 3: Evolution of the volumetric ow rate during the cessation of plane Poiseuille ow of a Bingham uid with M =2 and various Bingham numbers. 27
28 T f Q=.... Round Poiseuille ow Plane Poiseuille ow.... Bn Figure 4: Calculated times for Q= ;3 and ;5 in cessation of plane and round Poiseuille ows of Bingham uids with M =2. T f Q= ;5 Q= ;3.... Bn Figure 5: Calculated times for Q= ;3 and ;5 in cessation of plane Poiseuille ow of Bingham uids with M =2 (dashed) and M =5 (solid). 28
29 T f Theory. Numerical (Q= ;5 )... Bn Figure 6: Comparison of the computed stopping time (Q= ;5 ) in cessation of plane Poiseuille ow of a Bingham uid with the theoretical upper bound f = and M =5. 29
30 Q f = t.5 Q...5 f = t Figure 7: (a) Evolution of the volumetric ow rate for various values of the imposed pressure gradient f (b) Zoom of the same plot showing that a nite volumetric ow rate is reached when f> plane Poiseuille ow with Bn= and M =5. 3
31 .5 Q f f Figure 8: Volumetric ow rates reached with the regularized Papanastasiou model versus the imposed pressure gradient f plane Poiseuille ow, Bn= and M =5. 3 T f 2 Estimate for ideal Bingham uid Regularized Bingham uid (Q= ;3 ) f Figure 9: Comparison of the times required to reach Q= ;3 in cessation of plane Poiseuille ow of a regularized Bingham uid with the theoretical estimate of Huilgol et al. (22) for an ideal Bingham uid Bn= and M =5. 3
32 u x t= y Figure 2: Evolution of the velocity in cessation of round Poiseuille ow of a Newtonian uid. 2.5 u x t= y Figure 2: Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn= and M =2. 32
33 2.5 u x 2.5 t= y Figure 22: Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn=5 and M = u x t= y Figure 23: Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn=2 and M =2. 33
34 .25 u x t= y Figure 24: Evolution of the velocity in cessation of round Poiseuille ow of a Bingham uid with Bn=2 and M =2 (zoom near the wall). Q Bn= t Figure 25: Evolution of the volumetric ow rate during the cessation of round Poiseuille ow of a Bingham uid with M =2 and various Bingham numbers. 34
35 T f Theory. Numerical (Q= ;5 )... Bn Figure 26: Comparison of the computed stopping time (Q= ;5 ) in cessation of round Poiseuille ow of a Bingham uid with the theoretical upper bound f = and M =5. 35
Cessation 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 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 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 informationJournal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 165 (2010) 544 550 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm Short communication Numerical
More informationUnsteady circular Couette flow of a Bingham plastic with the Augmented Lagrangian Method
Unsteady circular Couette flow of a Bingham plastic with the Augmented Lagrangian Method Rheologica Acta ISSN 0035-4511 Volume 49 Combined 11-12 Rheol Acta (2010) 49:1197-1206 DOI 10.1007/s00397-010-0497-
More information6.1 Steady, One-Dimensional Rectilinear Flows Steady, Spherically Symmetric Radial Flows 42
Contents 6 UNIDIRECTIONAL FLOWS 1 6.1 Steady, One-Dimensional Rectilinear Flows 6. Steady, Axisymmetric Rectilinear Flows 19 6.3 Steady, Axisymmetric Torsional Flows 8 6.4 Steady, Axisymmetric Radial Flows
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 informationELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY. David O. Olagunju. University of Delaware. Newark, DE 19716
ELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY David O. Olagunju Department of Mathematical Sciences University of Delaware Newark, DE 19716 June 15, 1995 Abstract Consider the axisymmetric,
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 informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationBoundary-Layer Transition. and. NASA Langley Research Center, Hampton, VA Abstract
Eect of Far-Field Boundary Conditions on Boundary-Layer Transition Fabio P. Bertolotti y Institut fur Stromungsmechanik, DLR, 37073 Gottingen, Germany and Ronald D. Joslin Fluid Mechanics and Acoustics
More informationFlow instabilities of Herschel Bulkley fluids
J. Non-Newtonian Fluid Mech. 116 (2003) 19 32 Flow instabilities of Herschel Bulkley fluids Andreas N. Alexandrou a,b,, Philippe Le Menn b, Georgios Georgiou c, Vladimir Entov d a Department of Mechanical
More informationThe time-dependent extrudate-swell problem of an Oldroyd-B fluid with slip along the wall
The time-dependent extrudate-swell problem of an Oldroyd-B fluid with slip along the wall Eric Brasseur Unité de Mécanique Appliquée, Université Catholique de Louvain, Bâtiment Euler, 4 6 Avenue Georges
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationAnumerical and analytical study of the free convection thermal boundary layer or wall jet at
REVERSED FLOW CALCULATIONS OF HIGH PRANDTL NUMBER THERMAL BOUNDARY LAYER SEPARATION 1 J. T. Ratnanather and P. G. Daniels Department of Mathematics City University London, England, EC1V HB ABSTRACT Anumerical
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
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 informationKostas D. Housiadas. Teaching experience: University of Patras: Simulations of transport phenomena, Spring 2005.
Kostas D. Housiadas Personal: Born: in Athens, Greece. Present position: Professor, Department of Mathematics, University of the Aegean, Karlovassi, Samos, Greece. Phone number: +30-22730-82152, E-mail:
More informationIntroduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:
Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir
More informationAn improved MPS method for numerical simulations of convective heat transfer problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2006; 51:31 47 Published online 15 November 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/d.1106
More informationTime-dependent plane Poiseuille flow of a Johnson±Segalman fluid
J. Non-Newtonian Fluid Mech. 82 (1999) 105±123 Abstract Time-dependent plane Poiseuille flow of a Johnson±Segalman fluid Marios M. Fyrillas a,b, Georgios C. Georgiou a,*, Dimitris Vlassopoulos b a Department
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 informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
More informationDrag reduction in turbulent MHD pipe ows. By P. Orlandi 1. This is a preliminary study devoted to verifying whether or not direct simulations
Center for Turbulence Research Proceedings of the Summer Program 1996 447 Drag reduction in turbulent MHD pipe ows By P. Orlandi 1 This is a preliminary study devoted to verifying whether or not direct
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationViscosity * Desmond Schipper Andrew R. Barron. 1 Introduction
OpenStax-CNX module: m50215 1 Viscosity * Desmond Schipper Andrew R. Barron This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 4.0 Abstract This module discusses
More informationFUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL
Annual Journal of Hydraulic Engineering, JSCE, Vol.54, 2010, February FUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL How Tion PUAY1 and Takashi HOSODA2 1 Member of JSCE, Phd Student,
More informationTime (s)
TO: Joseph M. Powers FROM: Christian D. Hanson DATE: 18 April 1997 RE: ME 334 Project (Part II) 1 Problem Description AFortune 500 hundred company, Gooey Foods, Inc., wants to develop a piping system which
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 informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
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 informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationON THE EFFECTIVENESS OF HEAT GENERATION/ABSORPTION ON HEAT TRANSFER IN A STAGNATION POINT FLOW OF A MICROPOLAR FLUID OVER A STRETCHING SURFACE
5 Kragujevac J. Sci. 3 (29) 5-9. UDC 532.5:536.24 ON THE EFFECTIVENESS OF HEAT GENERATION/ABSORPTION ON HEAT TRANSFER IN A STAGNATION POINT FLOW OF A MICROPOLAR FLUID OVER A STRETCHING SURFACE Hazem A.
More informationTHE INITIAL-VALUE PROBLEM FOR VISCOUS CHANNEL FLOWS 1. W.O. Criminale. University ofwashington. Seattle, Washington T.L.
THE INITIAL-VALUE PROBLEM FOR VISCOUS CHANNEL FLOWS 1 W.O. Criminale Department of Applied Mathematics University ofwashington Seattle, Washington 98195 T.L. Jackson Institute for Computer Applications
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationRegularization of the Chapman-Enskog Expansion and Its Description of Shock Structure
NASA/CR-2001-211268 ICASE Report No. 2001-39 Regularization of the Chapman-Enskog Expansion and Its Description of Shock Structure Kun Xu Hong Kong University, Kowloon, Hong Kong ICASE NASA Langley Research
More informationNumerical solution of general boundary layer problems by the method of dierential quadrature
Scientia Iranica B (2013) 20(4), 1278{1301 Sharif University of Technology Scientia Iranica Transactions B: Mechanical Engineering www.scientiairanica.com Research Note Numerical solution of general boundary
More informationApplied Mathematics and Computation
Applied Mathematics and Computation 291 (2016) 98 114 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc Development lengths in
More informationUnsteady Hydromagnetic Couette Flow within a Porous Channel
Tamkang Journal of Science and Engineering, Vol. 14, No. 1, pp. 7 14 (2011) 7 Unsteady Hydromagnetic Couette Flow within a Porous Channel G. S. Seth*, Md. S. Ansari and R. Nandkeolyar Department of Applied
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationCRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th
ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationStokes bubble with volumetric change. Qing Nie. Institute for Mathematics and Its Applications. 514 Vincent Hall. 206 Church Street S.E.
The evolution of an axi-symmetric Stokes bubble with volumetric change Qing Nie Institute for Mathematics and Its Applications 514 Vincent Hall 206 Church Street S.E. Minneapolis, MN 55455-0436 nie@ima.umn.edu
More informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
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 information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationEnsemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3
Center for Turbulence Research Proceedings of the Summer Program 1996 237 Ensemble averaged dynamic modeling By D. Carati 1,A.Wray 2 AND W. Cabot 3 The possibility of using the information from simultaneous
More informationResearch Article Innovation: International Journal of Applied Research; ISSN: (Volume-2, Issue-2) ISSN: (Volume-1, Issue-1)
Free Convective Dusty Visco-Elastic Fluid Flow Through a Porous Medium in Presence of Inclined Magnetic Field and Heat Source/ Sink 1 Debasish Dey, 2 Paban Dhar 1 Department of Mathematics, Dibrugarh University,
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 informationTHE ACOUSTIC AND INSTABILITY WAVES OF JETS CONFINED INSIDE AN ACOUSTICALLY LINED RECTANGULAR DUCT FANG Q. HU. Department of Mathematics and Statistics
THE ACOUSTIC AND INSTABILITY WAVES OF JETS CONFINED INSIDE AN ACOUSTICALLY LINED RECTANGULAR DUCT FANG Q. HU Department of Mathematics and Statistics Old Dominion University Norfolk, VA 23529 ABSTRACT
More informationFlow Transition in Plane Couette Flow
Flow Transition in Plane Couette Flow Hua-Shu Dou 1,, Boo Cheong Khoo, and Khoon Seng Yeo 1 Temasek Laboratories, National University of Singapore, Singapore 11960 Fluid Mechanics Division, Department
More informationERIK STERNER. Discretize the equations in space utilizing a nite volume or nite dierence scheme. 3. Integrate the corresponding time-dependent problem
Convergence Acceleration for the Navier{Stokes Equations Using Optimal Semicirculant Approximations Sverker Holmgren y, Henrik Branden y and Erik Sterner y Abstract. The iterative solution of systems of
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationCHAPTER 6 Effect of slip and heat transfer on the Peristaltic flow of a Williamson fluid in an incliped channel
CHAPTER 6 Effect of slip and heat transfer on the Peristaltic flow of a Williamson fluid in an incliped channel 6.1. Introduction Peristalsis is a well-known mechanism for pumping biological and industrial
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More informationPARALLEL PSEUDO-SPECTRAL SIMULATIONS OF NONLINEAR VISCOUS FINGERING IN MIS- Center for Parallel Computations, COPPE / Federal University of Rio de
PARALLEL PSEUDO-SPECTRAL SIMULATIONS OF NONLINEAR VISCOUS FINGERING IN MIS- CIBLE DISPLACEMENTS N. Mangiavacchi, A.L.G.A. Coutinho, N.F.F. Ebecken Center for Parallel Computations, COPPE / Federal University
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 informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
More informationmodel and its application to channel ow By K. B. Shah AND J. H. Ferziger
Center for Turbulence Research Annual Research Briefs 1995 73 A new non-eddy viscosity subgrid-scale model and its application to channel ow 1. Motivation and objectives By K. B. Shah AND J. H. Ferziger
More informationA j = 0.1 cm 2 10 cm 10 cm 10 cm. W j Wj. W j W j. W j. 10 cm 10 cm 10 cm. r i
ME 131B Fluid Mechanics Solutions to Week Eight Problem Session: Angular Momentum Principle (3/2/98) 1. In control volume analysis, all governing principles share the same common structure: storage = inow
More informationSchiestel s Derivation of the Epsilon Equation and Two Equation Modeling of Rotating Turbulence
NASA/CR-21-2116 ICASE Report No. 21-24 Schiestel s Derivation of the Epsilon Equation and Two Equation Modeling of Rotating Turbulence Robert Rubinstein NASA Langley Research Center, Hampton, Virginia
More informationBy B. Perot 1 AND P. Moin 2. A new approach to Reynolds averaged turbulence modeling is proposed which has
Center for Turbulence Research Proceedings of the Summer Program 1996 35 A new approach to turbulence modeling By B. Perot 1 AND P. Moin 2 A new approach to Reynolds averaged turbulence modeling is proposed
More informationExpansions-contractions Flows
III Expansions-contractions Flows III.1 Introduction Flows of viscoplastic materials through internal passages of abruptly varying cross sections are found in a wide variety of industrial and natural processes.
More informationIterative procedure for multidimesional Euler equations Abstracts A numerical iterative scheme is suggested to solve the Euler equations in two and th
Iterative procedure for multidimensional Euler equations W. Dreyer, M. Kunik, K. Sabelfeld, N. Simonov, and K. Wilmanski Weierstra Institute for Applied Analysis and Stochastics Mohrenstra e 39, 07 Berlin,
More informationThe dissertation of Bin Jiang is approved: Chair Date Date Date Date University of California at Santa Barbara 999
Non-overlapping Domain Decomposition and Heterogeneous Modeling Used in Solving Free Boundary Problems by Bin Jiang A dissertation submitted in partial satisfaction of the requirements for the degree of
More information1. (a) +EA; (b) EA; (c) 0; (d) 0 2. (a) 2; (b) 3; (c) 1 3. (a) equal; (b) equal; (c) equal e; (b) 150e 5. 3 and 4 tie, then 2, 1
CHAPTER 24 GAUSS LAW 659 CHAPTER 24 Answer to Checkpoint Questions 1. (a) +EA; (b) EA; (c) ; (d) 2. (a) 2; (b) 3; (c) 1 3. (a) eual; (b) eual; (c) eual 4. +5e; (b) 15e 5. 3 and 4 tie, then 2, 1 Answer
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationNatural Convection over a Non-Isothermal Vertical Flat Plate in Supercritical Fluids
Transaction B: Mechanical Engineering Vol. 6, No. 6, pp. 470{478 c Sharif University of Technology, December 2009 Natural Convection over a Non-Isothermal Vertical Flat Plate in Supercritical Fluids Abstract.
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More informationΓιάννης E. Δημακόπουλος
ΔΗΜΟΣΙΕΥΣΕΙΣ Γιάννης E. Δημακόπουλος ΑΡΘΡΑ ΣΕ ΔΙΕΘΝΗ ΠΕΡΙΟΔΙΚΑ Α1. Kouris, Ch., Dimakopoulos, J., Georgiou, G. and Tsamopoulos, J., Comparison of spectral and finite element methods applied to the study
More informationI. INTRODUCTION It is well known that, of the three shear ows most commonly used to model transition to turbulence, plane Poiseuille ow is linearly un
STABILITY ANALYSIS OF PERTURBED PLANE COUETTE FLOW Dwight Barkley Mathematics Institute, University of Warwick, Coventry, CV4 7AL, United Kingdom barkley@maths.warwick.ac.uk Laurette S. Tuckerman LIMSI-CNRS,
More informationA Numerical Study of Chamber Size and Boundary Eects on CPT Tip Resistance in NC Sand
Scientia Iranica, Vol. 15, No. 5, pp 541{553 c Sharif University of Technology, October 2008 A Numerical Study of Chamber Size and Boundary Eects on CPT Tip Resistance in NC Sand M.M. Ahmadi 1; and P.K.
More informationFLUID MECHANICS PROF. DR. METİN GÜNER COMPILER
FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 5. FLOW IN PIPES 5.1.3. Pressure and Shear Stress
More informationUNIVERSITY OF MARYLAND Department of Physics College Park, Maryland. PHYSICS Ph.D. QUALIFYING EXAMINATION PART II
UNIVERSITY OF MARYLAND Department of Physics College Park, Maryland PHYSICS Ph.D. QUALIFYING EXAMINATION PART II August 23, 208 9:00 a.m. :00 p.m. Do any four problems. Each problem is worth 25 points.
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 informationVorticity and Potential Vorticity
Chapter 4 Vorticity and Potential Vorticity In this chapter we explore a way of solving the shallow water equations for motions with characteristic time scales longer than the rotation period of the earth.
More informationModelling of dispersed, multicomponent, multiphase flows in resource industries Section 4: Non-Newtonian fluids and rheometry (PART 1)
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 4: Non-Newtonian fluids and rheometry (PART 1) Globex Julmester 2017 Lecture #3 05 July 2017 Agenda Lecture #3 Section
More informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationstrain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip
Nonlocal damage models: Practical aspects and open issues Milan Jirasek LSC-DGC, Swiss Federal Institute of Technology at Lausanne (EPFL), Switzerland Milan.Jirasek@ep.ch Abstract: The purpose of this
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 informationNUMERICAL SIMULATIONS OF CONSTANT VELOCITY SQUEEZE FLOW
U.P.B. Sci. Bull., Series D, Vol. 75, Iss. 2, 2013 ISSN 1454-2358 NUMERICAL SIMULATIONS OF CONSTANT VELOCITY SQUEEZE FLOW Daniela COBLAŞ 1, Diana BROBOANĂ 2, Corneliu BĂLAN 3, Mohamed HAJJAM 4 The paper
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationCandidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2011 2012 FLUID DYNAMICS MTH-3D41 Time allowed: 3 hours Attempt FIVE questions. Candidates must show on each answer book the type
More informationWall turbulence with arbitrary mean velocity profiles
Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the
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 informationStability of two-layer viscoelastic plane Couette flow past a deformable solid layer
J. Non-Newtonian Fluid Mech. 117 (2004) 163 182 Stability of two-layer viscoelastic plane Couette flow past a deformable solid layer V. Shankar Department of Chemical Engineering, Indian Institute of Technology,
More informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationDynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary Conditions
Transaction B: Mechanical Engineering Vol. 16, No. 3, pp. 273{279 c Sharif University of Technology, June 2009 Research Note Dynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary
More informationOn fully developed mixed convection with viscous dissipation in a vertical channel and its stability
ZAMM Z. Angew. Math. Mech. 96, No. 12, 1457 1466 (2016) / DOI 10.1002/zamm.201500266 On fully developed mixed convection with viscous dissipation in a vertical channel and its stability A. Barletta 1,
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 informationRenormalization Group analysis of 2D Ising model
Renormalization Group analysis of D Ising model Amir Bar January 7, 013 1 Introduction In this tutorial we will see explicitly how RG can be used to probe the phase diagram of d > 1 systems, focusing as
More informationInertial effect on stability of cone-and-plate flow Part 2: Non-axisymmetric modes
J. Non-Newtonian Fluid Mech., 78 (1998) 27 45 Inertial effect on stability of cone-and-plate flow Part 2: Non-axisymmetric modes Yuriko Renardy *, David O. Olagunju 1 Department of Mathematics and ICAM,
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 informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More information6 VORTICITY DYNAMICS 41
6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some
More informationSOE3211/2 Fluid Mechanics lecture 2
SOE311/ Fluid Mechanics lecture example example Fluid ows governed by conservation of mass, momentum. We can use this to solve ow problems. Draw box (control volume) around region of interest, then equate
More information