1 files changed, 46 insertions, 0 deletions
diff --git a/3733/CH24/EX24.17/Ex24_17.sce b/3733/CH24/EX24.17/Ex24_17.sce
new file mode 100644
index 000000000..488919804
--- /dev/null
+++ b/3733/CH24/EX24.17/Ex24_17.sce
@@ -0,0 +1,46 @@
+// Example 24_17
+clc;funcprot(0);
+//Given data 
+P=1600;// Power plant capacity in kW
+T_1=300;// K
+T_3=1050;// K
+p_1=1;// bar
+T_5=1100;// K
+p_2=5;// bar
+e=0.7;// The effectiveness of regenerator
+n_c=0.8;// Isentropic efficiency of compressor
+n_t1=0.85;// Efficiency of compressor turbine
+n_t2=0.90;// Efficiency of power turbine
+n_com=0.95;// Combustion efficiency
+n_m=0.90;// Mechanical efficiency of both turbines
+n_g=1;// Generation efficiency
+CV=40000;// kJ/kg
+C_pa=1;// kJ/kg.K
+C_pg=1.1;// kJ/kg.K
+r_a=1.4;// Specific heat ratio for air
+r_g=1.35;// Specific heat ratio for gases
+
+//Calculation
+T_2a=T_1*(p_2/p_1)^((r_a-1)/r_a);// K
+T_2=((T_2a-T_1)/n_c)+T_1;// K
+T_4a=T_3*(p_1/p_2)^((r_g-1)/r_g);// K
+T_4=T_3-(n_t1*(T_3-T_4a));// K
+T_6a=T_5*(p_1/p_2)^((r_g-1)/r_g);// K
+T_6=T_5-(n_t2*(T_5-T_6a));// K
+m_a2=P/(C_pg*(T_5-T_6)*n_m*n_g);// kg/sec
+//Power developed by compressor turbine = Power absorbed by compressor
+//m_a1=y(1)
+function[X]=mass(y);
+    X(1)=((y(1)+m_a2)*C_pa*(T_2-T_1))-(y(1)*C_pg*(T_3-T_4)*n_m);
+endfunction
+y=[10];
+z=fsolve(y,mass);
+m_a1=z(1);// kg/sec
+T_y=((m_a1/(m_a1+m_a2))*T_4)+((m_a2/(m_a1+m_a2))*T_6);// The temperature after mixing in °C
+T_x=T_2+((e*C_pg*(T_y-T_2))/C_pa);// K
+m_f=((C_pg*m_a1*(T_3-T_x))+(C_pg*m_a2*(T_5-T_x)))/(CV*n_com);// kg/sec
+n_th=(P/(m_f*CV))*100;// Plant efficiency in percentage
+Sfc=(m_f*3600)/P;// kg/kWh
+Afr=(m_a1+m_a2)/m_f;// Air fuel ratio
+printf('\n(a)Plant efficiency=%0.1f percentage \n(b)Specific fuel consumption=%0.3f kg/kW-hr \n(c)Air fuel ratio=%0.0f:1',n_th,Sfc,Afr);
+// The answer vary due to round off error