diff options
author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
---|---|---|
committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3511/CH6/EX6.11/Ex6_11.sce | |
parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
download | Scilab-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.sce | 34 |
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 = ");
|