Published on *PNG 520: Phase Behavior of Natural Gas and Condensate Fluids* (https://www.e-education.psu.edu/png520)

Any function f(x) that possesses the characteristic mapping:

$$\begin{array}{l}x\to \lambda x\\ f(x)\to \lambda f(x)\end{array}$$

is said to be *homogeneous,* with respect to x, to degree 1. By the same token, if f(x) obeys the mapping:

$$\begin{array}{l}x\to \lambda x\\ f(x)\to {\lambda}^{k}f(x)\end{array}$$

then f(x) is homogeneous to degree “k”. In general, a multivariable function f(x_{1},x_{2},x_{3},…) is said to be homogeneous of degree “k” in variables x_{i}(i=1,2,3,…) if for any value of $\lambda $ ,

$$f(\lambda {x}_{1},\lambda {x}_{2},...)={\lambda}^{k}f({x}_{1},{x}_{2},...)$$

For example, let us consider the function:

$$f(x,y)=\frac{x}{{x}^{2}+{y}^{2}}$$

How do we find out if this particular function is homogeneous, and if it is, to what degree? We evaluate this function at $x=\lambda x\text{}and\text{}y=\text{}\lambda y$ to obtain:

$$f(\lambda x,\lambda y)=\frac{\lambda x}{{\lambda}^{2}{x}^{2}+{\lambda}^{2}{y}^{2}}={\lambda}^{-1}\frac{x}{{x}^{2}+{y}^{2}}{\lambda}^{-1}f(x,y)$$

hence, the function f(x,y) in (15.4) is homogeneous to degree -1.

In regard to thermodynamics, extensive variables are homogeneous with degree “1” with respect to the number of moles of each component. They are, in fact, proportional to the mass of the system to the power of one (k=1 in equation 15.2 or 15.3). This is, if we triple the amount of mass in the system, the value of any given extensive property will be tripled as well. Notice that this is not the case for intensive properties of the system (such as temperature or pressure), simply because they are independent of mass. Hence, intensive thermodynamic properties are homogeneous functions with degree “0” — in such a case, k=0 in equation (15.2) or (15.3).

From the previous section, it is clear that we are not only interested in looking at thermodynamic functions alone, but that it is also very important to compute how thermodynamic functions change and how that change is mathematically related to their partial derivatives $\frac{\partial f}{\partial x},\frac{\partial f}{\partial y},\text{and}\frac{\partial f}{\partial z}$. Hence, to complete the discussion on homogeneous functions, it is useful to study the mathematical theorem that establishes a relationship between a homogeneous function and its partial derivatives. This is *Euler’s theorem.*

Euler’s theorem states that if a function f(a_{i}, i = 1,2,…) is homogeneous to degree “k”, then such a function can be written in terms of its partial derivatives, as follows:

$$k{\lambda}^{k-1}f({a}_{i})={\displaystyle \sum _{i}^{}{a}_{i}}{\left(\frac{\partial f({a}_{i})}{\partial (\lambda {a}_{i})}\right)|}_{\lambda x}$$

Since (15.6a) is true for all values of $\lambda $ , it must be true for $\lambda -1$ . In this case, (15.6a) takes a special form:

$$kf({a}_{i})={\displaystyle \sum _{i}^{}{a}_{i}}{\left(\frac{\partial f({a}_{i})}{\partial ({a}_{i})}\right)|}_{x}$$

So far, so good. But…what is the application of all this? Well, first of all, we have to know something more about extensive thermodynamic properties. A very neat thing about them is that they can be written as a function of a sufficient number of independent variables to completely define the thermodynamic state of the system. Such a set is said to be a *complete* set. As it turns out, any thermodynamic system is completely defined when both the masses of all the substances within it are defined and two additional independent variables are fixed. This is *Duhem’s theorem*. From a real-life perspective, it is natural to choose pressure and temperature as those “independent variables” — physical quantities that we have a “feel” for and we think we can control — rather than specific volume or entropy. As we will see later, they are also convenient variables of choice because they are homogeneous of degree zero in mass.

Let “$\Im $ ” be a given extensive property of a multi-component system. From the previous section, we know that the value of “$\Im $ ” must be fixed and uniquely determined once we fix the pressure, temperature, and number of moles of each component in the system. This is,

$$\Im =\Im (P,T,{n}_{1},{n}_{2},...,{n}_{N})$$

Additionally, we recall that extensive properties are homogeneous of degree one *with respect to number of moles* and homogeneous of degree zero with respect to pressure and temperature. Thus, expression (15.6b) is readily applicable:

$$\Im ={\displaystyle \sum _{i}{n}_{i}}{\left(\frac{\partial \Im}{\partial {n}_{i}}\right)}_{P,T,{n}_{i=1}}={\displaystyle \sum _{i}{n}_{i}{\Im}_{i}}$$

where we have just defined:

$${\Im}_{i}{\left(\frac{\partial \Im}{\partial {n}_{i}}\right)}_{P,T,{n}_{i=1}}$$

Equation (15.7c) is a very important definition. It defines the concept of a *partial molar quantity*. A partial molar quantity $\overline{{\Im}_{i}}$ represents the change in the total quantity ($\Im $ ) due to the addition of an infinitesimal amount of species “i” to the system *at constant pressure and temperature*. Euler’s theorem gave birth to the concept of *partial molar quantity* and provides the functional link between it (calculated for each component) and the total quantity. The selection of *pressure and temperature* in (15.7c) was not trivial. First, they are convenient variables to work with because we can measure them in the lab. But most important, they are* intensive variables*, homogeneous functions of degree zero in number of moles (and mass). This allowed us to use Euler’s theorem and jump to (15.7b), where only a summation with respect to number of moles survived. The definition of the partial molar quantity followed.

The conventional notation we are going to follow throughout the following section is:

$\Im $ = Total quantity (e.g., total volume, total internal energy, etc),

$\tilde{\Im}$ = Molar quantity, i.e., total quantity per unit mole:

$$\tilde{\Im}=\frac{\Im}{n}(\text{foramixture,}n={\displaystyle \sum _{i}{n}_{i}})$$

$\overline{\Im}$ = Partial molar quantity,

$\widehat{\Im}$ = Mass or specific quantity, i.e., total quantity per unit mass:

$$\widehat{\Im}=\frac{\Im}{m}(\text{foramixture,}m={\displaystyle \sum _{i}{m}_{i}})$$

We can rewrite equation (15.7b) in terms of molar quantity using the definition in (15.8a),

$$\tilde{\Im}=\frac{1}{n}{\displaystyle \sum _{i}{n}_{i}}{\overline{\Im}}_{i}={\displaystyle \sum _{i}{x}_{i}}{\overline{\Im}}_{i}$$

where:

${x}_{i}={\text{molarfractionofspecies"i"=n}}_{i}/n$

Any molar quantity in thermodynamics can be written in terms of the partial molar quantity of its constituents. If we set $\Im =G$ , we end up with,

$$\tilde{G}={\displaystyle \sum _{i}{x}_{i}}{\overline{G}}_{i}$$

Equivalently, if we set $\Im =V$ (total volume),

$$\tilde{v}={\displaystyle \sum _{i}{x}_{i}}{\overline{v}}_{i}$$

Notice that for single component systems (x_{i}=1), *partial molar properties are just equal to the molar property*:

E$$\tilde{G}={\overline{G}}_{i}{\text{(ifx}}_{i}=1,i=1)$$

$$\tilde{v}={v}_{i}{\text{(ifx}}_{i}=1,i=1)$$

This is also a consequence of the definition in (15.7c),

$${\overline{\Im}}_{i}={\left(\frac{\partial \Im}{\partial {n}_{i}}\right)}_{P,T,{n}_{i=1}}$$

For a pure component, $n={n}_{i}\text{,i=1,}\Im \text{=n}\tilde{\Im}$ and:

$$\overline{\Im}={\left(\frac{\partial n\tilde{\Im}}{\partial n}\right)}_{PI}=\tilde{\Im}$$

The reason for the introduction of the concept of a *partial molar quantity* is that often times we deal with mixtures rather than pure-component systems. The way to characterize the state of the mixtures is via partial molar properties. This concept provides the bridge between the thermodynamics of systems of constant composition, which we have studied so far, and the thermodynamics of systems of variable composition, which we will deal with in the next section. Basically, the definition in (15.7c):

$${\overline{\Im}}_{i}={\left(\frac{\partial \Im}{\partial {n}_{i}}\right)}_{P,T,{n}_{i=1}}$$

allows us to quantify how the total, extensive property $\Im $
changes with additions of n_{i} at constant pressure and temperature. If you look at (15.7b) and (15.9), you will also realize that (15.7c) is just an *allocation formula* that allows assigning to each species “i” a share of the total mixture property, such that:

$$\Im ={\displaystyle \sum _{i}{n}_{i}{\overline{\Im}}_{i}}$$

$$\tilde{\Im}={\displaystyle \sum _{i}{x}_{i}{\overline{\Im}}_{i}}$$

We can play with “$\Im $
” a little more. Let us say that we are now interested in looking at the differential *changes* of $\Im $. Since we also know that $\Im $ is a state function, and given the functional relationship in (15.7a), the total differential for $\Im $ is written:

$$d\Im ={\left(\frac{\partial \Im}{\partial P}\right)}_{T{n}_{i}}dP+{\left(\frac{\partial \Im}{\partial T}\right)}_{P,{n}_{i}}dT+{\displaystyle \sum _{i}{\left(\frac{\partial \Im}{\partial {n}_{i}}\right)}_{P,T,n=1}d{n}_{i}}$$

or

$$d\Im ={\left(\frac{\partial \Im}{\partial P}\right)}_{T{n}_{i}}dP+{\left(\frac{\partial \Im}{\partial T}\right)}_{P,{n}_{i}}dT+{\displaystyle \sum _{i}{\overline{\Im}}_{i}}d{n}_{i}$$

Basically, equations (15.13) tell us that any change in P, T, or n_{i} will cause a corresponding change in the total property, $\Im $
. This is a reinforcement of what is explicitly declared in (15.7a). If we recall (15.7b), an alternate expression for the total differential in (15.13) is written:

$$d\Im =d\left({\displaystyle \sum _{i}{n}_{i}{\overline{\Im}}_{i}}\right)={\displaystyle \sum _{i}{\overline{\Im}}_{i}}d{n}_{i}+{\displaystyle \sum _{i}{n}_{i}d{\overline{\Im}}_{i}}$$

If we subtract (15.14) from (15.13b), we get:

$$d\Im -d\Im ={\left(\frac{\partial \Im}{\partial P}\right)}_{T,{n}_{i}}dP+{\left(\frac{\partial \Im}{\partial T}\right)}_{P,{n}_{i}}dT+{\displaystyle \sum _{i}{\overline{\Im}}_{i}d{n}_{i}-}{\displaystyle \sum _{i}{\overline{\Im}}_{i}d{n}_{i}-{\displaystyle \sum _{i}{n}_{i}d{\overline{\Im}}_{i}}}$$

Therefore,

$${\left(\frac{\partial \Im}{\partial P}\right)}_{T,{n}_{i}}dP+{\left(\frac{\partial \Im}{\partial T}\right)}_{P,{n}_{i}}dT-{\displaystyle \sum _{i}{n}_{i}d{\overline{\Im}}_{i}}=0$$

Equation (15.16) is the well-known Gibbs-Duhem equation. It can be applied to any extensive thermodynamic property: U, S, H, G, A, and it must hold true. It represents a thermodynamic constraint between the intensive variables P, T and $\tilde{\Im}$ . Pressure, temperature and partial molar properties cannot vary in just any fashion; any change taking place among them must satisfy (15.16). The change in any one of them can be calculated as a function of the change in the other two by means of the Gibbs-Duhem equation. This equation is the basis for thermodynamic consistency checks of experimental data.