1 files changed, 39 insertions, 0 deletions
diff --git a/3685/CH12/EX12.12/Ex12_12.sce b/3685/CH12/EX12.12/Ex12_12.sce
new file mode 100644
index 000000000..473dabb72
--- /dev/null
+++ b/3685/CH12/EX12.12/Ex12_12.sce
@@ -0,0 +1,39 @@
+clc
+// From table 
+h1 = 2792.2 // Enthalpy at state 1 in kJ/kg 
+h4 = 122.96// Enthalpy at state 4 in kJ/kg 
+hb = 254.88 // Enthalpy at state b in kJ/kg 
+hc = 29.98// Enthalpy at state c in kJ/kg 
+ha = 355.98 // Enthalpy at state a in kJ/kg 
+hd = hc // Isenthalpic process
+h2 = 1949.27 // // Enthalpy at state 2 in kJ/kg 
+//
+m = (h1-h4)/(hb-hc) // Amount of mercury circulating
+Q1t = m*(ha-hd) // Heat addition
+W1t = m*(ha-hb) + (h1-h2) // Turbine work
+n = W1t/Q1t // first law efficiency
+
+printf("\n Example 12.12 \n")
+printf("\n Overall efficiency of the cycle is %f percent",n*100)
+//The answers vary due to round off error
+
+S = 50000 // Stem flow rate through turbine in kg/h
+wm = S*m // mercury flow rate
+printf("\n Flow through the mercury turbine is %e kg/h",wm)
+
+Wt = W1t*S/3600 // Turbine work
+printf("\n Useful work done in binary vapor cycle is %f MW",Wt/1e3)
+nm = 0.85 // Internal efficiency of mercury turbine
+ns = 0.87 // Internal efficiency of steam turbine
+WTm = nm*(ha-hb) // turbine work of mercury based cycle
+hb_ = ha-WTm // Enthalpy at state b in kJ/kg
+m_ = (h1-h4)/(hb_-hc) // mass flow rate of mercury
+h1_ = 3037.3 // Enthalpy at state 1 in kJ/kg
+Q1t = m_*(ha-hd)+(h1_-h1) // Heat addition
+x2_ = (6.9160-0.4226)/(8.47-0.4226) // steam quality
+h2_ = 121+(0.806*2432.9) // Enthalpy at state 2 in kJ/kg 
+WTst = ns*(h1_-h2_) // Turbine work
+WTt = m_*(ha-hb_)+WTst // Total turbine work
+N = WTt/Q1t //Overall efficiency 
+printf("\n Overall efficiency is %f percent",N*100)
+// The answers vary due to round off error