1 files changed, 59 insertions, 0 deletions
diff --git a/3733/CH22/EX22.20/Ex22_20.sce b/3733/CH22/EX22.20/Ex22_20.sce
new file mode 100644
index 000000000..240a3a64c
--- /dev/null
+++ b/3733/CH22/EX22.20/Ex22_20.sce
@@ -0,0 +1,59 @@
+// Example 22_20
+clc;funcprot(0);
+//Given data
+T_1=556;// °C
+T_2=222;// °C
+m_s=20;// kg/sec
+n_m=80/100;// Mechanical efficiency
+n_t=95/100;// Transmission efficiency
+n_g=85/100;// Generator efficiency
+W_act=50/100;
+h_f1=76;// kJ/kg
+h_f2=29;// kJ/kg
+h_fg1=290;// kJ/kg
+h_fg2=302;// kJ/kg
+h_g1=366;// kJ/kg
+h_g2=331;// kJ/kg
+s_f1=0.152;// kJ/kg-K
+s_f2=0.08;// kJ/kg-K
+s_fg1=0.359;// kJ/kg-K
+s_fg2=0.626;// kJ/kg-K
+s_g1=0.511;// kJ/kg-K
+s_g2=0.706;// kJ/kg-K
+p_a=17;// bar
+p_b=0.035;// bar
+h_fa=874;// kJ/kg
+h_fb=111;// kJ/kg
+h_fga=1932;// kJ/kg
+h_fgb=2453;// kJ/kg
+h_ga=2806;// kJ/kg
+h_gb=2564;// kJ/kg
+s_fa=2.37;// kJ/kg-K
+s_fb=0.388;// kJ/kg-K
+s_ga=6.42;// kJ/kg-K
+s_gb=0.388;// kJ/kg-K
+
+// Calculation
+//(a)
+x_2=(s_g1-s_f2)/s_fg2;// The condition of the mercury vapour at point2
+m_hg=h_fga/(x_2*h_fg2);// kg
+//(b)
+W=h_g1-(h_f2+(x_2*h_fg2));// kJ/kg
+W_m=W*m_hg;//Work done per kg of Hg vapour in kJ
+//(c)
+// From steam tables,
+T_sup=380+273;// K
+T_sa=203.4+273;// K
+T_b=26.5+273;// K
+x_b=(((s_ga+(2*2.303*log(T_sup/T_sa)))*(T_b))-s_fb)/(T_b/h_fga);
+T_sup=383+273;// K
+x_2=0.72;
+W_s=(h_ga+(2*(T_sup-T_sa)))-(h_fb+(x_2*h_fgb));// Work done per kg of steam in  kJ/kg
+//(d)
+W=W_m+W_s;// Total work done in kJ
+Q_s=(m_hg*(h_g1-h_f2))+(1*(h_fa-h_fb))+(2*(T_sup-T_sa));// Heat supplied in kJ
+n_o=(W/Q_s)*100;// Overall efficiency of the cycle in %
+E=((m_s*W*W_act)*n_m*n_t*n_g)/1000;// Total energy generated per sec in MW
+printf('\n(a)Mass of Hg required per kg of steam used=%0.1f kg \n(b)Work done per kg of Hg vapour=%0.1f kJ/kg \n(c)Work done per kg of steam=%0.0f kJ/kg \n(d)Overall efficiency of the cycle=%0.1f percentage \n(e)Total energy generated per sec=%0.3f MW',m_hg,W_m,W_s,n_o,E);
+// The answer provided in the textbook is wrong
+