26.12.2014 Views

Chapter 9 Viscous flow

Chapter 9 Viscous flow

Chapter 9 Viscous flow

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

10 CHAPTER 9. VISCOUS FLOW<br />

x<br />

3<br />

x 1<br />

x 2<br />

g<br />

U<br />

Figure 9.3: Sphere settling in a fluid.<br />

condition (∂u i /∂x i ) = 0 can be simplified as<br />

∂u i<br />

∂x i<br />

(<br />

∂ 1 ∂<br />

= A 1 U i<br />

∂x i r + A δ ij<br />

2U j<br />

∂x i r − 3 ∂ )<br />

x i x j<br />

+ A 3U j<br />

3 ∂x i r 5 2µ<br />

= − A 1U i x i<br />

r 3 + A 2 U j<br />

(<br />

+ A 3U j<br />

2µ<br />

= − A 1U i x i<br />

r 3<br />

(<br />

δii x j<br />

r 3<br />

+ δ ijx i<br />

r 3<br />

− 3δ ijx i<br />

r 5<br />

− 3δ ijx i<br />

r 5<br />

− 3x2 i x )<br />

j<br />

r 5<br />

− 3δ iix j<br />

r 5<br />

∂ x i x j<br />

∂x i r 3<br />

)<br />

+ 15x2 ix j<br />

r 7<br />

+ A 3U j x j<br />

2µr 3 (9.30)<br />

In simplifying equation 9.30, we have used x 2 i = (x 2 1+x 2 2+x 2 3) = r 2 . Therefore,<br />

from the incompressibility condition, we get A 1 = (A 3 /2µ). Inserting this<br />

into equation 9.30, the velocity field is<br />

(<br />

δij<br />

u i = A 1 U j<br />

r + x ) (<br />

ix j δij<br />

+ A<br />

r 3 2 U j<br />

r − 3x )<br />

ix j<br />

(9.31)<br />

3 r 5<br />

The constants A 1 and A 2 are determined from the boundary condition at the<br />

surface of the sphere (r = 1),<br />

u i | r=1<br />

= A 1 U i + A 1 U j x i x j + A 2 U i − 3A 2 U j x i x j<br />

= U i (9.32)<br />

From this, we get two equations for A 1 and A 2 ,<br />

A 1 + A 2 = 1<br />

A 1 − 3A 2 = 0 (9.33)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!