PNG 520
Phase Relations in Reservoir Engineering

Thermodynamics of Systems of Variable Composition (Open Multicomponent Systems)

PrintPrint

To extend all of the previous concepts to systems of variable mass, we must now consider at least one new variable: number of moles, n. To account for this effect, you will see below that we will have to introduce a “new” thermodynamic property. Let us take the case of the internal energy. For a constant composition system, we wrote:

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

and,

dU= ( U S ) V dS+ ( U V ) S dV This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(14.22a)
 

If this system has a single component, and now we allow its mass (and hence, number of moles, n) to change (an open system), the change in U (dU) is no longer just a function of dS and dV. We now have to account for changes in ‘n’, thus:

dU= ( U S ) V,n dS+ ( U V ) S,n dV+ ( U n ) S,V dn This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.17)
 

If we had a binary system, we would have two new variables, ‘n1’ and ‘n2’ (n=n1+n2) and we would have to expand (15.17) accordingly,

dU= ( U S ) V,n dS+ ( U V ) S,n dV+ ( U n 1 ) S,V, n 2 d n 1 + ( U n 2 ) S,V, n 1 d n 2 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.18)
 

Hence, for a multi-component system, we just keep on adding terms ( n= i n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. ):

dU= ( dU dS ) V,n dS+ ( dU dV ) S,n dV+ i ( n i ) S,V, n i n 1 d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.19)
 

Thermodynamics defines this ‘coefficient’ which multiplies the change in the number of moles of each component (dni) as the ‘chemical potential’ of that component(μi). See how the chemical potential is a thermodynamic property that must be defined for the proper description of a system of variable composition — that is, an open system.

Then, we write:

μ i = ( dU d n i ) S,V, n i n 1 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.20)
 

and finally,

dU=TdSPdV+ i μ i d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.21)
 

The same can be done with the thermodynamic definitions for dH, dG, and dA (the rest of equations 14.22). In fact, the chemical potential may be defined in at least four different and equivalent ways. You can now show that for equations (14.22) to account for systems of variable composition, we would have to expand them into:

dH= ( H S ) P,n dS= ( H P ) S,n dP+ i ( H n i ) S,P,n i n 1 d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.22a)
 
dG= ( G P ) T,n dP= ( G T ) P,n dT+ i ( G n i ) P,T,n i n 1 d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.22b)
 
dA= ( A V ) T,n dS= ( A T ) V,n dT+ i ( A n i ) T,V,n i n 1 d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.22c)
 

Because thermodynamics defines the ‘coefficient’ which multiplies the change in the number of moles of each component (dni) as the ‘chemical potential’ of that component( μ i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. ), we have three new ways to express it:

μ i = ( H n i ) S,P,n i n j This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.23a)
 
μ i = ( G n i ) P,T,n i n j This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.23b)
 
μ i = ( A n i ) T,V,n i n j This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.23c)
 

Don’t any of these three equations ring a bell? Don’t forgot the definition of (15.7c). Of all the four available definitions for chemical potential, there is one that suits our definition of a partial molar quantity perfectly. Compare equation (15.23b) to (15.7c). What we can say is that the chemical potential of a component “i” is equal to the partial molar Gibbs energy of such component:

μ i = ( G n i ) P,T,n i n j = G ¯ l This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.24)
 

Notice that, for a pure component, the chemical potential is equal to the molar Gibbs energy of the substance (see equation 15.9),

E μ= G ˜ This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.25)
 

When we introduce the definition of chemical potential into each of equations (15.22), the fundamental thermodynamic expressions that apply to systems of variable composition become:

dU=TdSPdV+ i μ i d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.26a)
 
dH=TdS+VdP+ i μ i d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.26b)
 
dG=VdPSdT+ i μ i d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.26c)
 
dA=PdVSdT+ i μ i d n i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.26d)
 

It is clear that equations (14.23) still hold — evaluated at constant “n” (total number of moles) — as shown:

T= ( U S ) V,n P= ( U V ) S,n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(14.23a)
 
T= ( H S ) P,n V= ( H P ) S,n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(14.23b)
 
V= ( G P ) T,n S= ( G T ) P,n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(14.23c)
 
P= ( A V ) T,n S= ( A T ) V,n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(14.23d)
 

The same reasoning applies to Maxwell’s relationships. From the above expressions, equation (15.22b) is the only one that matches equation (15.13). Notice that the following additional identities can be also identified (see equations 14.20):

( μ i S ) V,n = ( T n i ) S,V, n i1 ( μ i V ) S,n = ( P n i ) S,V, n i1 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.27a)
 
( μ i S ) P,n = ( T n i ) S,P, n i1 ( μ i P ) S,n = ( V n i ) S,P, n i1 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.27b)
 
( μ i P ) T,n = ( V n i ) P,T, n i1 = V ¯ i ( μ i T ) P,n = ( S n i ) P,T, n i1 = S ¯ i This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.27c)
 
( μ i V ) T,n = ( P n i ) T,V, n i1 ( μ i T ) V,n = ( S n i ) T,V, n i1 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.27d)
 

We can see that chemical potential can be calculated by solving any of these differential expressions. To do this, experimental information is needed on how the other properties (T, V, S, P) change with additions of the given species (ni) under certain restraining conditions.

In some instances, we may have an open system of constant composition. This is the case of a system of a pure component exchanging mass with its surroundings. For such a system, nc = 1 and equations (15.26) are written as:

dU=TdSPdV+μdn This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.28a)
 
dH=TdS+VdP+μdn This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.28b)
 
dG=VdPSdT+μdn This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.28c)
 
dA=PdVSdT+μdn This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(15.28d)
 

Equations (14.23) and (15.27) will still hold with n i =n,  μ i =μ This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. .