1 files changed, 55 insertions, 0 deletions
diff --git a/3685/CH12/EX12.6/Ex12_6.sce b/3685/CH12/EX12.6/Ex12_6.sce
new file mode 100644
index 000000000..4aaa681d3
--- /dev/null
+++ b/3685/CH12/EX12.6/Ex12_6.sce
@@ -0,0 +1,55 @@
+clc
+h1 = 3023.5  // Enthalpy of steam at state 1 in kJ/kg
+s1 = 6.7664 // Enthalpy of steam at state 1 in kJ/kgK
+s2 = s1 // Isentropic process
+s3 = s1 //Isentropic process
+s4 = s1 //Isentropic process
+t_sat_20 = 212 // Saturation temperature at 20 bar in degree Celsius
+t_sat_1 = 46 // Saturation temperature at 1 bar in degree Celsius
+dt = t_sat_20-t_sat_1 // Change in temperature
+n =3 // number of heaters
+t = dt/n // temperature rise per heater
+t1 = t_sat_20-t // Operational temperature of first heater
+t2 = t1-t// Operational temperature of second heater
+// 0.1 bar
+hf = 191.83 // Enthalpy of fluid in kJ/kg
+hfg = 2392.8 // Latent heat of vaporization in kJ/kg
+sf = 0.6493// Entropy of fluid in kJ/kgK
+sg = 8.1502// Entropy of gas in kJ/kgK
+// At 100 degree
+hf100 = 419.04 // Enthalpy of fluid in kJ/kg  
+hfg100 = 2257.0// Latent heat of vaporization in kJ/kg 
+sf100 = 1.3069 // Entropy of fluid in kJ/kgK 
+sg100 = 7.3549 // Entropy of gas in kJ/kgK
+// At 150 degree
+hf150 = 632.20 // Enthalpy of fluid in kJ/kg  
+hfg150 = 2114.3// Latent heat of vaporization in kJ/kg  
+sf150 = 1.8418 // Entropy of fluid in kJ/kgK 
+sg150 = 6.8379// Entropy of gas in kJ/kgK
+x2 = (s1-sf150)/4.9961 // Steam quality
+h2 = hf150+(x2*hfg150) // Enthalpy at state 2 in kJ/kg
+x3 = (s1-sf100)/6.0480 // Steam quality
+h3 = hf100+(x3*hfg100) // Enthalpy at state 3 in kJ/kg 
+x4 = (s1-sf)/7.5010 // Steam quality
+h4 = hf+(x4*hfg)//Enthalpy at state 4 in kJ/kg
+h5 = hf // Enthalpy at state 5 in kJ/kg
+h6 = h5 //Enthalpy at state 6 in kJ/kg
+h7 = hf100 // Enthalpy at state 7 in kJ/kg
+h8 = h7 // Enthalpy at state 8 in kJ/kg
+h9 = 632.2  // Enthalpy at state 9 in kJ/kg
+h10 = h9 // Enthalpy at state 10 in kJ/kg
+m1 = (h9-h7)/(h2-h7) // regenerative mass 
+m2 = ((1-m1)*(h7-h6))/(h3-h6) // regenerative mass
+Wt = 1*(h1-h2)+(1-m1)*(h2-h3)+(1-m1-m2)*(h3-h4) // Turbine work
+Q1 = h1-h9 // Heat addition
+Wp = 0  // Pump work is neglected
+n_cycle = 100*(Wt-Wp)/Q1 // Cycle efficiency
+sr = 3600/(Wt-Wp)  // Steam rate
+
+printf("\n Example 12.6\n")
+printf("\n Steam quality at turbine exhaust is %f ",x3)
+printf("\n Net work per kg of stem is %f kJ/kg",Wt)
+printf("\n Cycle efficiency is %f percent",n_cycle)
+printf("\n Stream rate is %f kg/kW h",sr)
+//The answers vary due to round off error
+