blob: e4224037184ee2482f3e467a0c111e2dc03bfbc9 (
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
|
// scilab Code Exa 18.22 Isentropic or Stage Terminal Velocity for Turbines
T01=1273; // in Kelvin
funcprot(0);
p01=5; // Initial Pressure in bar
p02=3.5; // exit gas Pressure in bar
cp=1.005; // Specific Heat at Constant Pressure in kJ/(kgK)
gamma=1.4; // Specific Heat Ratio
m=28; // mass flow rate of the gas in kg/s
n_tt=0.84; // stage efficiency
shi=1.7; // stage loading coefficient
pr_0=p01/p02;
delh01_03ss=cp*T01*(1-(pr_0^((1-gamma)/gamma)));
//part(a)stage terminal velocity
c0=sqrt(2*delh01_03ss*1e3);
disp("m/s",c0,"(a)stage terminal velocity is")
// part(b)isentropic blade to gas speed ratio
sigma_s=sqrt(0.5*n_tt/shi);
disp(sigma_s,"(b)the isentropic blade to gas speed ratio is")
//part(c) peripheral speed of the rotor
u=sigma_s*c0;
disp("m/s",u,"(c)peripheral speed of the rotor is")
//part(d) the power developed
P=m*n_tt*delh01_03ss;
disp("MW",P*1e-3,"(d) the power developed is")
|