1. Introduction
In classical thermodynamics, the internal energy, , of a one-component ideal gas is a function of three parameters: entropy, S, volume V and the number of particles N. Each parameter of state represents one independent way of system’s energy exchange with the external world. For fixed N, there are two ways of energy change: heat and work. The infinitesimal change of the internal energy satisfies the equation , where the first term is the heat and the last is the work term. Until recently, no such description was available for systems in non-equilibrium states, subjected to energy flow.
The classical theory of irreversible (non-equilibrium) thermodynamics [
1] is based on three differential non-linear equations representing conservation of mass, momentum (Navier–Stokes equation) and energy. These conservation laws are supplemented by the assumption of the local equilibrium, corresponding local equations of state and constitutive relations between fluxes and thermodynamic forces. The solutions of these equations are given in terms of velocity,
, temperature,
and number density of particles,
profiles. In the stationary states, the profiles and fluxes depend only on the position in space
, not on time,
t. The assumption of the local equilibrium is crucial in the formulation of the local equations of state. For the ideal gas it is justified, but contested for interacting systems [
2]. The assumption of local equilibrium for the ideal gas is valid for a small temperature gradient, i.e., such that
for the mean free path of the molecules,
[
3]. At the pressure of 1 bar at room temperature, the mean free path is of the order of
100 nm. In such conditions, the assumption of local equilibrium breaks down only for temperature gradients higher than
K/cm. Thus, this assumption is even satisfied inside the Sun! This local description contains the first law of thermodynamics, but in a different form than that given in the equilibrium thermodynamics and based on a few
global parameters of state. In our recent paper [
4] we provided the latter, i.e., global thermodynamic description for the ideal gas in a heat flow.
We have previously [
4] studied the one-component ideal gas in a heat flow between two parallel walls at distance
L, kept at two different temperatures
. The local equilibrium gives the local pressure,
, and the internal energy per unit volume,
:
and
with Boltzmann constant
, particle number density
, and the temperature
at position
z. We have shown rigorously that when we integrate both equations over the volume of the system we can formulate the global thermodynamics with the internal energy as a function of a few parameters of state. After integration we obtain:
and
We identify the system’s temperature,
,
In this way we made a mapping of the non-uniform system into the uniform one. We observe that obtained equations have the same form as in equilibrium for the ideal gas at temperature
. We demanded such form because, after the mapping, we treat the system as the uniform one as in equilibrium. Moreover, the obtained equations of state must reduce to the equilibrium equations of state when the heat flux is zero. Now we define the internal energy as a function of three parameters of state
, with the thermodynamic relation:
This mapping gives us the same formal structure as we know from equilibrium. The entropy
is responsible for the
net heat that enters or leaves the system [
5,
6,
7] and changes the internal energy. In general, the heat flows through the system all the time without changing the internal energy. Upon any process we would like to know how much of the heat transferred to the system changes the internal energy. This heat is called net heat and is given in the differential form by
.
is only part of the total entropy of the system. The total entropy of the system,
, where
is the volumetric entropy density given by local equilibrium assumption within irreversible thermodynamics, can be written as a sum [
4]:
Only
governs net heat in the system (heat absorbed/released in the system).
controls the dissipative background and solely depends on the entropy production given by [
1]
, where
is the heat conductivity. The difference between the total entropy and
vanishes,
, when the system approaches the equilibrium state,
. Therefore,
becomes in this limit the equilibrium entropy. Nonetheless the formal dependence of
on
U and
V at non-equilibrium state in a heat flow is the same as at equilibrium.
In this paper we want to apply the same mapping procedure to the binary mixture of ideal gases. The big difference between previous work and the current one is the fact that apart from the number density profile and the temperature profile we have additional profiles of the number densities of each component. As we shall see, these additional profiles lead to new parameters of state in the non-equilibrium state, which have no direct counterpart at equilibrium. The purpose of this work is to formulate the first law of global thermodynamics for ideal gas binary mixture in the heat flow. The paper is organized as follows: In
Section 2, we recall the equilibrium properties of the ideal gas mixture. In
Section 3 we discuss this mixture enclosed between two parallel walls at different temperatures and solve the equations of irreversible thermodynamics.
Section 4 and
Section 5 contain the main results of the present study. In
Section 4, we define all parameters of state for the mixture and perform the mapping of non-uniform system into the uniform one. We introduce the first law of non-equilibrium thermodynamics, which follows from these parameters, and discuss some of its consequences in
Section 5. We discuss the results in
Section 6.
2. Preliminaries
We consider a binary mixture of ideal gases enclosed between two parallel walls separated from each other by a distance
L in the
z direction. The volume of a system is
, where
A is an area in the
plane. The components of the mixture have number densities
and
, such that
is the total density. The equation of state of ideal gas at pressure
p and temperature
T,
, where
is the Boltzmann constant, can be written as a sum of partial pressures
which is the Dalton’s law [
8].
The internal energy density (per volume) is the sum of the internal energy density
(
) of the two components separately considered
. From the classical equipartition theorem applied to ideal gas it follows that
where
are translational and rotational degrees of freedom of the
i component.
takes the value 3 for monoatomic, 5 for diatomic, and 6 for polyatomic gas. Often, the ideal gas law (
9) is written using dimensionless specific heat capacity at constant volume
.
The equilibrium entropy density of a binary mixture is
Here,
, where
is the entropy per particle of the two components separately considered. The second term in (
10) is the mean-field expression for the mixing entropy. For ideal gas, one has
where
is independent of the thermodynamic state of the gas and has dimension of
. For a monoatomic ideal gas, a quantum mechanical theory of Sacur–Tetrode predicts that the constant
depends only upon the mass of the gas particle [
9,
10].
The total energy and entropy are given by
which reduces to
and
for homogeneous systems. In the absence of external fields, the equilibrium ideal gas is homogeneous.
From fundamental thermodynamic relation in terms of entropy, the differential
has the following form
In the next section we will use the chemical potential difference defined as
, where
. It is given by
where
is the number fraction of
i component;
in the absence of chemical reactions.
is the chemical potential of the component
i individually considered; it depends only on temperature
T and pressure
p.
3. Ideal Gas Mixture in Heat Flow
Now, we introduce heat flow into the system by setting different temperatures on the walls, i.e.,
In the non-equilibrium state induced by this boundary condition, the system becomes inhomogeneous and one has to consider spatially varying temperature
, pressure
, density
and difference in number densities
. In the hydrodynamic limit, the time evolution of the binary mixture is given by the conservation laws for
n,
, momentum and energy supplemented by the assumption of local equilibrium, relations between thermodynamic forces and fluxes, and thermodynamic equations of state [
1]. Note that ideal gas satisfies the local equilibrium exactly. This fact describes the strict absence of spatial correlations between particles, which is, of course, not the case for non-ideal systems.
We focus here on a stationary state with vanishing gas velocity field and constant pressure across the system. The latter, together with
and
, follows from conservation of mass and momentum. The assumptions
and
are in agreement with our previous simulation studies [
11,
12,
13] (for gas–liquid evaporating systems [
11,
12] and for Lennard–Jones fluid volumetrically heated [
13]), which show that mechanical equilibrium is established very fast (in comparison to heat flow). Thus, the balance equations for
and energy simplify to
where
is the diffusion current and
is the heat current. We choose the following phenomenological expressions for
and
[
1]:
In steady state,
because the system has no contact with particles reservoir.
First, we ignore the thermodiffusion effect, i.e., we assume that
. We also neglect the other cross-term by setting
. As a result, Equation (
17) reduces to
and
The latter yields a linear temperature profile (the system is translationally invariant in
directions)
From Equation (
18) we determine the relationship between the local number fraction
of the component 1 and the temperature. The local equilibrium allows us to use the local form of Equation (
14) for
; it depends on
, and the pressure
p, which is constant throughout the system. Introducing
, Equation (
18) with the use of Equation (
14) gives
With
we find
Note that if
, temperature inhomogeneity does not induce spatial variation of the composition, unless we take into account the effect of thermodiffusion, i.e.,
. In that case we find generally
Onsager coefficients
and
are assumed to be independent on the concentration and temperature [
14], therefore the solution of Equation (
23) is
This solution clearly shows that in the non-equilibrium steady state the difference
is a relevant parameter. An analogy can be drawn here with the order parameter: for non-zero
, the profile
is
qualitatively different from that for the zero difference
. Specifically, for
, the concentration of the first component of the mixture (with more degrees of freedom) increases towards the hotter wall, while for
the opposite is true (see
Figure 1). This observation is important in identifying the steady-state parameters in constructing global thermodynamics for this system. (In
Appendix B, we mention other approaches to thermodiffusion involving different, measurable thermodiffusion coefficients and relate them to the Onsager ones.)
The density profile follows from the local equation of state and the condition of constant pressure at steady state,
. For a given number of particles,
, this determines pressure as
Finally, the total energy of a mixture is the sum of the energies of the both components of the mixture (see Equations (
9) and (
12))
Thus, the set of independent parameters controlling the stationary state is
,
,
A,
L,
N,
,
,
.
4. Non-Equilibrium Parameters and Functions of State
In this section we start constructing non-equilibrium thermodynamics. To this end we rewrite the total energy (per volume) as follows
where
. This has the same form as in equilibrium but with temperature
T replaced by
and
replaced by parameters
, which we call the “effective” degrees of freedom. From Equation (
26), it follows that
and
where we have used
. The parameter
can be determined from the requirement that the pressure
Then, Equation (
30) can be written as
which is to say that partial pressures in steady state are related to the non-equilibrium energies in the same way as in equilibrium (see Equation (
9)) but with
replaced by
. Solving Equation (
31) with the use of Equation (
25), we find
can be interpreted as an “average” temperature in the non-equilibrium steady state. Eliminating
p in expressions for
we find:
Using
we can see that the new parameters of state
are not independent but obey the following relation:
We note that even if both components of the mixture have the same number of degrees of freedom
, the effective parameters
are not equal in non-equilibrium steady states providing we do not neglect thermal diffusion (see Equation (
24)). Using Equations (
20) and (
22) for the profiles of the number fraction
of the component 1 and temperature
, we can express the steady-state variable
in terms of the control parameters. The explicit formulas can be obtained for the case of
in terms of the special functions as presented in
Appendix A. For small, reduced temperature gradients
and for
we find
If
, the expression for
in this limit is more complicated, but has a similar structure in terms of dependence on
r. The coefficients are functions of
and
but not of
as shown in
Appendix A. In the case of
, for which
given by Equation (
24), the integral cannot be expressed in a closed form.
In the next step in the construction of global thermodynamics, we note that because the Equations (
27) and (
31) have the same structure as the equilibrium equation of state, we may formally write:
with
Equation (
36) has the functional form of the equilibrium fundamental relation for a binary ideal gas mixture with
replaced by
(compare Equations (
10) and (
11)). We will treat Equation (
36) as a definition of the
non-equilibrium steady-state entropy
that also provides the
fundamental relation for the non-equilibrium steady state. We note that
differs from total entropy defined as
, where
is the volumetric entropy density given by local equilibrium assumption within irreversible thermodynamics. In comparison with the equilibrium entropy,
depends on additional two state parameters
and
. As a consequence, for the difference
we find:
where
and
Equations (
32)and (
35) provide the effective parameters of state,
and
for given control parameters
and
at fixed
.
5. First Law and Its Consequences
We now consider the change of the total internal energy
in our system. From Equation (
38), we have
As argued in Ref. [
4], the change in internal energy
in the very slow process of moving from one steady state to another by varying the control parameters (so that the pressure remains homogeneous) is given by
where
is the net heat entering the system during this transition, and for the
fixed number of mixture components the mechanical work done is
The Equation (
42) is therefore the first law of thermodynamics for non-equilibrium steady states where the net heat identified by Equations (
42) and (
38) (at constant
and
) is
We note that the net heat flow during the transition between two steady states is a combination of two exact differentials: the effective entropy
and the effective degrees of freedom
(from Equation (
34) it follows that
). This is contrary to the equilibrium thermodynamics, in which heat depends only on temperature and the change in entropy.
At equilibrium, the most experimentally available thermodynamic quantities are response functions such as heat capacities, compressibility, susceptibility or chemical response functions. They follow from the first law of thermodynamics and the equilibrium fundamental relationship. Using the symmetry of the second derivatives, i.e., Maxwell relations, one can express one response function in terms of others. This is very useful as, for example, the heat capacity can be determined by measurements of other quantities, such as isothermal compressibility [
8].
Once we have established the fundamental relationship for steady states of binary mixtures of ideal gases, we can generalize the equilibrium response functions to steady-state response functions. First, we consider thermal response function, i.e., the heat capacity
C, which is a measure of amount of heat needed to raise the temperature of a system by a given amount. Generally, it is defined as derivative,
, but depending on which independent variables are fixed during the measurements, one has different heat capacities. For example, the heat capacity at constant volume and number of components
is defined as
, whereas
is the heat capacity at constant pressure and
. Meanwhile, in the equilibrium state at a constant volume and fixed
and
, we can change only the temperature of the system, in the non-equilibrium steady state we have more possibilities: we can independently vary
and
, or equivalently
and the reduced difference
. If one changes
at fixed
r then the state parameter
remains constant (see
Appendix A and Equation (
35)) we can define steady-state heat capacities for constant volume and pressure as follows
and
Concerning mechanical response functions, we can generalize isothermal compressibility and thermal expansivity to the steady states as follows
and
Because the fundamental equation of steady state has the same form as in equilibrium, the Maxwell relations for the path with
and
provide connection between the thermal and mechanical response function, which is the same as in the equilibrium, i.e.,
In order to determine
one need to measure the excess heat due to the small change of
, which is in principle possible due to the recent development of the appropriate experimental techniques [
15]. The coefficient
should be measured by changing the pressure at fixed both temperatures
and
, while the coefficient
should be determined by varying
at fixed reduced gradient
r.
In order to control the state parameter
during measurements of response function, one has to be able to determine it experimentally as a function of control parameters
, and
V for given
. To this end, one needs an experimental procedure to determine
and
(
is not an independent state parameter). Because the fundamental relation has the same form as in equilibrium thermodynamics, by performing the Legendre transform of the non-equilibrium fundamental relation, we can move to the variables
and
. First, we measure pressure and use Equation (
31) to determine
. The net heat, with the use of Equations (
34), (
38), (
40) and (
44) can be written as
. With the non-equilibrium temperature
determined as described above and the heat differential determined during the change from one to another steady state by slight change of
or
V, we determine
. This situation is similar to equilibrium thermodynamics, where mechanical and caloric measurements are necessary to determine the fundamental relation. This aspect is the same. The difference between steady-state and equilibrium thermodynamics lies in the fact that out-of-equilibrium, the measurement should be performed for a larger space of parameters. In our case, we have the two-dimensional space
in equilibrium while the three-dimensional space in steady-state
. Having the non-equilibrium parameters of state determined,
allows one to plan the experiment and change
in such a way that
is kept constant and the relation can be experimentally checked (
49).
Let us consider a more general case where the system’s temperature profile is unknown. This situation appears when the heat conductivity depends on the gas density. It is straightforward to apply the integration of equations over volume and obtain the fundamental relation without any changes. As a result, the fundamental relations hold with other (not necessarily linear) temperature profiles. Moreover, the effective parameters
can be determined experimentally according to the above-described procedure. In addition, because the pressure in the system is homogeneous, the mapping procedure also holds for the container of any shape. The introduced thermodynamic description thus works for a broad class of a binary mixture of ideal gases in vessels of different shapes and practically any temperature profile. There may be more control parameters, not only the two, that determine the temperature profile in the system. However, even when there are many control parameters, the fundamental relation with
, and
with only five parameters of state always carries the whole of thermodynamic information in the system and leads to thermodynamic relations such as relation (
49).
6. Summary
The internal energy of the binary ideal gas mixture in a heat flow is the function of five parameters of state
, irrespective of the number of boundary conditions. These parameters determine different ways of changing the internal energy of the system. The parameter
is responsible for additional net heat, not included in
.
is not the differential form describing total net heat in the system. The total net heat is given not only by the changes in
but also in
. The same observation was made for the van der Waals gas in the heat flow. In this case, the total net heat was given not only by
, but also by the renormalization in the mapping procedure of two parameters of state
describing attractive interactions and the excluded volume in the van der Waals gas [
16]. These parameters are constant at equilibrium since they are material parameters that define interactions in a particular system, but in non-equilibrium van der Waals gas they are parameters of state, which change the energy due to the change of density profiles. In our case of binary mixture the new parameter of state emerged from concentration profiles. Due to existence of profiles of different physical quantities, we expect all material parameters to become state parameters in the non-equilibrium systems.
Our construction of global thermodynamics for stationary states can be naturally extrapolated from binary to multi-component ideal gas mixtures in a stationary heat flow. Obviously, the number of state parameters of the stationary state will increase to include all
and
, where
and
m is the number of components. For the constant number of molecules
, the effective parameters
will not be independent. Each new parameter of state will add contribution to
as given by Equation (
41).
The formulation of the first law of global non-equilibrium thermodynamics opens the possibility to formulate the second law. The second law defines the direction of processes that take place in out-of-equilibrium systems. In our recent contributions [
17] we showed that formulating the first law of stationary thermodynamics for a specific class of problems is a necessary step towards pursuing the second law. This class of problems can be extended to additional forms of energies residing in the system, e.g., Van der Waals gas (internal potential), ideal gas in gravitational field (external potential) and in this contribution, the mixture of gases (multiple components). Formulation of the second law for these system is under development.