Chapter 5 : Properties of Liquids and Gases

Example 5.1 Page No : 182

In [1]:
# given data
p = 0.6988; 			#Unit:psia 			#absolute pressure
vg = 467.7; 			#Unit:ft**3/lbm 			#Saturated vapour specific volume
ug = 1040.2; 			#Unit:Btu/lbm 			#Saturated vapour internal energy
J = 778; 			#J = Conversion factor
# 1 Btu  =  778 ft*LBf
#h = u+(p*v)/J 
hg = ug+((p*vg*144)/J); 			#The enthalpy of saturated steam 			#1 ft**2 = 144 in**2 			#Btu/lbm

# Results
print "The enthalpy of saturated steam at 90 F is %.2f Btu/lbm"%(hg); 			#The value is matched with the value in table 1
The enthalpy of saturated steam at 90 F is 1100.69 Btu/lbm

Example 5.2 Page No : 187

In [2]:
# given data
p = 4.246; 			#Unit:kPa 			#absolute pressure
vg = 32.894; 			#Unit:m**3/kg 			#specific volume
ug = 2416.6; 			#Unit:kJ/kg 			#internal energy
J = 778; 			#J = Conversion factor
# 1 Btu  =  778 ft*LBf
#h = u+(p*v)
hg = ug+(p*vg); 			#The enthalpy of saturated steam 			#1 ft**2 = 144 in**2 			#unit:kJ/kg

# Results
print "The enthalpy of saturated steam at 30 C is %.2f kJ/kg"%(hg); 			#The value is matched with the value in table 1
The enthalpy of saturated steam at 30 C is 2556.27 kJ/kg

Example 5.3 Page No : 188

In [3]:
#The necessary interpolations are best done in tabular forms as shown:
#  p    hg
# 115  1190.4  table 2
# 118  1190.8  (hg)118 = 1190.8
# 120  1191.1
hg = 1190.4+(3./5)*(1191.1-1190.4); 			#Btu/lbm 			#enthaply
print "The enthalpy of saturated steam at 118 psia is %.2f Btu/lbm"%(hg);

#  p    vg
# 115  3.884  table 2
# 118  3.792  (vg)118 = 3.790
# 120  3.730
vg = 3.884-(3./5)*(3.884-3.730); 			#ft**3/lbm 			#specific volume
print "The specific volume of saturated steam at 118 psia is %.2f ft**3/lbm"%(vg);

#  p    sg
# 115  1.5921  table 2
# 118  1.5900  (sg)118 = 1.5900
# 120  1.5886
sg = 1.5921-(3./5)*(1.5921-1.5886); 			#entropy
print "The entropy of saturated steam at 118 psia is %.2f"%(sg);

#  p    ug
# 115  1107.7  table 2
# 118  1108.06  (ug)118 = 1180.1
# 120  1108.3
ug = 1107.7-(3./5)*(1108.3-1107.7); 			#internal energy
print "The internal energy of saturated steam at 118 psia is %.2f"%(ug);
#The interpolation process that was done in tabular form for this problem can also be demonstated by refering to figure 5.8 for the specific volume.It will be 
#seen that the results of this problem and the tabulated values are essentially in exact agreement and that linear interpolation is satisfactory in these tables.
The enthalpy of saturated steam at 118 psia is 1190.82 Btu/lbm
The specific volume of saturated steam at 118 psia is 3.79 ft**3/lbm
The entropy of saturated steam at 118 psia is 1.59
The internal energy of saturated steam at 118 psia is 1107.34

Example 5.4 Page No : 189

In [4]:
#By defination,
#hg = ug+(p*vg)/J
#hf = uf+(p*vf)/J
#hfg  =  hg-hf  =  (ug-uf) + p*(vg-vf)/J  =  ufg + p*(vg-vf)/J
#From table 2 at 115 psia,
p = 115; 			#Unit:psia 			#absolute pressure
ufg = 798.8; 			#Unit:Btu/lbm 			#Evap. internal energy
ug = 3.884; 			#Unit:ft**3/lbm 			#Saturated vapour internal energy
vf = 0.017850; 			#Unit:ft**3/lbm 			#Saturated liquid specific volume
J = 778; 			#J = Conversion factor 			#Unit:ft*lbf/Btu
#1 ft**2 = 144 in**2
hfg = ufg+(p*144*(ug-vf))/J; 			#Evap. Enthalpy 			#Unit:Btu/lbm

# Results
print "hfg for saturated steam at 115 psia is %.2f Btu/lbm"%(hfg); 			#The tabulated values are matched
hfg for saturated steam at 115 psia is 881.09 Btu/lbm

Example 5.5 Page No : 190

In [5]:
#From table 2 at 1.0 MPa,
p = 1000; 			#Unit:kN/m**2 			#absolute pressure
ufg = 1822.0; 			#Unit:kJ/kg 			#Evap. internal energy
vf = 0.0011273; 			#Unit:m**3/kg 			#Saturated liquid specific volume
vg = 0.19444; 			#Unit:m**3/kg 			#Saturated vapour specific volume

# Calculations
vfg = vg-vf;  			#Evap. specific volume 			#m**3/kg
hfg = ufg+(p*vfg); 			#Evap. Enthalpy 			#Unit:kJ/kg

# Results
print "hfg for saturated steam at 1.0 MPa is %.2f kJ/kg"%(hfg); 			#The tabulated values are matched
hfg for saturated steam at 1.0 MPa is 2015.31 kJ/kg

Example 5.6 Page No : 190

In [6]:
#For constant-temperature,reversible vaporization, hfg = deltah = T*deltas = T*sfg
hfg = (388.12+460)*(1.1042); 			#Evap. Enthalpy 			#Unit:Btu/lbm

# Results
print "By considering the process to be a reversible, \
\nconstant-temperature, hfg for saturated steam at 115 psia is %.2f Btu/lbm"%(hfg); 			#ans is wrong in the book
#Values are matched with tabulated values.Use of -459.67 F for absolute zero,which is the value used in table,gives almost exact agreement.
By considering the process to be a reversible, 
constant-temperature, hfg for saturated steam at 115 psia is 936.49 Btu/lbm

Example 5.7 Page No : 192

In [8]:
#Umath.sing Table 2 ans a quality of 80%(x = 0.8),we have
#at 120 psia
x = 0.8;
sf = 0.49201; 			#saturated liquid entropy 			#Unit:Btu/lbm*R
sfg = 1.0966; 			#Evap. Entropy 			#Unit:Btu/lbm*R 
hf = 312.67; 			#saturated liquid enthalpy 			#Unit:Btu/lbm 
hfg = 878.5; 			#Evap. Enthalpy 			#Unit:Btu/lbm
uf = 312.27; 			#saturated liquid internal energy 			#Unit:Btu/lbm 
ufg = 796.0; 			#Unit:Btu/lbm 			#Evap. internal energy
vf = 0.017886; 			#Saturated liquid specific volume 			#Unit:ft**3/lbm 
vfg = (3.730-0.017886); 			#evap. specific volume 			#Unit:ft**3/lbm 
sx = sf+(x*sfg); 			#entropy 			#Btu/lbm*R
print "Entropy of a wet steam mixture at 120 psia is %.2f Btu/lbm*R"%(sx);
hx = hf+(x*hfg); 			#enthalpy  			#Btu/lbm*R
print "Enthalpy of a wet steam mixture at 120 psia is %.2f Btu/lbm"%(hx);
ux = uf+(x*ufg); 			#internal energy 			#Btu/lbm*R
print "Internal energy of a wet steam mixture at 120 psia is %.2f Btu/lbm"%(ux);
vx = vf+(x*vfg); 			#/specific volume 			#ft**3/lbm
print "Specific Volume of a wet steam mixture at 120 psia is %.2f ft**3/lbm"%(vx);
#As a check,
J = 778; 			#ft*lbf/Btu 			#Conversion factor
px = 120; 			#psia 			#pressure
ux = hx-((px*vx*144)/J); 			#1 ft**2 = 144 in**2 			#internal energy
print "As a check"
print "Internal energy of a wet steam mixture at 120 psia is %.2f Btu/lbm"%(ux);
print "Which agrees with the values obtained above";
Entropy of a wet steam mixture at 120 psia is 1.37 Btu/lbm*R
Enthalpy of a wet steam mixture at 120 psia is 1015.47 Btu/lbm
Internal energy of a wet steam mixture at 120 psia is 949.07 Btu/lbm
Specific Volume of a wet steam mixture at 120 psia is 2.99 ft**3/lbm
As a check
Internal energy of a wet steam mixture at 120 psia is 949.11 Btu/lbm
Which agrees with the values obtained above

Example 5.8 Page No : 193

In [9]:
#Using Table 2 ans a quality of 85%(x = 0.85),we have
#at 1.0 MPa
x = 0.85;
sf = 2.1387; 			#saturated liquid entropy 			#Unit:kJ/kg*K
sfg = 4.4487; 			#Evap. Entropy 			#Unit:kJ/kg*K 
hf = 762.81; 			#saturated liquid enthalpy 			#Unit:kJ/kg
hfg = 2015.3; 			#Evap. Enthalpy 			#Unit:kJ/kg
uf = 761.68; 			#saturated liquid internal energy 			#Unit:kJ/kg
ufg = 1822.0; 			#Unit:kJ/kg 			#Evap. internal energy
vf = 1.1273; 			#Saturated liquid specific volume 			#Unit:m**3/kg 
vfg = (194.44-1.1273); 			#evap. specific volume 			#Unit:m**3/kg 
sx = sf+(x*sfg); 			#entropy 			#kJ/kg*K
print "Entropy of a wet steam mixture at 1.0 MPa  is %.2f kJ/kg*K"%(sx);
hx = hf+(x*hfg); 			#enthalpy 			#kJ/kg*K
print "Enthalpy of a wet steam mixture at 1.0 MPa is %.2f kJ/kg"%(hx);
ux = uf+(x*ufg); 			#internal energy  			#kJ/kg*K
print "Internal energy of a wet steam mixture at 1.0 MPa is %.2f kJ/kg"%(ux);
vx = (vf+(x*vfg))*(0.001); 			#specific volume 			#m**3/kg
print "Specific Volume of a wet steam mixture at 1.0 MPa is %.2f m**3/kg"%(vx);
#As a check,
px = 10**6; 			#psia 			#pressure
ux = hx-((px*vx)/10**3); 			#1 ft**2 = 144 in**2 			#internal energy
print "As a check%"
print "Internal energy of a wet steam mixture at 120 psia is %.2f kJ/kg"%(ux);
print "Which agrees with the values obtained above";
Entropy of a wet steam mixture at 1.0 MPa  is 5.92 kJ/kg*K
Enthalpy of a wet steam mixture at 1.0 MPa is 2475.81 kJ/kg
Internal energy of a wet steam mixture at 1.0 MPa is 2310.38 kJ/kg
Specific Volume of a wet steam mixture at 1.0 MPa is 0.17 m**3/kg
As a check%
Internal energy of a wet steam mixture at 120 psia is 2310.37 kJ/kg
Which agrees with the values obtained above

Example 5.9 Page No : 193

In [10]:
#For the wet mixture,hx = hf+(x*hfg),solving for x gives us
#Using table 1,we have,
hx = 900; 			#Btu/lbm 			#Enthalpy of wet mixture at 90F
hf = 58.07; 			#Btu/lbm 			#saturated liquid enthalpy
hfg = 1042.7; 			#Btu/lbm 			#Evap. Enthalpy

# Calculations
x = (hx-hf)/hfg; 			#quality

# Results
print "The quality is %.2f percentage of a wet steam at 90F"%(x*100);
The quality is 80.75 percentage of a wet steam at 90F

Example 5.10 Page No : 194

In [11]:
#For the wet mixture,hx = hf+(x*hfg),solving for x gives us
#Using table 1,we have,
hx = 2000; 			#kJ/kg 			#Enthalpy of wet mixture at 30 C
hf = 125.79; 			#kJ/kg 			#saturated liquid enthalpy
hfg = 2430.5; 			# 			#Evap. Enthalpy 			#kJ/kg

# Calculations
x = (hx-hf)/hfg; 			#quality

# Results
print "The quality is %.2f percentage of a wet steam at 30 C"%(x*100);
The quality is 77.11 percentage of a wet steam at 30 C

Example 5.11 Page No : 197

In [12]:
#The values of temperature and pressure are listed in Table 3(Figure 5.10) and can be read directly.
print "Specific volume of superheated steam at 330 psia and 450F is v = 1.4691 ft**3/lbm";
print "Internal Energy of superheated steam at 330 psia and 450F is u = 1131.8 Btu/lbm";
print "Enthalpy of superheated steam at 330 psia and 450F is h = 1221.5 Btu/lbm";
print "Entropy of superheated steam at 330 psia and 450F is s = 1.5219 Btu/lbm*R";
Specific volume of superheated steam at 330 psia and 450F is v = 1.4691 ft**3/lbm
Internal Energy of superheated steam at 330 psia and 450F is u = 1131.8 Btu/lbm
Enthalpy of superheated steam at 330 psia and 450F is h = 1221.5 Btu/lbm
Entropy of superheated steam at 330 psia and 450F is s = 1.5219 Btu/lbm*R

Example 5.12 Page No : 197

In [13]:
#The values of temperature and pressure are listed in Table 3(Figure 5.10) and can be read directly.
print "Specific volume of superheated steam at 2.0 MPa and 240 C is v = 0.10845 m**3/lbm";
print "Internal Energy of superheated steam at 2.0 MPa and 240 C is u = 2659.6 kJ/kg";
print "Enthalpy of superheated steam at 2.0 MPa and 240 C is h = 2876.5 kJ/kg";
print "Entropy of superheated steam at 2.0 MPa and 240 C is s = 6.4952 kJ/kg*K";
Specific volume of superheated steam at 2.0 MPa and 240 C is v = 0.10845 m**3/lbm
Internal Energy of superheated steam at 2.0 MPa and 240 C is u = 2659.6 kJ/kg
Enthalpy of superheated steam at 2.0 MPa and 240 C is h = 2876.5 kJ/kg
Entropy of superheated steam at 2.0 MPa and 240 C is s = 6.4952 kJ/kg*K

Example 5.13 Page No : 197

In [14]:
#The necessary interpolations(between 450F and 460F at 330 psia) are best done in tabular forms as shown:
#  t    v
# 460  1.4945  
# 455  1.4818  
# 450  1.4691
v = 1.4691+(1./2)*(1.4945-1.4691); 			#ft**3/lbm 			#specific volume
print "The specific volume of saturated steam at 330 psia & 455F is %.2f ft**3/lbm"%(v);

#  t    u
# 460  1137.0  
# 455  1134.4  
# 450  1131.8 
u = 1131.8+(1./2)*(1137.0-1131.8); 			#Btu/lbm 			#internal energy
print "The internal energy of saturated steam at 330 psia & 455F is %.2f Btu/lbm"%(u);

#  t    h
# 460  1228.2  
# 455  1224.9  
# 450  1221.5
h = 1221.5+(1./2)*(1228.2-1221.5); 			#enthaply 			#Btu/lbm
print "The enthalpy of saturated steam at 330 psia & 455F is %.2f Btu/lbm"%(h);

#  t    s
# 460  1.5293  
# 455  1.5256  
# 450  1.5219
s = 1.5219+(1./2)*(1.5293-1.5219); 			#entropy 			#Btu/lbm*R
print "The entropy of saturated steam at 330 psia & 455F is %.2f Btu/lbm*R"%(s);
The specific volume of saturated steam at 330 psia & 455F is 1.48 ft**3/lbm
The internal energy of saturated steam at 330 psia & 455F is 1134.40 Btu/lbm
The enthalpy of saturated steam at 330 psia & 455F is 1224.85 Btu/lbm
The entropy of saturated steam at 330 psia & 455F is 1.53 Btu/lbm*R

Example 5.14 Page No : 198

In [15]:
#From Table3, we first obtain the properties at 337 psia and 460 F and then 337 psia and 470 F.
#The necessary interpolations are best done in tabular forms as shown:
#Proceeding with the calculation,at 460 F,
#  p    v                                     			#  p    h
# 340  1.4448                                 			# 340  1226.7
# 337  1.4595                                 			# 337  1227.2
# 335  1.4693                                 			# 335  1227.5
v = 1.4696-(2./5)*(1.4693-1.4448);                
h = 1227.5-(2./5)*(1227.5-1226.7);
#ft**3/lbm 			#specific volume                   			#Btu/lbm 			#enthaply

#And at 470 F,
#  p    v                                     			#  p    h
# 340  1.4693                                 			# 340  1233.4
# 337  1.4841                                 			# 337  1233.9
# 335  1.4940                                 			# 335  1234.2
v = 1.4640-(2./5)*(1.4640-1.4693);                
h = 1234.2-(2./5)*(1234.2-1233.4);
#ft**3/lbm 			#specific volume                   			#Btu/lbm 			#enthaply

#Therefore,at 337 psia and 465 F
#  t    v                                    			#  t    h
# 470  1.4841                                			# 470  1233.9
# 465  1.4718                                			# 465  1230.7
# 460  1.4595                                			# 460  1227.5
v = 1.4595+(1./2)*(1.4841-1.4595);               
h = 1227.5+(1./2)*(1233.9-1227.5);
#ft**3/lbm 			#specific volume                  			#Btu/lbm 			#enthaply
print "At 465 F and 337 psia, specific volume = %.2f ft**3/lbm and enthalpy = %.2f Btu/lbm"%(v,h);
At 465 F and 337 psia, specific volume = 1.47 ft**3/lbm and enthalpy = 1230.70 Btu/lbm

Example 5.15 Page No : 202

In [16]:
#The values of temperature and pressure are listed in Table 4(Figure 5.10) and can be read directly.
print "Specific volume of subcooled water at 1000 psia and 300F is v = 0.017379 ft**3/lbm";
print "Internal Energy of subcooled water at 1000 psia and 300F is u = 268.24 Btu/lbm";
print "Enthalpy of subcooled water at 1000 psia and 300F is h = 271.46 Btu/lbm";
print "Entropy of subcooled water at 1000 psia and 300F is s = 0.43552 Btu/lbm*R";
Specific volume of subcooled water at 1000 psia and 300F is v = 0.017379 ft**3/lbm
Internal Energy of subcooled water at 1000 psia and 300F is u = 268.24 Btu/lbm
Enthalpy of subcooled water at 1000 psia and 300F is h = 271.46 Btu/lbm
Entropy of subcooled water at 1000 psia and 300F is s = 0.43552 Btu/lbm*R

Example 5.16 Page No : 202

In [17]:
#It is necessary to ontain the saturation values corresponding to 300 F.This is done by reading Table A.1 in Appendix 3,which gives
pf = 66.98; 			#psia 			#pressure
vf = 0.017448; 			#ft**3/lbm 			#specific volume
hf = 269.73; 			#Btu/lbm 			#enthaply
#Now,
p = 1000; 			#psia 			#pressure
J = 778; 			#Conversion factor 			#ft*lbf/Btu
#From eq.5.5,
h = hf+((p-pf)*vf*144)/J; 			#1ft**2 = 144 in**2 			#The enthalpy of subcooled water 			#Btu/lbm
print "The enthalpy of subcooled water is %.2f Btu/lbm"%(h);
#The difference between this value and the value found in problem 5.15,expressed as a percentage is
percentoferror = (h-271.46)/271.46;
print "Percent of error is %.2f"%(percentoferror*100);
The enthalpy of subcooled water is 272.74 Btu/lbm
Percent of error is 0.47

Example 5.25 Page No : 211

In [18]:
#On a chart in Appendix 3,it is necessary to estimate the 90 F point on the saturation line.From the chart or the table in the upper left of the chart,we note that 90 F is between 1.4 and 1.5 in. of mercury.Estimating the intersection of this value with the saturation curve yields
print "Enthalpy of saturated steam hg = 1100 Btu/lbm";
#This is a good agreement with results of problem 5.1
Enthalpy of saturated steam hg = 1100 Btu/lbm

Example 5.26 Page No : 212

In [19]:
#The Mollier chart has lines of constant moisture in the wet region which correspond to (1-x).Therefore,we read at 20% moisture(80% Quality) and 120 psia,
print "The enthalpy of a wet steam mixture at 120 psia having quality 80 percent is 1015 Btu/lbm";
#Which also agrees well with the calculated value in problem 5.7
The enthalpy of a wet steam mixture at 120 psia having quality 80 percent is 1015 Btu/lbm

Example 5.27 Page No : 213

In [20]:
#Entering the Mollier chart at 900 Btu/lbm and estimating 90 F(near the 1.5-in. Hg dashed line) yields a constant moisture percent of 19.2%.
print "The quality is %.2f percent"%((1-0.192)*100);
#We show good agreement with the calculated value.
The quality is 80.80 percent

Example 5.28 Page No : 214

In [21]:
#From the chart,
print "The enthalpy of steam at 330 psia is h = 1220 Btu/lbm";
#Compared to 1221.5 Btu/lbm found in problem 5.11
The enthalpy of steam at 330 psia is h = 1220 Btu/lbm

Example 5.29 Page No : 214

In [22]:
#We note that the steam is superheated.From the Mollier chart in SI units,
print "The enthalpy h = 2876.5 kJ/kg and entropy s = 6.4952 kJ/kg*K";
#Values are matched with problem 5.12
The enthalpy h = 2876.5 kJ/kg and entropy s = 6.4952 kJ/kg*K

Example 5.30 Page No : 215

In [23]:
#Because neither pressure nor temperature is shown directly,it is necessary to estimate to obtain the desired value.
print "The enthalpy of steam is h = 1231 Btu/lbm";
#In problem 5.14,h = 1230.7 Btu/lbm,Which is matched here.
The enthalpy of steam is h = 1231 Btu/lbm

Example 5.31 Page No : 215

In [24]:
#Reading the chart at 30 C and saturation gives us,
print "The enthalpy of saturated steam is hg = 2556 kJ/kg";
#Which matches with value of problem 5.2
The enthalpy of saturated steam is hg = 2556 kJ/kg

Example 5.32 Page No : 215

In [25]:
#Reading the chart in wet region at 1.0 MPa and x = 0.85(moisture of 15%) gives us
print "hx = 2476 kJ/kg and sx = 5.92 kJ/kg*K";
#The chart does not give ux or vx directly
hx = 2476 kJ/kg and sx = 5.92 kJ/kg*K

Example 5.33 Page No : 215

In [26]:
#Locate 30 C on the saturation line.Now follow a line of constant pressure,which is also a line of constant temperature in wet region,until an enthalpy of 2000kJ/kg is reached.
print "The moisture content is 23 percent or x = 77 percent";
The moisture content is 23 percent or x = 77 percent

Example 5.34 Page No : 216

In [27]:
#We enter the chart in the superheat region at 2.0MPa and 240 C to read the enthalpy and entropy.This procedure gives 
print "Enthalpy h = 2877 kJ/kg and entropy s = 6.495 kJ/kg*K";
#The other properties cant be obtained directly from the chart
Enthalpy h = 2877 kJ/kg and entropy s = 6.495 kJ/kg*K

Example 5.35 Page No : 218

In [28]:
#As already noted,h1 = h2 for this process.On the Mollier chart,h2 is found to be 1170 Btu/lbm at 14.7 psia and 250 F.Proceeding to the left on the chart,the constant-enthalpy value of 1170 Btu/lbm to 150 psia yields a moisture of 3% or a quality of 97%.
#If we use the tables to obtain the solution to this problem,we would first obtain h2 from the superheated vapor tables as 1168.8 Btu/lbm.Because hx = hf+(x*hfg),we    obtain x as
hx = 1168.8; 			#Btu/lbm 
hf = 330.75; 			#Btu/lbm 			#values of 150 psia
hfg = 864.2; 			#Btulbm 			#values of 150 psia
x = (hx-hf)/hfg; 			#Quality
print "Moisture in the steam flowing in the pipe is %.2f percent"%((1-x)*100);
print "or quality of the steam is %.2f percent"%(x*100);
#very often,it is necessary to perform multiple interpolations if the tables are used,and the Mollier chart yields results within the rquired accuracy for most engineering problems and saves considerable time.
#We can also use the computerised programs to solve this program.We first enter the 250F and 14.7 psia to obtain h of 1168.7 Btu/lbm.We then continue by entering h of 1168.7 Btu/lbm and p of 150 psia.The printout gives us x of 0.9699 or 97%.While the computer solution is quick and easy to use,you should still sketch out the    problem on an h-s or T-s diagram to show the path of the process.

# Saturation Properties
#--------------------------
# T = 250.00 degF
# P = 29.814 psia
#       z         z1        zg
# v(ft**3/lbm)   0.01700   13.830
# h(Btu/lbm)     218.62   1164.1
# s(Btu/lbm*F)  0.3678    1.7001
# u(Btu/lbm)    218.52    1087.8

#Thermo Properties
#------------------------
# T =  250.00 degF
# P =  14.700 psia
# v =  28.417 ft**3/lbm
# h =  1168.7 Btu/lbm
# s =  1.7831 Btu/lbm*F
# u =  1091.4 Btu/lbm

# Saturation Properties
#--------------------------
# T = 340.06 degF
# P = 118.00 psia
#       z        z1        zg
# v(ft**3/lbm)   0.01787   3.7891
# h(Btu/lbm)    311.39    1190.7
# s(Btu/lbm*F)  0.4904    1.5899
# u(Btu/lbm)    311.00    1108.0

#Thermo Properties
#------------------------
# T =  358.49 degF
# P =  150.00 psia
# v =  2.9248 ft**3/lbm
# h =  1168.7 Btu/lbm
# s =  1.5384 Btu/lbm*F
# u =  1087.5 Btu/lbm
# x =  0.9699

#Region:Saturated
Moisture in the steam flowing in the pipe is 3.03 percent
or quality of the steam is 96.97 percent

Example 5.36 Page No : 219

In [29]:
#Because the math.tank volume is 10 ft**3,the final specific volume of the steam is 10 ft**3/lbm.Interpolations in Table A.2 yield a final pressure of 42 psia.The heat    added is simply difference in internal energy between the two states.
u2 = 1093.0; 			#internal energy 			#Btu/lbm
u1 = 117.95; 			#internal energy 			#Btu/lbm

# Calculations
q = u2-u1; 			#heat added 			#Btu/lbm

# Results
print "The final pressure is 42 psia and the heat added is %.2f Btu/lbm"%(q);
The final pressure is 42 psia and the heat added is 975.05 Btu/lbm

Example 5.37 Page No : 220

In [30]:
#The mass in the math.tank is constant,and the heat added will be the change in internal energy of the contents of the math.tank between the two states.The initial mass in     the math.tank is found as follows:
Vf = 45; 			#volume of water 			#ft**2
vf = 0.016715;
Vg = 15; 			#Volume of steam 			#ft**2
vg = 26.80;
mf = Vf/vf; 			#lbm
mg = Vg/vg; 			#lbm
total = mf+mg; 			#total mass
#The internal energy is the sum of the internal energy of the liquid plus vapor:
ug = 1077.6;
uf = 180.1;
Ug = mg*ug; 			#Btu
Uf = mf*uf; 			#Btu
Total = Ug+Uf; 			#total internal energy
print "The total internal energy is %.2f Btu"%(Total);
#Because the mass in the math.tank is constant,the final specific volume must equal the initial specific volume,or
vx = (Vf+Vg)/(mf+mg); 			#ft**3/lbm
#But vx = vf+(x*vfg).Therefore Using table A.2 at 800 psia,
vx = 0.022282;
vf = 0.02087;
vfg = 0.5691-0.02087;
x = (vx-vf)/vfg;
print "The final amount of vapor is %.2f lbm"%(x*total);  			#x*total mass 
mg = x*total;
print "The final amount of liquid is %.2f lbm"%(total-x*total); 			#total mass minus final amount of vapor 
mf = total-(x*total);
#The final internal energy is found as before:
ug = 1115.0;
uf = 506.6;
Ug = mg*ug; 			#Btu
Uf = mf*uf; 			#Btu
Total1 = Ug+Uf;
difference = Total1-Total; 			#final internal energy-initial internal energy
#per unit mass heat added is,
print "The heat added per unit is %.2f Btu/lbm"%(difference/total); 			#the difference of internal energy/total mass
The total internal energy is 485467.03 Btu
The final amount of vapor is 6.94 lbm
The final amount of liquid is 2685.82 lbm
The heat added per unit is 327.88 Btu/lbm

Example 5.38 Page No : 222

In [31]:
#As shown in Fig. 5.21b,the process dscribed in this problem is a vertical line on the Mollier Chart.For 800 psia and 600F,the Mollier chart yeilds h1 = 1270 Btu/lbm and s1 = 1.485.Proceeding vertically down the chart at constant s to 200 psia yields a final enthalpy h2 = 1148 Btu/lbm.The change in enthaly Using the process is        1270-1148 = 122 Btu/lbm.
#We may also solve this problem Using the steam tables in Appendix 3.Thus,the enthalpy at 800 psia and 600 F is 1270.4 Btu/lbm,and its entropy is 1.4861 Btu/lbm*R.
#Because the process is isentropic,the final entropy at 200psia must be 1.4861.From the saaturation table,the entropy of saturated steam at 200 psia is 1.5464,which indicates the final steam condition must be wet because the entropy of the final steam is less than the entropy of saturation.Using the wet steam relation yields,
#sx = sf+(x*sfg)
h1 = 1270.4; 
sx = 1.4861; 
sf = 0.5440; 
sfg = 1.0025 ;
hf = 355.6; 
hfg = 843.7;

# Calculations
x = (sx-sf)/sfg; 			#Quality
#Therefore,the final enthalpy is
hx = hf+(x*hfg); 			#Btu/lbm

# Results
print "The final enthalpy is %.2f Btu/lbm"%(hx);
print "The change in enthalpy is %.2f Btu/lbm"%(h1-hx); 			#Note the agreement with the Mollier chart solution
#we can also use the computer program to solve this problem.For 600F and 800 psia, h = 1270. Btu/lbm and s = 1.4857 Btu/lbm*R.Now Using p = 200 psia and s = 1.4857,we obtain
#h = 1148.1 Btu/lbm.The change in enthalpy is 1270.0-1148.1 = 121.9 Btu/lbm.Note the effort saved Using either the Mollier chart or the computer program. 

# Saturation Properties
#--------------------------
# T = 600.00 degF
# P = 1541.7 psia
#       z        z1         zg
# v(ft**3/lbm)   0.02362   0.2675
# h(Btu/lbm)    616.59    1166.2
# s(Btu/lbm*F)  0.8129    1.3316
# u(Btu/lbm)    609.85    1089.9

#Thermo Properties
#------------------------
# T =  600.00 degF
# P =  800.00 psia
# v =  0. ft**3/lbm
# h =  1168.7 Btu/lbm
# s =  1.5384 Btu/lbm*F
# u =  1087.5 Btu/lbm
#Region:Superheated

# Saturation Properties
#--------------------------
# T = 381.87 degF
# P = 200.00 psia
#       z        z1        zg
# v(ft**3/lbm)   0.01839   2.2883
# h(Btu/lbm)    355.60    1199.0
# s(Btu/lbm*F)  0.5440    1.5462
# u(Btu/lbm)    354.92    1114.3

#Thermo Properties
#------------------------
# T =  381.87 degF
# P =  200.00 psia
# v =  2.1512 ft**3/lbm
# h =  1148.1 Btu/lbm
# s =  1.4857 Btu/lbm*F
# u =  1068.5 Btu/lbm
# x =  0.9396

#Region:Saturated
The final enthalpy is 1148.47 Btu/lbm
The change in enthalpy is 121.93 Btu/lbm

Example 5.39 Page No : 226

In [33]:
#As refering to figure 5.21,it will be seen that the final temperature and enthalpy will both be higher than for the isentropic case.
#80% of the isentropic enthalpy difference
deltah = 0.8*122; 			#change in enthalpy 			#Btu/lbm
h1 = 1270; 			#Btu/lbm 			#initial enthalpy
h2 = h1-deltah; 			#the final enthalpy 			#Btu/lbm
print "The final enthalpy is %.2f Btu/lbm"%(h2);
print "and the final pressure is 200 psia";
print "The Mollier chart indicates the final state to be in the wet region"
print "with 3.1percent moisture content and an entropy of 1.514 Btu/lbm*R";
The final enthalpy is 1172.40 Btu/lbm
and the final pressure is 200 psia
The Mollier chart indicates the final state to be in the wet region
with 3.1percent moisture content and an entropy of 1.514 Btu/lbm*R

Example 5.40 Page No : 226

In [35]:
#Using the Mollier chart,
h1 = 2942; 			#kJ/kg 			#initial enthalpy
#Proceeding as shown in figure 5.21b,that is,vertically at constant entropy to a pressure of 0.1 MPa,gives us
h2 = 2512; 			#kJ/kg 			#final enthalpy

# Results
print "Neglecting kinetic & potential energy, The change in enthalpy of the steam is %.2f kJ/kg"%(h1-h2);
Neglecting kinetic & potential energy, The change in enthalpy of the steam is 430.00 kJ/kg

Example 5.41 Page No : 226

In [38]:
import math

#From the conditions given in problem 5.38,the isentropic change in enthalpy is 122 Btu/lbm
#So,
h1minush2 = 122; 			#Btu/lbm 			#change in enthalpy
J = 778; 			#Conversion factor
gc = 32.17; 			#lbm*ft/lbf*s**2 			#constant of proportionality

# Calculations
V2 = math.sqrt(2*gc*J*(h1minush2)); 			#final velocity 			#ft/s

# Results
print "As the steam leaves the nozzle, The final velocity is %.2f ft/s"%(V2);
As the steam leaves the nozzle, The final velocity is 2471.21 ft/s

Example 5.42 Page No : 227

In [40]:
import math

#Because the process is irreversible,we cannot show it on the Mollier diagram.However,the analysis of problem 3.22 for the nozzle is still valid,and all that is      needed is the enthalpy at the beginning and the end of the expansion.From the problem 5.38,
h1 = 1270; 			#Btu/lbm 			#initial enthalpy
#For h2 we locate the state point on the Mollier diagram as being saturated vapor at 200 psia.This gives us
h2 = 1199; 			#Btu/lbm 			#final enthalpy
J = 778; 			#Conversion factor
gc = 32.17; 			#lbm*ft/lbf*s**2 			#constant of proportionality
V2 = math.sqrt(2*gc*J*(h1-h2)); 			#final velocity 			#Ft/s
print "As the steam leaves the nozzle, The final velocity is %.2f ft/s"%(V2);
As the steam leaves the nozzle, The final velocity is 1885.21 ft/s

Example 5.43 Page No : 229

In [41]:
#From the saturation table,500 psia corresponds to a temperature of 467.13F,and the saturated vapor has an enthalpy of 1205.3 Btu/lbm.At 500 psia and 800 F,the       saturated vapor has an enthalpy of 1412.1 Btu/lbm.Because this process is a steady-flow process at constant pressure,the energy equation becomes q = h2-h1,assuming    that differences in the kinetic energy and potential energy terms are negligible.Therefore,
h2 = 1412.1; 			#Btu/lbm 			#final enthalpy
h1 = 1205.3; 			#Btu/lbm 			#initial enthalpy

# Calculations
q = h2-h1; 			#heat added 			#Btu/lbm

# Results
print "%.2f Btu/lbm heat per pound of steam was added"%(q);
206.80 Btu/lbm heat per pound of steam was added

Example 5.44 Page No : 229

In [44]:
#From the saturation table at 1 psia,
hf = 69.74; 			#Btu/lbm 			#saturated liquid enthalpy
hfg = 1036.0; 			#Btu/lbm 			#Evap. Enthalpy
hg = 1105.8; 			#Btu/lbm 			#The enthalpy of saturated steam
x = 0.97; 			#Quality

# Calculations and Results
#Because the condensation process is carried out at constant pressure,the energy equation is q = deltah.
hx = hf+(x*hfg); 			#the initial enthalpy 			#Btu/lbm
print "The initial enthalpy is %.2f Btu/lbm"%(hx);
#The final enthalpy is hf = 69.74.So,
deltah = hx-hf; 			#The enthalpy difference 			#Btu/lbm
print "At 1 psia, The enthalpy difference is %.2f Btu/lbm"%(deltah);
print "By the computer solution,the enthalpy difference is 1004.6 Btu/lbm";
# Saturation Properties
#--------------------------
# T = 101.71 degF
# P = 1.0000 psia
#       z        z1         zg
# v(ft**3/lbm)   0.01614   333.55
# h(Btu/lbm)    69.725    1105.4
# s(Btu/lbm*F)  0.1326    1.9774
# u(Btu/lbm)    69.722    1043.6

#Thermo Properties
#------------------------
# T =  101.71 degF
# P =  1.0000 psia
# v =  323.55 ft**3/lbm
# h =  1074.3 Btu/lbm
# s =  1.9221 Btu/lbm*F
# u =  1014.4 Btu/lbm
# x =  0.9700

#Region:Saturated
The initial enthalpy is 1074.66 Btu/lbm
At 1 psia, The enthalpy difference is 1004.92 Btu/lbm
By the computer solution,the enthalpy difference is 1004.6 Btu/lbm