1 files changed, 32 insertions, 0 deletions
diff --git a/555/CH5/EX5.5/5.sce b/555/CH5/EX5.5/5.sce
new file mode 100644
index 000000000..dd75dbfc6
--- /dev/null
+++ b/555/CH5/EX5.5/5.sce
@@ -0,0 +1,32 @@
+// Implementation of example 5.5
+// Basic and Applied Thermodynamics by P.K.Nag
+
+clc
+clear
+
+//initial velocity 'V1', initial temperature of air 't1',Temperature after passing through heat exchanger 't2', Temperature after expansion 't3',Velocity after leaving turbine 'V3', Temperature after leaving turbine 't4', air flow rate 'w', enthalpy of air 'h', specific heat 'Cp' 
+t1 = 15; //degree C
+t2 = 800; //degree C
+V1 = 30; //m/s
+V2 = 30; //m/s
+t3 = 650; //degree C
+V3 = 60; //m/s
+t4 = 500; //degree C
+g = 9.8; //m/s2
+w = 2; //kg/s 
+cp = 1.005; //kJ/kg K 
+h2 = cp*t2;
+h1 = cp*t1;
+h3 = cp*t3;
+h4 = cp*t4;
+//Rate of heat transfer 'Q12'
+Q12 = w*(h2-h1); //kJ/s
+mprintf("Rate of heat transfer to air in heat exchanger, Q12 = %d kJ/s\n\n",round(Q12));
+//V = V2^2 - V3^2
+//Power output from turbine 'Wt'
+Wt = w*((V2^2-V3^2)*10^(-3)/2 + (h2-h3)); //kW
+mprintf("Power output from turbine, Wt = %0.1f kW\n\n",Wt);
+//Velocity at exit
+V4 = sqrt(2*(h3-h4)*1000 + V3^2);
+mprintf("Velocity at exit from the nozzle, V4 = %d m/s",round(V4));
+// end  
+\ No newline at end of file