initial commit / add all books

1 files changed, 32 insertions, 0 deletions

diff --git a/3733/CH24/EX24.5/Ex24_5.sce b/3733/CH24/EX24.5/Ex24_5.sce
new file mode 100644
index 000000000..c836bf522
--- /dev/null
+++ b/3733/CH24/EX24.5/Ex24_5.sce
@@ -0,0 +1,32 @@
+// Example 24_5

+clc;funcprot(0);

+//Given data

+T_1=19+273;// K

+p_1=100;//kN/m^2

+p_2=800;// kN/m^2

+n_c=0.85;// The isentropic efficiency of compressor

+n_t=0.88;// The isentropic efficiency of turbine

+n_pt=0.86;// The isentropic efficiency of power turbine

+m=7;//Air flow rate in kg/s

+T_3=980+273;// K

+C_p=1.006;// kJ/kg.K

+r=1.4;// Specific heat ratio

+

+//Calculation

+T_2a=T_1*(p_2/p_1)^((r-1)/r);// K

+T_2=((T_2a-T_1)/n_c)+T_1;// K

+//(1)For the first turbine

+// Compressor work= Turbine work

+T_4=T_3-(T_2-T_1);// Turbine exit temperature in K

+T_4a=T_3-((T_3-T_4)/(n_t));// K

+p_3=p_2;// bar

+p_4a=(p_3)/((T_3/T_4a)^(r/(r-1)));// kN/m^2

+p_4=p_4a;//kN/m^2

+//(2)For the power turbine

+p_5=p_1;// bar

+T_5a=T_4*(p_5/p_4)^((r-1)/r);// K

+T_5=T_4-(n_pt*(T_4-T_5a));// K

+P=(m*C_p*(T_4-T_5));// kW

+n_th=(C_p*(T_4-T_5))/(C_p*(T_3-T_2));// Thermal efficiency

+printf('\n1.The condition of air at the exit of the first turbine:T_4=%0.0f K & p_4=%0.0f kN/m^2 \n2.The power output of the turbine=%0.0f kW\nThe thermal efficiency of the plant=%0.3f or %0.1f percentage',T_4,p_4,P,n_th,n_th*100 );

+// The answer vary due to round off error