1 files changed, 49 insertions, 0 deletions
diff --git a/3733/CH24/EX24.36/Ex24_36.sce b/3733/CH24/EX24.36/Ex24_36.sce
new file mode 100644
index 000000000..884254949
--- /dev/null
+++ b/3733/CH24/EX24.36/Ex24_36.sce
@@ -0,0 +1,49 @@
+// Example 24_36
+clc;funcprot(0);
+//Given data 
+P=25;// MW
+p_r=8;// Pressure ratio
+T_1=300;// K
+p_1=1;// bar
+T_8=700;// K
+AF_1=80;// Air fuel ratio
+AF_2=70;// Air fuel ratio
+e=0.7;// Effectiveness of heat exchanger 
+CV=40*10^3;// kJ/kg
+C_pa=1;// kJ/kg.K
+C_pg=1;// kJ/kg.K
+r=1.4;// Specific heat ratio for air and gases
+
+//Calculation
+p_2=p_1*p_r;// bar
+T_2=T_1*(p_r)^((r-1)/r);// K
+function[Y]=temperature(x)
+    Y(1)=(e*(x(1)-T_2))-(x(1)-T_8);
+endfunction
+x=[100];
+T=fsolve(x,temperature);
+T_7=T(1);// K
+T_2a=(T_7-T_8)+T_2;// K
+//Assume m_f1=y(1);T_3=y(2);m_f2=y(3);T_5=y(4);T_6=y(5);T_4=y(6)
+function[X]=massoffuel(y)
+    X(1)=(((80*y(1))+y(1))*C_pg*(y(2)-T_2a))-(y(1)*CV);
+    X(2)=(((70*y(3))+y(3))*C_pg*(y(4)-T_2a))-(y(3)*CV);
+    X(3)=(((70*y(3))+y(3))*C_pg*(y(4)-y(5)))-(P*10^3);
+    X(4)=y(5)-((y(4))/((p_r)^((r-1)/r)));
+    X(5)=(((80*y(1))+y(1))*((y(2))-y(6)))-(((80*y(1))+(70*y(3)))*(T_2-T_1));
+    X(6)=y(6)-((y(2))/((p_r)^((r-1)/r)));
+endfunction
+y=[0.1 1000 0.1 1000 100 100 ];
+z=fsolve(y,massoffuel);
+m_f1=z(1);// kg/sec
+m_f2=z(3);// kg/sec
+T_3=z(2);// K
+T_4=z(6);// K
+T_5=z(4);// K
+T_6=z(5);// K
+m_f=(m_f1+m_f2)*3600;// Total mass of fuel consumed per hour in kg/hr
+m_a=((m_f1*AF_1)+(m_f2*AF_2))*60;// Mass of air compressed per minute in kg/hr
+Sfc=(m_f)/(P*10^3);// Specific fuel consumption in kg/kW-hr
+n_th=((P*10^3)/((m_f1+m_f2)*CV))*100;// Thermal efficiency in %
+printf('\n(a)Total mass of fuel consumed per hour=%0.0f kg/hr \n(b)Mass of air compressed per minute=%0.0f kg/hr \n(c)Specific fuel consumption=%0.3f kg/kW-hr \n(d)Thermal efficiency=%0.0f percentage',m_f,m_a,Sfc,n_th);
+// The answers provided in the textbook is wrong