Maxwell construction explained
In thermodynamic equilibrium, a necessary condition for stability is that pressure,
, does not increase with volume, or molar volume,
; this is expressed mathematically as
, where
is the temperature.
[1] This basic stability requirement, and similar ones for other
conjugate pairs of variables, is violated in analytic models of first order
phase transitions. The most famous case is the
van der Waals equation,
[2] [3] where
are dimensional constants. This violation is not a defect, rather it is the origin of the observed discontinuity in properties that distinguish liquid from vapor, and defines a first order phase transition.Figure 1 shows an isotherm drawn, for
, as a continuously differentiable solid black, dotted black, and dashed gray curve. The decreasing part of the curve to the right of point C in Fig. 1 describes a gas, while the decreasing part to the left of point E describes a liquid. These two parts are separated by a region between the local minimum and local maximum on the curve with positive slope that violates the stability criterion. This mathematical criterion expresses a physical condition which Epstein
[4] described as follows:
"It is obvious that this middle part, dotted in our curves [dashed in Fig.1 here], can have no physical reality. In fact, let us imagine the fluid in a state corresponding to this part of the curve contained in a heat conducting vertical cylinder whose top is formed by a piston. The piston can slide up and down in the cylinder, and we put on it a load exactly balancing the pressure of the gas. If we take a little weight off the piston, there will no longer be equilibrium and it will begin to move upward. However, as it moves the volume of the gas increases and with it its pressure. The resultant force on the piston gets larger, retaining its upward direction. The piston will, therefore, continue to move and the gas to expand until it reaches the state represented by the maximum of the isotherm. Vice versa, if we add ever so little to the load of the balanced piston, the gas will collapse to the state corresponding to the minimum of the isotherm." This situation is similar to a body exactly balanced at the top of a smooth surface that, with the slightest disturbance will depart from its equilibrium position and continue until it reaches a local minimum. As they are described such states are dynamically unstable, and consequently they are not observed. The gap
is a precursor of the actual phase change from liquid to vapor. The points E
and C
, where
, that delimit the largest possible liquid and smallest possible vapor states are called
spinodal points. Their locus forms a spinodal curve which bounds a region where no homogeneous stable states can exist.
Experiments show that if the volume of a vessel containing a fixed amount of liquid is heated and expands at constant temperature, at a certain pressure,
, vapor, (denoted by dots at points
and
in Fig. 1) bubbles nucleate so the fluid is no longer homogeneous, but rather it has become a heterogeneous mixture of boiling liquid and condensing vapor. Gravity separates the boiling (saturated) liquid,
, from the less dense condensing (saturated) vapor,
that coexist at the same saturation temperature and pressure. As the heating continues the amount of vapor,
increases and that of the liquid,
, decreases. All the while the pressure,
and temperature,
, remain constant and the volume
increases. In this situation the molar volume of the mixture is a weighted average of its components
where
, the mole fraction of the vapor,
, increases continuously; however, the molar volume of the substance itself has only the largest possible stable value for its liquid state, and smallest possible stable value for its vapor state at the given
. To repeat, although the mixture molar volume passes continuously from
to
(denoted by the dashed line in Fig. 1), the underlying fluid has a discontinuity in this property, and others as well. This equation of state of the mixture is called the
lever rule.
[5] [6] [7] The dotted parts of the curve in Fig. 1 are metastable states. For many years such states were an academic curiosity; Callen[8] gave as an example,
water that has been cooled below 0°C at a pressure of 1 atm. A tap on a beaker of water in this condition precipitates a sudden dramatic crystallization of the system.
However, studies of boiling heat transfer have made clear that metastable states occur routinely as an integral part of this process. In it the heating surface temperature is higher than the saturation temperature, often significantly so, hence the adjacent liquid must be superheated.
[9] Further the advent of devices that operate with very high heat fluxes has created interest in the metastable states, and the thermodynamic properties associated with them, in particular the superheated liquid states.
[10] Moreover, the fact that they are predicted by the van der Waals equation, and cubic equations in general, is compelling evidence of its efficacy in describing phase transitions; Sommerfeld described this as follows:
[11] It is very remarkable that the theory due to van der Waals is in a position to predict, at least qualitatively, the existence of the unstable [called metastable here] states along the branches AA` or BB` [BC and FE in Fig. 1 here].
Equal area rule
The discontinuity in
, and other properties,
e.g. internal energy,
, and entropy,
, of the substance, is called a first order phase transition.
[12] [13] In order to specify the unique experimentally observed pressure,
, at which it occurs another thermodynamic condition is required, for from Fig.1 it could clearly occur for any pressure in the range
. Such a condition was first enunciated in a clever thermodynamic argument by Maxwell at a lecture he delivered to the British Chemical Society on Feb 18, 1875
[14] (Fig.1, including the letters B C D E F, is the curve he described):"The portion of the curve from
C to
E represents points which are essentially unstable, and which cannot therefore be realized.
Now let us suppose the medium to pass from B to F along the hypothetical curve B C D E F in a state always homogeneous, and to return along the straight line path F B in the form of a mixture of liquid and vapor. Since the temperature has been constant throughout, no heat can have been transformed into work. Now the heat transformed into work is represented by the excess of the area F D E over B C D. Hence the condition which determines the maximum pressure of the vapor at given temperature is that the line B F cuts off equal areas from the curve above and below."
The easiest way to understand Maxwell's argument is to consider the cycle he suggested on a temperature—molar entropy plane. Every introductory thermodynamics text[15] [16] presents the fact that on such a plane the area under any curve is the heat transfer to the substance per mole, positive going from left to right and negative from right to left; moreover, in a cyclic process the net heat transfer to the substance is the area enclosed by the cycle's closed curve. Since the cycle he considered is composed of the two gray dashed isothermals at the same temperature, one proceeding from B to F (through C D and E), and the other directly back from F to B, the two lines are identical, just traversed in reverse; there is zero area enclosed, and hence
.
Further the same texts describe the area under these curves when plotted on a pressure—molar volume, see Fig. 1, as being the work done by the substance, positive going from left to right, and negative from right to left. Likewise the net work done in a cycle is the area enclosed by the closed curve. Since the first law of thermodynamics yields in the special case of a cycle
, for the cycle envisioned by Maxwell
; then since the area enclosed is I+II=0, see Fig.1, with I positive and II negative, the transition pressure must be such that the two areas are equal.
Written as a mathematical equation in terms of the work done in each process this isThis equation together with the equation of state written for each of the states
and
are three equations for the four variables,
, so given any one of them, say
, the other three are determined. In other words, there is a unique value of
, as well as
and
, at which the phase transition can occur.
Gibbs criterion
At the end of his lecture, after complimenting van der Waals by referring to his work as "an exceedingly ingenious thesis", Maxwell finished it by saying:
I must not, however, omit to mention a most important American contribution to this part of thermodynamics by Prof. Willard Gibbs of Yale College U.S., who has given us a remarkably simple and thoroughly satisfactory method of representing the relations of the different states of matter by means of a model. By means of this model, problems which had long resisted the efforts of myself and others may be solved at once.
This remark proved prescient because in 1876-1878
Gibbs published his definitive work on thermodynamics
[17] in which he showed that thermodynamic equilibrium of a
heterogeneous substance requires that, in addition to
mechanical equilibrium (the same pressure for each component) and
thermal equilibrium (the same temperature for each component), there must also be
material equilibrium (the same chemical potential for each component). In the present instance of one substance and two phases in addition to
and
, material equilibrium requires
(for the special case of one substance its chemical potential is the molar Gibbs function,
where
).
[18] This condition can be deduced by a simple physical argument as follows: the energy required to vaporize a mole is from the second law at constant temperature
, and from the first law at constant pressure
, then equating these two and rearranging produces the result since
. The conditions of material equilibrium lead to the famous Gibbs
phase rule,
, where
is the number of substances,
the number of phases, and
the number of independent intensive variables required to specify the state.
[19] [20] In the case of one substance and two phases discussed here this gives,
, the experimentally observed number.
Now
is a
thermodynamic potential function, its differential is
[21] Integrating this at constant temperature produceshere
is a constant of integration, but the constant is different for each isotherm hence it is written as a function of
.
[22] In order to evaluate
one must invert
to obtain
. However, it is the nature of the phase transition phenomenon that this inversion is not unique; for example the van der Waals quation written for
is,
a cubic with either 1 or, in this case, 3 real roots. Thus there are three curves, as seen in Fig. 2, consisting of stable (shown solid black), metastable (shown dotted black), and unstable (shown dashed gray) states.
Actually, the figure was not produced by solving the cubic and integrating, rather
was obtained from its definition by first obtaining
and
, which is easily done analytically for the van der Waals equation, and plotting it
parametrically with
, using
as the parameter. Considering only its stable states
is continuous with discontinuous
partial derivatives,
and
, at the phase transition point. In the
Ehrenfest classification, a first order phase transition refers to the discontinuity of the first partial derivatives of
while a second order phase transition would involve discontinuities of the second partial derivatives.
[23] Relationship between the Gibbs and Maxwell criteria
Evaluating the integral expression for
given previously between the saturated liquid and vapor states and applying the Gibbs criterion of material equilibrium to this phase change process requires writing it as
Here the integral has been split into three parts using the three real roots of the cubic corresponding to the liquid,
, unstable,
, and vapor,
, states respectively. These integrals can best be visualized by viewing Fig. 1 rotated
counterclockwise in the paper plane then
about the
axis so that
appears on the leftside ordinate of the curve as shown in the accompanying graph. In this view the function
clearly is multi-valued; this is the reason it requires three real functions to describe its behavior between
and
. Now on splitting the middle integral into two
The first two integrals here are area I while the second two are the negative of area II. The two areas add to zero hence their magnitudes are equal according to this Gibbs criterion. This is again the equal area rule of Maxwell, the Maxwell construction, and it can also be shown analytically. Since
,
Integrating this for constant temperature from state
to
with the Gibbs condition produces
which is Maxwell's result. This equal area rule can also be derived by making use of the Helmholtz free energy.
[24] In any event the Maxwell construction derives from the Gibbs condition of material equilibrium. However, even though
is more fundamental it is more abstract than the equal area rule, which is understood geometrically.
Common tangent construction
Another method to determine the coexistence points is based on the Helmholtz potential minimum principle, which states that in a system in diathermal contact with a heat reservoir
,
and
, namely at equilibrium the Helmholtz potential is a minimum.
[25] Since, like
, the molar Helmholtz function
is also a potential function whose differential is,
[26] this minimum principle leads to the stability condition
\partial2f/\partial
p/\partialv|T>0
.
[27] This condition requires that at any stable state of the system the function
is strictly
convex, namely that in its vicinity the curve lies on or above its tangent.
[28] Moreover, for those states the previous stability condition for the pressure is necessarily satisfied as well.
A plot of this function for the same subcritical isotherm of the vdW equation as Figs. 1 and 2 is shown in Fig. 3. Included in this figure is the (dashed/solid) straight line that has a double (common) tangent with the curve of the function
at B and F. This straight line is,
, with
constant, which can be written as
. The last equality follows from the relation
,
[29] together with
.
[18] All this means that every point on the line has the same values of
, in particular the points B and F, which produces the Gibbs condition for material equilibrium
as well as eqality of temperature and pressure.
[30] Therefore this construction is equivalent to both the Gibbs conditions and the Maxwell construction.
This construction, based on
defined earlier by Gibbs,
[31] [32] was originally used by van der Waals (he called it both a double and common tangent),
[33] because it could be easily extended to include binary fluid mixtures for which an isotherm of
, with
a composition variable, forms a surface that can have a common tangent plane. It has subsequently become a popular way to treat phase change problems.
[34] [35] Application to the van der Waals equation
From the van der Waals equation applied to the saturated liquid,
, and vapor,
, states
These two equations specify 4 variables so they can be solved for
in terms of
. This results in
where
,
, and
are a characteristic pressure, molar volume, and temperature defined by the constants (note that
). Applying the Maxwell construction to the van der Waals equation gives
These three equations can be solved numerically. This has been done given a value for either
or
, and tabular results presented;
[36] [37] however, the equations also admit an analytic parametric solution that, according to Lenkner,
[38] was obtained by Gibbs. Lenkner himself devised a simple, elegant, method to obtain this solution, by eliminating
and
from the equations, and writing them in terms of a stretched dimensionless number density,
, that varies between
and 0 as
varies from
to
; this produces
Although transcendental, this equation has a simple analytic, parametric solution obtained by writing the left side of the equation, which is just
, as
Then
and when used to eliminate
from the right hand side a linear equation for
is obtained, whose solution is
Accordingly, the fundamental variable that specifies all the others in this phase transition process is
. This solution to the saturation problem is easily extended to encompass all its variables
where
g(y)=[1+\varrhog(y)][1+\varrhof(y)]=1+2f(y)\coshy+f(y)2
.The values of all other property discontinuities across the saturation curve also follow from this solution.
[39] These functions define the coexistence curve which is the locus of the saturated liquid and saturated vapor states of the van der Waals fluid. In Fig. 4 this curve is plotted in blue together with the spinodal curve in black, calculated from where
is a parameter. The variables used in making these plots are the reduced (dimensionless) variables,
,
, and
where the
subscripted quantities are the critical point values. They are defined by,
, and
at the critical point,
[40] and are measurable quantities. The relations
,
,
are used to convert the star quantities in the solution to the
quantities used in the figures. The curve agrees completely with the numerical results referenced earlier. In the region inside the spinodal curve there are two states at each point, one stable and one metastable, either superheated liquid to the right of the blue curve, or subcooled vapor to the left, while outside the spinodal curve there is one stable state at each point. In Fig. 5 the region under the (dot dash black) spinodal curve contains no homogeneous stable states while between the (dot dash red) coexistence and spinodal curves there is one metastable state at every point, and outside the coexistence curve there is one stable state at each point. The two blue and two green circles denote the saturated liquid and vapor states on their respective isotherms. There are also observed heterogeneous states everywhere under the coexistence curve that satisfy the lever rule; however, they are not homogeneous states of the van der Waals equation, so their existence, indicated by horizontal lines connecting the saturation points on each subcritical isotherm, is not displayed. Also the abscissa in this figure is logarithmic, not linear, in order to show more of the vapor region at large
without excessively compressing the liquid and unstable regions at small
; however, this device distorts areas, so the two areas I and II in Fig.1 would not appear equal here.
Over the parameter range
,
decreases monotonically from
and approaches 0 as
in the limit
. Therefore
\varrhof(0)=\varrhog(0)=1/2
and in the limit
,
and
\varrhog(y)\sim2(y-1)e-2y
. The behavior of
and
follow from the equations. Both these properties also decrease monotonically from
and
, and approach 0 as
and
in the limit
. Note from these that
ps\sim(p*/T
g\sim(p*v*/T
g=\rhogRTs
; the van der Waals saturated vapor is an ideal gas in this limit. To paraphrase Sommerfeld, it is remarkable that the theory due to van der Waals is able to predict that when
the saturated vapor behaves like an ideal gas; the saturated vapor of real gases behave exactly this way.
In addition for
the liquid spinodal point occurs at a negative pressure, and the isotherm
is included in Fig. 4 to illustrate this point. This means that some part of those liquid metastable states are in tension, and the lower the temperature the greater the tensile stress. Although this seems counterintuitive it is known that under some circumstances liquids can support tension.
Tien and
Lienhard[41] noted this and wrote:
The van der Waals equation predicts that at low temperatures liquids sustain enormous tension---a fact that has led some authors to take the equation lightly. In recent years measurements have been made that reveal this to be entirely correct.[42] Liquids that are clean and free of dissolved gas can be subjected to tensions greater in magnitude than
.
This is another interesting feature of the van der Waals theory.
References
- Barrufet. M.A.. Eubank. P.T.. 1989. Generalized Saturation Properties of Pure Fluids Via Cubic Equations of State. Chemical Engineering Education. 23. 3. 168–175.
- Book: Callen, H.B.. 1960. Herbert Callen. Thermodynamics. NY. John Wiley and Sons.
- DeBoer. J.. 1974. Van der Waals in his time and the present revival opening address. Physica. 73. 1–27. 10.1016/0031-8914(74)90223-7.
- Book: Epstein, P.S. . 1937 . Paul Sophus Epstein. Textbook of Thermodynamics . NY . John Wiley and Sons.
- Book: Gibbs, J.W.. 1928 . 1876, 1878. Josiah Willard Gibbs. Longley. W.R.. Van Name. R.G.. The Collected Works of J. Willard Gibbs Volume I Thermodynamics. On the Equilibrium of heterogeneous Substances. 1948. New Haven. Yale University Press. 55–353.
- Book: Goodstein, D.L.. David Goodstein . 1985 . 1975. States of Matter . NY . Dover.
- Book: Johnston, D.C.. Advances in Thermodynamics of the van der Waals Fluid . 2014 . 10.1088/978-1-627-05532-1 . 1402.1205. 978-1-627-05532-1 .
- Lekner. J.. 1982. Parametric solution of the van der Waals liquid–vapor coexistence curve. Am. J. Phys.. 50. 2. 161–163. 10.1119/1.12877.
- Book: Lienhard IV. J.H.. Lienhard V. 2019. John H. Lienhard. John H. Lienhard V. A Heat Transfer Textbook. Fifth . Mineola, NY. Dover Publications Inc..
- Lienhard . J.H. . 1986 . The Properties and Behavior of Superheated Liquids . Lat. Am. J. Heat and Mass Transfer . 10 . 169–187.
- Maxwell. J.C.. 1875. James Clerk Maxwell. On the Dynamical Evidence of the Molecular Constitution of Bodies. Nature. 11. 279. 357–359. 10.1038/011357a0 .
- Book: Moran. M.J.. Shapiro. H.N.. 2000. Fundamentals of Engineering Thermodynamics. 4th. NY. McGraw-Hill.
- Book: Rectorys. K.. 1969. Rectories. K.. Survey of Applicable Mathematics. Differential Calculus of Functions of a Real Variable. 397–439. Cambridge, MA. The M.I.T. Press.
- Shamsundar. N.. Lienhard. J.H.. 1983. Saturation and Metastable Properties of the van der Waals Fluid. Canad J Chem Eng. 61. 6. 876–880. 10.1002/cjce.5450610617.
- Book: Sommerfeld . A. . Arnold Sommerfeld. Bopp . F. . Meixner . J. . Kestin . J. . 1956 . Thermodynamics and Statistical Mechanics - Lectures on Theoretical Physics Volume V . NY . Academic Press.
- Book: Tien . C.L. . Lienhard . J.H. . Chang-Lin Tien. 1979 . Statistical Thermodynamics Revised Printing . NY . Hemisphere Publishing.
- Book: Van Wylen. G.J.. Sonntag. R.E.. 1973. Fundamentals of Classical Thermodynamics. Second . NY. John Wiley and Sons. Gordon Van Wylen.
- Book: van der Waals . J.D. . Johannes Diderik van der Waals. Rowlinson . J.S. . John Shipley Rowlinson. 1984 . On the Continuity of the Gaseous and Liquid States, edited and with an Introduction by J.S. Rowlinson . NY . Dover Phoenix Editions.
Notes and References
- Callen, pp 131-135, (1960}
- van der Waals, p. 174, (1988)
- Epstein, p 9
- Epstein, p 10
- Callen, pp 146-163
- Goodstein, pp 443-452
- Book: Modern Thermodynamics: From Heat Engines to Dissipative Structures. Kondepudi. Dilip. Prigogine. Ilya. 2014-12-31. John Wiley & Sons. 9781118371817. en.
- Callen, p 163
- Lienhard IV and Lienhard V, pp 467-469
- Lienhard, 169-187
- Sommerfeld, p. 66
- Callen, p 146
- Goodstein, p 443
- Maxwell, J.C., pp. 357-359
- Van Wylen and Sonntag, p. 211
- Moran and Shapiro, H.N. p. 251
- Gibbs,pp. 62-65
- Callen, pp. 99
- Gibbs, pp. 96-100
- Callen, pp 163-167
- Callen, pp. 98-100
- Callen, p 150
- Callen, pp. 172-173
- Sommerfeld, A., pp 67-68, (1956)
- Callen, p. 105
- Callen, p. 120
- Callen, p. 135
- Rectorys, p. 429
- Callen, p.98
- van der Waals, p. 246
- De Boer, p.8
- Gibbs, pp. 87-90
- Van der Waals, pp. 246-247
- Book: Energy Landscapes: Applications to Clusters, Biomolecules and Glasses. Wales. David. Wales. 2003. Cambridge University Press. 9780521814157. en. 444.
- Web site: First order phase transitions and the dynamics of spinodal decomposition. www.mhkoepf.de. 2019-11-12.
- Shamsundar and Lienhard, 876-880
- Barrufet and Eubank, pp. 168-175
- Lekner, p. 161
- Johnston, pp. 16-18
- Sommerfeld, A., pp 56-57
- Tien and Lienhard, p.254
- Temperley, H.N.V., The Behavior of Water under Hydrostatic Tension: III, Proc. Phys. Soc. (London), 59, p 199, (1947)