1 files changed, 50 insertions, 0 deletions
diff --git a/2223/CH18/EX18.9/Ex18_9.sce b/2223/CH18/EX18.9/Ex18_9.sce
new file mode 100755
index 000000000..4e07f3e37
--- /dev/null
+++ b/2223/CH18/EX18.9/Ex18_9.sce
@@ -0,0 +1,50 @@
+// scilab Code Exa 18.9 Constant Pressure Gas Turbine Plant
+
+T1=298; // Minimum Temperature in Kelvin
+beeta=4.5; // Maximum to Minimum Temperature ratio(T_max/T_min)
+m=115; // mass flow rate through the turbine and compressor in kg/s
+n_C=0.79; // Compressor Efficiency
+n_T=0.83; // Turbine Efficiency
+gamma_g=1.33;
+R=0.287;
+cp=(gamma_g/(gamma_g-1))*R; // Specific Heat at Constant Pressure in kJ/(kgK)
+alpha=beeta*n_C*n_T;
+t_opt=sqrt(alpha); // For maximum power output, the temperature ratios in the turbine and compressor
+
+// part(a) Determining optimum pressure ratio of the plant
+r=t_opt^(gamma_g/(gamma_g-1));
+disp(r,"(a)optimum pressure ratio of the plant is")
+
+// part(b)Carnot's efficiency
+n_Carnot=1-(1/beeta);
+disp("%",n_Carnot*100,"(b)Carnot efficiency of the plant is")
+
+// part(c) Determining Joule's cycle efficiency
+n_Joule=1-(1/t_opt);
+disp("%",n_Joule*100,"(c)efficiency of the Joule cycle is")
+
+// part(d) Determining thermal efficiency of the plant for maximum power output
+n_th=(t_opt-1)^2/((beeta-1)*n_C-(t_opt-1));
+disp("%",n_th*100,"(d)thermal efficiency of the plant for maximum power output is")
+
+// part(e) Determining power output
+wp_max=cp*T1*((t_opt-1)^2)/n_C; // maximum work output
+P_max=m*wp_max;
+disp ("MW",P_max/1e3,"(e)Power output is")
+
+// part(f) Determining power generated by the turbine required to drive the compressor
+T3=beeta*T1;  // Maximum Temperature in degree K
+T4s=T3*(r^(-((gamma_g-1)/gamma_g)));
+T4=T3-((T3-T4s)*n_T);
+P_T=m*cp*(T3-T4);
+disp ("MW",P_T/1e3,"(f)Power generated by the turbine is")
+
+// part(g) Determining power absorbed by the compressor
+T2s=T1*(r^((gamma_g-1)/gamma_g));
+T2=T1+((T2s-T1)/n_C);
+P_C=m*cp*(T2-T1);
+disp ("MW",P_C/1e3,"(g)Power absorbed by the compressor is")
+
+//part(h)heat supplied in the combustion chamber
+Qs=m*cp*(T3-T2);
+disp("MW",Qs/1e3,"(h)heat supplied in the combustion chamber is")