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 = ");