blob: 42407227c560bd90a52663ade9848267c890a3b8 (
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
25
26
27
28
29
30
31
32
|
clc
//initialization of new variables
clear
r=0.5 //m average radius
Cz=140 //m/s Axial velocity
w=4000 //RPM turn rate
T01=300 //K Stagnation temperature ahead of rotor
alpha1=20 //degrees Incoming velocity angle
Cp=0.45 //Pressure rise coefficient
eta=0.98 //efficiency
cp=0.24 //specific heat
gama=1.4
//calculations
temp=alpha1 //just to store it
cp=cp*4200
U=r*w*2*%pi/60
alpha1=alpha1*%pi/180
beta1=atan(U/Cz-tan(alpha1))
beta2=acos(sqrt(cos(beta1)^2/(1-Cp)))
DCth=U-Cz*(tan(alpha1)+tan(beta2))
Pr=(1+eta*U*DCth/(cp*T01))^(gama/(gama-1))
printf('part (1)')
printf('\n stagnation pressure rise (ratio)is %.2f',Pr)
// part (2)
Cz=100 //m/s
alpha1=atan(U/Cz-tan(beta1))
Dalpha=alpha1*180/%pi-temp
printf('\n Stator ahead of this stage must be rotated by %.1f degrees',Dalpha)
DCth=U-Cz*(tan(alpha1)+tan(beta2))
Pr=(1+eta*U*DCth/(cp*T01))^(gama/(gama-1))
printf('\n part (2)')
printf('\n stagnation pressure rise (ratio)is %.2f',Pr)
|
