summaryrefslogtreecommitdiff
path: root/2015/CH9/EX9.17/9_17.sce
blob: 097441e0e210e717a84fa6179b3087f488afeba8 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

clc
//initialisation of variables
p2=6 //prressure in bar
p1=1 //pressure in bar
t1=313 //temp in k
a1=45 //angle in degrees
a2=10 //angle in degrees
a3=55 //angle in degrees
r=1.4
cp=1.005 //kj/kg
ieff=0.85 //isentropic efficiency
c=200 //m/s
//CALCULATIONS
t2s=(t1*(p2/p1)^((r-1)/r))
t2=(((t2s-t1)/ieff)+t1)
w=cp*(t2-t1)
cro=(c*(sin(45*(%pi/180))/sin(55*(%pi/180))))
cv=c-cro
n=w/cv
//RESULTS
printf('actual work is %2fkj/kg',w)
printf('\nchange in whirl velocities is %2fkj/kg/stage',cv)
printf('\nnumber of stages is %2fstages',n)