Pure Substances

Chapter 3
Properties of a pure substance

Read BS, Chapter 3
3.1 The pure substance
We dene a
Pure substance: a material with homogeneous and invariable composition.
To elaborate,
Pure substances can have multiple phases: an ice-water mixture is still a pure sub-
stance.
An air-steam mixture is not a pure substance.
Air, being composed of a mixture of N
2
, O
2
, and other gases, is formally not a pure
substance. However, experience shows that we can often treat air as a pure substance
with little error.
3.2 Vapor-liquid-solid phase equilibrium
Often we nd that dierent phases of pure substances can exist in equilibrium with one
another. Let us consider an important gedankenexperiment (Latin-German for thought
experiment) in which we boil water. Ordinary water boiling is shown in Fig. 3.1. However,
this ordinary experiment has constraints which are too loose. Most importantly, the mass
of water leaks into the atmosphere; thus, the water vapor and the air become a mixture and
no longer a pure substance.
Let us instead consider a more controlled piston-cylinder arrangement. Inside the cylin-
der, we begin with pure liquid water at T = 20

C. The piston is free to move in the cylinder,
but it is tightly sealed, so no water can escape. On the other side of the piston is a constant
pressure atmosphere, which we take to be at P = 100 kPa = 0.1 MPa = 10
5
Pa = 1 bar.
43
44 CHAPTER 3. PROPERTIES OF A PURE SUBSTANCE
Figure 3.1: Water boiling isobarically in an open environment.
We slowly add heat to the cylinder, and observe a variety of interesting phenomena. A sketch
of what we observe is given in Fig. 3.2. We notice the following behavior:
P = 100 kPa
liquid water
T = 20 C
P = 100 kPa
liquid water
T > 20 C
P = 100 kPa
saturated
liquid water
T = 99.62 C
P = 100 kPa
liquid,T = 99.62 C
P = 100 kPa
saturated
water vapor
T = 99.62 C
Q
P = 100 kPa
water vapor
T > 99.62 C
Q Q Q Q
vapor, T=99.62 C
Figure 3.2: Sketch of experiment in which heat is added isobarically to water in a closed
piston-cylinder arrangement.
The pressure remains at a constant value of 100 kPa. This is an isobaric process.
The total volume increases slightly as heat is added to the liquid.
The temperature of the liquid increases signicantly as heat is added to the liquid.
At a special value of temperature, observed to be T = 99.62

C, we have all liquid, but
cannot add any more heat and retain all liquid. We will call this state the saturated
liquid state. We call T = 99.62

C the saturation temperature at P = 100 kPa. As we
continue to add heat,
CC BY-NC-ND. 04 May 2012, J. M. Powers.
3.2. VAPOR-LIQUID-SOLID PHASE EQUILIBRIUM 45
The temperature remains constant (this is isothermal now as well as isobaric).
The total volume continues to increase.
We notice two phases present: liquid and vapor, with a distinct phase boundary.
The liquid is dense relative to the vapor. That is
f
>
g
, where f denotes uid
or liquid and g denotes gas or vapor. Thus, v
g
> v
f
.
As more heat is added, more vapor appears, all while P = 100 kPa and T =
99.62

C.
At a certain volume, we have all vapor and no liquid, still at P = 100 kPa,
T = 99.62

C. We call this state the saturated vapor state.
As heat is added, we nd both the temperature and the volume rise, with the pressure
remaining constant. The water remains in the all vapor state.
We have just boiled water! We sketch this process in the temperature-specic volume plane,
that is, the T v plane, in Fig. 3.3. Note that the mass m of the water is constant in this
v
T
two-phase
mixture
v
f
v
g
P

=

1
0
0

k
P
a

(
i
s
o
b
a
r
)
saturated
liquid saturated
vapor
Figure 3.3: Isobar in the T v plane for our thought experiment in which heat is added
isobarically to water in a piston-cylinder arrangement.
problem, so the extensive V is strictly proportional to specic volume, v = V/m.
We next repeat this experiment at lower pressure (such as might exist on a mountain
top) and at a higher pressure (such as might exist in a valley below sea level). For moderate
pressures, we nd qualitatively the exact same type of behavior. The liquid gets hotter,
turns into vapor isothermally, and then the vapor gets hotter as the heat is added. However,
we note the following important facts:
The saturation temperature (that is the boiling point) increases as pressure increases,
as long as the pressure increase is not too high.
As pressure increases v
f
becomes closer to v
g
.
Above a critical pressure, P = P
c
= 22.089 MPa, there is no phase change observed!
1
At the critical pressure, the temperature takes on a critical value of T
c
= 374.14

C.
At the critical pressure and temperature, the specic volume takes the value v
f
= v
g
=
v
c
= 0.003155 m
3
/kg.
We see how the boiling point changes with pressure by plotting the saturation pressure
as a function of saturation temperature in the T P plane in Fig. 3.4. This is the so-called
vapor pressure curve. Here, we focus on liquid-vapor mixtures and keep T high enough to
T
P
superheated
vapor
compressed
liquid
critical
point
T
c
P
c
v
a
p
o
r
p
r
e
s
s
u
r
e
c
u
r
v
e
Figure 3.4: Saturation pressure versus saturation temperature sketch.
prevent freezing. Note the curve terminates abruptly at the critical point.
We adopt the following nomenclature:
saturated liquid: the material is at T
sat
and is all liquid.
saturated vapor: the material is at T
sat
and is all vapor.
compressed (subcooled) liquid: the material is liquid with T < T
sat
.
superheated vapor: the material is vapor with T > T
sat
.
two-phase mixture: the material is composed of co-existing liquid and vapor with
both at T
sat
.
1
This behavior may have rst been carefully documented by T. Andrews, 1869, The Bakerian lecture:
on the continuity of the gaseous and liquid states of matter, Philosophical Transactions of the Royal Society
of London, 159: 575-590.
3.2. VAPOR-LIQUID-SOLID PHASE EQUILIBRIUM 47
For two-phase mixtures, we dene a new property to characterize the relative concentra-
tions of liquid and vapor. We dene the
quality= x: as the ratio of the mass of the mixture which is vapor (vap) to the total
mixture mass:
x =
m
vap
m
total
. (3.1)
We can also take the total mass to be the sum of the liquid and vapor masses:
m
total
= m
liq
+ m
vap
. (3.2)
So
x =
m
vap
m
liq
+ m
vap
. (3.3)
There are two important limits to remember:
x = 0: corresponds to m
vap
= 0. This is the all liquid limit.
x = 1: corresponds to m
vap
= m
total
. This is the all gas limit.
We must have
0 x 1. (3.4)
We sketch waters T v plane again for a wide variety of isobars in Fig. 3.5. We sketch
v
T
saturated
vapor line saturated
liquid line
superheated
vapor
two-phase
mixture
c
o
m
p
r
e
s
s
e
d
l
i
q
u
i
d
critical
point
P = 0.1 MPa
P
c
= 22.089 MPa
P = 10 MPa
P = 1 MPa
P = 40 MPa
T
c
= 374.14C
T= 99.62C
v
c
= 0.003155 m
3
/kg
Figure 3.5: Sketch of T v plane for water for a variety of isobars.
waters P v plane for a wide variety of isotherms in Fig. 3.6. We can perform a similar
v
P
saturated
vapor line
saturated
liquid line
superheated
vapor
two-phase
mixture
c
o
m
p
r
e
s
s
e
d
l
i
q
u
i
d
critical
point
T
T
low
high
T
c
= 374.14C
P
c
= 22. 089 MPa
v
c
= 0.003155 m
3
/kg
Figure 3.6: Sketch of P v plane for water for a variety of isotherms.
thought experiment for ice. We can start with ice at P = 100 kPa and add heat to it. We
observe the ices temperature rise until T = T
sat
0

C. At that temperature, the ice begins
to melt and the temperature remains constant until all the ice is melted. At this point the
liquid temperature begins to rise. If we continued to add heat, we would boil the water.
We note if we perform this experiment for P < 0.6113 kPa the ice in fact goes directly
to vapor. It is said to have undergone sublimation. There exists a second important point
where ice being heated isobarically can transform into either liquid or gas. This is the so-
called triple point. At the triple point we nd the saturation pressure and temperature are
P
tp
= 0.6113 kPa and T
tp
= 0.01

C, respectively. It is better described as a triple line,
because in the P v T space we will study, it appears as a line with constant P and T, but
variable v. In the P T projected plane of the P v T volume, it projects as a point. We
sketch waters P T plane again for a wider range to include the vapor-liquid-solid phase
behavior in Fig. 3.7.
These characteristics apply to all pure substances. For example, nitrogen has a triple
point and a critical point. Table A.2 in BS lists critical constants for many materials. Note
also that phase transitions can occur within solid phases. This involves a re-arrangement of
the crystal structure. This has important implications for material science, but will not be
considered in detail in this course.
3.3 Independent properties
Let us dene a
3.4. THERMAL EQUATIONS OF STATE 49
T (C)
P (kPa)
superheated
vapor
compressed
liquid
critical
point
solid
triple
point
f
u
s
i
o
n

l
i
n
e
v
a
p
o
r
i
z
a
t
i
o
n

l
i
n
e
s
u
b
lim
a
t
io
n
lin
e
supercritical
uid
0.01
0.6113
22089
374.14
Figure 3.7: Sketch of P T plane for water.
Simple compressible substance: a material that can be worked upon by pressure
forces.
While this is indeed restrictive, note we neglect electric, magnetic, and chemical work modes,
it will be important for mechanical engineering applications. The following important state-
ment can be proved (but will not be so here):
For a simple compressible substance, two independent intensive thermodynamic prop-
erties dene the state of the system.
Consider the implications for
superheated vapor: If we consider P, T, and v, this states that we must allow one of the
variables to be functions of the other two. We could have P = P(T, v), v = v(T, P),
or T = T(P, v). All are acceptable.
two-phase mixture: If we have a two-phase mixture, our experiments show that P and
T are not independent. In this case, we need another property to characterize the
system. That property is the quality, x. So for two-phase mixtures, we might have
v = v(T, x), for example.
3.4 Thermal equations of state
Here, we will describe some of the many dierent ways to capture the relation between two
independent properties and a third dependent property for a simple compressible substance.
We will focus on a so-called
Thermal equation of state: an equation which gives the pressure as a function
of two independent state variables. An example is the general form:
P = P(T, v). (3.5)
We will progress from simple thermal equations of state to more complex.
3.4.1 Ideal gas law
For many gases, especially at low density and far from the critical point, it is possible to
write a simple thermal equation of state which accurately describes the relation between
pressure, volume, and temperature. Such equations were developed in the 1600s and early
1800s based entirely on macroscopic empirical observation. In the late 1800s, statistical
mechanics provided a stronger theoretical foundation for them, but we will not consider that
here.
Let us start with the most important equation of state:
Ideal gas law: This equation, which is a combination of Boyles law,
2
Charles law,
3
and Avogadros law,
4
is most fundamentally stated as
PV = nRT. (3.6)
On the continent, Boyles law is sometimes known as Mariottes law after Edme Mariotte
(1620-1684), but Boyle published it fourteen years earlier.
5
A reproduction of Boyles data
is given in Fig. 3.8.
6
The data in (V, 1/P) space is t well by a straight line with intercept
at the origin; that is 1/P = KV , where K is a constant. Thus, PV = C, where C = 1/K is
a constant.
2
R. Boyle, 1662, New Experiments Physico-Mechanical, Touching the Air: Whereunto is Added A Defence
of the Authors Explication of the Experiments, Against the Obiections of Franciscus Linus and Thomas
Hobbes, H. Hall, Oxford. There exists a 1725 text modied from the 1662 original in The Philosophical
Works of Robert Boyle, Vol. 2, P. Shaw, ed., Innis, Osborn, and Longman, London. Boyles law holds that
PV is constant in an isothermal process.
3
attributed to J. A. C. Charles by J. L. Gay-Lussac, 1802, Recherches sur la dilatation des gaz des
vapeurs, Annales de Chimie, 43(1): 137-175. Charles law holds that V/T is constant for ideal gases
undergoing isobaric processes. Additionally, Guillaume Amontons (1663-1705) performed some of the early
experimentation which led to Charles law. John Dalton (1766-1844) is said to have also written on a version
of Charles Law in 1801.
4
A. Avogadro, 1811, Essai dune mani`ere de determiner les masses relatives des molecules elementaires
des corps, et les proportions selon lesquelles elles entrent dans ces combinaisons, Journal de Physique,
de Chimie et dHistorie Naturelle, 73:58-76. Here, Avogadro hypothesized that equal volumes of ideal or
perfect gases, at the same temperature and pressure, contain the same number of particles, or molecules.
5
C. Webster, 1965, The discovery of Boyles law, and the concept of the elasticity of air in the seventeenth
century, Archive for History of Exact Sciences, 2(6): 441-502. Also described here is how Henry Power
(1623-1668) and Richard Towneley (1629-1707) did important preliminary work which helped Boyle formulate
his law.
6
see J. B. West, 1999, The original presentation of Boyles law, Journal of Applied Physiology, 98(1):
31-39.
0 10 20 30 40 50
0.01
0.02
0.03
0.04
V (inches
3
)
1
/
P

(
i
n
c
h
e
s

H
g
)
-
1
a)
b)
Figure 3.8: a) Boyles 1662 data to validate his law (PV is constant for an isothermal
process), b) plot of Boyles data: V (column A) versus reciprocal of P (reciprocal of column
D), demonstrating its near linearity.
Depictions of Boyle, Charles, and Avogadro are given in Fig. 3.9. The ideal gas law was
a) b) c)
Figure 3.9: b) Robert Boyle (1627-1691), Irish scientist who developed an important special
case of the ideal gas law. Image from http://en.wikipedia.org/wiki/Robert Boyle, c)
Jacques Alexandre Cesar Charles (1746-1823), French scientist credited in 1802 by Joseph
Louis Gay-Lussac (1778-1850) for developing an important special case of the ideal gas
law in the 1780s. Image from http://en.wikipedia.org/wiki/Jacques Charles,
a) Lorenzo Romano Amedeo Carlo Avogadro di Quarengna e di Cerroto
(1776-1856), Italian physicist, nobleman, and revolutionary. Image from
http://en.wikipedia.org/wiki/Amedeo Avogadro.
rst stated in the form roughly equivalent to Eq. (3.6) by Clapeyron,
7
depicted in Fig. 3.10.
It is critical that the temperature here be the absolute temperature. For the original
argument, see Thomson.
8
Here, n is the number of moles. Recall there are N = 6.02214179
10
23
molecules in a mole, where Nis Avogadros number. Also R is the universal gas constant.
From experiment, it is determined to be
R = 8.314472
kJ
kmole K
. (3.7)
In this class the over bar notation will denote an intensive property on a per mole basis.
Intensive properties without over bars will be on a per mass basis. Recall the mass-basis
7
E. Clapeyron, 1834, Memoire sur la puissance motrice de la chaleur, Journal de l
Ecole Polytechnique,
14: 153-190.
8
W. Thomson, 1848, On an absolute thermometric scale founded on Carnots theory of the motive power
of heat and calculated from Regnaults observations, Proceedings of the Cambridge Philosophical Society,
1(5): 66-71.
Figure 3.10: Benot Paul

Emile Clapeyron (1799-1824), French engineer and
physicist who furthered the development of thermodynamics. Image from
http://en.wikipedia.org/wiki/Emile Clapeyron.
specic volume is v = V/m. Let us dene the mole-based specic volume as
v =
V
n
. (3.8)
Thus, the ideal gas law can be represented in terms of intensive properties as
P
V
n
..
=v
= RT, (3.9)
Pv = RT. (3.10)
There are other ways to write the ideal gas law. Recall the molecular mass M is the
mass in g of a mole of substance. More common in engineering, it is the mass in kg of a
kmole of substance. These numbers are the same! From chemistry, for example, we know
the molecular mass of O
2
is 32 g/mole = 32 kg/kmole. Symbolically, we can say that
M =
m
n
. (3.11)
Now, take the ideal gas law and divide by m:
PV = nRT, (3.12)
P
V
m
..
=v
=
n
m
..
=1/M
RT, (3.13)
Pv =
R
M
..
R
T. (3.14)
Now, let us dene
R
R
M
. (3.15)
Lets check the units:
[R] =
kJ
kmole K
kmole
kg
=
kJ
kg K
. (3.16)
We have actually just lost some universality. Recall R is independent of material. But since
each dierent gas has a dierent M, then each gas will have its own R. These values for
various gases are tabulated in Table A.5 of BS.
With this denition, the ideal gas law becomes
Pv = RT. (3.17)
This is the form we will use most often in this class. Note the useful fact that
Pv
T
= R. (3.18)
Thus, if an ideal gas undergoes a process going from state 1 to state 2, we can safely say
P
1
v
1
T
1
=
P
2
v
2
T
2
. (3.19)
Example 3.1
Find R for air.
We can model air as a mixture of N
2
and O
2
. Its average molecular mass is known from Table A.5
of BS to be M = 28.97 kg/kmole. So R for air is
R =
R
M
=
8.3145
kJ
kmole K
28.97
kg
kmole
= 0.287
kJ
kg K
. (3.20)
Consider some notions from algebra and geometry. The function f(x, y) = 0 describes a
curve in the x y plane. In special cases, we can solve for y to get the form y = y(x). The
function f(x, y, z) = 0 describes a surface in the x y z volume. In special cases, we can
solve for z to get z = z(x, y) to describe the surface in the x y z volume.
Example 3.2
Analyze the surface described by f(x, y, z) = z
2
x
2
y
2
= 0.
Here, we can solve for z exactly to get
z =
_
x
2
+ y
2
. (3.21)
This surface is plotted in Fig. 3.11. We can also get three two-dimensional projections of this surface
in the x y plane, the y z plane, and the x z plane. Orthographic projections of this surface are
plotted in Fig. 3.12.
-2
-1
0
1
2
x
-2
0
1
2
y
-2
-1
0
1
2
z
-1
Figure 3.11: The surface z
2
x
2
y
2
= 0.
-1 0 1 2
-2
-1
0
1
2
x
y
-2 -1 0 1 2
-2
-1
0
1
2
-2 -1 0 1 2
-2
-1
0
1
2
x y
z z
iso-z contours
iso-y
contours
iso-x
contours
-2
Figure 3.12: Contours of constant x, y and z in orthographic projection planes of the surface
z
2
x
2
y
2
= 0.
For the x y plane we consider
z
o
=
_
x
2
+ y
2
. (3.22)
for various values of z
o
. This yields a family of circles in this plane. For the y z plane, we consider
x
o
=
_
z
2
y
2
. (3.23)
This gives a family of hyperbolas. For real x
o
, we require z
2
y
2
. For the z x plane, we consider
y
o
=
_
z
2
x
2
. (3.24)
This gives a similar family of hyperbolas. For real y
o
, we require z
2
x
2
.
Similarly, the ideal gas equation P(v, T) = RT/v describes a surface in the P v T
volume. A surface for air is shown in Fig. 3.13. Often, it is easier to understand the behavior
100
200
300
400
500
T (K)
1
2
3
4
5
v (m /kg)
0
100
200
P (kPa)
3
Figure 3.13: Thermodynamic surface for air modeled as an ideal gas.
of the thermodynamic surfaces by projection into various thermodynamic planes and plotting
various iso-contours. Let us do this for an ideal gas.
isobars:
Consider curves in the T v plane on which P is constant. Thus, for the ideal
gas, we consider
T =
_
P
R
_
. .
slope
v. (3.25)
If we insist that P is constant, this gives the equation of an isobar in the T v
plane. Moreover, for the ideal gas, we see that in the T v plane isobars are
straight lines with slope P/R. The slope is always positive since P > 0 and
R > 0. So if the pressure is high, the slope is positive and steep. If the pressure
is low, the slope is positive and shallow.
Consider curves in the P v plane in which P is constant. Thus, we consider
P = constant, (3.26)
which are straight horizontal lines in the P v plane.
Consider curves in the P T plane in which P is a constant. Thus, we consider
P = constant, (3.27)
which are straight horizontal lines in the P v plane.
Isobars in various planes are shown in Fig. 3.14.
v
v
T
T
P P
P
high
P
low
P
high
P
low
P
high
P
low
Figure 3.14: Isobars for an ideal gas in T v, P v, and P T planes.
isotherms
Consider curves in the T v plane on which T is constant. Thus, for the ideal
gas, we have
T = constant (3.28)
These are straight horizontal lines in the T v plane.
Consider curves in the P v plane on which T is a constant. Thus, for the ideal
gas, we have
P = (RT)
1
v
. (3.29)
These are hyperbolas in the P v plane.
Consider curves in the P T plane on which T is a constant. Thus, for the ideal
gas, we have
T = constant. (3.30)
These are straight vertical lines in the P T plane.
Isotherms in various planes are shown in Fig. 3.15.
v
v
T
T
P P
T
high
T
low
T
high
T
low
T
high
T
low
Figure 3.15: Isotherms for an ideal gas in T v, P v, and P T planes.
isochores
Consider curves in the T v plane on which v is constant. Thus, for the ideal
gas, we have
v = constant (3.31)
These are straight vertical lines in the T v plane.
Consider curves in the P v plane on which v is a constant. Thus, for the ideal
gas, we have
v = constant. (3.32)
These are straight vertical lines in the P v plane.
Consider curves in the P T plane on which v is a constant. Thus, for the ideal
gas, we have
P =
_
R
v
_
. .
constant
T. (3.33)
These are straight lines in the P T plane with slope R/v. Since R > 0 and
v > 0, the slope is always positive. For large v, the slope is shallow. For small v,
the slope is steep.
Isochores in various planes are shown in Fig. 3.16.
Example 3.3
Given air in a cylinder with stops and a frictionless piston with area A = 0.2 m
2
, stop height of
1 m, and total height of 2 m, at initial state P
1
= 200 kPa and T
1
= 500

C with cooling, nd
the temperature when the piston reaches the stops, and
the pressure if the cooling continues to T = 20

C.
v
v
T
T
P P v
high
v
low
v
high
v
low
v
high
v
low
Figure 3.16: Isochores for an ideal gas in T v, P v, and P T planes.
1 m
1 m
air
A = 0. 2 m
2
P
atm
A
P A
free body diagram
at the initial state
Figure 3.17: Sketch for example problem of cooling air.
The initial state along with a free body diagram is sketched in Fig. 3.17.
We have three distinct states:
state 1: initial state
state 2: piston reaches the stops
state 3: nal state, where T = 20
C.
At the initial state, the total volume is
V
1
= A((1 m) + (1 m)) = (0.2 m
2
)(2 m) = 0.4 m
3
. (3.34)
We also know that P
1
= 200 kPa. For use of the ideal gas law, we must use absolute temperature. So
T
1
= 500 + 273.15 = 773.15 K. (3.35)
Now, use the ideal gas law to get v
1
:
v
1
=
RT
1
P
1
=
_
0.287
kJ
kg K
_
(773.15 K)
200 kPa
= 1.10947
m
3
kg
. (3.36)
Let us check the units:
kJ
kg K
K
kPa
kPa
kN
m
2
kN m
kJ
=
m
3
kg
.
As long as we employ kJ for energy and kPa for pressure, we will have few problem with units.
Now, the mass of the air, m, is constant in this problem. Since we have V
1
and v
1
, we can get m:
m =
V
1
v
1
=
0.4 m
3
1.10947
m
3
kg
= 0.360532 kg. (3.37)
Now, as long as the piston does not touch the stops, it will be in a force balance giving P
atm
A = PA.
So our atmosphere must be at P
atm
= 200 kPa. As the air cools, its temperature will go down. Since
Pv = RT, as temperature goes down with constant P, we expect the volume to decrease. Just when
the piston hits the stops, the stops still exert no force on the piston, so
P
2
= P
1
= 200 kPa. (3.38)
Now
V
2
= A(1 m) = (0.2 m
2
)(1 m) = 0.2 m
3
. (3.39)
So
v
2
=
V
2
m
=
0.2 m
3
0.360532 kg
= 0.554735
m
3
kg
. (3.40)
Use the ideal gas law to get T
2
:
T
2
=
P
2
v
2
R
=
(200 kPa)
_
0.55735
m
3
kg
_
0.287
kJ
kg K
= 386.575 K. (3.41)
Now, after the piston reaches the stops, the volume is constant. So the process from 2 to 3 is
isochoric, and
V
3
= V
2
= 0.2 m
3
. (3.42)
Thus
v
3
= v
2
= 0.554735
m
3
kg
. (3.43)
So
P
3
=
RT
3
v
3
=
_
0.287
kJ
kg K
_
(20 + 273.15) K)
0.554735
m
3
kg
= 151.665 kPa. (3.44)
We generate Table 3.1 to summarize the problem. It is usually useful to include sketches of the
process in the various thermodynamic planes. This process is sketched in each of the relevant planes
in Fig. 3.18.
Table 3.1: Numerical values for ideal gas cooling example
variable units state 1 state 2 state 3
T K 773.15 386.575 293.15
P kPa 200 200 151.665
v
m
3
kg
1.10947 0.554735 0.554735
V m
3
0.4 0.2 0.2
v (m
3
/kg) v (m
3
/kg) T (K)
T (K) P (kPa) P (kPa)
200 - 200 -
151.6 - 151.6 -
773 -
386 -
293 -
1
2
3
1 2
3
1 2
3
0.555 1.11
0.555 1.11 386 773
T
=
2
9
3
K
T
=
3
8
6
K
T

=

7
7
3
K
v

=

1
.
1
1

m
3 /
k
g

P

=

2
0
0

k
P
a
P
=
1
5
1
.6
k
P
a
v

=

0
.
5
5
5

m
3
/
k
g
Figure 3.18: Sketch of T v, P v, and P T planes for problem of air cooling example
problem.
3.4.2 Non-ideal thermal equations of state
The ideal gas law is not a good predictor of the P v T behavior of gases when
the gas has high enough density that molecular interaction forces become large and
the molecules occupy a signicant portion of the volume; this happens near the vapor
dome typically, or
the temperature is high enough to induce molecular dissociation (e.g. N
2
+ N
2

2N + N
2
)
One alternative is a corrected thermal equation of state.
3.4.2.1 van der Waals
For the van der Waals
9
equation of state, one has
P =
RT
v b

a
v
2
, (3.45)
9
J. D. van der Waals, 1873, Over de Continuiteit van den Gas -en Vloeistoftoestand, Ph.D. Dissertation,
U. Leiden; see also J. D. van der Waals, 1910, Nobel Lecture.
with
a =
27
64
R
2
T
2
c
P
c
, b =
1
8
R
T
c
P
c
. (3.46)
A depiction of van der Waals is given in Fig. 3.19.
Figure 3.19: Johannes Diderik van der Waals (1837-1923), Dutch physicist
and Nobel laureate who developed a corrected state equation. Image from
http://en.wikipedia.org/wiki/Johannes Diderik van der Waals.
3.4.2.2 Redlich-Kwong
For the Redlich-Kwong
10
equation of state, one has
P =
RT
v b

a
v(v + b)
T
, (3.47)
with
a = (0.42748)
R
2
T
5/2
c
P
c
, b = (0.08664)
RT
c
P
c
. (3.48)
3.4.3 Compressibility factor
In some cases, more detail is needed to capture the behavior of the gas, especially near the
vapor dome. Another commonly used approach to capturing this behavior is to dene the
Compressibility factor: the deviation from ideality of a gas as measured by
Z =
Pv
RT
. (3.49)
P (MPa)
2
1
0
1
10
300 K
1
3
0

K
2
0
0
K
saturated liquid
saturated vapor
critical
point
Ideal gas, Z = 1
Z

=

P

v

/

R

T
Figure 3.20: Sketch of compressibility chart for N
2
.
For ideal gases, Pv = RT, so Z = 1. Experiments show the behavior of real gases, and this
can be presented in graphical form, as shown for N
2
in Fig. 3.20. Note
for all T, Z 1 as P 0. Thus, one has ideal gas behavior at low pressure
for T > 300 K, Z 1 for P < 10 MPa.
Hold at P = 4 MPa and decrease temperature from 300 K; we see Z decrease below
unity. Now
Z =
Pv
RT
=
P
RT
, =
P
ZRT
. (3.50)
Since Z < 1, the density is higher than we would nd for an ideal gas with Z = 1.
Thus, in this region, there is an attractive force between molecules.
For P > 30 MPa, we nd Z > 1. Thus, a repulsive force exists in this regime. The
forces are complicated.
Note that generalized compressibility charts have been developed for general gases. These
are based on the so-called reduced pressures and temperatures, P
r
and T
r
, where
P
r
=
P
P
c
, T
r
=
T
T
c
. (3.51)
10
O. Redlich and J. N. S. Kwong, 1949, On the thermodynamics of solutions. V. an equation of state.
fugacities of gaseous solutions, Chemical Reviews, 44(1): 233-244.
The reduced pressure and temperature are dimensionless. Values with the c subscript are the
critical constants for the individual gases. Appendix D of BS gives generalized compressibility
charts.
3.4.4 Tabular thermal equations of state
Often equations are too inaccurate for engineering purposes. This is generally because we
may be interested in behavior under a vapor dome. Consider that the surface for steam is
well represented by that shown in Fig. 3.21.
Figure 3.21: P v T surface for H
2
O, showing solid, liquid, and vapor phases.
In such cases, one should use tables to nd a third property, given two independent
properties. We can say that the thermal equation of state is actually embodied in the
tabular data.
We lay down some rules of thumb for this class:
If steam, use the tables.
If air or most other gas, use the ideal gas law, but check if the pressure is high or
the properties are near the vapor dome, in which case use compressibility charts or
non-ideal state equations.
Let us look at how the tables are organized.
3.4.4.1 Saturated liquid-vapor water, temperature tables
For water, the most important table is the saturated steam table. One should go to such
tables rst. If the water is a two-phase mixture, tables of this type must be used as the
equation of state. Recall, for two-phase mixtures, pressure and temperature are not inde-
pendent thermodynamic variables. Two properties still determine the state, but quality x is
now important. So for two-phase mixtures we allow
T = T(v, x),
P = P(v, x), or
v = v(T, x),
for example. But P = P(T, v) as for ideal gases.
Consider the structure of saturation tables, as shown in Table 3.2, extracted from BSs
Table B.1.1. Data from the steam tables is sketched in Fig. 3.22. We have the notation:
Specic Volume,
m
3
kg
Temp. Press. Sat. Liquid Evap. Sat. Vapor
C kPa v
f
v
fg
v
g
0.01 0.6113 0.001000 206.131 206.132
5 0.8721 0.001000 147.117 147.118
10 1.2276 0.001000 106.376 106.377
15 1.705 0.001001 77.924 77.925
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
35 5.628 0.001006 25.2148 25.2158
40 7.384 0.001008 19.5219 19.5229
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
374.1 22089 0.003155 0 0.00315
Table 3.2: Saturated liquid-vapor water tables, temperature entry, from BS, Table B.1.1.
f: saturated liquid,
g: saturated vapor,
v
f
: specic volume of saturated liquid, and
v
g
: specic volume of saturated vapor.
v ( m
3
/kg )
T ( C )
compressed
liquid
superheated
vapor
15 C
v = 0.001000 m
3
/kg v = 77.925 m
3
/kg
f
g
v = 77.924 m
3
/kg
fg
f
g
Figure 3.22: Vapor dome for H
2
O with data for v
f
, v
g
, and v
fg
at T = 15
C.
Note for liquid-vapor mixtures, this table begins at the triple point temperature 0.01

C and
ends at the critical temperature 374.1

C. At P = P
c
and T = T
c
, we have v
f
= v
g
. Note
that
v
f
constant
v
g
decreases with increasing T
We dene v
fg
as
v
fg
v
g
v
f
. (3.52)
Recall the quality x is
x =
m
vap
m
total
.
Consider a mass of uid m in total volume V . We must have
V = V
liq
+ V
vap
, (3.53)
m = m
liq
+ m
vap
. (3.54)
Now, use the denition of specic volume and analyze to get
mv = m
liq
v
f
+ m
vap
v
g
, (3.55)
v =
m
liq
m
v
f
+
m
vap
m
v
g
, (3.56)
v =
mm
vap
m
v
f
+
m
vap
m
v
g
, (3.57)
v = (1 x)v
f
+ xv
g
, (3.58)
v = v
f
+ x(v
g
v
f
)
. .
=v
fg
. (3.59)
We get the nal important results:
v = v
f
+ xv
fg
. (3.60)
x =
v v
f
v
fg
. (3.61)
3.4.4.2 Saturated liquid-vapor water, pressure tables
Sometimes we are given the pressure of the mixture, and a saturation table based on the
pressure is more useful. An example of a portion of such a table is shown in Table 3.3.
Specic Volume,
m
3
kg
Press. Temp. Sat. Liquid Evap. Sat. Vapor
kPa

C v
f
v
fg
v
g
0.6113 0.01 0.001000 206.131 206.132
1.0 6.98 0.001000 129.20702 129.20802
1.5 13.03 0.001001 87.97913 87.98013
2.0 17.50 0.001001 67.00285 67.00385
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
22089 374.1 0.003155 0 0.00315
Table 3.3: Saturated water tables, pressure entry from BS, Table B.1.2.
Example 3.4
Given a vessel with V = 0.4 m
3
lled with m = 2 kg of H
2
O at P = 600 kPa, nd
the volume and mass of liquid, and
the volume and mass of vapor.
liquid water
water vapor
P=600 kPa
m = 2 kg
V = 0.4 m
Figure 3.23: Schematic for liquid-vapor mixture example problem.
The problem is sketched in Fig. 3.23. While the problem statement suggests we have a two-phase
mixture, that is not certain until one examines the tables. First, calculate the specic volume of the
water:
v =
V
m
=
0.4 m
3
2 kg
= 0.2
m
3
kg
. (3.62)
Next go to the saturated water tables with pressure entry to see if the water is a two-phase mixture.
We nd at P = 600 kPa that
v
f
= 0.001101
m
3
kg
, v
g
= 0.31567
m
3
kg
. (3.63)
Now, for our mixture, we see that v
f
< v < v
g
, so we have a two-phase mixture. Now, apply Eq. (3.61)
to nd the quality.
x =
v v
f
v
fg
=
v v
f
v
g
v
f
=
0.2
m
3
kg
0.001101
m
3
kg
0.31567
m
3
kg
0.001101
m
3
kg
= 0.632291. (3.64)
Now, from Eq. (3.1), we have x = m
vap
/m
total
, so
m
vap
= xm
tot
= 0.632291(2 kg) = 1.26458 kg. (3.65)
Now, for the liquid mass we have
m
liq
= m
total
m
vap
= (2 kg) (1.26458 kg) = 0.735419 kg. (3.66)
Most of the mass is vapor, but the fraction that is liquid is large.
Now, let us calculate the volumes.
V
vap
= m
vap
v
g
= (1.26458 kg)
_
0.31567
m
3
kg
_
= 0.39919 m
3
, (3.67)
V
liq
= m
liq
v
f
= (0.735419 kg)
_
0.001101
m
3
kg
_
= 0.000809696 m
3
. (3.68)
The volume is nearly entirely vapor.
3.4.4.3 Superheated water tables
The superheat regime is topologically similar to an ideal gas. For a superheated vapor,
the quality x is meaningless, and we can once again allow pressure and temperature to be
independent. Thus, we can have v = v(T, P). And the tables are in fact structured to give
v(T, P) most directly. An example of a portion of such a table is shown in Table 3.4. This
Temp. v u h s
C
m
3
kg
kJ
kg
kJ
kg
kJ
kg K
P = 10 kPa (45.81

C)
Sat. 14.67355 2437.89 2584.63 8.1501
50 14.86920 2443.87 2592.56 8.1749
100 17.19561 2515.50 2687.46 8.4479
150 19.51251 2587.86 2782.99 8.6881
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Table 3.4: Superheated water tables, from BS, Table B.1.3.
portion of the superheated tables focuses on a single isobar, P = 10 kPa. At that pressure,
the saturation temperature is 45.81

C, indicated in parentheses. As long as T > 45.81

C,
we can use this table for P = 10 kPa water. And for various values of T > 45.81

C, we
nd other properties, such as specic volume v, and properties we have not yet focused on,
internal energy u, enthalpy h, and entropy s.
3.4.4.4 Compressed liquid water tables
Liquids truly have properties which vary with both T and P. To capture such variation,
we can use compressed liquid tables as an equation of state. An example for water is given
in Table 3.5. If compressed liquid tables do not exist, it is usually safe enough to assume
properties are those for x = 0 saturated liquid at the appropriate temperature.
3.4.4.5 Saturated water, solid-vapor
Other types of saturation can exist. For example, below the triple point temperature, one
can have solid water in equilibrium with water vapor. The process where ice transforms
directly to water vapor is known as sublimation. Saturation tables for ice-vapor equilibrium
exist as well. For example, consider the structure of saturation tables, as shown in Table
3.6, extracted from BSs Table B.1.5.
Temp. v u h s
C
m
3
kg
kJ
kg
kJ
kg
kJ
kg K
P = 500 kPa (151.86

C)
Sat. 0.001093 639.66 640.21 1.8606
0.01 0.000999 0.01 0.51 0.0000
20 0.001002 83.91 84.41 0.2965
40 0.001008 167.47 167.98 0.5722
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Table 3.5: Compressed liquid water tables, from BS, Table B.1.4.
Specic Volume,
m
3
kg
Temp. Press. Sat. Solid Evap. Sat. Vapor
C kPa v
i
v
ig
v
g
0.01 0.6113 0.0010908 206.152 206.153
0 0.6108 0.0010908 206.314 206.315
-2 0.5177 0.0010905 241.662 241.663
-4 0.4376 0.0010901 283.798 283.799
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Table 3.6: Saturated solid-vapor water tables, temperature entry, from BS, Table B.1.5.
3.4.4.6 Tables for other materials
For many materials similar tables exist, e.g., ammonia, NH
3
. Consider the ammonia satura-
tion tables, as shown in Table 3.7, extracted from BSs Table B.2.1. One also has tables for
Specic Volume,
m
3
kg
C kPa v
f
v
fg
v
g
-50 40.9 0.001424 2.62557 2.62700
-45 54.5 0.001437 2.00489 2.00632
-40 71.7 0.001450 1.55111 1.55256
-35 93.2 0.001463 1.21466 1.21613
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
132.3 11333.2 0.004255 0 0.004255
Table 3.7: Saturated liquid-vapor ammonia tables, temperature entry, from BS, Table B.2.1.
superheated ammonia vapor. An example of a portion of such a table is shown in Table 3.8.
Other tables in BS, include those for carbon dioxide, CO
2
, a modern refrigerant, R-410a,
11
Temp. v u h s
C
m
3
kg
kJ
kg
kJ
kg
kJ
kg K
P = 50 kPa (46.53

C)
Sat. 2.1752 1269.6 1378.3 6.0839
-30 2.3448 1296.2 1413.4 6.2333
-20 2.4463 1312.3 1434.6 6.3187
-10 2.5471 1328.4 1455.7 6.4006
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Table 3.8: Superheated ammonia tables, from BS, Table B.2.2.
another common refrigerant, R-134a,
12
diatomic nitrogen, N
2
, and methane, CH
4
.
3.4.4.7 Linear interpolation of tabular data
interpolation is often required when exact values are not tabulated.
11
a common cooling uid invented in 1991, a near-azeotropic mixture of diuoromethane and pentauo-
roethane.
12
a cooling uid which became common in the 1990s, 1,1,1,2-tetrauoroethane.
in this course we will primarily use linear interpolations.
use extrapolations only if there is no other choice.
occasionally double interpolations will be necessary.
3.4.4.7.1 Single interpolation The most common interpolation is the single interpo-
lation of variables. We give an example here.
Example 3.5
Given steam at T = 36.7

C, with v = 10 m
3
/kg, nd the pressure and the quality if a two-phase
mixture.
A wise rst step is to go to the saturated tables. We check Table B.1.1 from BS and nd there are
no values at T = 36.7

C. So we must create our own personal steam table values at this temperature,
just to determine if where we are on the thermodynamic surface. We list the important part of the
saturated water liquid-vapor tables in Table 3.9.
Specic Volume,
m
3
kg
C kPa v
f
v
fg
v
g
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
35 5.628 0.001006 25.2148 25.2158
36.7 ? ? ? ?
40 7.384 0.001008 19.5219 19.5229
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Table 3.9: Relevant portion of saturated liquid-vapor water tables, temperature entry, from
BS, Table B.1.1.
We seek to get appropriate values for P, v
f
, v
fg
, and v
g
at T = 36.7

C. Let us nd P rst. The
essence of linear interpolation is to t known data to a straight line, then use the formula of that line to
predict intermediate values of variables of interest. We know values of P at T = 35

C and T = 40

C.
In fact we have two points: (T
1
, P
1
) = (35

C, 5.628 kPa), and (T
2
, P
2
) = (40

C, 7.384 kPa). This lets
us t a line using the familiar point-slope formula:
P P
1
=
_
P
2
P
1
T
2
T
1
_
. .
slope
(T T
1
). (3.69)
We could have used the other point. Note when T = T
1
, that P = P
1
. Also, when T = T
2
, P = P
2
.
Substituting numbers, we get
P =
_
(7.384 kPa) (5.628 kPa)
(40

C) (35

C)
_
(T (35

C)) + (5.628 kPa), (3.70)
=
_
0.3512
kPa
C
_
(T (35

C)) + (5.628 kPa), (3.71)
=
_
0.3512
kPa
C
_
((36.7

C) (35

C)) + (5.628 kPa), (3.72)
= 6.225 kPa. (3.73)
The interpolation is sketched is sketched in Fig. 3.24.
T (C )
P (kPa)
35 36.7 40
7.384
6.225
5.628
1
2
Figure 3.24: Sketch of linear interpolation to nd P when T = 36.7

C, v = 10 m
3
/kg for
water.
Now, we need to interpolate for v
f
and v
g
as well. Let us apply the same technique. For v
f
, we
have
v
f

_
0.001006
m
3
kg
_
=
_
_
_
0.001008
m
3
kg
_
_
0.001006
m
3
kg
_
(40

C) (35

C)
_
_
(T (35

C)). (3.74)
When T = 36.7

C, we get
v
f
= 0.00100668
m
3
kg
. (3.75)
For v
g
, we get
v
g
_
25.2158
m
3
kg
_
=
_
_
_
19.5229
m
3
kg
_
_
25.2158
m
3
kg
_
(40

C) (35

C)
_
_
(T (35

C)). (3.76)
When T = 36.7

C, we get
v
g
= 23.2802
m
3
kg
. (3.77)
Knowing v
f
and v
g
, we do not need to interpolate for v
fg
. We can simply use the denition:
v
fg
= v
g
v
f
=
_
23.2802
m
3
kg
_
_
0.00100668
m
3
kg
_
= 23.2792
m
3
kg
. (3.78)
Now, v = 10 m
3
/kg. Since at T = 36.7
C, we have v
f
< v < v
g
, we have a two-phase mixture.
Let us get the quality. From Eq. (3.61), we have
x =
v v
f
v
fg
=
_
10
m
3
kg
_
_
0.00100668
m
3
kg
_
23.2792
m
3
kg
= 0.429525. (3.79)
Thus
x =
m
vap
m
tot
= 0.429525. (3.80)
3.4.4.7.2 Double interpolation Sometimes, we need to do extra linear interpolations.
Say we are given superheated water with v
o
and T
o
and we are asked to nd P
o
. But neither
v
o
nor T
o
are listed in the tables. Then we need to do a multi-step procedure.
Go to the tables and for the given T
o
and v
o
, estimate approximately the value of P
o
by visual examination.
For a nearby value of P = P
1
, get a linear interpolation of the form T = T(v, P
1
). Use
this to get T
1
= T(v
o
, P
1
).
For a dierent nearby value of P = P
2
, get another linear interpolation of the form
T = T(v, P
2
). Use this to get T
2
= T(v
o
, P
2
). We now have two points (T
1
, P
1
) and
(T
2
, P
2
), both valid at v = v
o
.
Use the two points (T
1
, P
1
), (T
2
, P
2
) to develop a third interpolation P = P(T, v
o
).
Estimate P
o
by P
o
= P(T
o
, v
o
).
Example 3.6
Consider m = 1 kg of H
2
O initially at T
1
= 110

C, x
1
= 0.9. The H
2
O is heated until T
2
= 200

C.
As sketched in Fig. 3.25, the H
2
O is conned in a piston-cylinder arrangement, where the piston
is constrained by a linear spring with dP/dv = 40 kPa/m
3
/kg. At the initial state, the spring is
unstretched. Find the nal pressure.
While this problem seems straightforward, there are many challenges. Let us rst consider what
we know about the initial state. Since we have a numerical value for x
1
, we know state 1 is a two-phase
mixture. From the tables, we nd that
P
1
= 143.3 kPa, v
f1
= 0.001052
m
3
kg
, v
g1
= 1.2101
m
3
kg
, v
fg1
= 1.20909
m
3
kg
. (3.81)
P A
P
atm
A
free body diagram
water
F
s
= k ( y - y )
Figure 3.25: Sketch of piston-cylinder arrangement.
We can then calculate v
1
for the mixture:
v
1
= v
f1
+ x
1
v
fg1
=
_
0.001052
m
3
kg
_
+ (0.9)
_
1.20909
m
3
kg
_
= 1.08923
m
3
kg
. (3.82)
We now know everything we need about state 1.
At state 2, we only know one intensive thermodynamic property, the temperature, T
2
= 200

C.
To get a second, and thus dene the nal state, we will need to bring in information about the process.
Now, we will need to consider a force balance on the piston. Newtons second law for the piston says
m
piston
d
2
y
dt
2
=
F
y
. (3.83)
From our free body diagram, we note three major forces:
force due to the interior pressure from the water,
force due to the exterior pressure from the atmosphere,
force due to the linear spring, which we call F
s
= k(y y
1
) where k is the spring constant, y is the
position of the piston, and y
1
is the initial position of the piston. Note that F
s
= 0 when y = y
1
.
We write this as
m
piston
d
2
y
dt
2
= PAP
atm
Ak(y y
1
). (3.84)
Now, in classical thermodynamics, we make the assumption that the inertia of the piston is so small
that we can neglect its eect. We are really requiring that a force balance exist for all time. Thus,
even though the piston will move, and perhaps accelerate, its acceleration will be so small that it can
be neglected relative to the forces in play. We thus take
m
piston
d
2
y
dt
2
|F
s
|, |P
atm
A|, |PA|. (3.85)
With this assumption, we have
0 PAP
atm
A k(y y
1
). (3.86)
Solve for P, the water pressure, to get
P = P
atm
+
k
A
(y y
1
). (3.87)
Now, V = Ay and V
1
= Ay
1
, so we can say
P = P
atm
+
k
A
2
(V V
1
). (3.88)
Let us use the fact that V = mv and V
1
= mv
1
to rewrite as
P = P
atm
+
km
A
2
(v v
1
). (3.89)
This equation is highlighted because it provides an algebraic relationship between two intensive thermo-
dynamic properties, P and v, and such a tactic will be useful for many future problems. Using numbers
from our problem, with dP/dv = km/A
2
, we can say
P = (143.3 kPa) +
_
40
kPa
m
3
kg
_
_
v
_
1.08923
m
3
kg
__
, (3.90)
P = (99.7308 kPa) +
_
40
kPa
m
3
kg
_
v.
. .
linear spring rule
(3.91)
Now, at state 1, we have V = V
1
and P = P
1
= 143.3 kPa, so we must have P
atm
= 143.3 kPa for this
problem.
Let us now consider the possibilities for state 2. We are constrained to be on the line in P v space
given by our force balance, Eq. (3.91). We are also constrained to be on the T = 200

C isotherm,
which is also a curve in P v space. So let us consider the P v plane, as sketched in Fig. 3.26. The
v
P
T = 110C
T =200C
v
P
T =110C
T =200C
km/A ~ large
2
km/A ~ small
2
1
2
1
2
1
2
1
2
Figure 3.26: Sketch of P v plane for piston-cylinder-linear spring problem for water.
isotherms for T
1
= 110

C and T
2
= 200

C are set in both parts of Fig. 3.26. Since both T
1
and T
2
are
well below T
c
, both isotherms pierce the vapor dome. Our nal state has a line in P v space from the
force balance intersecting the state 2 isotherm. There are two distinct possibilities for the nal state:
for a sti spring, i.e. large km/A
2
, our line will intersect the isotherm within the vapor dome, or
for a loose spring, i.e. small km/A
2
, our line will intersect the isotherm in the superheated vapor
region.
Let us consider the rst possibility: state 2 is under the vapor dome. If that is the case, then the
tables tell us that P
2
= 1553.8 kPa. At this pressure, Eq. (3.91) gives us v = (1553.8 99.7308)/40 =
36.3517 m
3
/kg. However, at this pressure v
g
= 0.12736 m
3
/kg. Since we just found v > v
g
, our
assumption that the nal state was under the dome must be incorrect!
Therefore, let us go to the more dicult case posed by the second possibility: state 2 is a superheated
vapor. In general the intersection of the straight line with the isotherm is dicult. We can use
linearization to assist us. Let us choose a small region of the tables, and locally t the 200

C isotherm
to a straight line. This will give us a second independent equation in P v space. We will then solve
two equations in two unknowns for the nal state.
Our initial pressure, P
1
= 143.3 kPa lies between 100 kPa and 200 kPa. We have values from the
superheat tables at these pressures for v at 200

C. So, we approximate the isotherm by the line
P (100 kPa) =
(200 kPa) (100 kPa)
_
1.08034
m
3
kg
_
_
2.17226
m
3
kg
_
_
v
_
2.17226
m
3
kg
__
. (3.92)
P = (298.939 kPa)
_
91.5818
kPa
m
3
kg
_
v.
. .
linear approximation of isotherm from tabular thermal EOS
(3.93)
We simultaneously solve the two linear equations, Eqs. (3.91, 3.93), and get the unique solution
P
2
= 160.289 kPa, v
2
= 1.151395
m
3
kg
. (3.94)
Since we found 100 kPa < P
2
= 160.289 kPa < 200 kPa, we made a good assumption on the nal
pressure, and our interpolation values from the tables are acceptable. Lastly, we sketch the process in
the T v and P T planes in Fig. 3.27.
v
T
T
P
1
2
1
2
P1
=
1
4
3
.
3

k
P
a
P
2

=

1
6
0
.
2
8
9

k
P
a
liquid
solid
vapor
Figure 3.27: Sketch of T v and P T planes for piston-cylinder-linear spring problem for
water.

Pure Substances

Uploaded by

Copyright:

Available Formats

Pure Substances

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Pure Substances

Uploaded by

Copyright:

Available Formats

Chapter 3

Properties of a pure substance

E. Clapeyron, 1834, Memoire sur la puissance motrice de la chaleur, Journal de l

You might also like