1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
|
// determine flow rate,theoretical head required power,pressure across the impeller
clc
b1=44
r1=21e-3 // mm
B=11e-3 // mm
r2=66e-3 //mm
b2=5e-3//mm
N=2500 //rpm
h1=0
g=9.8
alpha=90 // degree
D1=2
D2=2
u1=(2*%pi*N*r1)/60
gamm1=9800
p1=1000
mprintf('\n peripherial velocity at inlet u1=wR1 =%f m/s',u1)
u2=(2*%pi*N*r2)/60
mprintf('\n peripherial velocity at exit u2= wR2=%f m/s',u2)
V1=tand(b1)*u1
mprintf('\n V1f = %f m/s',V1)
Q=%pi*2*r1*B*V1
mprintf('\n Q = %f m3/s',Q)
V2f=Q/(2*%pi*r2*b2)
mprintf('\n V2f =%f m/s',V2f)
V2w=V2f/(tand(30))
mprintf('\n u2-V2w = %f ',V2w)
v2w=u2-V2w
mprintf('\n V2w = %f m/s',v2w)
alpha2=atand(V2f/v2w)
mprintf('\n alpha2 = %f degree',alpha2)
v2=v2w/cosd(18.9)
mprintf('\n V2= %f m/s',v2)
H1=(u2*v2w)/g
mprintf('\n H1 = %f m',H1)
p=gamm1*Q*H1
mprintf('\n H1 = %f watt',p)
P2=(p1*g*H1)-((p1/2)*(v2^2-V1^2))
mprintf('\n p2-p1 = %e Pa',P2)
mprintf('\n p2-p1 = %f bar',P2/10^5)
|
