1 files changed, 30 insertions, 0 deletions
diff --git a/2223/CH18/EX18.6/Ex18_6.sce b/2223/CH18/EX18.6/Ex18_6.sce
new file mode 100755
index 000000000..fae4d033b
--- /dev/null
+++ b/2223/CH18/EX18.6/Ex18_6.sce
@@ -0,0 +1,30 @@
+// scilab Code Exa 18.6 Calculations on a Gas Turbine 
+
+m=472; // flow rate of hot gases in kg/s
+T01=1335;  // Turbine inlet temp in Kelvin
+p01=10; // Turbine Inlet Pressure in bar
+c2=150; // exit velocity in m/s
+pr0=10; // Turbine pressure ratio
+gamma_g=1.67;
+T2=560;  // Temperature of gases at exit in Kelvin
+cp_g=1.157; // Specific Heat of gas at Constant Pressure in kJ/(kgK)
+
+// part(a) Determining total to total efficiency
+T02=T2+(0.5*(c2^2)/(cp_g*1e3));
+T02s=T01/(pr0^((gamma_g-1)/gamma_g));
+n_tt=(T01-T02)/(T01-T02s);
+disp("%",n_tt*100,"(a)total to total efficiency is")
+
+
+// part(b) Determining total to static efficiency
+T2s=T02s-(0.5*(c2^2)/(cp_g*1e3));
+n_ts=(T01-T02)/(T01-T2s);
+disp("%",n_ts*100,"(b)total to static efficiency is")
+
+// part(c) Determining the polytropic efficiency
+n_p=((gamma_g)/(gamma_g-1))*((log(T01/T02))/(log(pr0)));
+disp("%",n_p*100,"(c)polytropic efficiency is")
+
+// part(d) Determining power developed by the turbine
+P=m*cp_g*(T01-T02);
+disp("MW",P/1e3,"(d)Power developed by the turbine is")