blob: d76388239947fcbf80d5db9f560ae6905863279f (
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
33
34
35
36
37
38
|
// scilab Code Exa 3.2 Gas Turbine Plant with an exhaust HE
T1=300; // Minimum cycle Temperature in Kelvin
funcprot(0);
pr=10; // pressure ratio of the turbine and compressor
T3=1500; // Maximum cycle Temperature in Kelvin
m=10; // mass flow rate through the turbine and compressor in kg/s
e(1)=0.8; // thermal ratio of the heat exchanger
e(2)=1;
n_c=0.82; // Compressor Efficiency
n_t=0.85; // Turbine Efficiency
gamma=1.4; // Specific Heat Ratio
cp=1.005; // Specific Heat at Constant Pressure in kJ/(kgK)
beeta=T3/T1;
T2s=T1*(pr^((gamma-1)/gamma));
T2=T1+((T2s-T1)/n_c);
T4s=T3*(pr^(-((gamma-1)/gamma)));
T4=T3-((T3-T4s)*n_t);
for i=1:2
T5=T2+e(i)*(T4-T2);
T6=T4-(T5-T2);
Q_s=cp*(T3-T5);
Q_r=cp*(T6-T1);
// part(a) Determining power developed
w_p=Q_s-Q_r;
P=m*w_p;
printf("for effectiveness=%f, \n (a)the power developed is %f kW",e(i),P)
// part(b) Determining thermal efficiency of the plant
n_th=1-(Q_r/Q_s);
disp ("%",n_th*100,"(b)thermal efficiency of the plant is")
end
// part(c) Determining efficiencies of the ideal Joules cycle
n_Joule=1-(pr^((gamma-1)/gamma)/beeta);
disp("%",n_Joule*100,"(c)efficiency of the ideal Joules cycle with perfect heat exchange is")
n_Carnot=1-(T1/T3);
disp("%",n_Carnot*100,"and the Carnot cycle efficiency is")
|
