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)