Next Article in Journal
Incremental Learning of Goal-Directed Actions in a Dynamic Environment by a Robot Using Active Inference
Next Article in Special Issue
An Entropy Generation Rate Model for Tropospheric Behavior That Includes Cloud Evolution
Previous Article in Journal
Graph Partitions in Chemistry
Previous Article in Special Issue
Complex Networks and Interacting Particle Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameters of State in the Global Thermodynamics of Binary Ideal Gas Mixtures in a Stationary Heat Flow

1
Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warszawa, Poland
2
Max-Planck-Institut für Intelligente Systeme Stuttgart, Heisenbergstr. 3, D-70569 Stuttgart, Germany
3
Department of Physics, Lancaster University, Lancaster LA1 4YB, UK
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2023, 25(11), 1505; https://doi.org/10.3390/e25111505
Submission received: 21 September 2023 / Revised: 26 October 2023 / Accepted: 30 October 2023 / Published: 31 October 2023
(This article belongs to the Special Issue Entropy, Time and Evolution II)

Abstract

:
In this paper, we formulate the first law of global thermodynamics for stationary states of the binary ideal gas mixture subjected to heat flow. We map the non-uniform system onto the uniform one and show that the internal energy U ( S * , V , N 1 , N 2 , f 1 * , f 2 * ) is the function of the following parameters of state: a non-equilibrium entropy S * , volume V, number of particles of the first component, N 1 , number of particles of the second component N 2 and the renormalized degrees of freedom. The parameters f 1 * , f 2 * , N 1 , N 2 satisfy the relation ( N 1 / ( N 1 + N 2 ) ) f 1 * / f 1 + ( N 2 / ( N 1 + N 2 ) ) f 2 * / f 2 = 1 ( f 1 and f 2 are the degrees of freedom for each component respectively). Thus, only 5 parameters of state describe the non-equilibrium state of the binary mixture in the heat flow. We calculate the non-equilibrium entropy S * and new thermodynamic parameters of state f 1 * , f 2 * explicitly. The latter are responsible for heat generation due to the concentration gradients. The theory reduces to equilibrium thermodynamics, when the heat flux goes to zero. As in equilibrium thermodynamics, the steady-state fundamental equation also leads to the thermodynamic Maxwell relations for measurable steady-state properties.

1. Introduction

In classical thermodynamics, the internal energy, U ( S , V , N ) , 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 d U = T d S p d V , 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, v ( r , t ) , temperature, T ( r , t ) and number density of particles, n ( r , t ) profiles. In the stationary states, the profiles and fluxes depend only on the position in space r , 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 l f p T / T 1 , for the mean free path of the molecules, l f p [3]. At the pressure of 1 bar at room temperature, the mean free path is of the order of l f p 100 nm. In such conditions, the assumption of local equilibrium breaks down only for temperature gradients higher than 10 7 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 T 1 ( z = 0 ) > T 2 ( z = L ) . The local equilibrium gives the local pressure, p ( z ) , and the internal energy per unit volume, u ( z ) :
p z = k B T z n z ,
and
u z = 3 2 n z k B T z ,
with Boltzmann constant k B , particle number density n z , and the temperature T z 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:
p V = N k B T 2 T 1 log T 2 T 1 .
and
U = 3 2 V N k B T 2 T 1 log T 2 T 1 .
We identify the system’s temperature, T * ,
T * = T 2 T 1 log T 2 T 1 .
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 T * . 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 U ( S * , V , N ) , with the thermodynamic relation:
S * U V , N = 1 T * , S * V U , N = p T * .
This mapping gives us the same formal structure as we know from equilibrium. The entropy S * 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 T * d S * .
S * is only part of the total entropy of the system. The total entropy of the system, S t o t = A d z s ( z ) , where s ( z ) is the volumetric entropy density given by local equilibrium assumption within irreversible thermodynamics, can be written as a sum [4]:
S tot U , V , T 2 T 1 = S * U , V + Δ S U , V , T 2 T 1 , Δ S U , V , T 2 / T 1 = N k B log T 2 T 1 5 / 4 log T 2 T 1 T 2 T 1 1 5 / 2 .
Only S * governs net heat in the system (heat absorbed/released in the system). Δ S controls the dissipative background and solely depends on the entropy production given by [1] σ = A 0 L d z κ T z · 1 T z = A κ L T 2 T 1 + T 1 T 2 2 , where κ is the heat conductivity. The difference between the total entropy and S * vanishes, Δ S U , A , L , T 2 / T 1 0 , when the system approaches the equilibrium state, T 2 / T 1 1 . Therefore, S * becomes in this limit the equilibrium entropy. Nonetheless the formal dependence of S * 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 V = A L , where A is an area in the x y plane. The components of the mixture have number densities n 1 and n 2 , such that n = N / V = n 1 + n 2 is the total density. The equation of state of ideal gas at pressure p and temperature T, p = n k B T , where k B is the Boltzmann constant, can be written as a sum of partial pressures p i
p = p 1 + p 2 = n 1 k B T + n 2 k B T
which is the Dalton’s law [8].
The internal energy density (per volume) is the sum of the internal energy density u i ( i = 1 , 2 ) of the two components separately considered u = u 1 + u 2 . From the classical equipartition theorem applied to ideal gas it follows that
u i = f i 2 n i k B T = f i 2 p i , i = 1 , 2
where f i , i = 1 , 2 are translational and rotational degrees of freedom of the i component. f i 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 c v ( i ) = f i 2 .
The equilibrium entropy density of a binary mixture is
n · s = n s ( 0 ) + k B n ln n n 1 ln n 1 n 2 ln n 2 .
Here, n s ( 0 ) = n 1 s 1 ( 0 ) + n 2 s 2 ( 0 ) , where s i ( 0 ) ( i = 1 , 2 ) 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
s i ( 0 ) k B = f i 2 + 1 + f i 2 ln 2 Φ i u i f i k B n i 1 + 2 / f i ,
where Φ i is independent of the thermodynamic state of the gas and has dimension of n 2 / f i / T . For a monoatomic ideal gas, a quantum mechanical theory of Sacur–Tetrode predicts that the constant Φ i depends only upon the mass of the gas particle [9,10].
The total energy and entropy are given by
U = V u d 3 r and S = V n s d 3 r ,
which reduces to U = V u and S = N s 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 d S has the following form
d S = 1 T d U + p T d V μ 1 T d N 1 μ 2 T d N 2 .
In the next section we will use the chemical potential difference defined as μ T = S φ u , n , where φ = n 1 n 2 . It is given by
μ T = μ 1 ( 0 ) μ 2 ( 0 ) 2 T + k B 2 ln x 1 1 x 1 ,
where x i = n i / n is the number fraction of i component; x 1 + x 2 = 1 in the absence of chemical reactions. μ i ( 0 ) T = n i s i ( 0 ) n i u i , n j i 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.,
T z = 0 = T 1 , T z = L = T 2 .
In the non-equilibrium state induced by this boundary condition, the system becomes inhomogeneous and one has to consider spatially varying temperature T r , pressure p r , density n r and difference in number densities φ r . 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 n r and φ r , follows from conservation of mass and momentum. The assumptions v = 0 and p r = p 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
φ t = 2 i J i d , e t = i J i Q ,
where J d J d , 1 = J d , 2 is the diffusion current and J q is the heat current. We choose the following phenomenological expressions for J d and J q [1]:
J d = L 11 μ T + L 12 1 T , J q = L 21 μ T + L 22 1 T
In steady state, J d = 0 because the system has no contact with particles reservoir.
First, we ignore the thermodiffusion effect, i.e., we assume that L 12 = 0 . We also neglect the other cross-term by setting L 21 = 0 . As a result, Equation (17) reduces to
μ T = 0
and
T = c o n s t 2 T = 0 .
The latter yields a linear temperature profile (the system is translationally invariant in x , y directions)
T ( z ) = T 1 + ( T 2 T 1 ) z L .
From Equation (18) we determine the relationship between the local number fraction x 1 ( z ) of the component 1 and the temperature. The local equilibrium allows us to use the local form of Equation (14) for μ / T ; it depends on T ( z ) , x 1 ( z ) , x 2 ( z ) , and the pressure p, which is constant throughout the system. Introducing T ˜ ( z ) = T ( z ) / T 1 , Equation (18) with the use of Equation (14) gives
ln x 1 ( z ) x 2 ( z ) = ln T ˜ ( z ) f 1 f 2 2 .
With x 1 ( z ) = 1 x 2 ( z ) we find
x 1 ( z ) = N 1 T ˜ ( z ) ( f 1 f 2 ) / 2 N 2 + N 1 T ˜ ( z ) ( f 1 f 2 ) / 2 .
Note that if f 1 = f 2 , temperature inhomogeneity does not induce spatial variation of the composition, unless we take into account the effect of thermodiffusion, i.e., L 12 0 . In that case we find generally
ln x 1 ( z ) x 2 ( z ) = f 1 f 2 2 2 L 12 k B L 11 T ( z ) ln T ˜ ( z ) .
Onsager coefficients L 12 and L 11 are assumed to be independent on the concentration and temperature [14], therefore the solution of Equation (23) is
x 1 ( z ) = N 1 T ˜ ( z ) ( f 1 f 2 ) / 2 N 2 e 2 L 12 k B T 1 L 11 1 1 T ˜ ( z ) + N 1 T ˜ ( z ) ( f 2 f 2 ) / 2 .
This solution clearly shows that in the non-equilibrium steady state the difference f 1 f 2 is a relevant parameter. An analogy can be drawn here with the order parameter: for non-zero f 1 f 2 , the profile x 1 ( z ) is qualitatively different from that for the zero difference f 1 f 2 . Specifically, for f 1 f 2 > 0 , the concentration of the first component of the mixture (with more degrees of freedom) increases towards the hotter wall, while for f 1 = f 2 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, n z = p / k B T z . For a given number of particles, N = A 0 L d z n z , this determines pressure as
p = N V k B T 2 T 1 log T 2 T 1 .
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))
U = U 1 + U 2 = A f 1 2 k B 0 L n ( z ) x 1 ( z ) T ( z ) d z + A f 2 2 k B 0 L n ( z ) x 2 ( z ) T ( z ) d z .
Thus, the set of independent parameters controlling the stationary state is T 1 , T 2 , A, L, N, N 1 , f 1 , f 2 .

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
u = U V = f 1 * 2 n 1 k B T * + f 2 * 2 n 2 k B T * ,
where n i = N i / V = ( 1 / L ) 0 L n i ( z ) d z . This has the same form as in equilibrium but with temperature T replaced by T * and f 1 , f 2 replaced by parameters f 1 * , f 2 * , which we call the “effective” degrees of freedom. From Equation (26), it follows that
f 1 * = f 1 1 L 0 L n ( z ) x 1 ( z ) k B T ( z ) 2 d z n 1 k B T * 2 = f 1 p 1 L 0 L x 1 ( z ) d z n 1 k B T *
and
f 2 * = f 2 1 L 0 L n ( z ) x 2 ( z ) k B T ( z ) 2 d z n 2 k B T * 2 = f 2 p 1 L 0 L x 2 ( z ) d z n 2 k B T *
where we have used n ( z ) T ( z ) = p / k B = c o n s t . The parameter T * can be determined from the requirement that the pressure
p = k B 1 L 0 L n 1 ( z ) T ( z ) + n 2 ( z ) T ( z ) d z = p 1 + p 2 .
Then, Equation (30) can be written as
p = n 1 k B T * + n 2 k B T * ,
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 f i replaced by f i * . Solving Equation (31) with the use of Equation (25), we find
T * = p k B 1 n 1 + n 2 = p k B V N = T 2 T 1 ln T 2 T 1
T * can be interpreted as an “average” temperature in the non-equilibrium steady state. Eliminating p in expressions for f 1 * , f 2 * we find:
f 1 * = f 1 x 1 1 L 0 L x 1 ( z ) d z , f 2 * = f 2 x 2 1 L 0 L x 2 ( z ) d z .
Using x 2 ( z ) = 1 x 1 ( z ) we can see that the new parameters of state f 1 * , f 2 * are not independent but obey the following relation:
f 1 * f 1 x 1 + f 2 * f 2 x 2 = 1 .
We note that even if both components of the mixture have the same number of degrees of freedom f 1 = f 2 , the effective parameters f 1 * , 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 x 1 ( z ) of the component 1 and temperature T ( z ) , we can express the steady-state variable f 1 * in terms of the control parameters. The explicit formulas can be obtained for the case of f 1 f 2 in terms of the special functions as presented in Appendix A. For small, reduced temperature gradients r = ( T 2 T 1 ) / T 1 and for N 1 = N 2 we find
f 1 * f 1 x 1 1 2 + f 1 f 2 16 r + O ( r 2 ) .
If N 1 N 2 , the expression for f 1 * in this limit is more complicated, but has a similar structure in terms of dependence on r. The coefficients are functions of N 1 / N 2 and f 1 , f 2 but not of T 1 as shown in Appendix A. In the case of f 1 = f 2 , for which x 1 ( z ) 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:
S * = ( S 1 * ) ( 0 ) + ( S 2 * ) ( 0 ) N 1 ln N 1 N N 2 ln N 2 N ,
with
( S i * ) ( 0 ) N i k B = f i * 2 + 1 + f i * 2 ln 2 Φ i U k B ( N 1 f 1 * + N 2 f 2 * ) V N i 2 / f i * .
Equation (36) has the functional form of the equilibrium fundamental relation for a binary ideal gas mixture with f i replaced by f i * (compare Equations (10) and (11)). We will treat Equation (36) as a definition of the non-equilibrium steady-state entropy S * ( U , V , N 1 , N 2 , f 1 * , f 2 * ) that also provides the fundamental relation for the non-equilibrium steady state. We note that S * differs from total entropy defined as S t o t = A d z s ( z ) , where s ( z ) is the volumetric entropy density given by local equilibrium assumption within irreversible thermodynamics. In comparison with the equilibrium entropy, S * depends on additional two state parameters f 1 * and f 2 * . As a consequence, for the difference d S * we find:
d S * = d U T * + p T * d V μ 1 * T * d N 1 μ 2 * T * d N 2 F 1 T * d f 1 * F 2 T * d f 2 * ,
where
μ i * k B T * = S * N i U , V , N j i , f 1 * , f 2 * = f i * 2 ln Φ i k B T * n i 2 / f 1 * + ln x i
and
F i T * = S * f i * U , V , N 1 , N 2 , f j i * = N i f i * ln x i μ i * T * .
Equations (32)and (35) provide the effective parameters of state, T * and f 1 * for given control parameters T 1 , T 2 and N 1 , N 2 at fixed N 1 + N 2 .

5. First Law and Its Consequences

We now consider the change of the total internal energy d U in our system. From Equation (38), we have
d U = T * d S * p d V + μ 1 * d N 1 + μ 2 * d N 2 + F 1 d f 1 * + F 2 d f 2 * .
As argued in Ref. [4], the change in internal energy d U 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
d U = đ Q + đ W ,
where đ Q is the net heat entering the system during this transition, and for the fixed number of mixture components the mechanical work done is
đ W = p d V .
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 N 1 and N 2 ) is
đ Q = d U + p d V = T * d S * + F 1 d f 1 * + F 2 d f 2 * .
We note that the net heat flow during the transition between two steady states is a combination of two exact differentials: the effective entropy d S * and the effective degrees of freedom d f 1 * (from Equation (34) it follows that d f 2 * = x 2 f 1 x 1 f 2 d f 1 * ). 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, C = đ Q / d T , 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 { N j } is defined as C V = đ Q / d T V , { N j } , whereas C p = đ Q / d T p , { N j } is the heat capacity at constant pressure and { N j } . Meanwhile, in the equilibrium state at a constant volume and fixed { N j } and { f j } , we can change only the temperature of the system, in the non-equilibrium steady state we have more possibilities: we can independently vary T 1 and T 2 , or equivalently T 2 and the reduced difference r = ( T 2 T 1 ) / T 1 . If one changes T 2 at fixed r then the state parameter f 1 * remains constant (see Appendix A and Equation (35)) we can define steady-state heat capacities for constant volume and pressure as follows
C V * = đ Q d T * V , { N j } , f 1 *
and
C p * = đ Q d T * p , { N j } , f 1 * .
Concerning mechanical response functions, we can generalize isothermal compressibility and thermal expansivity to the steady states as follows
κ T , N j * = 1 V V p T * , { N j } , f 1 *
and
α p , N j * = 1 V V T * p , { N j } , f 1 * .
Because the fundamental equation of steady state has the same form as in equilibrium, the Maxwell relations for the path with r = c o n s t and f 1 * = c o n s t provide connection between the thermal and mechanical response function, which is the same as in the equilibrium, i.e.,
κ T , N j * C p * C V * = T * V α p , N j * 2 .
In order to determine C p * one need to measure the excess heat due to the small change of T 2 , which is in principle possible due to the recent development of the appropriate experimental techniques [15]. The coefficient κ T , N j * should be measured by changing the pressure at fixed both temperatures T 1 and T 2 , while the coefficient α p , N j * should be determined by varying T 2 at fixed reduced gradient r.
In order to control the state parameter f 1 * during measurements of response function, one has to be able to determine it experimentally as a function of control parameters T 1 , T 2 , and V for given N 1 , N 2 . To this end, one needs an experimental procedure to determine S * and f 1 * ( f 2 * 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 T * and f 2 * . First, we measure pressure and use Equation (31) to determine T * ( T 1 , T 2 , V ) . The net heat, with the use of Equations (34), (38), (40) and (44) can be written as d Q = α ( T * , f 1 * , V ) d T * + β ( T * , f 1 * , V ) d f 1 * . With the non-equilibrium temperature T * determined as described above and the heat differential determined during the change from one to another steady state by slight change of T 1 , T 2 or V, we determine f 1 * ( T 1 , T 2 , V ) . 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 T 1 ( = T 2 ) , V in equilibrium while the three-dimensional space in steady-state T 1 , T 2 , V . Having the non-equilibrium parameters of state determined, T * ( T 1 , T 2 , V ) , f 1 * ( T 1 , T 2 , V ) , f 2 * ( T 1 , T 2 , V ) allows one to plan the experiment and change T 1 , T 2 , V in such a way that f 1 * 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 S * , f 1 * , f 2 * 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 S * , V , N 1 , N 2 , f 1 * , and f 2 * 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 U ( S * , V , N 1 , N 2 , f * ) , irrespective of the number of boundary conditions. These parameters determine different ways of changing the internal energy of the system. The parameter f * is responsible for additional net heat, not included in T * d S * . T * d S * is not the differential form describing total net heat in the system. The total net heat is given not only by the changes in S * but also in f * . 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 S * , but also by the renormalization in the mapping procedure of two parameters of state a * , b * 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 N i and f i * , where i = 1 , , m and m is the number of components. For the constant number of molecules N = i = 1 m N i , the effective parameters f i * will not be independent. Each new parameter of state will add contribution to d U 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.

Author Contributions

Conceptualization, R.H., A.M. and K.M.; methodology, R.H., K.M., A.M. and P.J.Ż.; formal analysis, R.H., A.M. and K.M.; investigation, R.H., A.M., K.M., K.G. and P.J.Ż.; writing—original draft preparation, A.M.; writing—review and editing, R.H., K.M., A.M., K.G. and P.J.Ż. All authors have read and agreed to the published version of the manuscript.

Funding

A.M. was supported by the Polish National Science Center (Opus Grant No. 2022/45/B/ST3/00936). P.J.Ż. would like to acknowledge the support of a project that has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie Grant Agreement No. 847413 and was a part of an international cofinanced project founded from the program of the Minister of Science and Higher Education entitled “PMW” in the years 2020–2024; Agreement No. 5005/H2020-MSCA-COFUND/2019/2.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Variables of Steady State in Terms of Control Parameters

The integral of the concentration profile x 1 ( z ) in Equation (28) for the integrand given by Equation (22) with N 1 = N 2 can be obtained analytically in terms of the Gaussian hypergeometric function 2 F 1 ( a , b ; c ; w ) and polygamma function ψ ( m ) ( w ) [18]:
1 L 0 L x 1 ( z ) d z = 1 r δ ψ ( 0 ) ( d ) ψ ( 0 ) ( d + 1 / 2 ) + 2 r ( 2 + δ ) 1 + r 1 + δ 2 2 F 1 ( 1 , 2 d ; 2 d + 1 ; ( 1 + r ) δ 2 ) ,
where d = 1 / δ + 1 / 2 with δ = f 1 f 2 > 0 and r = ( T 2 T 1 ) / T 1 > 0 is the reduced slope of the temperature profile. The state variable f 1 * is thus
f 1 * = f 1 x 1 B 1 ( δ , r ) + B 2 ( δ , r ) ,
where
B 1 ( δ , r ) = 1 r δ ψ ( 0 ) ( d ) ψ ( 0 ) ( d + 1 / 2 ) B 2 ( δ , r ) = 2 r ( 2 + δ ) 1 + r 1 + δ 2 2 F 1 ( 1 , 2 d ; 2 d + 1 ; ( 1 + r ) δ 2 ) .
For N 1 N 2 we obtain
f 1 * = f 1 x 1 D 1 δ , r , N 1 N 2 + D 2 δ , r , N 1 N 2 ,
where
D 1 δ , r , N 1 N 2 = 2 N 1 N 2 r ( 2 + δ ) 2 F 1 1 , 2 d ; 2 d + 1 ; N 1 N 2 D 2 δ , r , N 1 N 2 = 2 N 1 1 + r 1 + δ 2 N 2 r ( 2 + δ ) 2 F 1 1 , 2 d ; 2 d + 1 ; N 1 N 2 ( 1 + r ) δ 2 .
For small reduced temperature gradients across the system we find for N 1 N 2
f 1 * f 1 x 2 2 F 1 1 , 2 d ; 2 d + 1 ; x 1 x 2 f 1 x 1 x 1 x 2 2 δ B x 1 x 2 ; 2 d + 1 ; 1 + f 1 x 2 f 1 f 2 4 ( 1 + x 1 / x 2 ) 2 r + O ( r 2 ) ,
where B ( r ; p , q ) is incomplete Euler beta function [18]. For N 1 = N 2 the above expression reduces to Equation (35).

Appendix B. Other Approaches to Thermodiffusion

It is worth stressing that many numerical results as well as measurements of thermodiffusion coefficients are reported in the literature (see, e.g., [19,20,21,22] and references therein). They indicate that, unlike Onsager coefficients, these coefficients are functions of composition and temperature. For example, within the framework of the generalized Stefan–Maxwell thermodiffusion equations for f 1 = f 2 , one obtains exactly the same equation as Equation (23), but with 2 L 12 / ( L 11 k B T ) replaced by α 12 = A 12 / D 12 . A i j is the Newman–Soret thermal diffusion coefficient and D i j is the Stefan–Maxwell diffusivity (see e.g., Equation (6.8) from Ref. [22]). It is argued that coefficients A i j = D i T / ρ i D j T / ρ j , where D i T is the thermal diffusivity coefficient of component i and ρ i is its mass density, are more approximately constant with respect to composition in binary systems [23]. Assuming that the ratio α 12 does not depend on both temperature and composition, the solution of Equation (23) for the profile x 1 ( z ) is
x 1 ( z ) = N 1 T ˜ ( z ) ( f 1 f 2 ) / 2 α 12 N 2 + N 1 T ˜ ( z ) ( f 1 f 2 ) / 2 α 12 .
For noble gas pairs at T = 300 K, the value of α 12 is equal to 0.3432 for a Ne-He pair, 0.1741 for an Ar-Ne pair, or 0.0262 for a Xe-Kr pair [24,25].

References

  1. De Groot, S.R.; Mazur, P. Non-Equilibrium Thermodynamics; Dover Publications Inc.: New York, NY, USA, 2011. [Google Scholar]
  2. Vilar, J.M.G.; Rubi, J.M. Thermodynamics “beyond” local equilibrium. Proc. Natl. Acad. Sci. USA 2001, 98, 11081–11084. [Google Scholar] [CrossRef] [PubMed]
  3. Kreuzer, H.J. Nonequilibrium Thermodynamics and Its Statistical Foundations; Clarendon Press: Oxford, UK, 1984. [Google Scholar]
  4. Hołyst, R.; Makuch, K.; Maciołek, A.; Żuk, P.J. Thermodynamics of stationary states of the ideal gas in a heat flow. J. Chem. Phys. 2022, 157, 194108. [Google Scholar] [CrossRef] [PubMed]
  5. Oono, Y.; Paniconi, M. Steady state thermodynamics. Prog. Theor. Phys. Suppl. 1998, 130, 29–44. [Google Scholar] [CrossRef]
  6. Nakagawa, N.; Sasa, S.I. Global thermodynamics for heat conduction systems. J. Stat. Phys. 2019, 177, 825–888. [Google Scholar] [CrossRef]
  7. Chiba, Y.; Nakagawa, N. Numerical determination of entropy associated with excess heat in steady-state thermodynamics. Phys. Rev. E 2016, 94, 022115. [Google Scholar] [CrossRef] [PubMed]
  8. Callen, H.B. Thermodynamics and an Introduction to Thermostatistics; John Wiley & Sons: New York, NY, USA, 1960. [Google Scholar]
  9. Sackur, O. Die Anwendung der kinetischen Theorie der Gase auf chemische Probleme. Ann. Phys. 1911, 341, 958–980. [Google Scholar] [CrossRef]
  10. Tetrode, H. Die chemische Konstante der Gase und das elementare Wirkungsquantumt. Ann. Phys. 1912, 343, 434–442. [Google Scholar] [CrossRef]
  11. Babin, V.; Holyst, R. Evaporation of a Sub-Micrometer Droplet. J. Phys. Chem. B 2005, 109, 11367. [Google Scholar] [CrossRef] [PubMed]
  12. Holyst, R.; Litniewski, M. Heat Transfer at the Nanoscale: Evaporation of Nanodroplets. Phys. Rev. Lett. 2008, 100, 055701. [Google Scholar] [CrossRef] [PubMed]
  13. Holyst, R.; Maciołek, A.; Zhang, Y.; Litniewski, M.; Knychała, P.; Kasprzak, M.; Banaszak, M. Flux and storage of energy in nonequilibrium stationary states. Phys. Rev. E 2019, 99, 042118. [Google Scholar] [CrossRef] [PubMed]
  14. Onsager, L. Reciprocal Relations in Irreversible Processes. I. Phys. Rev. E 1931, 37, 405. [Google Scholar] [CrossRef]
  15. Yamamoto, T.; Nagae, Y.; Wakabayashi, T.; Kamiyama, T.; Suzuki, H. Calorimetry of phase transitions in liquid crystal 8CB under shear flow. Soft Matter 2023, 19, 1492–1498. [Google Scholar] [CrossRef] [PubMed]
  16. Hołyst, R.; Makuch, K.; Giżyński, K.; Maciołek, A.; Żuk, P.J. Fundamental Relation for Gas of Interacting Particles in a Heat Flow. Entropy 2023, 25, 1295. [Google Scholar] [CrossRef] [PubMed]
  17. Makuch, K.; Giżynski, K.; Hołyst, R.; Maciołek, A.; Żuk, P.J. Steady-state thermodynamics of a system with heat and mass flow coupling. J. Chem. Phys. 2023, in press. [Google Scholar]
  18. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Dover Publications: New York, NY, USA, 1972. [Google Scholar]
  19. Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. Transport Phenomena; Wiley: New York, NY, USA, 2001. [Google Scholar]
  20. Venerus, D.C.; Öttinger, H.C. A Modern Course in Transport Phenomena; Cambridge University Press: Cambridge, UK, 2018. [Google Scholar]
  21. Chapman, S. The characteristics of thermal diffusion. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1940, 177, 38–62. [Google Scholar]
  22. Van-Brunt, A.; Farrell, P.E.; Monroe, C.W. Consolidated theory of fluid thermodiffusion. AIChE J. 2022, 68, e17599. [Google Scholar] [CrossRef]
  23. Newman, J. Thermoelectric Effects in Electrochemical Systems. Ind. Eng. Chem. Res. 1995, 34, 3208–3216. [Google Scholar] [CrossRef]
  24. Marrero, T.R.; Mason, E.A. Gaseous diffusion coefficients. J. Phys. Chem. Ref. Data 1972, 1, 3–118. [Google Scholar] [CrossRef]
  25. Taylor, W.L. Thermal diffusion factors for the binary noble gas mixtures. J. Chem. Phys. 1980, 72, 4973–4981. [Google Scholar] [CrossRef]
Figure 1. Schematic illustration of a binary ideal gas mixture between two parallel walls held at different temperatures T 2 and T 1 . The resulting temperature gradient induces a gradient of fractions x i of the mixture components. Solid (green) lines indicate profiles of x 1 ( z / L ) . Two cases are shown: (i) f 1 > f 2 for which the profile is given by Equation (22) and (2) f 1 = f 2 for which the profile is given by Equation (A3). For the latter case, we took α 12 = 0.3432 , which corresponds to the Ne-He pair at T 1 = 300 K. In both cases the temperature profile is given by (20) with the reduced temperature gradient r = ( T 2 T 1 ) / T 1 chosen to be 0.5 and N 1 = N 2 .
Figure 1. Schematic illustration of a binary ideal gas mixture between two parallel walls held at different temperatures T 2 and T 1 . The resulting temperature gradient induces a gradient of fractions x i of the mixture components. Solid (green) lines indicate profiles of x 1 ( z / L ) . Two cases are shown: (i) f 1 > f 2 for which the profile is given by Equation (22) and (2) f 1 = f 2 for which the profile is given by Equation (A3). For the latter case, we took α 12 = 0.3432 , which corresponds to the Ne-He pair at T 1 = 300 K. In both cases the temperature profile is given by (20) with the reduced temperature gradient r = ( T 2 T 1 ) / T 1 chosen to be 0.5 and N 1 = N 2 .
Entropy 25 01505 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maciołek, A.; Hołyst, R.; Makuch, K.; Giżyński, K.; Żuk, P.J. Parameters of State in the Global Thermodynamics of Binary Ideal Gas Mixtures in a Stationary Heat Flow. Entropy 2023, 25, 1505. https://doi.org/10.3390/e25111505

AMA Style

Maciołek A, Hołyst R, Makuch K, Giżyński K, Żuk PJ. Parameters of State in the Global Thermodynamics of Binary Ideal Gas Mixtures in a Stationary Heat Flow. Entropy. 2023; 25(11):1505. https://doi.org/10.3390/e25111505

Chicago/Turabian Style

Maciołek, Anna, Robert Hołyst, Karol Makuch, Konrad Giżyński, and Paweł J. Żuk. 2023. "Parameters of State in the Global Thermodynamics of Binary Ideal Gas Mixtures in a Stationary Heat Flow" Entropy 25, no. 11: 1505. https://doi.org/10.3390/e25111505

APA Style

Maciołek, A., Hołyst, R., Makuch, K., Giżyński, K., & Żuk, P. J. (2023). Parameters of State in the Global Thermodynamics of Binary Ideal Gas Mixtures in a Stationary Heat Flow. Entropy, 25(11), 1505. https://doi.org/10.3390/e25111505

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop