fluid_mechanics
312 Chapter 6 ■ Differential Analysis of Fluid Flow Lubricating oil ω r o Rotating shaft r i Housing F I G U R E 6.33 journal bearing. Flow in the narrow gap of a Flow between parallel plates with one plate fixed and the other moving is called Couette flow. The actual velocity profile will depend on the dimensionless parameter P b2 2mU a 0p 0x b Several profiles are shown in Fig. 6.32b. This type of flow is called Couette flow. The simplest type of Couette flow is one for which the pressure gradient is zero; that is, the fluid motion is caused by the fluid being dragged along by the moving boundary. In this case, with 0p0x 0, Eq. 6.140 simply reduces to u U y b (6.142) which indicates that the velocity varies linearly between the two plates as shown in Fig. 6.31b for P 0. This situation would be approximated by the flow between closely spaced concentric cylinders in which one cylinder is fixed and the other cylinder rotates with a constant angular velocity, v. As illustrated in Fig. 6.33, the flow in an unloaded journal bearing might be approximated by this simple Couette flow if the gap width is very small 1i.e., r o r i r i 2. In this case U r i v, b r o r i , and the shearing stress resisting the rotation of the shaft can be simply calculated as t mr i v1r o r i 2. When the bearing is loaded 1i.e., a force applied normal to the axis of rotation2, the shaft will no longer remain concentric with the housing and the flow cannot be treated as flow between parallel boundaries. Such problems are dealt with in lubrication theory 1see, for example, Ref. 92. E XAMPLE 6.9 Plane Couette Flow GIVEN A wide moving belt passes through a container of a viscous liquid. The belt moves vertically upward with a constant velocity, V 0 , as illustrated in Fig. E6.9a. Because of viscous forces the belt picks up a film of fluid of thickness h. Gravity tends to make the fluid drain down the belt. Assume that the flow is laminar, steady, and fully developed. V 0 y h Fluid layer FIND Use the Navier–Stokes equations to determine an expression for the average velocity of the fluid film as it is dragged up the belt. x g SOLUTION Since the flow is assumed to be fully developed, the only velocity component is in the y direction 1the v component2 so that u w 0. It follows from the continuity equation that 0v0y 0, and for steady flow 0v0t 0, so that v v1x2. Under these conditions the Navier–Stokes equations for the x direction 1Eq. 6.127a2 and the z direction 1perpendicular to the paper2 1Eq. 6.127c2 simply reduce to 0p 0x 0 0p 0z 0 F I G U R E E6.9a This result indicates that the pressure does not vary over a horizontal plane, and since the pressure on the surface of the film 1x h2 is atmospheric, the pressure throughout the film must be
6.9 Some Simple Solutions for Laminar, Viscous, Incompressible Fluids 313 atmospheric 1or zero gage pressure2. The equation of motion in the y direction 1Eq. 6.127b2 thus reduces to or Integration of Eq. 1 yields On the film surface 1x h2 we assume the shearing stress is zero—that is, the drag of the air on the film is negligible. The shearing stress at the free surface 1or any interior parallel surface2 is designated as , where from Eq. 6.125d Thus, if t xy 0 at x h, it follows from Eq. 2 that A second integration of Eq. 2 gives the velocity distribution in the film as At the belt 1x 02 the fluid velocity must match the belt velocity, V 0 , so that and the velocity distribution is therefore With the velocity distribution known we can determine the flowrate per unit width, q, from the relationship h h q v dx 0 t xy 0 rg m d 2 v dx 2 d 2 v dx 2 g m dv dx g m x c 1 t xy m a dv dx b c 1 gh m v g 2m x2 gh m x c 2 v g 2m x2 gh m x V 0 0 c 2 V 0 a g 2m x2 gh m x V 0b dx (1) (2) (3) and thus The average film velocity, V 1where q Vh2, is therefore COMMENT form as q V 0 h gh3 3m V V 0 gh2 3m (Ans) Equation (3) can be written in dimensionless v c a x 2 V 0 h b 2c a x h b 1 where c h 2 /2V 0 . This velocity profile is shown in Fig. E6.9b. Note that even though the belt is moving upward, for c 1 (e.g., for fluids with small enough viscosity or with a small enough belt speed) there are portions of the fluid that flow downward (as indicated by v/V 0 0). It is interesting to note from this result that there will be a net upward flow of liquid (positive V) only if V 0 h 2 /3. It takes a relatively large belt speed to lift a small viscosity fluid. v/V 0 1 0.8 0.6 0.4 0.2 c = 1.0 0 0 0.2 0.4 0.6 0.8 1 –0.2 –0.4 c = 1.5 –0.6 –0.8 –1 F I G U R E E6.9b x/h c = 0 c = 0.5 c = 2.0 An exact solution can be obtained for steady, incompressible, laminar flow in circular tubes. 6.9.3 Steady, Laminar Flow in Circular Tubes Probably the best known exact solution to the Navier–Stokes equations is for steady, incompressible, laminar flow through a straight circular tube of constant cross section. This type of flow is commonly called Hagen–Poiseuille flow, or simply Poiseuille flow. It is named in honor of J. L. Poiseuille 11799– 18692, a French physician, and G. H. L. Hagen 11797–18842, a German hydraulic engineer. Poiseuille was interested in blood flow through capillaries and deduced experimentally the resistance laws for laminar flow through circular tubes. Hagen’s investigation of flow in tubes was also experimental. It was actually after the work of Hagen and Poiseuille that the theoretical results presented in this section were determined, but their names are commonly associated with the solution of this problem. Consider the flow through a horizontal circular tube of radius R as is shown in Fig. 6.34a. Because of the cylindrical geometry it is convenient to use cylindrical coordinates. We assume that the flow is parallel to the walls so that v r 0 and v u 0, and from the continuity equation 16.342 0v z0z 0. Also, for steady, axisymmetric flow, is not a function of t or u so the velocity, v z , v z
- Page 286 and 287: 262 Chapter 5 ■ Finite Control Vo
- Page 288 and 289: 264 Chapter 6 ■ Differential Anal
- Page 290 and 291: ) 266 Chapter 6 ■ Differential An
- Page 292 and 293: 268 Chapter 6 ■ Differential Anal
- Page 294 and 295: 270 Chapter 6 ■ Differential Anal
- Page 296 and 297: 272 Chapter 6 ■ Differential Anal
- Page 298 and 299: 274 Chapter 6 ■ Differential Anal
- Page 300 and 301: 276 Chapter 6 ■ Differential Anal
- Page 302 and 303: 278 Chapter 6 ■ Differential Anal
- Page 304 and 305: 280 Chapter 6 ■ Differential Anal
- Page 306 and 307: 282 Chapter 6 ■ Differential Anal
- Page 308 and 309: 284 Chapter 6 ■ Differential Anal
- Page 310 and 311: 286 Chapter 6 ■ Differential Anal
- Page 312 and 313: 288 Chapter 6 ■ Differential Anal
- Page 314 and 315: 290 Chapter 6 ■ Differential Anal
- Page 316 and 317: 292 Chapter 6 ■ Differential Anal
- Page 318 and 319: 294 Chapter 6 ■ Differential Anal
- Page 320 and 321: 296 Chapter 6 ■ Differential Anal
- Page 322 and 323: 298 Chapter 6 ■ Differential Anal
- Page 324 and 325: 300 Chapter 6 ■ Differential Anal
- Page 326 and 327: 302 Chapter 6 ■ Differential Anal
- Page 328 and 329: 304 Chapter 6 ■ Differential Anal
- Page 330 and 331: 306 Chapter 6 ■ Differential Anal
- Page 332 and 333: 308 Chapter 6 ■ Differential Anal
- Page 334 and 335: 310 Chapter 6 ■ Differential Anal
- Page 338 and 339: 314 Chapter 6 ■ Differential Anal
- Page 340 and 341: 316 Chapter 6 ■ Differential Anal
- Page 342 and 343: 318 Chapter 6 ■ Differential Anal
- Page 344 and 345: 320 Chapter 6 ■ Differential Anal
- Page 346 and 347: 322 Chapter 6 ■ Differential Anal
- Page 348 and 349: 324 Chapter 6 ■ Differential Anal
- Page 350 and 351: 326 Chapter 6 ■ Differential Anal
- Page 352 and 353: 328 Chapter 6 ■ Differential Anal
- Page 354 and 355: 330 Chapter 6 ■ Differential Anal
- Page 356 and 357: 7 Dimensional Analysis, Similitude,
- Page 358 and 359: 334 Chapter 7 ■ Dimensional Analy
- Page 360 and 361: 336 Chapter 7 ■ Dimensional Analy
- Page 362 and 363: 338 Chapter 7 ■ Dimensional Analy
- Page 364 and 365: 340 Chapter 7 ■ Dimensional Analy
- Page 366 and 367: 342 Chapter 7 ■ Dimensional Analy
- Page 368 and 369: 344 Chapter 7 ■ Dimensional Analy
- Page 370 and 371: 346 Chapter 7 ■ Dimensional Analy
- Page 372 and 373: 348 Chapter 7 ■ Dimensional Analy
- Page 374 and 375: 350 Chapter 7 ■ Dimensional Analy
- Page 376 and 377: 352 Chapter 7 ■ Dimensional Analy
- Page 378 and 379: 354 Chapter 7 ■ Dimensional Analy
- Page 380 and 381: 356 Chapter 7 ■ Dimensional Analy
- Page 382 and 383: 358 Chapter 7 ■ Dimensional Analy
- Page 384 and 385: 360 Chapter 7 ■ Dimensional Analy
6.9 Some Simple Solutions for Laminar, Viscous, Incompressible Fluids 313<br />
atmospheric 1or zero gage pressure2. The equation of motion in<br />
the y direction 1Eq. 6.127b2 thus reduces to<br />
or<br />
Integration of Eq. 1 yields<br />
On the film surface 1x h2 we assume the shearing stress is<br />
zero—that is, the drag of the air on the film is negligible. The<br />
shearing stress at the free surface 1or any interior parallel surface2<br />
is designated as , where from Eq. 6.125d<br />
Thus, if t xy 0 at x h, it follows from Eq. 2 that<br />
A second integration of Eq. 2 gives the velocity distribution in<br />
the film as<br />
At the belt 1x 02 the <strong>fluid</strong> velocity must match the belt velocity,<br />
V 0 , so that<br />
and the velocity distribution is therefore<br />
With the velocity distribution known we can determine the<br />
flowrate per unit width, q, from the relationship<br />
h<br />
h<br />
q v dx <br />
0<br />
t xy<br />
0 rg m d 2 v<br />
dx 2<br />
d 2 v<br />
dx 2 g m<br />
dv<br />
dx g m x c 1<br />
t xy m a dv<br />
dx b<br />
c 1 gh<br />
m<br />
v g<br />
2m x2 gh<br />
m x c 2<br />
v g<br />
2m x2 gh<br />
m x V 0<br />
0<br />
c 2 V 0<br />
a g<br />
2m x2 gh<br />
m x V 0b dx<br />
(1)<br />
(2)<br />
(3)<br />
and thus<br />
The average film velocity, V 1where q Vh2, is therefore<br />
COMMENT<br />
form as<br />
q V 0 h gh3<br />
3m<br />
V V 0 gh2<br />
3m<br />
(Ans)<br />
Equation (3) can be written in dimensionless<br />
v<br />
c a x 2<br />
V 0 h b 2c a x h b 1<br />
where c h 2 /2V 0 . This velocity profile is shown in Fig. E6.9b.<br />
Note that even though the belt is moving upward, for c 1 (e.g.,<br />
for <strong>fluid</strong>s with small enough viscosity or with a small enough belt<br />
speed) there are portions of the <strong>fluid</strong> that flow downward (as indicated<br />
by v/V 0 0).<br />
It is interesting to note from this result that there will be a net<br />
upward flow of liquid (positive V) only if V 0 h 2 /3. It takes a<br />
relatively large belt speed to lift a small viscosity <strong>fluid</strong>.<br />
v/V 0<br />
1<br />
0.8<br />
0.6<br />
0.4<br />
0.2<br />
c = 1.0<br />
0<br />
0 0.2 0.4 0.6 0.8 1<br />
–0.2<br />
–0.4<br />
c = 1.5<br />
–0.6<br />
–0.8<br />
–1<br />
F I G U R E E6.9b<br />
x/h<br />
c = 0<br />
c = 0.5<br />
c = 2.0<br />
An exact solution<br />
can be obtained for<br />
steady, incompressible,<br />
laminar flow in<br />
circular tubes.<br />
6.9.3 Steady, Laminar Flow in Circular Tubes<br />
Probably the best known exact solution to the Navier–Stokes equations is for steady, incompressible,<br />
laminar flow through a straight circular tube of constant cross section. This type of flow is commonly<br />
called Hagen–Poiseuille flow, or simply Poiseuille flow. It is named in honor of J. L. Poiseuille 11799–<br />
18692, a French physician, and G. H. L. Hagen 11797–18842, a German hydraulic engineer. Poiseuille<br />
was interested in blood flow through capillaries and deduced experimentally the resistance laws<br />
for laminar flow through circular tubes. Hagen’s investigation of flow in tubes was also<br />
experimental. It was actually after the work of Hagen and Poiseuille that the theoretical results<br />
presented in this section were determined, but their names are commonly associated with the<br />
solution of this problem.<br />
Consider the flow through a horizontal circular tube of radius R as is shown in Fig. 6.34a.<br />
Because of the cylindrical geometry it is convenient to use cylindrical coordinates. We assume that<br />
the flow is parallel to the walls so that v r 0 and v u 0, and from the continuity equation 16.342<br />
0v z0z 0. Also, for steady, axisymmetric flow, is not a function of t or u so the velocity, v z ,<br />
v z