PNG 520
Phase Relations in Reservoir Engineering

Fugacity

PrintPrint

We have seen that, for a closed system, the Gibbs energy is related to pressure and temperature as follows:

dG=VdPSdT This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.15)

For a constant temperature process,

dG=VdP @ constant T This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.16)

For an ideal gas,

dG= RT P dP This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.17)
dG=RTd In P @ constant T This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.18)

This expression by itself is strictly applicable to ideal gases. However, Lewis, in 1905, suggested extending the applicability of this expression to all substances by defining a new thermodynamic property called fugacity, f, such that:

dG=RTd In f @ constant T This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.19)

This definition implies that for ideal gases, ‘f’ must be equal to ‘P’. For mixtures, this expression is written as:

d G i ¯ =RTd  In f i  @ constant T This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.20)

where G i ¯  and f i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. are the partial molar Gibbs energy and fugacity of the i-th component, respectively. Fugacity can be readily related to chemical potential because of the one-to-one relationship of Gibbs energy to chemical potential, which we have discussed previously. Therefore, the definition of fugacity in terms of chemical potential becomes:

For a pure substance,

d In f= dμ RT @ constant T This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.21a)
lim P0  f=P (ideal gas limit) This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.21b)

For a component in a mixture,

d In   f i =  d μ i RT @ constant T This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.22c)
lim P0   f i = y i P = partial pressure (ideal gas limit) This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.22d)

The fugacity coefficient ( ϕ i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. ) is defined as the ratio of fugacity to its value at the ideal state. Hence, for pure substances:

ϕ= f P This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.23a)

and for a component in a mixture,

ϕ i = f i y i P This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.23b)

The fugacity coefficient takes a value of unity when the substance behaves like an ideal gas. Therefore, the fugacity coefficient is also regarded as a measure of non-ideality; the closer the value of the fugacity coefficient is to unity, the closer we are to the ideal state.

Fugacity turns out to be an auxiliary function to chemical potential. Even though the concept of thermodynamic equilibrium which we discussed in the previous section is given in terms of chemical potentials, above definitions allow us to restate the same principle in terms of fugacity. To do this, previous expressions can be integrated for the change of state from liquid to vapor at saturation conditions to obtain:

l v d In  f i = 1 RT l v d μ i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.24a)
In  f i (v)  In  f i (v) = 1 RT ( μ i (v) μ i (l) ) This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.24b)

For equilibrium, μ i (l) = μ i (v) This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. ,hence,

In( f i (v) f i (l) )=0 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.24c)

Therefore:

f i (l) = f i (v); i=1, 2, ...  n c This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(16.25)

For equilibrium, fugacities must be the same as well! This is, for a system to be in equilibrium, both the fugacity and the chemical potential of each component in each of the phases must be equal. Conditions (16.14) and (16.25) are equivalent. Once one of them is satisfied, the other is satisfied immediately. Using μ i (l) = μ i (v)  or   f i (l) = f i (v) This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. to describe equilibrium is a matter of choice, but generally the fugacity approach is preferred.