diff options
Diffstat (limited to '2223/CH18/EX18.10/Ex18_10.sce')
-rwxr-xr-x | 2223/CH18/EX18.10/Ex18_10.sce | 54 |
1 files changed, 54 insertions, 0 deletions
diff --git a/2223/CH18/EX18.10/Ex18_10.sce b/2223/CH18/EX18.10/Ex18_10.sce new file mode 100755 index 000000000..b3df9356c --- /dev/null +++ b/2223/CH18/EX18.10/Ex18_10.sce @@ -0,0 +1,54 @@ +// scilab Code Exa 18.10 Calculation on combined cycle power plant
+
+P_gt=25.845; // Power Output of gas turbine plant in MW
+P_st=21; // Power Output of steam turbine plant in MW
+m_gt=115; // mass flow rate of the exhaust gas in kg/s
+n_T=0.86; // Turbine Efficiency
+gamma_g=1.33;
+R=0.287;
+cp=(gamma_g/(gamma_g-1))*R; // Specific Heat at Constant Pressure in kJ/(kgK)
+T3=1341; // Maximum Temperature in gas turbine in degree K from Ex18.9
+p1=84; // steam Pressure at the entry of steam turbine in bar
+// from steam tables
+t_6s=298.4; // saturation temperature at 84 bar in degree C
+t_5s=t_6s;
+h_6s=1336.1; // from steam table liquid vapour enthalpy at 84 bar
+t6=535; // steam temperature at the entry of steam turbine in degree C
+T6=t6+273; // in Kelvin
+h_4s=3460; // from mollier diagram at t=535 degree C
+h_7=2050;
+p_c=0.07; // Condenser pressure in bar
+r=8.8502464; //optimum pressure ratio from Ex18.9
+T4=875.92974; //from Ex 18.9
+t4=T4-273; // in degree C
+h_7s=163.4; // Specific Enthalpy of water in kJ/kg
+m_st=P_st*1e3/((h_4s-h_7)*n_T); // mass flow rate of the steam in kg/s
+
+// part(a)Exhaust gas temperature at stack
+t_7=t4-((m_st*(h_4s-h_7s))/(m_gt*cp)); // energy balance for the economiser entry(7') to the superheater exit(4')
+disp("degree celsius",t_7,"(a)Exhaust gas temperature at stack is")
+
+// part(b)mass of steam per kg of gas
+disp("kg",m_st/m_gt,"(b)mass of steam per kg of gas is")
+
+// part(c) Pinch Point(PP)
+t_6=t_7+((m_st*(h_6s-h_7s))/(m_gt*cp)); // energy balance for the economiser
+PP=t_6-t_6s;
+disp("degree celsius",PP,"(c)Pinch Point(PP) is")
+
+// part(d)thermal efficiency of steam turbine plant
+delh4s_7ss=(h_4s-h_7)*n_T;
+n_st=delh4s_7ss/(h_4s-h_7s);
+disp("%",n_st*100,"(d)thermal Efficiency of steam turbine plant is")
+
+// part(e) thermal efficiency of the combined cycle plant
+n_B=0.978; // Assuming Combustion chamber Efficiency
+Qs=102.72554; // heat supplied in the combustion chamber from Ex 18.9
+Qss=Qs/n_B; // power supplied to the combined cycle
+n_gst=(P_gt+P_st)/Qss;
+disp ("%" ,n_gst*100,"(e)thermal Efficiency of combined gas and steam power plant is")
+
+// part(f)the dryness fraction of steam at the turbine exhaust
+x=0.875; // from Mollier diagram at p=0.07 bar
+disp(x,"(f)the dryness fraction of steam at the turbine exhaust is")
+
|