Residual property (physics) explained

In thermodynamics a residual property is defined as the difference between a real fluid property and an ideal gas property, both considered at the same density, temperature, and composition, typically expressed as

X(T,V,n)=Xid(T,V,n)+Xres(T,V,n)

where

X

is some thermodynamic property at given temperature, volume and mole numbers,

Xid

is value of the property for an ideal gas, and

Xres

is the residual property. The reference state is typically incorporated into the ideal gas contribution to the value, as

Xid(T,V,n)=X\circ,(T,n)+\DeltaidX(T,V,n)

where

X\circ,

is the value of

X

at the reference state (commonly pure, ideal gas species at 1 bar), and

\DeltaidX

is the departure of the property for an ideal gas at

(T,V,n)

from this reference state.

Residual properties should not be confused with excess properties, which are defined as the deviation of a thermodynamic property from some reference system, that is typically not an ideal gas system. Whereas excess properties and excess models (also known as activity coefficient models) typically concern themselves with strictly liquid-phase systems, such as smelts, polymer blends or electrolytes, residual properties are intimately linked to equations of state which are commonly used to model systems in which vapour-liquid equilibria are prevalent, or systems where both gases and liquids are of interest. For some applications, activity coefficient models and equations of state are combined in what are known as "

\gamma

-

\phi

models" (read: Gamma-Phi) referring to the symbols commonly used to denote activity coefficients and fugacities.

Significance

In the development and implementation of Equations of State, the concept of residual properties is valuable, as it allows one to separate the behaviour of a fluid that stems from non-ideality from that stemming from the properties of an ideal gas. For example, the isochoric heat capacity is given by

CV=\left(

\partialU
\partialT

\right)V,n=T\left(

\partialS
\partialT

\right)V,n=T\left[\left(

\partialSid
\partialT

\right)V,n+\left(

\partialSres
\partialT

\right)V,n\right] =

id
C
V

+

res
C
V

Where the ideal gas heat capacity,

id
C
V
, can be measured experimentally, by measuring the heat capacity at very low pressure. After measurement it is typically represented using a polynomial fit such as the Shomate equation. The residual heat capacity is given by

T\left(

\partialSres
\partialT

\right)V,=-T\left(

\partial2Ares
\partialT2

\right)V,

,

and the accuracy of a given equation of state in predicting or correlating the heat capacity can be assessed by regarding only the residual contribution, as the ideal contribution is independent of the equation of state.

In Equilibrium Calculations

In fluid phase equilibria (i.e. liquid-vapour or liquid-liquid equilibria), the notion of the fugacity coefficient is crucial, as it can be shown that the equilibrium condition for a system consisting of phases

\alpha

,

\beta

,

\gamma

, ... the condition for chemical equilibrium is
\alpha
x
i
\alpha
\Phi
i

=

\beta
x
i
\beta
\Phi
i

=

\gamma
x
i
\gamma
\Phi
i

=...

for all species

i

, where
j
x
i
denotes the mole fraction of species

i

in phase

j

, and
j
\Phi
i
is the fugacity coefficient of species

i

in phase

j

. The fugacity coefficient, being defined by

\mui=

\circ
\mu
i

+RTln

\Phiixip
p\circ

is directly related to the residual chemical potential, as

\mui=

id
\mu
i

+

res
\mu
i

=

\circ
\mu
i

+RTln

xip
p\circ

+

res
\mu
i

\implies

res
\mu
i

=RTln\Phii

,

thus, because

res
\mu
i

=\left(

\partialAres
\partialni

\right)T,

, we can see that an accurate description of the residual Helmholtz energy, rather than the total Helmholtz energy, is the key to accurately computing the equilibrium state of a system.

Residual Entropy Scaling

The residual entropy of a fluid has some special significance. In 1976, Yasha Rosenfeld published a landmark paper, showing that the transport coefficients of pure liquids, when expressed as functions of the residual entropy, can be treated as monovariate functions, rather than as functions of two variables (i.e. temperature and pressure, or temperature and density).[1] This discovery lead to the concept of residual entropy scaling, which has spurred a large amount of research, up until the modern day, in which various approaches for modelling transport coefficients as functions of the residual entropy have been explored.[2] Residual entropy scaling is still very much an area of active research.

Dependence on variable set

While any real state variable

X

, in a real state (

T,V,p,n

), is independent of whether one evaluates

X(T,p,n)

or

X(T,V,n)

, one should be aware that the residual property is in general dependent on the variable set, i.e.

Xres(T,p,n)Xres(T,V,n)

This arises from the fact that the real state

(T,V,p,n)

is in general not a valid ideal gas state, such that the ideal part of the property will be different depending on variable set. Take for example the chemical potential of a pure fluid: In a state

(T,V,p,n)

that does not satisfy the ideal gas law, but may be a real state for some real fluid. The ideal gas chemical potential computed as a function of temperature, pressure and mole number is

\muid(T,p,n)=\mu\circ+RTln

p
p\circ
,

while computing it as a function of concentration (

c=n/V

), we have

\muid(T,V,n)=\mu\circ+RTln

c
c\circ
,

such that

\muid(T,p,n)-\muid(T,V,n)=RTln

p
p\circ

-RTln

c
c\circ

=RTln

pV
nRT

=RTlnZ

,

where we have used

p\circ=c\circRT

, and

Z

denotes the compressibility factor. This leads to the result

\mui(T,p,n)-\mui(T,V,n)=0\implies

res
\mu
i

(T,V,n)-

res
\mu
i

(T,p,n)=RTlnZ

.

Practical Calculation

In practice, the most significant residual property is the residual Helmholtz energy. The reason for this is that other residual properties can be computed from the residual Helmholtz energy as various derivatives (see: Maxwell relations). We note that

\left(\partialA
\partialV

\right)T,=\left(

\partialAid
\partialV

\right)T,+\left(

\partialAres
\partialV

\right)T,\iff\left(

\partialAres
\partialV

\right)T,=\left(

\partialA
\partialV

\right)T,-\left(

\partialAid
\partialV

\right)T,=-p(T,V,n)-(-pid(T,V,n))

such thatA^(T, V', n) - A^(T, V = \infty, n)= \int_\infty ^ \left(\frac \right) dV = \int_\infty ^ p^(T, V, n) - p(T, V, n) d Vfurther, because any fluid reduces to an ideal gas in the limit of infinite volume,

A(T,V=infty,n)=Aid(T,V=infty,n)\iffAres(T,V=infty,n)=0

.

Thus, for any Equation of State that is explicit in pressure, such as the van der Waals Equation of State, we may compute

Ares(T,V,n)=

V
\int
infty
nRT
V'

-p(T,V',n)dV'

.

However, in modern approaches to developing Equations of State, such as SAFT, it is found that it can be simpler to develop the equation of state by directly developing an equation for

Ares

, rather than developing an equation that is explicit in pressure.

Correlated terms

References

Notes and References

  1. Rosenfeld . Yaakov . 1977-06-01 . Relation between the transport coefficients and the internal entropy of simple systems . Physical Review A . 15 . 6 . 2545–2549 . 10.1103/PhysRevA.15.2545.
  2. Bell . Ian H. . Hellmann . Robert . Harvey . Allan H. . 2020-03-12 . Zero-Density Limit of the Residual Entropy Scaling of Transport Properties . Journal of Chemical & Engineering Data . en . 65 . 3 . 1038–1050 . 10.1021/acs.jced.9b00455 . 0021-9568.