blob: 473dabb7232d8b11c59673810ad19847505ee13f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
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
|
