blob: 129cef6dd68dd9f46e5b74625d0eee62740ae1bf (
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
|
clc;
eff_pe=0.88; // Compressor and turbine polytropic efficiencies
T01=310; // Temperature at LP compressor inlet in kelvin
p01=14; // Pressure at LP compressor inlet in bar
rp=2; // Compressor pressure ratio
T03=300;// Temperature at HP compressor inlet in kelvin
m=180; // Mass flow of Helium in kg/s
Q=500; // Heat input to gas turbine in MW
T07=700; // Helium Temperature at entry to reactor channels in kelvin
P_precoller=0.34; // Pressure loss in pre-cooler and intercooler in bar
P_loss_HE=0.27; // Pressure loss in heat exchanger in bar
P_loss_RC=1.03; // Pressure loss in reactor channel in bar
eff_pc=0.88; // Polytropiic efficiency
Cp=5.19;// Specific heat at constant pressure in kJ/kg K
r=1.66; // Specific heat ratio
n_1_n=((r-1)/r)*(1/eff_pc);
T02=T01*rp^n_1_n;
T04=T03*rp^n_1_n;
T05=((Q*10^3)/(m*Cp))+T07;
T_press_loss=P_precoller+P_loss_HE+P_loss_RC; // Total pressure loss
p05=56-T_press_loss;
p06=p01+P_precoller+P_loss_HE;
n__1_n=eff_pc*((r-1)/r);
T06=T05/(p05/p06)^n__1_n;
WC=m*Cp*((T02-T01)+(T04-T03)); // Work of compressor
WT=m*Cp*(T05-T06); // Work of Turbine
WN=WT-WC; // Net work output
eff_th=WN/(Q*10^3); // Efficiency
eff=(T07-T04)/(T06-T04); // Effectiveness
disp ("MW (roundoff error)",WN/1000,"Power output = ");
disp ("% (roundoff error)",eff_th*100,"Thermal efficiency = ");
disp ("% (roundoff error)",eff*100,"Effectiveness = ");
|
