2.1. Two Extremes: Differential Change of Total Internal Energy in a Closed System and Differential Change of Total Entropy in an Isolated System—De Donder’s Equation
Applying the chain rule to the differential change of the total internal energy and the differential change of the total entropy of the composite system (
Figure 1) [
7,
9]:
and using identity
, the following expression is obtained (
xj is an extensive quantity of the state of the system that changes during the thermodynamic transformation and can be the quantity of the state from any compartment of the composite system; in thermodynamic potentials, it is the primary system, i.e., a system for which the thermodynamic parameters that define a given thermodynamic potential apply – or a system in which a chemical reaction takes place):
Multiplied expression (3) by the infinitesimal change of the extensive parameter d
xj of the primary system due to the infinitesimal progression of a thermodynamic process (chemical reaction) towards the maximum entropy (equilibrium state) – example for an irreversible process in an isolated composite system (
UT=const.); or progression of a quasi-static reversible process towards minimal internal energy (final equilibrium state) – if the rate of change of extensive parameter is negligible small, succession of states of internal equilibrium [
4] – in a closed composite system (
ST=const.):
If the change in the system’s state starts from the
xj (
i ) initial value of an extensive variable, regardless of whether an isolated composite system with the constant total internal energy is observed (total entropy of the composite system increases to equilibrium) or an isentropically closed composite system (the total internal energy of a composite system decreases to the final equilibrium); the system is reached the same equilibrium state
xj (
e) (
Figure 1).
Equation (5) is actually a more general form of De Donder’s Equation (6) [
1,
16], which describes the production of entropy in an isolated system during a chemical reaction:
where
μi and
νi represent respectively the chemical potential and stoichiometric coefficient of the
i-th particle in the observed reaction
. In Equation (6)
A is the De Donder affinity of the reaction and d
ξ corresponds to the infinitesimal extent of the chemical reaction [
17]. As the chemical potential is the partial molar Gibbs energy, it can then be expressed with the partial molar enthalpy (
hi) and entropy (
si):
In Equation (7), rxn means that the thermodynamic functions change during the progress of the chemical reaction. In an isolated composite system,
represents the differential change of entropy in the pressure-temperature reservoir, while d
Srxn corresponds to the differential change of entropy in the primary system where the chemical reaction takes place, which means that De Donder’s differential change of entropy corresponds to the differential change of total entropy (
) of the isolated composite system. According to Calen, the equality
holds, from which Equation (8) follows. Equality
from De Donder’s equation can also be proven under
V = constant conditions when the system is without a reservoir of an intensive state parameter, i.e., the system is not composite. Namely, as the chemical potential is, in addition to the partial molar Gibbs energy, also the partial derivative of the internal energy by the amount of the
i-th component, it follows that the right-hand side of the equation in expression (7) is also the differential change of the total internal energy (
Figure 2):
2.2. Irreversibility: as a Loss of Working Capacity of the System in Which the Chemical Reaction Takes Place
Let the observed system be divided into two parts (
Figure 3). One part of the system is where the quasi-static reversible (isentropic) process always occurs (at each point of the compartment, the intensive quantities of the state have the same values). In contrast, the other part of the system is where both the quasi-static reversible and irreversible processes can occur (if an irreversible process also occurs, then the intensive state quantities have different values in the distinct places of the compartment). As internal energy and entropy are additives over the elements of an system, where the elements can be real compartments of the composite system or imaginary [
7], therefore, at some point in the thermodynamic process the total internal energy of the composite system can be represented as the sum of:
In Equation (10),
represents the fraction of internal energy of the composite (real or imaginary) system where a quasi-static reversible (isentropic process) process takes. Let it be assumed that the fractions
and
remain constant during the thermodynamic process (
+
= 1). Then the differential change of the total internal energy of the composite system due to the infinitesimal progress of a quasi-static reversible thermodynamic process (without an irreversible process in the
part of CS –
Figure 2 path
a) is:
Given Equality (5), Equation (11) is:
Generally, let’s consider a closed composite system (Cs) that becomes isolated. In this scenario, the total internal energy change (11) due to the thermodynamic process remains in the system as thermal energy increases the total entropy (12) - spontaneous, irreversible process.
Now, let’s return to the closed (imaginary) composite system shown in
Figure 3. For the part of the composite system where the quasi-static reversible process always occurs, the relation:
remains valid. Suppose in the part of the imaginary composite system where both the quasi-static reversible and the irreversible process can occur, the quasi-static reversible process takes place. In that case, it is (
Figure 2):
Presume that a spontaneous irreversible process occurs in
part of the system CS. In that case, part of the differential change in the total internal energy (
) of the hypothetical quasi-static reversible process returns to the system as differential heat (
refers to a quasi-static thermodynamic process with
, i.e., formally, a quasi-static reversible process first happened and then a spontaneous irreversible process,
Figure 3 path
b):
Since the system is further closed from the part where the quasi-static reversible process exclusively takes place (
), the differential change in internal energy (
) leaves CS without contributing to the differential entropy change of CS (13) (
Figure 3). While from part
of CS, where it now, in addition to a quasi-static reversible thermodynamic process (which first happened), also takes place spontaneous irreversible process – part of the total internal energy change
of a quasi-static reversible thermodynamic process – remains in the imaginary part of the system generating entropy (
). Therefore, Equation (12) in the case of a system transition from a quasi-static reversible regime to a system where, in addition to a quasi-static reversible process, there is also a spontaneous irreversible process (the system is closed,
Figure 3 path
b):
In Equation (17),
refers to the differential change of the total internal energy in CS if a quasi-statically reversible thermodynamic process takes place in each real or imaginary compartment (formally, which occurred before the irreversible process). The actual differential change in total internal energy of CS is
and leaves the composite system in the form of useful differential work (
Figure 2 path
b):
It follows from Equation (17):
and taking into account Equation (16):
In the last expression (20),
corresponds to the useful differential work, which would be obtained if the thermodynamic process in the composite system were entirely quasi-statically reversible
(
Figure 3 path
a);
represents the real useful differential work that is obtained if there is also a parallel irreversible process:
(
Figure 3 path
b). So the difference on the left side of Equation (20) is the differential part of the changes in the total internal energy that is wasted in the form of thermal energy during an irreversible process (wasted work potential – dissipated energy) and leads to the generation of entropy:
In Equation (21) are
and
since these are differential changes in the total internal energy of the composite system that leaves the system in the form of useful work. However, infinitesimal (differential) dissipated energy – differential irreversibility (
) is
[
18]. Namely, irreversibility is part of the change in the total internal energy of the composite system (i.e., the useful work of a quasi-static reversible process when the total entropy change would be zero), which remains in the given composite system – i.e., formally, returns to the composite system during an irreversible process as
(
Figure 3 path
b). If, at every moment of the thermodynamic process in the composite system, there is an internal equilibrium (hypothetical limit case), then a quasi-static reversible process takes place in the entire composite system. Hence, the irreversibility value is equal to zero:
Expression (21) for the final thermodynamic process is:
2.3. Change of Gibbs Free Energy and Irreversibility of a Thermodynamic Process
In quasi-static reversible processes, when the temperature and pressure remain constant (isothermal-isobaric conditions), the differential change in Gibbs free energy is equivalent to the maximum differential non-expansion work [
3,
4,
19,
20,
21]. Additionally, the differential change in the total internal energy of the composite system is also equal to the differential change in Gibbs free energy. This equality is represented by the Equation (25), which holds true when quasi-static reversible (isentropic) processes occur (
Figure 4) [
7]:
In Equation (25), represents the differential change of internal energy in the primary system (where the main thermodynamic process occurs, for example, a chemical reaction), while represents the differential change of internal energy in the reservoir that provides isothermal-isobaric conditions.
Calen introduced the concept of a reversible work (source) reservoir (RWR), which formally represents a hypothetical auxiliary system where the energy of the change in the total internal energy of the composite system (CS) is "stored", i.e., maximum non-expansion work (in the case of quasi-static reversible processes in the primary system) [
7]. Maximum non-expansion work from the CS in the RWR transforms into the potential energy (compression or stretching of the spring). RWR is a fictitious system with a conservative force. In the view of the observer (O) from the RWR, the energy of the change in the internal energy of the composite system, which is the maximum non-expansion work, leaves the CS, resulting in
and
. However, the same energy enters the RWR. Therefore, it is
(
Figure 4):
If a chemical reaction takes place in the primary system of the composite system CS with
T,
p reservoir, then the differential change in the internal energy of the primary system Sy is [
22]:
Differentiating Equation (27) by the chemical reaction extent at
T, p = const., the following equation is obtained:
Multiplying Equation (28) by the infinitesimal change in the chemical reaction extent gives:
The thermodynamic parameters of the primary system change due to the chemical reaction (
,
and
). The change in the internal energy of the
p,
T reservoir (R) of the composite system is possible due to volumetric work and thermal energy exchange, so Equation (25) is:
However, as the total volume of CS is constant, and the chemical reaction is reversible (quasi-static isentropic process), i.e., the CS’s total entropy change is zero, therefore Equation (31) at
T, p = const. has the following form:
which proves that the change in total internal energy of the composite system is equal to the change in Gibbs free energy resulting from the chemical reaction in the primary system (
Figure 4).
If, in a closed composite system, the total entropy change is not zero (
and
, i.e.,
), then a spontaneous irreversible process occurs during the thermodynamic process. Therefore part of the maximum non-expansion work turns into thermal energy (into a dissipated form of energy) – irreversibility. Considering Equation (21), the difference between the maximum non-expansion work of a quasi-static reversible process and the real useful non-expansion work whenever a spontaneous irreversible process occurs is (
Figure 5):
The difference (32) represents part of the change in the total internal energy of CS (due to the quasi-static reversible process) that is dissipated in the CS as heat (32) and increases the total entropy of the composite system – entropy is generated in the CS due to an irreversible (secondary) process. The absolute value of the differential non-expansion work
is smaller than the absolute value of the maximum differential non-expansion work
(Figure 5).
The criteria for quasi-static reversible (isentropic process) and irreversible spontaneous process in a closed system follow from the Equation (32) (
Table 1):
Clausius’s inequality [
23,
24] is hidden in the above inequality. Namely, the difference
represents part of the change in the total internal energy of the CS due to the primary thermodynamic process (i.e., chemical reaction), which is dissipated in the form of thermal energy in the CS, i.e., according to Clausius, it represents uncompensated entropy (i.e., entropy production due to an irreversible process). Multiplying the inequality by
yields Claussiu’s equation in its most generalized form:
If a quasi-static reversible process occurs in CS, then
so that the inequality (33) becomes equality (
Figure 5A and
Table 1):
Within a closed system, a spontaneous irreversible process generates entropy alongside a quasi-static reversible process (33). This entropy, in turn, resulting in less absolute value of negative work. Conversely, a hypothetical quasi-static reversible process within the same system (CS) would result in a higher absolute value of negative work. As such, the negative inequality value signifies that a spontaneous irreversible process has occurred within a closed composite system (36), thereby generating entropy (as depicted in
Figure 5B and
Table 1).
Suppose the value of
in Equation (33) equals zero. In that case, it signifies that the spontaneous irreversible process in the composite system utilizes the complete alteration in the internal energy of the quasi-static reversible process. As a result, the entire change in the Gibbs free energy of the quasi-static process transforms into heat energy – dissipates into the environment (i.e., thermoreservoir). In a closed CS, the maximum irreversibility is attained, which is the highest value of the produced entropy – the closed composite system behaves like an isolated CS, as depicted in
Figure 5C (37 and 38). It is crucial to note that no useful work can be extracted from the composite system (CS) in this scenario (
Table 1).
For CS in a natural equilibrium state (without a quasi-static reversible process), and hold, so inequality (33) is the equality 0 = 0.
Suppose CS is the real isolated system (constant total internal energy). In that case, the total enthalpy change of the spontaneous irreversible thermodynamic process remains in the system (
) so that inequality (33) is the equality (37 and 38) or:
Bearing in mind Equation (12) for an isolated system with a spontaneous irreversible process, it follows that the coefficient , i.e., for a specific change in the CS state the generated entropy have a maximum value (39).
The change in Gibbs free energy (ΔG and dG) represents the thermodynamic potential for the maximum non-expansion work of a quasi-static reversible process in a closed system (ΔST = 0 and dST = 0) or the change in a thermodynamic quantity -G/T (i.e., -ΔG/T or -dG/T) is the thermodynamic potential for the generated maximum total entropy of a spontaneous irreversible process in an isolated system (ΔUT = 0 and dUT = 0).
2.4. Black Box Method Application in the Volume Element of the Primary System of CS: the Gibbs Free Energy and the Irreversibility
Consider the primary system (Sy) - where the chemical reaction takes place intending to obtain useful work - and which (i.e., Sy) is located in the
T,
p reservoir forming together a composite system (
Figure 4). If only a quasi-static chemical reaction takes place in the primary system, then maximum work is obtained from the composite system (isentropic condition,
); while in the primary system, in addition to a quasi-static chemical reaction (isentropic conditions), a non-quasi-static reaction also takes place (or a non-quasi-static thermodynamic process is coupled with a quasi-static chemical reaction) then
(
) work is obtained from the composite system. Based on the above, the
parameter can be defined:
As the internal energy of the primary system (
) is an additive quantity over real or imaginary elements, Equations (10) and (11) apply, i.e.,
and
. The internal energy
in the equilibrium state is a function of the first order:
Therefore,
can be represented as the product of the parameter
and the element of the primary system’s internal energy (where the reversible process always takes place)
, for which the volume element
corresponds, where
holds. The volume element
corresponds to the volume for which the law of local equilibrium is valid [
1], i.e., that the internal energy and entropy densities are independent of the system’s coordinates and time. Similarly, for the
(part of the primary system where, in addition to the reversible thermodynamic process, an irreversible process can also take place), in the equilibrium state, it is valid that it is a homogeneous function of the first order.
However,
or
(which depends on the irreversible process), i.e., the volume element
, can consist of several volume elements
, between which an irreversible thermodynamic process takes place. The volumes
and
do not have to be compact, but their elements
and
are randomly distributed over the volume of the primary system (
Figure 6).
If the reaction
(isothermal-isobaric conditions) takes place in the volume element
– black box [
25] then after time
t from the beginning of the reaction as well as at time
t = 0, the internal energy of the observed volume element of the primary system is:
Where
corresponds to the amount of the
i-th reactive component at time
t,
corresponds to the concentration of the
i-th component in the volume element
at time
t, and
represents the partial molar internal energy of component
i at time
t. Let there be the same state (concentrations,
p,
T =
const.) in each element of the volume (
) of the primary system at time
t = 0, and if the reaction
proceeds infinitely slowly. So that the concentrations of the reactive components are equal in each
element, then each element of the volume can be viewed as a closed reactive system, i.e., there is no diffusion between volume elements. This means that in each volume element (
), the amount of the
i-th component changes solely due to the reaction:
is the amount of the
i-th component (reactant or product) at
t = 0, and
is the stoichiometric coefficient of component
i in the reaction
. The differential change of internal energy in the volume element
for the d
t time interval is:
As the volume element
corresponds to a closed system, it follows that the differential change in internal energy in d
t time corresponds to the differential amount of heat exchange (
) and the differential useful work (
). During the differential change of time (d
t) based on Equation (44), the change in the amount (concentrate) of the
i-th reactive component is:
Where
r represents the chemical reaction velocity (
[
26],
V corresponds to the volume where the reaction takes place). Inserting Equation (47) into Equation (45) yields:
The Gibbs-Duham equation also appears in equation (49) [
25,
27]:
Since the reaction
in the volume element
takes place at
p,
T =
const. conditions, the Gibbs-Duham equation is:
By introducing the Gibbs-Duham Equation (51) into Equation (49), we get:
If the reaction takes place in a liquid aggregate state, then there is
(
h = partial molar enthalpy) [
25], therefore Equation (52) is:
The entropy change in the volume element
due to the reaction
is (
s = partial molar entropy and the volume elements
of the primary system are isentropic):
By introducing Equality (55) into Equation (53), we get:
Since the volume contains the number of the element of the volume and since Equation (57) holds for each element of the volume, it follows that the useful (non-expansion) work is obtained from the volume .
Let a quasi-static isentropic process (chemical reaction) take place formally in the volume element
, and then an irreversible thermodynamic process, i.e., Equation (49) is valid with
:
In the quasi-static isentropic phase of the thermodynamic process, equality (55) applies, so Equation (58) is:
However, after an isentropic chemical reaction, an irreversible process formally follows (in reality, both thermodynamic processes occur in parallel), for example due to the friction of a piston, parts of the volume element
are heated, i.e., useful work is converted into thermal energy of friction (
), which is in agreement with Lucia’s statement according to which: "
at the beginning of the process, the system could be able to use all the energy, while during the process the irreversibility generates a loss of energy" [
28]. Therefore, the Gibbs-Duham Equation (50) at constant pressure is:
By introducing the Gibbs-Duham Equation (60) and
into Equation (59), we get:
When defining the volume element of the primary system
, it was assumed that with these volume elements, useful work is completely dissipated in an irreversible process, i.e., the equation
is valid for each volume element
– volume element
for the useful work formally becomes an isolated system:
For the part of the volume of the primary system (
) in which, in addition to the quasi-static reaction (with isentropic conditions in the composite system), an irreversible process also takes place,
applies. In the case of the dissipation of useful work into thermal energy, the equality between the volume elements can be taken (
) so that
(
Figure 6). In this way, the maximum useful non-expansion)work (
) that is obtained from the composite system (CS) when a quasi-static chemical reaction (with isentropic conditions) flows in each volume element
of the primary system is:
Where
corresponds to the differential change of Gibbs free energy per volume element
, while
is the total number of volume elements
. If in the
number of the volume element
, an irreversible process also flows, i.e., Equation (62) applies, then the differential useful non-expansion work from the composite system is:
In Equation (64)
represents irreversibility, i.e., dissipated useful work per volume element of the primary system (
). If the elements of the volume are not equal
, the equality
applies, i.e.,
(
Figure 6), then Equation (63) is:
that is, Equation (64) is:
Exergy is the maximum useful energy obtained from a system that reaches equilibrium with the environment during a hypothetical reversible (quasi-static isentropic) thermodynamic process. In the case of a chemical reaction (for example, in galvanic cells), the equilibrium of chemical potentials with the environment can be replaced by a zero value of De Donder’s affinity (
A = 0, Equation (6)). Therefore, Equations (63) and (65) correspond to the maximum non-expansion useful work (
), i.e., exergy that can be obtained from the composite system. Gouy (1889) and Stodola (1905) independently confirmed that the loss of exergy (maximum useful work) is proportional to the generated entropy [
29], which follows Equations (64) and (66).
Equation (62) can be generalized, namely since irreversibility is the dissipated useful work then it is:
, i.e.,:
Considering De Donder’s Equations (6) and (8), expression (67) is:
Where
represents the differential production of entropy in an irreversible process per volume element
(
, where
is the direct entropy change in the volume element
due to the chemical reaction and
represents the entropy change in the volume element
due to reaction heat). The differential entropy change in an irreversible thermodynamic process can be represented as the product of the thermodynamic force (
Fk) – proportional to the gradient of some intensive quantity – and an extensive physical quantity (corresponding to the intensive quantity) that is transported during the d
t time interval (
) [
1,
30]:
Alternatively, taking into account the flux of the extensive variable (
), i.e., entropy generation per unit of time (
) [
30]:
In this case, Equations (69) and (70) shed light on the volume element
of the primary system of the composite system, revealing its dynamic nature. For instance, when the volume elements of the primary system experience fluctuations in the chemical reaction velocity, the system’s behavior becomes more complex. In some elements of the volume, the reaction runs faster compared to the hypothetical infinitesimally slow chemical reaction, leading to a deviation from the state with uniform velocities for each element of the volume (
Figure 6). This dynamic behavior results in the concentrations of reactive components becoming non-uniform across the volume elements of the primary system. In addition to the thermodynamic force of the chemical reaction
A/
T (i.e.,
, Equation (70)), the dynamic system also gives rise to the thermodynamic force of diffusion of reactive components (diffusion between the volume elements of the primary system). (
-differential amount of the
k-th reactive component that diffuses in the d
t time interval and
is the difference of the
k-th component chemical potential between two points of the volume element
; for example, if
, i.e.,
then the difference
is between the two volume elements
):
The emergence of the thermodynamic force of diffusion can reduce the chemical potential of a reactive component in some volume element of the primary system. This, in turn, significantly influences the De Donder affinity (A) of the chemical reaction. The result is the absence of in some volume elements of the primary system, leading to the absence of useful non-expansion work. However, the most significant impact is the generation of entropy due to the irreversible diffusion process (Equation (71)), a crucial aspect of the system’s behavior.
2.5. Change of the Gibbs Free Energy whiout Non-Expansion Work
Suppose non-expansion work is not obtained from the closed CS, i.e., in that case, no chemical reaction takes place in the primary system Sy, but a quasi-static reversible expansion of the Sy’s volume, during which the Sy’s pressure also changes while
T = const. (
Figure 7). The change in Gibbs free energy for the isothermal quasi-static reversible process, i.e., an ideal gas’s expansion (no friction between cylinder and piston), is (
):
During the isothermal volume expansion process (pressure-volume work) of an ideal gas in the primary system of CS, the change in internal energy is zero (
). Hence, Equation (72) has the following form:
Let isothermal pressure-volume work first occur (for example, expansion of the primary system from CS), whereby a differential change in the volume of the primary system also means a differential change in the total volume of CS. The transfer of heat from the reservoir toward the primary system follow the expansion of the primary system – the primary system is located in the thermoreservoir, according to which the wall of the primary system is fixed, the wall of the primary system changes according to the RWR (
Figure 7). As d
SSy = -d
SR is valid, Equation (73) is the quasi-static reversible (maximum) differential pressure-volume work:
During expansion, the system performs work on RWR, so that in terms of CS, the maximum differential pressure-volume work has a negative value:
It follows from Equations (73) and (74) that the maximum differential pressure-volume work (expansion) comes from the change in the internal energy of the reservoir. A similar result is obtained starting from the Helmholtz free energy. If friction occurs during the movement of the piston – spontaneous dissipation of pressure-volume work into thermal energy, then a smaller absolute value of the pressure-volume work is obtained, i.e., Equation (32) applies:
Of course, even in this case, it is valid that at maximum irreversibility, the pressure-volume work of expansion of the primary system is zero, i.e., the change in Gibbs free energy is dissipated as thermal energy in CS (76). In the state of thermal and mechanical equilibrium of the primary system, i.e., without quasi-static pressure-volume work, dG = 0.
However, it should be emphasized that during the isothermal expansion of the primary system Sy (
Figure 7), at every moment of the infinitely slow thermodynamic process, the primary Sy is in mechanical equilibrium with the environment, i.e., in this case with RWR, so the quasi-static process is reversible with a total entropy change that is zero. However, suppose a chemical reaction occurs in the primary system where De Donder’s affinity is not zero (
). In that case, the chemical reaction (and, therefore, the primary system) is not in equilibrium – therefore, it is not a reversible process. According to Calen, if every infinitesimal step of a chemical reaction is described with an expression that refers to a quasi-static process:
=
(
p,
T =
const.) and, at the same time, the process is isentropic (i.e., the total entropy does not change in the composite system) [
7]. This hypothetical reaction in obtaining the maximum non-expansion work is identical to a reversible thermodynamic process.