summaryrefslogtreecommitdiff
path: root/3511/CH6/EX6.11/Ex6_11.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3511/CH6/EX6.11/Ex6_11.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3511/CH6/EX6.11/Ex6_11.sce')
-rw-r--r--3511/CH6/EX6.11/Ex6_11.sce34
1 files changed, 34 insertions, 0 deletions
diff --git a/3511/CH6/EX6.11/Ex6_11.sce b/3511/CH6/EX6.11/Ex6_11.sce
new file mode 100644
index 000000000..129cef6dd
--- /dev/null
+++ b/3511/CH6/EX6.11/Ex6_11.sce
@@ -0,0 +1,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 = ");