Transcript Document 7360896

Lecture 16. The van der Waals Gas (Ch. 5)
Nobel 1910
The simplest model of a liquid-gas phase
transition - the van der Waals model of
“real” gases – grasps some essential
features of this phase transformation.
(Note that there is no such transformation
in the ideal gas model). This will be our
attempt
to
take
intermolecular
interactions into account.
Outline
• Thermodynamics of vdW gas
• critical point
• liquid-gas phase transition
van der Waals gas = ideal gas + interaction
short-distance
repulsion
The main reason for the transformation of
gas into liquid at decreasing T and (or)
increasing P - interaction between the
molecules.
4
U(r)
3
Energy
2
1
Two ingredients of the model:
0
the weak long-range attraction:
-1
-2
N
N 2a
U   N   a  
V
V
-3
1.5
2.0
2.5
3.0
r
3.5
4.0
distance
r=
U vdW  U IG
long-distance
Lennard-Jones
12
6
attraction
 
 
U r      
r
r
-
N 2a
N 2a

 Peff  P  2
V
V
the constant a is a measure of the longrange attraction
the strong short-range repulsion: the molecules are rigid: P   as soon as the
molecules “touch” each other.
Veff  V  Nb
The vdW equation
of state
-
the constant b (~ 43/3) is a measure of the short-range
repulsion, the “excluded volume” per particle

N 2a 
 P  2 V  Nb  NkBT
V 

NkBT
N 2a
P
 2
V  Nb V
The van der Waals Parameters
b – roughly the volume of a molecule, (3.5·10-29 – 1.7 ·10-28) m3 ~(few Å)3
a – varies a lot [~ (8·10-51 – 3 ·10-48) J · m3] depending on the intermolecular
interactions (strongest – between polar molecules, weakest – for inert gases).
a’
b’
(J. m3/mol2)
(x10-5 m3/mol)
Air
.1358
3.64
Carbon Dioxide (CO2)
.3643
4.27
Nitrogen (N2)
.1361
3.85
Hydrogen (H2)
.0247
2.65
Water (H2O)
.5507
3.04
Ammonia (NH3)
.4233
Helium (He)
Freon (CCl2F2)
Substance
Pc
(MPa)
3.77 2
NA
7.39
Tc
(K)
133 K
a 304.2
a
K ' N Ab  b '
3.39
126.2 K
3.73
 J  m3 
a' 
3
mol 2 


J

m
1.30
33.2 K

a 

2
2 
22.09 molecule
647.3 K
molecules 
 
 NA

11.28
406 K
mole 

.00341
2.34
0.23
5.2 K
1.078
9.98
4.12
385 K

N 2a 
When can  P  2 V  Nb  NkBT be reduced to PV  NkBT ?


V


- low densities
Nb  V
Na 

PV  N  k BT 

V


k BT 
Na
V
- high temperatures
(kinetic energy >> interaction energy)
Problem
The vdW constants for N2: NA2a = 0.136 Pa·m6 ·mol-2, NAb = 3.85·10-5 m3 ·mol-1.
How accurate is the assumption that Nitrogen can be considered as an ideal gas at
normal P and T?

N 2a 
 P  2 V  Nb  NkBT
V 

1 mole of N2 at T = 300K occupies V1 mol  RT/P  2.5 ·10-2 m3 ·mol-1
NAb = 3.9·10-5 m3 ·mol-1
NAb / V1 mol ~1.6%
NA2a / V2 = 0.135 Pa·m6 ·mol-2 /(2.5 ·10-2 m3 ·mol-1) 2 = 216 Pa
NA2a / V2P = 0.2%
Entropy of monatomic van der Waals gas
1
V N  2 mU 
Multiplicity of ideal gas: IG  N ,V ,U  


N !  3 N / 2 !  h 2 
van der Waals gas:
Veff  V  Nb
U vdW  U IG
N 2a

V
Multiplicity of monatomic van der Waals gas:

1 V  Nb   2 m U 

 vdW  N ,V ,U  
N !  3N / 2 ! 
h2

N
N 2a
V
 


3N / 2
2p
Entropy of monatomic van der Waals gas:
 V  Nb  4 m U  N 2 a 3/ 2  5 


V

SvdW ( N ,V ,U )  Nk B ln 


 2

N
3
h
N
2

 


 

Is it correct?
3N / 2
2p
Equations of state (vdW)
 V  Nb  4 m U  N 2 a 3/ 2  5 


V

SvdW ( N ,V ,U )  Nk B ln 


 2

N   2
  N  3h


Equations of state:
1  S 


T  U V , N
 S 
P  T 

 V U , N
U vdW
3
N 2a
 Nk BT 
2
V
NkBT
N 2a
P
 2
V  Nb V
 S 

 N U ,V
  T  
V  Nb  2 mk BT 3/ 2  2 Na
Nb
  k BT ln 


 
2
V  Nb
  V
 N  h
(take home exercise: derive these results)
S, F, and G for the monatomic van der Waals gas
S
C
f
dT
NkB
 S 
 S 
 P 
dS  

dV
 dT  
 dV  V dT  
 dV  NkB
T
2
T
V  Nb
 T V
 V T
 T V
(see Pr. 5.12)
SvdW
S vdW 
f
Nk B ln T  Nk B ln V  Nb   const
2
3/ 2
  V  Nb  2  m
 Nb 
  5 
 Sideal  NkB ln 1    N k B ln 
 2 k BT    
V
N


 h
  2 
 
F
FvdW
3
N 2a
 U  TS  NkBT 
 N k BT
2
V
- the same “volume” in the
momentum space, smaller
accessible volume in the
coordinate space.
3/ 2
  V  Nb  2  m
  5 
 2 k BT    
ln 
 h
  2 
  N
3/ 2
 V  Nb  2  m
aN 2
 
  N k BT ln 
 2 k BT    NkBT 
V
 h
 
 N
 Nb  N a
 Fideal  N k BT ln 1 

V  V

2
FvdW
G
GvdW
PvdW
Nk BT aN 2
 F 
 
 2
 

V
V

Nb
V

T
3/ 2
2
 V  Nb   2 m
N 2b
  2aN
 F  PV   Nk BT ln 
 k BT
 N
 2 k BT   
N
V
V  Nb
 h
 

(take home exercise: derive these results)
Quasi-static Processes in a vdW Gas (low n, high T)
U  Q  W
dU  TdS  PdV
isochoric ( V = const )
P
W12  0
2
T2
T1
1
V1,2
V
3
Q12  Nk B T2  T1   0  CV T 
2
same as ideal gas
dU  Q12
isobaric
P
( P = const )
2
W12    P(V , T )dV   PV2  V1   0
2
1
V1
1
T2
T1
V2
Q12  CP T
V
 U 
 dV 
CP  CV  
  P 

 V T
 dT  P
dU  W12  Q12
complicate
Isothermal Process in a vdW Gas (low n, high T)
isothermal ( T = const ) :
P
PV= NkBT
V2
W
V2
 U T
N 2a N 2a


V1
V2
V1
V
W12    P(V , T )dV
V1
 Nk BT
N 2a 
  
 2  dV
V  Nb V 
V1 
V2
Wi-f > 0 if Vi >Vf (compression)
Wi-f < 0 if Vi <Vf (expansion)
V2  Nb  N 2 a N 2 a 
  Nk BT ln



V1  Nb  V1
V2 
Q12  U  W12
V2  Nb
 Nk BT ln
V1  Nb
Isothermal Process in a vdW Gas (low n, high T)
For N2, the vdW coefficients are N2a = 0.138 kJ·liter/mol2 and Nb = 0.0385
liter/mol. Evaluate the work of isothermal and reversible compression of N2
(assuming it is a vdW gas) for n=3 mol, T=310 K, V1 =3.4 liter, V2 =0.17 liter.
Compare this value to that calculated for an ideal gas. Comment on why it is
easier (or harder, depending on your result) to compress a vdW gas relative to an
ideal gas under these conditions.
WvdW T  9mol 2  0.138 kJ  l2 
1
1 
J
 0.17l  3mol  0.0385l / mol 

 310 K ln 
  3mol  8.3

mol  3.4l 0.17l 
K  mol
3
.
4
l

3
mol

0
.
0385
l
/
mol


 6.94kJ  31.64kJ  24.7kJ
 
Wideal T   NkBT ln  V2   3mol  8.3 J  310K  ln  0.17l   23.12kJ
mol  K
 3.4l 
 V1 
Upot
r
Depending on the interplay between the 1st and 2nd terms, it’s
either harder or easier to compress the vdW gas in comparison
with an ideal gas. If both V1 and V2 >> Nb, the interactions between
molecules are attractive, and WvdW < Wideal However, as in this
problem, if the final volume is comparable to Nb , the work against
repulsive forces at short distances overweighs that of the attractive
forces at large distances. Under these conditions, it is harder to
compress the vdW gas rather than an ideal gas.
Adiabatic Process in a vdW Gas
adiabatic (thermally isolated system)
3/ 2

 V  Nb  2 mk BT   5 

SvdW ( N ,V , T )  Nk B ln 



 
2
N
h
2





 

P
2
1
V2
Q12  0  S  0 dU  W12
V1
 V  Nb  T
T2
T1
V

V2
W12    P(V , T )dV
3/ 2
 const.  T V  Nb 
 1
 const.

N 2a 

P

V

Nb
 const



2 
V 

V1
recall:
2
2

3
N
a
N
a
3
N a

 U  NkB T2  T1   
 Nk BT 

2
V2 
2
V
 V1
2
U vdW
take home exercise: Carnot engine efficiency using vdW gas
Problem (vdW + heat engine)
The working substance in a heat engine is the vdW gas with a known constant b and
a temperature-independent heat capacity cV (the same as for an ideal gas). The gas
goes through the cycle that consists of two isochors (V1 and V2) and two adiabats.
Find the efficiency of the heat engine.
P
e  1
A
D
A-D
B
B-C
 QC
 QH
QH  cV TA  TD 
e  1
QC  cV TB  TC 
TB  TC 
TA  TD 
C
V1
S vdW 
V2
V
The relationship between TA, TB, TC, TD – from
the adiabatic processes B-C and D-A
f
Nk B ln T  Nk B ln V  Nb   const
2
Tf
V f  Nb
f
NkB ln
 NkB ln
0
2
Ti
Vi  Nb
T  TC 
QC
e 1
1 B
QH
TA  TD 
 V  Nb 
TA  1
V2  Nb 

1
SvdW 
f
2
Tf
V f  Nb
f
Nk B ln
 Nk B ln
2
Ti
Vi  Nb
T V  Nb  const
2/ f
 V  Nb 
 TD  1

 V2  Nb 
TA  TD
2/ f
adiabatic process
for the vdW gas
 V  Nb 
1  1

 V2  Nb 
2/ f
 V  Nb 
1  1

 V2  Nb 
 1
The van der Waals Isotherms
NkBT
aN 2
P
 2
V  Nb V
2
3
Nk
T
aN
abN


V 3   Nb  B  V 2 
V
0
P 
P
P

unstable
0
N·b
The Critical Point
The critical point is the unique point where both
(dP/dV)T = 0 and (d2P/dV 2)T = 0 (see Pr. 5.48)
Critical parameters:
1 a
8 a
VC  3Nb PC 
k BTC 
2
27 b
27 b
P
Pˆ 
PC
V
Vˆ 
VC
T
Tˆ 
TC
- in terms of P,T,V normalized by the critical parameters:
 ˆ 3  ˆ 1  8kBTˆ
 P  ˆ 2 V   
3
3
V 

- the materials parameters vanish if we introduce the proper scales.
RTC 8
KC 
  2.67
PCVC 3
- the critical coefficient
substance
H2
He
N2
CO2
H20
KC
3.0
3.1
3.4
3.5
4.5
TC (K)
33.2
5.2
126
304
647
PC (MPa)
1.3
0.23
3.4
7.4
22.1
Problems
For Argon, the critical point occurs at a pressure PC =
4.83 MPa and temperature TC = 151 K. Determine values
for the vdW constants a and b for Ar and estimate the
diameter of an Ar atom.
1 a
PC 
27 b 2
8 a
k BTC 
27 b

k T
b B C
8 PC
27 k BTC 
27 1.38 1023 J/K 151 K
a

64 PC
64
4.83 106 Pa
2
k BTC 1.38 10 23 J/K 151 K
 29
3
b


5
.
4

10
m
8PC
8  4.83 106 Pa
27 k BTC 
a
64 PC
2

2
 3.8 1049 J 2 /Pa
b1/ 3  3.86 1010 m  3.86 A
Per mole: a=0.138 Pa m6 mol-2; b=3.25x10-5 m3 mol-1
Problem
One mole of Nitrogen (N2) has been compressed at T0=273 K to the volume
V0=1liter. The critical parameters for N2 are: VC = 3Nb = 0.12 liter/mol, TC =
(8/27)(a/kBb) = 126K. The gas goes through the free expansion process (Q = 0,
W = 0), in which the pressure drops down to the atmospheric pressure Patm=1 bar.
Assume that the gas obeys the van der Waals equation of state in the compressed
state, and that it behaves as an ideal gas at the atmospheric pressure. Find the
change in the gas entropy.
f
R ln T  R ln V  N Ab   g ( N , m )
2
f
Sideal  R ln T  R lnV  g ( N , m )
2
Tf
Tf
Vf
 Vf  5


5




S  R ln
 R ln 

R
ln

R
ln


2
T0
T0
 V0  Nb  2
 V0  VC / 3 
S vdW 
2 aN A2
9 TCVC
T f  T0 
 T0 
 266.2 K
5 RV0
20 V0
Vf 
RT f
Patm
 2.2  102 m3


5
266.2
2.2 102
  5.24 101  26  25.5 J/K
S  R ln
 R ln 
3
3 
2
273
 110  0.04 10 
P
T = const (< TC)
Phase Separation in the vdW Model
The phase transformation in the vdW model is easier to
analyze by minimizing F(V) rather than G(P) (dramatic
changes in the term PV makes the dependence G(P) very
complicated, see below).
V1
V2
V
F
V
F1
(liquid)
At T< TC, there is a region on the F(V) curve in which F makes
a concave protrusion (2F/V 2<0) – unlike its ideal gas
counterpart. Due to this protrusion, it is possible to draw a
common tangent line so that it touches the bottom of the left
dip at V = V1 and the right dip at V = V2. Since the common
tangent line lies below the free energy curve, molecules can
minimize their free energy by refusing to be in a single
homogeneous phase in the region between V1 and V2, and by
preferring to be in two coexisting phases, gas and liquid:
V  V1 N liquid V2  V
N  N gas  N liquid


N
V2  V1
N
V2  V1
F
F
V  V2
V V
F  1 N liquid  2 N gas  F1
 F2 1
N
N
V1  V2
V1  V2
N gas
F2
(gas)
- we recognize this as the common tangent line.
V < V1
V1 < V < V 2
V2 < V
As usual, the minimum free energy principle controls
the way molecules are assembled together.
P
Pvap
7
T < TC
Phase Separation in the vdW Model (cont.)
3
Since the tangent line F(V) maintains the same slope
between V1 and V2, the pressure remains constant
between V1 and V2:
 F 
6
4
2
1
5
V
F
V
Fliq

  P

V

T , N
In other words, the line connecting points on the PV
plot is horizontal and the two coexisting phases are in a
mechanical equilibrium. For each temperature below
TC, the phase transformation occurs at a well-defined
pressure Pvap, the so-called vapor pressure.
Fgas
P
P
critical
point
Pvap(T1)
triple
point
Vliq
Vgas
V
T1
T
Two stable branches 1-2-3 and
5-6-7 correspond to different
phases. Along branch 1-2-3 V
is large, P is small, the density
is also small – gas. Along
branch 5-6-7 V is small, P is
large, the density is large –
liquid. Between the branches –
the
gas-liquid
phase
transformation, which starts
even before we reach 3 moving
along branch 1-2-3.
Phase Separation in the vdW Model (cont.)


T / TC
P / PC
n / nC
n / nC
P / PC
  3n 
 9
 n
1
 3nC 
 F 
C 

     k BT ln 
 ln 
 1 
 k BTC  

 N T ,V
 n
 3nC / n  1 4
  nQ 
 nC 
For T<TC, there are three values of n with the same  . The outer two values of n
correspond to two stable phases which are in equilibrium with each other.
The kink on the G(V) curve is a signature of the 1st order transition. When we move
along the gas-liquid coexistence curve towards the critical point, the transition
becomes less and less abrupt, and at the critical point, the abruptness disappears.
P
The Maxwell Construction
T < TC
7
3
Pvap
6
[finding the position of line 2-6 without analyzing F(V)]
2
4
1
5
V
F
V
On the one hand, using the dashed line on the F-V
plot:
 F 
Fgas  Fliquid  
 V2  V6    Pvap V2  V6 
 V T , N
Fliq
On the other hand, the area under the vdW isoterm
2-6 on the P-V plot:
Fgas
 F 
PdV


V
V  V T , N dV   Fgas  Fliquid 
6
6
V2
G
4
5
3
V2
Thus, V2
7
 P V dV P V
2,6
vdW
2
 V6 
V6
1
the lowest branch
represents the stable
phase, the other branches
are unstable
vap
P
the areas 2-3-4-2 and 4-5-6-4
must be equal !
-Maxwell construction
Problem
P
The total mass of water and its saturated vapor
(gas) mtotal = mliq+ mgas = 12 kg. What are the
masses of water, mliq, and the gas, mgas, in the
state of the system shown in the Figure?
V
Vliq
Vgas
Vliq increases from 0 to V1 while the total volume
decreases from V2 to V1. Vgas decreases from V2
to 0 while the total volume decreases from V2 to
V1. When V = V1, mtotal = mliq. Thus, in the state
shown in the Figure, mliq  2 kg and mgas  10
kg.
V1
V2
V
P
Phase Diagram in T-V Plane
T < TC
Pvap
V
F
V
Fliq
Fgas
T
Single Phase
TC
At T >TC, the N molecules can exist in a single phase in
any volume V, with any density n = N/V. Below TC, they
can exist in a homogeneous phase either in volume V < V1
or in volume V > V2. There is a gap in the density allowed
for a homogeneous phase.
There are two regions within the two-phase “dome”:
metastable (P/V< 0) and unstable (P/V>0). In the
unstable region with negative compressibility, nothing can
prevent phase separation. In two metastable regions,
though the system would decrease the free energy by
phase separation, it should overcome the potential
barrier first. Indeed, when small droplets with radius R are
initially formed, an associated with the surface energy term
tends to increase F. The F loss (gain) per droplet:
condensation:
unstable
interface:
V1(T)
P
2F
 0,
0
V
V 2
V2(T)
V
N1  4

3
  R / V1 
N 3

 4 R 2
Total balance:
VC
 F
 F   4 R 2  F
F
N1  4

3
  R / V1 
N 3

RC
R
heating
Joule-Thomson Process for the vdW Gas
cooling
The JT process corresponds to an isenthalpic expansion:
 H 
 H 
H  
 T  
 P  0

T

P

P

T
 T 


 P  H
 H 

  CP
 T  P
 H 
 V 

T

 V

P

T
T  
P
 
CP
CP
 H 
CP T  
 P
 P T
H  TS  VP
 S 
 V 
   

 P T
 T  P
 H 
 S 

  T   V
 P T
 P T
(see Pr. 5.12)
This is a pretty general (model-independent) result. By applying this result to the vdW
equation, one can qualitatively describe the shape of the inversion curve (requires solving
cubic equations...).
We’ll consider the vdW gas at low densities:

N 2a 
 P  2 V  Nb  NkBT
V 

N 2 a  V 
 V 
P
  2 
  Nk B

T
V

T

P

P
N 2a
P 
V2
V  Nb
N 2a

...P
PV 
 PNb  NkBT
V
T
Nk B
 V 

 
2
 T  P P  N a
V2
 T 

 

P

H
2
Na
 Nb
k BT
CP
Joule-Thomson Process for the vdW Gas (cont.)
Cooling:
2a
b  0
k BT
Heating:
2a
b  0
k BT
If b = 0, T always decreases in the JT process: an increase of Upot at the expense of K.
If a = 0, T always increases in the JT process (despite the work of molecular forces is 0):
PV  Nb  RT
PV  RT  PNb
H  U  PV  CV  R T  PNb  CPT  PNb
The upper inversion temperature:
(at low densities)
heating
cooling
2a
b  0
k BTINV
TINV 
2a
27

TC
k Bb
4
(TC – the critical temperature
of the vdW gas, see below)
Substance
TINV
(P=1 bar)
CO2
(2050)
CH4
(1290)
O2
893
N2
621
H2
205
4He
51
3He
(23)
Thus, the vdW gas can be liquefied by
compression only if its T < 27/4TC.
Problem
The vdW gas undergoes an isothermal expansion from volume V1 to volume V2.
Calculate the change in the Helmholtz free energy.
In the isothermal process, the change of the Helmholtz free energy is
dF T , N   SdT  PdV  dN T , N
  PdV
 V2  Nb 
 RT
N 2a 
1
2  1





F    PdV    
 2 dV N a     RT ln 
V  Nb V 
 V1 V2 
 V1  Nb 
V1
V1 
V2
compare with
V2
F  U  TS

U vdW T  N 2a  1 
 V1
1

V2 
TSvdW T