Problem Set 5 Solution
Problem Set 5 Solution
We seek the position of the junction point J(x, y), or the values of the branching angles
1
,
and
2
, for which H is a minimum. Differentiating the above equation first with respect to
x and then with respect to y we get:
2
1
2
2
2
2
2 2
2
1
2
1
2
1
1 1
0
} ) ( ) {(
) (
} ) ( ) {(
) (
y y x x
x x h
y y x x
x x h
h
x
H
+
+
+
+ =
c
c
and
2
1
2
2
2
2
2 2
2
1
2
1
2
1
1 1
} ) ( ) {(
) (
} ) ( ) {(
) (
y y x x
y y h
y y x x
y y h
y
H
+
+
+
=
c
c
which can be written as:
0 cos cos
2 2 1 1 0
= + + =
c
c
u u h h h
x
H
and 0 sin sin
2 2 1 1
= + =
c
c
u u h h
y
H
Therefore:
1 0
2
2
2
1
2
0
1
2
cos
h h
h h h +
= u and
2 0
2
1
2
2
2
0
2
2
cos
h h
h h h +
= u
In a fully developed Poiseuille flow through a vessel segment of radius r and length l, if q
denotes the flow and denotes the viscosity of blood, the pumping power H, can be
written as:
H=h.l
Where h=8q
2
/r
4
.
For a given branching mode, therefore, the optimum branching angles for minimum
power, can be obtained simply by substituting for the general property h in the
expressions obtained in the previous section. Thus, we obtain the following result:
4
1
2
1
4
0
2
0
8
2
4
2
8
1
4
1
8
0
4
0
1
2
cos
r
q
r
q
r
q
r
q
r
q
+
= u and
4
2
2
2
4
0
2
0
8
1
4
1
8
2
4
2
8
0
4
0
2
2
cos
r
q
r
q
r
q
r
q
r
q
+
= u
Considering the fact that
3
r q for a circular cross-sections when the power is
minimized (cube law), and substituting
2
1
2
2
r
r
= o :
3
2
3
4
3
4
3
1
) 1 ( 2
1 ) 1 (
o
o o
u
+
+ +
= Cos and
3
2
3 2
4
3
4
3
2
) 1 ( 2
1 ) 1 (
o o
o o
u
+
+ +
= Cos
2) Consider pulsatile flow in a rigid tube with frequency =3.0 and driving pressure
gradient of the form KsCost. Using the equation:
t i
s
e
J
J
u
t r u
e
|
|
.
|
\
|
=
) (
) (
1
4
) , (
0
0
2
and value of Bessel functions (for =3.0, =-2.1213+2.1213i and J
0
()=-
0.2214+1.9376i), find the velocity at the center of the tube (a) at the beginning of the
cycle and (b) At a quarter way through the cycle (50 points).
At the center of the tube (r=0 ), =0. Knowing that J
0
(0)=1, we then have after
simplifying:
) sin )(cos 4703 . 0 2264 . 0 ( ) 4703 . 0 2264 . 0 (
) , 0 (
t i t i e i
u
t u
t i
s
e e
e
+ = =
The velocity corresponding to the given pressure gradient in given by the real part of the
above expression, since K
s
cost is the real part of the pressure gradient. Therefore:
a. At the beginning of the cycle where t=0, the required velocity is 0.2264
b. At a quarter away through the cycle where t=90 degree, the required velocity is
0.4703