Engineering
Vol.10 No.05(2018), Article ID:84802,21 pages
10.4236/eng.2018.105019

Absolute Reference Values of the Real Gas

Albrecht Elsner

Am Mühlbach 14, D-85748 Garching, Germany

Copyright © 2018 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: April 2, 2018; Accepted: May 22, 2018; Published: May 25, 2018

ABSTRACT

With his publication in 1873 [1] J. W. Gibbs formulated the thermodynamic theory. It describes almost all macroscopically observed properties of matter and could also describe all phenomena if only the free energy U S T were explicitly known numerically. The thermodynamic uniqueness of the free energy obviously depends on that of the internal energy U and the entropy S, which in both cases Gibbs had been unable to specify. This uncertainty, lasting more than 100 years, was not eliminated either by Nernst’s hypothesis S = 0 at T = 0 . This was not achieved till the advent of additional proof of the thermodynamic relation U = 0 at T = T c . It is noteworthy that from purely thermodynamic consideration of intensive and extensive quantities it is possible to derive both Gibbs’s formulations of entropy and internal energy and their now established absolute reference values. Further proofs of the vanishing value of the internal energy at the critical point emanate from the fact that in the case of the saturated fluid both the internal energy and its phase-specific components can be represented as functions of the evaporation energy. Combining the differential expressions in Gibbs’s equation for the internal energy, d ( μ / T ) / d ( 1 / T ) and d ( p / T ) / d ( 1 / T ) , to a new variable d ( μ / T ) / d ( p / T ) leads to a volume equation with the lower limit vc as boundary condition. By means of a variable transformation one obtains a functional equation for the sum of two dimensionless variables, each of them being related to an identical form of local interaction forces between fluid particles, but the different particle densities in the vapor and liquid spaces produce different interaction effects. The same functional equation also appears in another context relating to the internal energy. The solution of this equation can be given in analytic form and has been published [2] [3] . Using the solutions emerging in different sets of problems, one can calculate absolutely the internal energy as a function of temperature-dependent, phase-specific volumes and vapor pressure.

Keywords:

Entropy Reference Value S ( M , V , 0 ) = 0 , Internal Energy Reference Value U ( M , V , T c ) = 0 , Critical Volume V = M [ d ( μ / T ) / d ( p / T ) ] c , Two-Phase Equilibrium M ( u v u l ) ( 0 ) S T + M μ = M [ v v + v l ( v v v l ) / ln ( v v / v l ) ] d( p/T )/ d( 1/T ) + Vdp/ dT T=U +VpVpV p c , Chemical Potential Negative, M μ ( T ) M μ ( 0 ) = U ( M , V , 0 ) = M [ ( u v u l ) ( 0 ) ] < 0 .

1. Introduction

The purpose of this paper is to show that Gibbs’s theory [1] contains intrinsically “natural” reference values of entropy and internal energy that are reciprocal and thus represent thermodynamic reference values [4] . Nevertheless, it was and is still accepted till the recent past that the value of the internal energy of the real gas cannot be given absolutely, e.g. [5] [6] . Instead, the calculation of entropy and internal energy could be based on a so-called fiducial reference value [7] . To put an end to the diversity of individually, arbitrarily chosen fiducial values, it was decided at conferences in the 1950s to assign the values of entropy and internal energy (or enthalpy) of a liquid at the triple point the reference values zero, which in any case the Nernst hypothesis contradicts. On this basis mathematically complex equations were put forward, but an analysis (see Appendix) shows that they contain thermodynamic inconsistencies. These then yielded incorrect thermodynamic data published for a large number of gases, e.g. [8] - [15] . None of these skeleton tables presents data on the chemical potential, which could have been derived from the differences h v , l s v , l T . A paper that specifically investigates calculation of the chemical potential of the generally accepted conference agreements comes to the surprising conclusion that the chemical potential increases as the temperature [16] . This result cannot, however, account for daily observation that flow of freely-moving matter occurs from cold to warm regions [17] . One has to take the consequence from the numerous thermodynamic discrepancies ensuing from calculating entropy and internal energy from the assumptions mentioned, s l ( T t ) = 0 and u l ( T t ) = 0 (or h l ( T t ) = 0 ). Here, too, as so often in the history of physics, one has to abandon a trusted hypothesis. Here it is the assumption of the naive addition of fiducial constants to thermodynamic fundamental quantities.

The paper shows four possible ways of finding directly the thermodynamic reference values of entropy and internal energy and thus comply with the correct treatment of real properties of matter called for by Gibbs. First there is the possibility of studying the characteristic features of intensive and extensive quantities in order to describe thermodynamically the real properties of matter. A second investigation deals with the question what follows from representation of the internal energy of the fluid as a function of the evaporation energy. The third possibility is concerned with solution of a functional equation for the sum of two dimensionless variables, the one referring to the local interaction potential of fluid particles in the vapor space and the other to that in the liquid. The functional equation emerges from the original equation for the internal energy and has been solved [3] . The functional equation is, on the other hand, also encountered when one represents the particular phase-specific internal energy as a funtion of the evaporation or condensation energy. It has of course the same (physically unique) solution and allows the internal energy to be explicitly calculated as an absolutely determined temperature function of the measurable quantities: phase-specific volumes and vapor pressure [2] .

2. Intensive and Extensive Quantities of the Saturated Fluid

The homogeneity of the fluid allows its macroscopic properties to be described by intensive and extensive quantities.

The intensive quantities are the temperature T, vapor pressure p and chemical potential μ; vapor pressure and chemical potential are pure temperature functions below the critical point of the fluid. With the finite critical values pc and d n 1 p / d T c n 1 > d n p / d T c n ( n = 1 , 2 , 3 , ) the ò-expansion of p to second order yields

0 T T c , 0 ϵ = 1 T / T c 1 ,

0 p = p c d p d T | c T c ϵ + 1 2 d 2 p d T 2 | c T c 2 ϵ 2 p c ,

0 d p d T = d p d T | c d 2 p d T 2 | c T c ϵ + 1 2 d 3 p d T 3 | c T c 2 ϵ 2 d p d T | c ,

0 d 2 p d T 2 = d 2 p d T 2 | c d 3 p d T 3 | c T c ϵ + d 2 p d T 2 | c ,

d ( p / T ) d ( 1 / T ) = p d p d T T 0 , d d T d ( p / T ) d ( 1 / T ) = d 2 p d T 2 T 0. (1.1)

For T 0 the vapor pressure is a positive and convexly curved temperature function increasing with T from 0 to pc. Also the temperature derivatives increase with T from d n p / d T n ( 0 ) = 0 to finite values d n p / d T n ( T c ) (where n = 1 , 2 , ). In contrast, the chemical potential is a negative and concavely curved temperature function decreasing with T. With the finite critical values μ c and d μ / d T c , and, on the other hand, the divergent terms d n μ / d T c n ( n = 2 , 3 , ) an ò-expansion of μ is not possible:

0 > μ μ c , 0 d μ d T d μ d T | c , 0 d 2 μ d T 2 d 2 μ d T 2 | c , (1.2)

0 > d ( μ / T ) d ( 1 / T ) d ( μ / T ) d ( 1 / T ) | c , 0 d d T d ( μ / T ) d ( 1 / T ) = d 2 μ d T 2 T < d 2 μ d T 2 | c T c + ,

0 > μ p > μ p | c , 0 > d μ d p > d μ d p | c , 0 < d 2 μ d p 2 | c < d 2 μ d p 2 , 0 < d ( μ / T ) d ( p / T ) | c d ( μ / T ) d ( p / T ) .

Thermodynamics treats the quantity X as an extensive quantity, which means that X is proportional to the fluid mass M. The validity of the relation X = x M leads to x having the property of additivity and ensures its uniqueness. At temperatures below the critical point the fluid mass M in the volume V is additively composed of the vapor mass Mv and the condensed mass Ml in the sub-volumes Vv and Vl: M = M v + M l in V = V v + V l . The same applies to other extensive quantities such as the entropy S, internal energy U, enthalpy H, free energy F, and heat capacity C. Denoting such quantities by X and X v , l and the corresponding mass-specific quantities by x = X / M and x v , l = X v , l / M v , l , one obtains the following definitions:

M = M v + M l , X = X v + X l , X = { V , S , U , H , F , C } , x = X / M , x v , l = X v , l / M v , l , x = { v , s , u , h , f , c } . (1.3)

It is worth mentioning that the temperature variation of the ratio of the differences ( x x l ) and ( x v x ) of quantities such as the volume, entropy, internal energy, enthalpy, free energy, and specific heat is the same and equal to that of the ratio of vapor to condensed masses, viz.

0 M v M l = v v l v v v = s s l s v s = u u l u v u = h h l h v h = f f l f v f = c c l c v c 1 , (1.4)

where the equality signs are valid for T = 0 and T = T c , respectively. From Equation (1.4) one arrives at the interdependence of two extensive quantities x and y and their phase-specific values x v , l and y v , l in the following form:

x ( y l y v ) + y ( x v x l ) x v y l + x l y v = 0. (1.5)

This relationship can also be deduced from the correlations x = ( M v / M ) x v + ( M l / M ) x l and y = ( M v / M ) y v + ( M l / M ) y l with ( M v / M ) + ( M l / M ) = 1 .

The decomposition of mass M into M v and M l below the critical point occurs within limits and is given by

0 M v M = v v l v v v l 1 2 M l M = v v v v v v l 1. (1.6)

While the mass-specific quantity x constitutes an average of the quantity X in V in relation to the total mass M and is thus a function of T and V/M, the mass- and phase-specific quantity xv describes the quantity Xv in the volume Vv and is related to the vapor mass Mv, and the quantity xl describes Xl in Vl and is related to Ml. The quantities x v , l then give thermodynamic information on the masses M v , l in V v , l , which are subject to equal values of temperature T, vapor pressure p = p ( T ) and chemical potential μ = μ ( T ) in V, and different density values v v 1 in Vv and v l 1 in Vl. Since the densities of the vapor and condensate in Vv and Vl are functions of the temperature alone, the quantities x v , l in V v , l are likewise pure temperature functions. As vv is different to vl for T < T c , the value xv is different to xl for T < T c , and as vv and vl are equal to the critical value vc at the critical point ( T c , p c , μ c ) , the values xv and xl are equal to the critical value xc at ( T c , v c ) . The information from xv and xl implicitly contains all particle interactions that can be expressed in terms of various imaginable types of descriptions of thermodynamic properties.

If a thermodynamic quantity is represented in its domain of definition [ T , v ] by a thermodynamic function x ( T , v ) , a thermodynamic quantity is always an absolute quantity. For example, the physics of the real gas operates in the temperature range [ 0, T c ] in the limits 0 = S ( 0 , M , V ) < S ( T , M , V ) < S ( T c , M , V ) = [ V p c M μ c ] / T c and M ( u v u l ) ( 0 ) = U ( 0 , M , V ) < U ( T , M , V ) < U ( T c , M , V ) = 0 .

The critical value xc of the quantity x = { v , s , u , f } is finite. In fact, when xv approaches the finite value xc from below, then xl approaches xc from above and vice versa; in any case, one has x v ( T c ) x l ( T c ) = 0 . The approach is determined by lim T T c ( d x v ( T ) / d T ) = lim T T c ( d x l ( T ) / d T ) . The ratio value of x v ( T ) / x l ( T ) at the critical point is then different for the two possible cases of equal or opposite sign of x v ( T ) and x l ( T ) . One has

x v ( T c ) = x l ( T c ) 0 , x v ( T c ) / x l ( T c ) = + 1 for x v ( T ) / x l ( T ) > 0 , x v ( T c ) = x l ( T c ) = 0 , x v ( T c ) / x l ( T c ) = 1 for x v ( T ) / x l ( T ) < 0. (1.7)

The relations of the first line of (1.7) are valid for the quantities x = { v , s , f } .

The consequence from relations (1.7) for the quantity x = u shall be investigated. It can be stated that the vapor energy U v at low temperatures is positive since vapor particles are so far apart that their (negative) interaction potentials are vanishingly small in comparison with their (positive) thermal energies. At low temperatures one thus has u v = U v / M v > 0 . Under these conditions the evaporation energy ( u v u l ) M is very much larger than U v . From ( u v u l ) M U v one obtains ( u v u l ) U v / M = u v M v / M or, with M v / M 1 / 2 , ( u v u l ) > u v > u v / 2 > u v M v / M 0 . This gives an estimate of the mass-specific energy ul in the form of u l = u v ( u v u l ) 0 , i.e. the vapor energy uv is not negative and the condensate energy ul is not positive. The second line of relations (1.7) then states u c = 0 . This yields the important results,

u l ( T ) 0 u v ( T ) , u l ( T c ) = 0 = u v ( T c ) , u v ( T c ) / u l ( T c ) = 1. (1.8)

From Equations (1.3) one obtains the thermodynamic relation of the mean fluid quantity x (which is a function of T and v) to the phase-specific quantities x v , l (which are pure temperature functions) as follows:

x ( T , v ) = M v M ( T , v ) x v ( T ) + M l M ( T , v ) x l ( T ) = x v M l M ( x v x l ) = x l + M v M ( x v x l ) = v v l v v v l x v + v v v v v v l x l = x v v v v v v v l ( x v x l ) = x l + v v l v v v l ( x v x l ) . (1.9)

Differentiation of x ( T , v ) with respect to v at fixed T yields the coefficient of isothermal phase transition,

( x v ) T = x v x l v v v l . (1.10)

Thus the fluid quantity x ( T , v c ) can be expressed in terms of x v , l , v v , l , v c , and ( x / v ) T :

x ( T , v c ) = x v , l ( T ) ( v v , l ( T ) v c ) ( x / v ) T . (1.11)

If the saturated fluid does not have the critical volume vc, but the volume v, the following conversion has to be made:

x ( T , v ) = x ( T , v c ) ( v c v ) ( x / v ) T = x v , l ( v v , l v ) ( x / v ) T . (1.12)

At T = T c the fluid takes the critical volume vc. From Equation (1.9) it is immediately obvious that the critical values x ( T c , v c ) , x v ( T c ) , x l ( T c ) , and the zero-point values x ( 0, v ) and x l ( 0 ) are respectively equal,

x ( T c , v c ) = x v ( T c ) = x l ( T c ) = x c , x ( 0 , v ) = x l ( 0 ) , (1.13)

and from Equation (1.10) it follows that

x l ( T ) < x ( T , v ) < x v ( T ) if ( x / v ) T > 0 ,

x l ( T ) > x ( T , v ) > x v ( T ) if ( x / v ) T < 0. (1.14)

In the theory it is not only the difference of the phase-specific quantities xv and xl, i.e.

x v x l = ( v v v l ) ( x ( T , v ) / v ) T , (1.15)

that is of importance, but also their sum ( x v + x l ) = 2 x ( T , v ) + ( v v + v l 2 v ) ( x ( T , v ) / v ) T . Since ( v v + v l 2 v ) / ( v v v l ) = ( [ v v v ] [ v v l ] ) / ( v v v l ) = M l / M M v / M 0 , one has

x v + x l = 2 x + ( x v x l ) ( M l / M M v / M ) . (1.16)

The sign of the function ( x v + x l ) will subsequently be of interest. It is the same as that of the function ( x v x l ) if the product function ( x v + x l ) ( x v x l ) has a positive sign, while the signs of ( x v + x l ) and ( x v x l ) are opposite if ( x v + x l ) ( x v x l ) = 2 x ( x v x l ) + ( x v x l ) 2 ( M l / M M v / M ) is negative. The latter can be the case if either the product function x ( x v x l ) or the difference ( x v 2 x l 2 ) is negative (see Equations (1.25) and (1.26) below). At the critical point one has ( x v + x l ) ( T c ) = 2 x ( T c , v c ) and at absolute zero ( x v + x l ) ( 0 ) = 2 x ( 0, v ) + ( x v x l ) ( 0 ) and x v ( 0 ) = x ( 0 , v ) + ( x v x l ) ( 0 ) ; for the vapor-phase quantity x v one then obtains

x v ( T c ) = 2 x ( T c , v c ) x l ( T c ) , x v ( 0 ) = x ( 0 , v ) + ( x v x l ) ( 0 ) . (1.17)

The critical mean fluid value x ( T c , v c ) is thus equal to the mean of the

phase-specific critical values 1 2 [ x v ( T c ) + x l ( T c ) ] and, if the fluid value

x ( 0, v ) at absolute zero is given by the condensation energy value ( x v x l ) ( 0 ) , the vapor value x v ( 0 ) vanishes, i.e. it holds, for example, that u v ( 0 ) = 0 .

At this place the quantum state of the Bose-Einstein condensation should be noticed. In contrast to the thermodynamic temperature absolute zero the lowest temperature available is the transition temperature T t r = ( 2 / m k ) ( 2 π / e 5 / 3 ) n 2 / 3 > 0 , where = h / 2 π is the Planck constant, k the Boltzmann constant, m the particle mass and n the particle density in the condensate. The atomic densities achieved in experiments range from 10-14 to 10-15 cm-3 and transition temperatures from 100 nK to a few μK [18] . The internal energy of the dilute gas is positive for T > T t r and vanishes at T = T t r , whereas that of the condensate is negative [19] .

2.1. Interdependence of Extensive and Intensive Quantities

Since the two-phase equilibrium can be described by extensive as well by intensive quantities, an interdependence between these quantities exists. Thermodynamics yields for the quotients ( x v x l ) / ( v v v l ) of the particular quantities x = { f , s , u , c } the well-known vapor pressure relations,

( f v ) T = f v f l v v v l = p 0 , ( s v ) T = s v s l v v v l = d p d T 0 , ( u v ) T = u v u l v v v l = d ( p / T ) d ( 1 / T ) 0 , ( c v ) T = c v c l v v v l = d 2 p d T 2 T 0. (1.18)

Correspondingly, for the quotients ( x v v l x l v v ) / ( v v v l ) the chemical potential relations are

f v v l f l v v v v v l = μ > 0 , s v v l s l v v v v v l = d μ d T 0 , (1.19)

u v v l u l v v v v v l = d ( μ / T ) d ( 1 / T ) > 0 , c v v l c l v v v v v l = d 2 μ d T 2 T 0.

Equations (1.18) and (1.19) allow one to define volume functions ( x v v l x l v v ) / ( x v x l ) , which can be represented in different ways:

x v v l x l v v x v x l = v v + v l x v v v x l v l x v x l = v l + v v v l 1 x v / x l = v v v v v l 1 x l / x v = v v + v l 2 v v v l 2 x v + x l x v x l . (1.20)

In evaluating the critical value of a volume function (1.20) one should give heed, in respect of Equations (1.7) in the case x v ( T c ) / x l ( T c ) = 1 , to whether a finite limiting value exists; in the case x v ( T c ) / x l ( T c ) = 1 the limiting value is vc. Hence the result is:

u v v l u l v v u v u l ( T c ) = d ( μ / T ) d ( p / T ) ( T c ) = v c . (1.21)

With x v , l = { f v , l , s v , l , u v , l , c v , l } one obtains

μ p = f v v l f l v v f v f l < 0 , d μ d p = s v v l s l v v s v s l < 0 , (1.22)

d ( μ / T ) d ( p / T ) = u v v l u l v v u v u l > 0 , d 2 μ d T 2 / d 2 p d T 2 = c v v l c l v v c v c l < 0.

From Equations (1.22) it immediately follows that the critical value is finite for x = { f , s , u } and divergent for x = c .

2.2. Entropy and Internal Energy Relations

According to relation (1.5), the interdependence of volume and entropy is v ( s l s v ) + s ( v v v l ) v v s l + v l s v = 0 . Taking relations (1.18) and (1.19) into account, viz. ( s v s l ) = ( v v v l ) d p / d T and ( v v s l v l s v ) / ( v v v l ) = d μ / d T , this can be transformed to 0 d μ d T = s l v v v l s v v v v l = s v d p d T . Thus one gets Gibbs’s entropy relations, including the thermodynamic reference value 0:

0 v d p d T d μ d T + v d p d T = s ( v , T ) , (1.23)

0 s l ( T ) , s v , l ( T ) = d μ d T + v v , l d p d T = s ( v , T ) + [ v v , l ( T ) v ] d p d T .

The value s ( v , T ) is the sum of the terms v d p / d T 0 and ( d μ / d T ) 0 and is positive for T > 0 and vanishes for T = 0 . The same is valid for the phase-specific entropy s l ( T ) ; and since 0 < v l ( T ) v ( T c ) = v c v v ( T ) one obtains the following sequences

0 s l ( T ) s ( v c , T ) s ( v c , T c ) = s c s v ( T ) . (1.24)

The interdependence of volume and internal energy is v ( u l u v ) + u ( v v v l ) v v u l + v l u v = 0 , which with ( u v u l ) = ( v v v l ) d ( p / T ) / d ( 1 / T ) and ( v v u l v l u v ) / ( v v v l ) = d ( μ / T ) / d ( 1 / T ) gives Gibbs’s internal energy relations and with respect to relation (1.21) the thermodynamic reference value 0:

u ( v , T ) = d ( μ / T ) d ( 1 / T ) v d ( p / T ) d ( 1 / T ) d ( μ / T ) d ( 1 / T ) ( T c ) v c d ( p / T ) d ( 1 / T ) ( T c ) 0 , u v , l ( T ) = d ( μ / T ) d ( 1 / T ) v v , l ( T ) d ( p / T ) d ( 1 / T ) = u ( v , T ) [ v v , l ( T ) v ] d ( p / T ) d ( 1 / T ) . (1.25)

For 0 T T c the value u ( v , T ) is the sum of the negative term d ( μ / T ) / d ( 1 / T ) and the positive term ( v d ( p / T ) / d ( 1 / T ) ) , where it holds that d ( μ / T ) / d ( 1 / T ) v d ( p / T ) / d ( 1 / T ) 0 and therefore u ( v , T ) is not positive. Furthermore, the phase-specific internal energies obey the relations ( u v u l ) = ( v v v l ) d ( p / T ) / d ( 1 / T ) 0 and ( u v + u l ) 0 ; the last relation

follows from ( u v + u l ) ( u v u l ) = ( u v 2 u l 2 ) = 2 ( v v v l ) [ d ( p / T ) / d ( 1 / T ) ] 2

[ d ( μ / T ) / d ( p / T ) ( v v + v l ) / 2 ] 0 since d ( μ / T ) / d ( p / T ) ( v v + v l ) / 2 0 . Hence one gets

( u v u l ) 0 and ( u v + u l ) 0, (1.26)

( u v u l ) u l u ( u v + u l ) u c = 0 u v ( u v u l ) . (1.27)

The relations ( u v + u l ) = 2 u v ( u v u l ) 0 and ( u v + u l ) = 2 u l ( u v u l ) 0 lead to the limits of the energies u v , l in relation to the transient energies ± ( u v u l ) :

0 u v u v u l 1 2 u l ( u v u l ) 1. (1.28)

In other words: For 0 < T < T c the vapor internal energy uv is positive and

always lower than half the evaporation energy, u v 1 2 ( u v u l ) , and the liquid internal energy ul is negative and lower than half the condensation energy, u l 1 2 [ ( u v u l ) ] . The critical values of the entropy are obtained from

the relations s c T c = v c p c μ c = ( v c d p / d T d μ / d T ) | c T c , yielding

s ( v c , T c ) = v c d p d T | c d μ d T | c = v c p c T c μ c T c = 1 2 [ v c ( p T + d p d T ) c ( μ T + d μ d T ) c ] = f c T c > 0. (1.29)

Because 0 c ( v , T ) = d u ( v , T ) / d T = T d s ( v , T ) / d T , the functions u ( v , T ) and s ( v , T ) increase monotonically with increasing T, and so the critical values present the maximum internal energy and entropy of the saturated fluid.

The interdependence of entropy and internal energy is calculated from the equation s ( u l u v ) + u ( s v s l ) s v u l + s l u v = 0 and leads to the following identities and estimates:

s = u s v s l u v u l + u v s l u l s v u v u l 0 , u = s u v u l s v s l + s v u l s l u v s v s l 0. (1.30)

The estimates follow from [ u v M v / M + u l M l / M ] [ s v s l ] + [ M v / M + M l / M ] [ u v s l u l s v ] 0 (where the equality sign is valid for T = 0 ) and [ s v M v / M + s l M l / M ] [ u v u l ] + [ M v / M + M l / M ] [ s v u l s l u v ] 0 (where the equality sign is valid for T = T c ).

The two obviously equivalent Equations (1.30) present the opportunity for proving the correctness of the reference data mentioned in Equations (1.8), (1.13), (1.17), (1.21) and (1.27). For example, it follows from s ( v c , 0 ) = 0 that u ( v c , 0 ) = μ ( 0 ) = u l ( 0 ) = ( u v u l ) ( 0 ) with u v ( 0 ) = 0 and from s ( v c , T c ) = d μ / d T | c + v c d p / d T | c that u ( v c , T c ) = 0 . The data mentioned are thus reciprocal to one another. These data are thermodynamic reference values.

2.3. Heat Capacity Relations

The measurable heat capacity is defined by

c ( v , T ) = d u ( v , T ) d T = d s ( v , T ) d T T = s ( v , T ) + d s ( v , T ) d T T + d f ( v , T ) d T . (1.31)

Calculation of the specific heat capacities requires the temperature derivatives of the quantities given in Equations (1.23) and (1.25). Taking into consideration Equations (3.1)-(3.3) below, the heat capacity relations are

c ( T , v ) = d d T [ [ v v + v l v v v v l ln ( v v / v l ) ] d ( p / T ) d ( 1 / T ) ] = d 2 μ d T 2 T + v d 2 p d T 2 T 0 , c v , l ( T ) = c ( T , v ) + ( v v , l v ) ( c ( T , v ) / v ) T = d 2 μ d T 2 T + v v , l d 2 p d T 2 T 0. (1.32)

As the values of c ( T , v ) , s ( T , v ) , and d 2 p / d T 2 vanish at absolute zero, those of ( d μ / d T ) and ( d 2 μ / d T 2 ) also vanish there according to Equations (1.18) and (1.19). And as the value d 2 μ / d T 2 diverges at the critical point, the fluid heat capacity and the specific heat capacities also diverge there: c ( T c , v ) = c l ( T c ) = c v ( T c ) + .

Further expressions for the heat capacity and temperature derivatives of internal energies can be given as follows:

d u ( v , T ) d T = c ( v , T ) = c v v v l v v v l + c l v v v v v v l = [ d u v d T + d v v d T d ( p / T ) d ( 1 / T ) ] v v l v v v l + [ d u l d T + d v l d T d ( p / T ) d ( 1 / T ) ] v v v v v v l , (1.33)

d u v , l d T = T d 2 μ d T 2 + v v , l T d 2 p d T 2 d v v , l d T d ( p / T ) d ( 1 / T ) = c ( v , T ) + ( v v , l v ) c v c l v v v l + d v v , l d T u v u l v v v l = c v v v l v v v l + c l v v v v v v l + ( v v , l v ) c v c l v v v l + d v v , l d T u v u l v v v l . = c v , l + d v v , l d T u v u l v v v l = c v , l d v v , l d T d ( p / T ) d ( 1 / T ) , (1.34)

d ( u v u l ) d T = ( c v c l ) + d ( v v v l ) d T u v u l v v v l = ( v v v l ) d 2 p d T 2 T d ( v v v l ) d T d ( p / T ) d ( 1 / T ) ,

d ( u v + u l ) d T = 2 d 2 μ d T 2 T + ( v v + v l ) d 2 p d T 2 T d ( v v + v l ) d T d ( p / T ) d ( 1 / T ) > 0. (1.35)

It holds that ( u v u l ) ( 0 ) = u l ( 0 ) = u ( 0 ) = ( u v + u l ) ( 0 ) < u l ( T ) < u ( T ) < ( u v + u l ) ( T ) < u ( T c ) = ( u v + u l ) ( T c ) = 0 and for 0 < T < T c that

( u v + u l ) u = [ v v v l ln ( v v / v l ) v c ] [ d ( p / T ) d ( 1 / T ) ] > 0. (1.36)

In order to show that ( u v + u l ) is a convexly curved temperature function, one has to prove that d 2 ( u v + u l ) / d T 2 > 0 or that the straight line u ( 0 ) ( 1 T / T c ) is above ( u v + u l ) . The condition of convexity for ( u v + u l ) then reads ( u v + u l ) u ( 0 ) ( 1 T / T c ) , which can be transformed to

0 u v ( T ) + u l ( T ) u l ( 0 ) u l ( 0 ) T / T c . (1.37)

Indeed, for T > 0 the terms in Equation (1.37) are positive since u v ( T ) > 0 and u l ( T ) > u l ( 0 ) , which meets the condition mentioned. Likewise, the function ( u v u l ) ( T ) is concavely curved when its values are above the straight line u ( 0 ) ( 1 T / T c ) , which means that ( u v u l ) ( T ) u ( 0 ) ( 1 T / T c ) . Division by u ( 0 ) = ( u v u l ) ( 0 ) > 0 gives the correct relations 1 ( u v u l ) ( T ) / ( u v u l ) ( 0 ) ( 1 T / T c ) 0 and thus confirms that the condition for concavity of ( u v u l ) is met.

From c v , l = d u v , l / d T ( u v u l ) / ( v v v l ) d v v , l / d T and ( u v u l ) / ( v v v l ) = d ( p / T ) / d ( 1 / T ) it follows that

c v c l u v u l = d ln [ ( u v u l ) / ( v v v l ) ] d T = d ln [ d ( p / T ) / d ( 1 / T ) ] d T = d ( 1 / T ) d ( p / T ) d 2 p d T 2 T . (1.38)

It is of interest to take T = 0 in entropy and heat capacity relations. This is immediately possible in the case of the entropy because Nernst’s theorem states that the entropy vanishes at absolute zero, where only the condensed phase exists, and increases with the temperature:

s ( 0 , v ) = s l ( 0 ) = 0 , ( s ( T , v ) / T ) v = c ( T , v ) / T > 0. (1.39)

For 0 T T c one thus gets

0 = s ( 0 , v ) = s l ( 0 ) < s l ( T ) < s ( T , v ) < s ( T c , v ) = s l ( T c ) = s v ( T c ) < s v ( T ) . (1.40)

Similarly, in the case of the heat capacity one gets for 0 T < T c

0 = c ( 0 , v ) = c l ( 0 ) < c l ( T ) < c ( T , v ) < c v ( T ) . (1.41)

The temperature derivatives of the chemical potential function ( d μ / d T ) and the phase-specific entropies s v , l can be determined from measurements of c ( T , v ) , v v , l ( T ) and p ( T ) since

d 2 μ d T 2 = c ( T , v ) T v d 2 p d T 2 > 0 ,

d s v , l d T = c ( T , v ) T + ( v v , l v ) d 2 p d T 2 + d v v , l d T d p d T . (1.42)

With the result ( d μ / d T ) > 0 one obtains from Equations (1.23) the relations

s l v l d p d T = s v v v d p d T = s v d p d T = d μ d T > 0. (1.43)

These state that the entropy values are always positive and greater than the product of the volumes and vapor pressure coefficient. From relations (1.43), in turn, one can derive the following relations:

s l v l > s v > s v v v > d p d T > 0. (1.44)

Temperature properties of the phase-specific heats are similarly derived. With the result ( d 2 μ / d T 2 ) > 0 one obtains from Equations (1.32) the relations

c l v l d 2 p d T 2 T = c v d 2 p d T 2 T = c v v v d 2 p d T 2 T = d 2 μ d T 2 T > 0 , (1.45)

from which in turn the relations

c l v l > c v > c v v v > d 2 p d T 2 > 0 (1.46)

can be derived. Experimental verification of C V / V = c / v > d 2 p / d T 2 T also proves confirmation of ( d 2 μ / d T 2 ) > 0 .

2.4. Chemical Potential Relations

The identities

x v y l x l y v y v y l = x v , l + y v , l ( x y ) T , (1.47)

d d T x v y l x l y v y v y l = d x v , l d T + d y v , l d T ( x y ) T + y v , l d d T ( x y ) T

are now used to put the chemical potential functions in explicit form as energy functions:

μ = f v , l + v v , l p = u v , l s v , l T + v v , l p < 0 ,

d μ d T = s v , l + v v , l d p d T = d f v , l d T + d v v , l d T p + v v , l d p d T < 0 , (1.48)

d ( μ / T ) d ( 1 / T ) = u v , l + v v , l d ( p / T ) d ( 1 / T ) = d ( f v , l / T ) d ( 1 / T ) d v v , l d T T p + v v , l d ( p / T ) d ( 1 / T ) < 0 ,

d 2 μ d T 2 = c v , l T + v v , l d 2 p d T 2 = d s v , l d T + d v v , l d T d p d T + v v , l d 2 p d T 2 < 0.

The relations state that μ is a negative, concavely curved function, decreasing with increasing T. For μ as a function of measurable quantities see Equations (3.15) and (3.16) below.

3. Internal Fluid Energy as a Function of the Condensation Energy

It can be shown that the fluid energy u ( T , v ) can be expressed in terms of the condensation energy [ ( u v u l ) ] :

u = d ( μ / T ) d ( 1 / T ) v d ( p / T ) d ( 1 / T ) = d ( p / T ) d ( 1 / T ) [ d ( μ / T ) d ( p / T ) v ] = ( u v u l ) v v v l [ u v v l u l v v u v u l v ] = ( u v u l ) [ v l v v v v l u l u v u l ] = ( u v u l ) [ v v v v v v l u v u v u l ] . (2.1)

Since the temperature coefficients ( v v l ) / ( v v v l ) and ( v v v ) / ( v v v l ) are represented by M v / M and M l / M , respectively, the mass distributions M v / M = 0 and M l / M = 1 at T = 0 give the three relations u ( 0 ) = u l ( 0 ) , u ( 0 ) = ( u v u l ) ( 0 ) + u v ( 0 ) , u ( 0 ) = ( u v u l ) ( 0 ) , and the distributions M v / M = M l / M = 1 / 2 at T = T c give the three relations u ( T c ) = ( u v u l ) ( T c ) / 2 + u l ( T c ) , u ( T c ) = ( u v u l ) ( T c ) / 2 + u v ( T c ) , u ( T c ) = u v ( T c ) / 2 + u l ( T c ) / 2 . One thus obtains the following characteristic thermodynamic reference values, which are valid for every gas:

u ( 0 ) = u l ( 0 ) = ( u v u l ) ( 0 ) , u v ( 0 ) = 0 , u ( T c ) = u l ( T c ) = u v ( T c ) = 0. (2.2)

4. Internal Fluid Energy as an Expression of Measurable Quantities

Endeavors to publish data of the functions u ( v c , T ) and u v , l ( T ) are prominent in the current literature. The energy equations can be written in the form

u ( v c , T ) = [ d ( μ / T ) d ( p / T ) v c ] d ( p / T ) d ( 1 / T ) 0 , u v , l ( T ) = [ d ( μ / T ) d ( p / T ) v v , l ] d ( p / T ) d ( 1 / T ) . (3.1)

According to Equations (1.20)-(1.22) the diffential quotient d ( μ / T ) / d ( p / T ) is a positive volume quantity monotonically decreasing with increasing T from high values v v ( T ) near absolute zero to the lowest value vc at the critical point. The expression of measurable volumes v v ( T ) and v l ( T ) instead of d ( μ / T ) / d ( p / T ) as published in (2015) [3] reads:

d ( μ / T ) d ( p / T ) = v v + v l v v v l l n ( v v / v l ) . (3.2)

It is symmetric in the variables and linear in both vv and vl, and at the critical point it yields vc. Then inserting the solution (3.2) in the internal energy Equations (3.1) yields the results

u ( T , v c ) = [ v v v l ln ( v v / v l ) v v v v + v c ] [ d ( p / T ) d ( 1 / T ) ] 0 , u v , l ( T ) = [ v v v l ln ( v v / v l ) v l , v ] [ d ( p / T ) d ( 1 / T ) ] . (3.3)

These equations state that the internal energies can be expressed in terms of the measurable quantities, phase-specific volumes and vapor pressure, and are given by absolute figures; in particular, it holds that u v ( T c ) = u l ( T c ) = 0 since ( v v v l ) / ln ( v v / v l ) = v c and v v = v l = v c at T = T c .

If the functions uv and ul are expressed as dependent on the volume ratio z = v v / v l , one obtains for the energy ratios u v / ( u v u l ) and u l / ( u v u l ) according to Equations (1.28) the following relations:

0 u v u v u l = 1 ln z 1 z 1 1 2 for 0 T T c and z = v v ( T ) v l ( T ) 1 , 1 2 u l u v u l = 1 ln z + z z 1 1. (3.4)

The sum of them is, of course, equal to 1. Equations (3.4) allow one to calculate the relation between the energy ratio η = u v / u l and volume ratio z of fluid particles in the vapor and liquid spaces, viz.

0 η = u v u l = z 1 ln z z 1 z ln z 1 for 0 T T c and z = v v ( T ) v l ( T ) 1. (3.5)

According to Equations (1.17) and (1.8), the energy ratio η assumes the value 0 at absolute zero (where z ) and, respectively, the value -1 at the critical point (where z = 1 ). The relation η ( z ) represents a universal law of the two-phase equilibrium of real gases.

Calculation of the energy ratio ( u v + u l ) / ( u v u l ) as a function of z starts from Equation (1.20) or from ( η + 1 ) / ( η 1 ) and ends in any case with the result

1 u v + u l u v u l = 2 ln z z + 1 z 1 0 for 0 T T c and z = v v v l 1. (3.6)

Rigorous thermodynamic calculations combine Equations (1.8) and (1.44) and yield the following relations for the internal energy and entropy:

1 u v / u l < 0 , 1 s v / s l v v / v l for T T c . (3.7)

The constraints (3.7) state that the ratios of absolute energy and entropy for vapor and liquid are restricted within certain limits for temperatures in the two-phase region. And combining Equations (1.20), (1.22) and (3.2) gives

v v + v l u v v v u l v l u v u l = u v v l u l v v u v u l = v v + v l v v v l ln ( v v / v l ) , u v v v u l v l u v u l = v v v l ln ( v v / v l ) . (3.8)

Algebraic rearrangement of Equation (3.8) leads to a data criterion for consistent thermodynamic values v v ( T ) , v l ( T ) , u v ( T ) and u l ( T ) , which is valid for T T c and reads:

[ u v u v u l + v l v v v l ] ln ( v v v l ) = [ u l u v u l + v v v v v l ] ln ( v v v l ) = 1. (3.9)

In turn, from Equations (3.9) one obtains

u v = u v u l ln ( v v / v l ) v l u v u l v v v l 0 , u l = u v u l ln ( v v / v l ) v v u v u l v v v l 0 , ( u v u l ) = ( h v h l ) ( v v v l ) p = ( s v s l ) T ( v v v l ) p = ( v v v l ) d ( p / T ) d ( 1 / T ) 0. (3.10)

Equations (3.9) and (3.10) clearly state that, if table data T, p, vv, vl, ( h v h l ) , or ( s v s l ) and ( u v u l ) are thermodynamically consistent, then the internal energies, uv, ul, u = u v ( v c v l ) / ( v v v l ) + u l ( v v v c ) / ( v v v l ) , can be given absolutely. An entropy data criterion of the same kind as for the internal energy can also be formulated. To this end one has to confirm the validity of the relation.

s v , l + μ / T [ v v + v l ( v v v l ) / ln ( v v / v l ) ] p / T s v s l = u v , l u v u l (3.11)

to obtain with u v , l / ( u v u l ) = 1 / ln ( v v / v l ) v l , v / ( v v v l ) the criterion desired:

[ s v , l s v s l + v l , v v v v l + μ ( s v s l ) T [ v v + v l v v v l 1 ln ( v v / v l ) ] p / T d p / d T ] ln ( v v v l ) = 1. (3.12)

This relation allows, in principle, to give the chemical potential in terms of p, v v , l and s v , l as follows:

μ = s v , l T v l , v d p d T T + v v v l ln ( v v / v l ) d p d T T + [ v v + v l v v v l ln ( v v / v l ) ] p . (3.13)

The entropy expressions read

0 s ( v c , T ) = d μ d T + v c d p d T = s l + ( v c v l ) d p d T = s v ( v v v c ) d p d T . (3.14)

It holds that s v , l T v l , v d p d T T = s T ( v v + v l v c ) d p d T T and

μ = s T [ v v + v l v v v l ln ( v v / v l ) ] d p d T T + v c d p d T T + [ v v + v l v v v l ln ( v v / v l ) ] p .

According to relations (1.39) one has s T = T 0 T d s d T d T = T 0 T c T d T . The

chemical potential can thus be determined, on the one hand, by measuring the two-phase heat capacity, phase-specific volumes and vapor pressure:

μ ( T ) = T 0 T c ( T , v c ) T d T + v c d p d T T [ v v + v l v v v l l n ( v v / v l ) ] [ d ( p / T ) d ( 1 / T ) ] , (3.15)

or, on the other, by measuring the phase-specific volumes and vapor pressure only:

μ ( T ) = T [ μ ( T c ) T c p / T p ( T c ) / T c [ v v + v l v v v l l n ( v v / v l ) ] d p T ] < 0 , μ c T c = s c + v c p c T c . (3.16)

The energy sum ( μ + s T ) is also measurable and calculable and it holds that

μ + s T = [ v v v l l n ( v v / v l ) v v v l ] [ d ( p / T ) d ( 1 / T ) ] + v c d p d T T = u + v c p . (3.17)

From d ( μ + s T ) / d T = c + v c d p / d T > 0 and d 2 ( μ + s T ) / d T 2 = d c / d T + v c d 2 p / d T 2 > 0 it is concluded that the energy sum μ + s T = u + v c p is a convex temperature function strongly increasing from μ ( 0 ) = ( u v u l ) ( 0 ) < 0 at absolute zero to μ c + s c T c = v c p c > 0 at the critical point.

5. Phase-Specific Energy as a Function of the Evaporation Energy

The Carnot-Clapeyron-Clausius equation

u v u l = ( v v v l ) [ d ( p / T ) / d ( 1 / T ) ] (4.1)

suggests a unique relation between the energy density and particle density in the coexisting bulk phases, vapor and liquid. As the local interaction potentials in the partial volumens Vv and Vl are unique functions of the local particle density and determine the values uv and ul, respectively, it should be possible by means of a density coefficient and a temperature function to find suitable ansatzes for the functions u v , l depending on the density coefficient and temperature function. As density coefficient, a function ρ of the density variable z = v v / v l is now chosen, and as temperature function the evaporation energy ( u v u l ) . The ansatz proposed for the energy in the vapor phase is then u v = ρ ( 1 / z ) ( u v u l ) . The energy in the condensed phase must be u l = ρ ( z ) ( u v u l ) since the transfer of a particle from Vv to Vl causes the value uv to change to ul, the density variable from 1/z to z and the phase transition energy from ( u v u l ) to ( u l u v ) = ( u v u l ) , while the form of ρ is preserved because the functional density dependence of the microscopic interaction forces is of course phase-invariant. The conjectures [2]

u v = ρ ( 1 / z ) ( u v u l ) , u l = ρ ( z ) ( u v u l ) (4.2)

then yield, on the one hand, the reference values

0 = u v ( T c ) = u l ( T c ) = 1 2 [ u v ( T c ) + u l ( T c ) ] = u ( v c , T c ) and, on the other, for the density coefficient the equation

ρ ( 1 / z ) + ρ ( z ) = 1 (4.3)

and because u v ( 0 ) = 0 and u v ( T c ) / u l ( T c ) = 1 the boundary conditions

0 ρ ( 1 / z ) / ρ ( z ) 1. (4.4)

The functional Equation (4.3) for ρ under condition (4.4) is satisfied by

0 ρ ( 1 / z ) = 1 ln z 1 z 1 1 2 ρ ( z ) = 1 ln z + z z 1 1. (4.5)

From Equation (4.2) it follows that

0 u v u v u l = 1 ln z 1 z 1 1 2 1 ln z + z z 1 = u l u v u l 1. (4.6)

The relations (4.6) and (3.4) are identical and valid for 0 T T c and z = v v ( T ) / v l ( T ) 1 .

6. Results and Discussion

The paper treats thermodynamic properties of the saturated fluid. It is shown that the fluid state is completely determined by the internal energy and entropy in the vapor and liquid spaces. The ratios of the absolute phase-specific internal energies and entropies are restricted within certain limits. If for temperatures T [ T 1 , T 2 ] measured saturation data of the vapor pressure, reciprocal phase-specific densities and isothermal transient energy obey the data criterion Equation (3.9), then the internal energy as a function of T calculated according to Equations (3.10) is an absolute thermodynamic quantity. This fundamental procedure in gaining thermodynamic data excludes any application of so-called fiducial reference data since they cannot yield correct values. As the state of saturation is maintained, it is not possible to distinguish between a constant pressure and a constant volume condition. This is the characteristic difference between the heat capacities of a two-phase and a one-phase fluid. There is only one phase-specific heat capacity in the vapor space, c v ( T ) , and one in the liquid space, c l ( T ) . The measurable heat capacity of the fluid, C ( v c , T ) = M c ( v c , T ) = M v c v + M l c l , can also be calculated. Further results are the concavity of the positive, measurable temperature function ( u v u l ) and the convexity of the negative temperature functions ( u v + u l ) . The chemical potential is a negative, measurable and calculable temperature function. The ratios of phase-specific energies to isothermal transient

energies obey the relations 0 u v / ( u v u l ) 1 2 u l ( u l u v ) 1 , which should

be heeded when a state chart of the fluid under consideration is constructed.

Acknowledgements

The author is grateful to A. M. Nicol for the English translation.

Cite this paper

Elsner, A. (2018) Absolute Reference Values of the Real Gas. Engineering, 10, 270-290. https://doi.org/10.4236/eng.2018.105019

References

  1. 1. Gibbs, J.W. (1873) A Method of Geometrical Representation of the Thermodynamic Properties of Substances by Means of Surfaces. Transactions of the Connecticut Academy, 2, 382-404. http://www3.nd.edu/-powers/ame.20231/gibbs1873b.pdf

  2. 2. Elsner, A. (1988) Thermodynamic Representation of the State of a Saturated Fluid. Physics Letters A, 130, 225-233. https://doi.org/10.1016/0375-9601(88)90598-1

  3. 3. Elsner, A. (2015) Thermodynamic Equilibrium of the Satu-rated Fluid with a Free Surface Area and the Internal Energy as a Function of the Phase-Specific Volumes and Vapor Pressure. Engineering, 7, 577-596. https://doi.org/10.4236/eng.2015.79053

  4. 4. Elsner, A. (2012) Applied Thermodynamics of the Real Gas with Respect to the Thermodynamic Zeros of the Entropy and Internal Energy. Physica B: Physics of Condensed Matter, 407, 1055-1067. https://doi.org/10.1016/j.physb.2011.12.118

  5. 5. Strunk, Ch. (2015) Moderne Thermodynamik: Von einfachen Systemen zu Nanostrukturen. 32. https://doi.org/10.1515/9783110371062

  6. 6. Lauth, J.G. and Kowalczyk, J. (2015) Thermodynamik. Springer-Verlag, Berlin Heidelberg, 67. https://doi.org/10.1007/978-3-662-46229-4

  7. 7. Callen, H.B. (1960) Thermodynamics. John Wiley & Sons, Hoboken, 52.

  8. 8. Johnson, V.J. (1961) Properties of Materials at Low Temperature (Phase 1). National Bureau of Standards, Cryogenic Engineering Laboratory, Pergamon Press, Oxford, London, New York, Paris.

  9. 9. Journal of Physical and Chemical Reference Data, Annual Series of Properties of Gases, published since 1972.

  10. 10. IUPAC Commission on Thermodynamics (1972) International Thermodynamic Tables of the Fluid State. Blackwell Scientific Publications, Hoboken, Vol. 1.

  11. 11. Selover, T.B. (1987) A Series of Property Tables. Hemisphere Publishing Corporation, New York.

  12. 12. Reynolds, W.C. (1979) Thermodynamic Properties in SI. Graphs, Tables and Computational Equations for 40 Substances. The Department of Mechanical Engineering, Stanford University, Stanford.

  13. 13. Grigull, U. and Schmidt, E. (1989) Properties of Water and Steam in SI-Units. Springer-Verlag, Berlin, Heidelberg, New York.

  14. 14. Wagner, W. and Pruss, A. (1993) International Equations for the Saturation Properties of Ordinary Water Substance. Revised According to the International Temperature Scale of 1990. Addendum to J. Phys. Chem. Ref. Data 16, 893 (1987). Journal of Physical and Chemical Reference Data, 22, 783. https://doi.org/10.1063/1.555926

  15. 15. Kohlrausch, F. (1996) Praktische Physik. Band 3, B. G. Teubner Stuttgart.

  16. 16. Feistel, R. and Wagner, W. (2006) A New Equation of State for H2O Ice Ih. Journal of Physical and Chemical Reference Data, 35, 1021-1047. https://doi.org/10.1063/1.2183324

  17. 17. Elsner, A. (2014) The Dominant Role of the Chemical Potential for Driving Currents in Oceans and Air. Journal of Geoscience and Environment Protection, 2, 117-125.

  18. 18. Pethick, C.J. and Smith, H. (2002) Bose-Einstein Condensation in Dilute Gases. Cambridge University Press, Cambridge.

  19. 19. Elsner, A. (2017) Absolute Internal Energy of the Real Gas. Engineering, 9, 361-375.

Appendix: Comments on the Internationally Accepted Equations for the Saturation Properties of Water

The International Association for the Properties of Water and Steam (IAPWS) [13] [14] provides internationally accepted formulations for the properties of water. There are special correlation equations for the vapor-liquid saturation properties of water. Formulas are given for the vapor pressure p, phase-specific volumes v v , l , internal energies u v , l , and entropies s v , l as functions of the saturation temperature T. This affords a unique description of the temperature dependence of every property of the saturated water.

This study treats the IAPWS equations in the framework of thermodynamics. The IAPWS equations read, on the one hand [13] [14] ,

u v , l = α v v , l d ( p / T ) d ( 1 / T ) , s v , l = ϕ + v v , l d p d T , (A1)

where the so-called auxiliary quantities α and ϕ are given as functions of the temperature T; and, on the other, in terms of the thermodynamic fundamental equation,

u v , l s v , l T = μ v v , l p , (A2)

where μ is the chemical potential of the saturated fluid. Equations (A1) and (A2) are related as follows:

α ϕ T = μ , α = d ( μ / T ) d ( 1 / T ) , ϕ = d μ d T . (A3)

Let us now investigate the thermodynamic conditions that have to be satisfied by the temperature functions α and ϕ if they are to define the two-phase chemical potential μ according to Equations (A3). First a few thermodynamic relations are taken and then transformed into IAPWS parlance. The Carnot-Clapeyron-Clausius equations read

u v u l v v v l = d ( p / T ) d ( 1 / T ) , s v s l v v v l = d p d T . (A4)

These are satisfied by Equations (A1). The Gibbs-Duhem equations read

d u v , l d T + d v v , l d T p = d s v , l d T T . (A5)

They are satisfied by Equations (A1) if the following relation between α and ϕ is valid:

d α d T = d ϕ d T T . (A6)

Note that Equation (A6) is implicitly contained in relations (A3) since d μ / d T = d α / d T d ϕ / d T T ϕ = ϕ and d ( μ / T ) / d ( 1 / T ) = μ d μ / d T T = ( α ϕ T ) + ϕ T = α . Relations (A4) and (A5) therefore afford nothing new, but merely confirm relations (A3). The mathematical structures of Equations (A1) thus conform to the thermodynamic internal energy and entropy expressions in respect of Equations (A3).

Condition (A6) is numerically satisfied by the formulas in [13] [14] for α ( T ) and ϕ ( T ) , whose temperature dependences are expressed as follows:

α [ J / g ] = d α + d 1 θ 19 + d 2 θ + d 3 θ 4.5 + d 4 θ 5 + d 5 θ 54.5 , where θ T T c , (A7)

ϕ [ J / g K ] = [ d ϕ + 19 20 d 1 θ 20 + d 2 ln θ + 9 7 d 3 θ 3.5 + 5 4 d 4 θ 4 + 109 107 d 5 θ 53.5 ] 1 T c , (A8)

where [ d α , d ϕ ] = [ 1.135905627715 E 3 , 2.3195246 E 3 ] , T c = 647.096 K ,

[ d 1 , d 2 , d 3 , d 4 , d 5 ] = [ 5.65134998 E 8 , 2.69066631 E 3 , 1.27287297 E 2 , 1.35003439 E 2 , 9.81825814 E 1 ] .

The vapor pressure p and chemical potential μ obey for temperatures T T c the following relations:

p > 0 , d p d T > 0 , d ( p / T ) d ( 1 / T ) < 0 , d 2 p d T 2 > 0 , (A9)

μ < 0 , d μ d T < 0 , d ( μ / T ) d ( 1 / T ) < 0 , d 2 μ d T 2 < 0. (A10)

Relations (A10) then lead to further conditions for α and ϕ:

α < 0 , d α d T > 0 , ϕ > 0 , d ϕ d T > 0. (A11)

Accordingly, α must be a negative function increasing as T, and ϕ must be a positive function likewise increasing as T. It is found that d α / d T > 0 , but not that α < 0 for all T in the range [ T t , T c ] ; similarly, it does hold that d ϕ / d T > 0 , but not that ϕ > 0 for all T [ T t , T c ] . The auxiliary Equations (A1), despite their correct formal structure, are therefore not thermodynamically appropriate for justifying the chemical potential according to Equations (A3). Further consideration of the temperature dependence of u v , l shows that the sum ( u v + u l ) yields a concavely, but not a convexly curved function, as should be.

The fit function for measured vapor pressure data reads [13] [14]

ln ( p p c ) = T c T [ a 1 τ + a 2 τ 1.5 + a 3 τ 3 + a 4 τ 3.5 + a 5 τ 4 + a 6 τ 7.5 ] , where τ 1 θ (A12)

and p c = 22.064 [ MPa ] . The terms τ 1.5 , τ 3.5 , τ 7.5 produce divergent temperature derivatives d n p / d T n at T c for n 2 , whereas all vapor pressure derivatives have finite values. In contrast, the derivatives of the chemical potential d n μ / d T n diverge at T c for n 2 and determine the divergence of the heat capacity C ( M , V , T ) = M c ( V / M , T ) . In addition, it should be noted here that, if the saturation state is maintained, it is not possible to distinguish between a constant pressure and a constant volume condition. There is thus only one phase-specific heat capacity in the vapor space, c ( V v / M v , T ) , and one in the liquid space, c ( V l / M l , T ) , and it holds that C ( M , V , T ) = M v c ( V v / M v , T ) + M l c ( V l / M l , T ) = C ( M v , V v , T ) + C ( M l , V l , T ) .

Numerous thermodynamic deficiencies have been mentioned, viz. violation of the conditions α ( T ) < 0 and ϕ ( T ) > 0 for every temperature T [ T t , T c ] , the incorrect temperature dependence of ( u v + u l ) , the proposed vapor pressure fit formula with a divergent term d 2 p / d T 2 at T c and, finally, tables which list two different phase-specific heat capacities for vapor as well for liquid. From all this it is concluded that the published data [13] [14] are in need of basic correction.