1 files changed, 25 insertions, 0 deletions
diff --git a/2223/CH3/EX3.1/Ex3_1.sce b/2223/CH3/EX3.1/Ex3_1.sce
new file mode 100755
index 000000000..adcbf9c38
--- /dev/null
+++ b/2223/CH3/EX3.1/Ex3_1.sce
@@ -0,0 +1,25 @@
+// scilab Code Exa 3.1 Constant Pressure Gas Turbine Plant
+
+t1=50;  // Minimum Temperature in degree C
+T1=t1+273; // in Kelvin
+t3=950;  // Maximum Temperature in degree C
+T3=t3+273; // in Kelvin
+n_c=0.82; // Compressor Efficiency
+n_t=0.87; // Turbine Efficiency
+gamma=1.4; // Specific Heat Ratio
+cp=1.005; // Specific Heat at Constant Pressure in kJ/(kgK)
+beeta=T3/T1;
+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 pressure ratio of the turbine and compressor
+pr=T_opt^(gamma/(gamma-1));
+disp(pr,"(a)Pressure Ratio is")
+
+// part(b) Determining maximum power output per unit flow rate
+wp_max=cp*T1*((T_opt-1)^2)/n_c;
+disp("kW/(kg/s)",wp_max,"(b)maximum power output per unit flow rate is")
+
+// part(c) 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,"(c)thermal efficiency of the plant for maximum power output is")