blob: 0d2488c092d34a7ed65403b68a2b6a33bf6bd0a9 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
clc;
rp=4; // Pressure ratio
WN=1500; // Net work output in kW
T01=25+273; // Inlet temperature in kelvin
p01=1; // Inlet pressure in bar
p03=4; // Turbine inlet pressure in bar
T03=700+273;// turbine inlet temperature in kelvin
eff_c=0.85; // Compressor efficiency
eff_over=0.21; // Overall efficiency
Cp=1.005;// Specific heat of air at constant pressure in kJ/kg K
r=1.4; // Specific heat ratio of air
T02=T01+T01*(rp^((r-1)/r)-1)/eff_c;
Q=WN/eff_over;
m=Q/(Cp*(T03-T02));
Wn=WN/m; // Net work per kg
T04=T03-T02+T01-(Wn/Cp);
T_04=T03/rp^((r-1)/r);
eff_T=(T03-T04)/(T03-T_04);
disp ("kg/s",m,"Mass flow rate = ");
disp ("%",eff_T*100,"Isentropic efficiency of the Turbine = ");
|