Next Article in Journal
A van der Waals Model of Solvation Thermodynamics
Next Article in Special Issue
The Application of Pinch Technology to a Novel Closed-Loop Spray Drying System with a Condenser and Reheater
Previous Article in Journal
APCSMA: Adaptive Personalized Client-Selection and Model-Aggregation Algorithm for Federated Learning in Edge Computing Scenarios
Previous Article in Special Issue
Energy, Exergetic, and Thermoeconomic Analyses of Hydrogen-Fueled 1-kW Proton-Exchange Membrane Fuel Cell
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Direction of Spontaneous Processes in Non-Equilibrium Systems with Movable/Permeable Internal Walls

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
*
Authors to whom correspondence should be addressed.
Entropy 2024, 26(8), 713; https://doi.org/10.3390/e26080713
Submission received: 13 July 2024 / Revised: 2 August 2024 / Accepted: 20 August 2024 / Published: 22 August 2024
(This article belongs to the Special Issue Thermal Science and Engineering Applications)

Abstract

:
We consider three different systems in a heat flow: an ideal gas, a van der Waals gas, and a binary mixture of ideal gases. We divide each system internally into two subsystems by a movable wall. We show that the direction of the motion of the wall, after release, under constant boundary conditions, is determined by the same inequality as in equilibrium thermodynamics d U đ Q 0 . The only difference between the equilibrium and non-equilibrium laws is the dependence of the net heat change, đ Q , on the state parameters of the system. We show that the same inequality is valid when introducing the gravitational field in the case of both the ideal gas and the van der Waals gas in the heat flow. It remains true when we consider a thick wall permeable to gas particles and derive Archimedes’ principle in the heat flow. Finally, we consider the Couette (shear) flow of the ideal gas. In this system, the direction of the motion of the internal wall follows from the inequality d E đ Q đ W s 0 , where d E is the infinitesimal change in total energy (internal plus kinetic) and đ W s is the infinitesimal work exchanged with the environment due to the shear force imposed on the flowing gas. Ultimately, we synthesize all these cases within a general framework of the second law of non-equilibrium thermodynamics.

1. Introduction

The direction of spontaneous processes follows the second law of equilibrium thermodynamics [1]. It states that isolated systems reaching equilibrium increase their entropy. The second law is practical. It allows us to predict which chemical reactions occur spontaneously, under what thermodynamic conditions liquids evaporate or freeze, and how to design efficient engines. However, such a law, although needed, has yet to be discovered for out-of-equilibrium systems characterized by the continuous flux of energy flowing across them. This remains true despite large efforts [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28], and sucesses are limited to either isothermal situations or small temperature differences [21,22,23,24,25,26,27,28]. Promising results have been obtained in small-scale systems allowing for a discrete description [29,30], and they pave the way for a trajectory-based approach to macroscopic stochastic thermodynamics [31].
Non-equilibrium thermodynamics today (called the thermodynamics of irreversible processes) is a set of non-linear differential equations. Depending on the specific system under consideration, these equations include the conservation of mass, momentum, and energy [32,33,34], as well as charge conservation with Maxwell’s equations in the case of magnetohydrodynamics or the Poisson–Nernst–Planck equations for electrokinetics, among others. The conservation laws are typically accompanied by closing relations like the ideal gas equation of state. This paper will present preliminary observations for a few systems without macroscopic motion and for a single system with macroscopic motion. These results may provide insights toward formulating a general second law of non-equilibrium thermodynamics.
Non-equilibrium states involve macroscopic energy fluxes flowing across the system [33]. These fluxes are sustained within the system by non-vanishing gradients of temperature (in the case of heat flow), pressure (in the case of mass flow), or chemical potential (in the case of particle diffusion). Thus, a non-equilibrium state has non-uniform temperature, pressure, chemical potential, and velocity. The local equations of state, relating internal energy and pressure to density and temperature, and the spatial profiles of density, temperature, and concentrations fully characterize non-equilibrium states. In our recent papers, we formulated the first law of global thermodynamics for various non-equilibrium systems [35,36,37], including cases involving gravity [36] or Couette flow [38]. We represented internal energy as a function of a few global state parameters. These parameters were obtained by mapping a non-equilibrium, and by definition, non-uniform system into a uniform one. We averaged local equations of state over the system’s volume. This averaging resulted in global equations of state, which we wrote in the same form as at equilibrium. These global non-equilibrium equations of state included new state parameters. For example, the internal energy of the van der Waals gas subjected to continuous heat flux U ( S , V , N , a , b ) is a function of five state parameters, where S is the non-equilibrium entropy, V is the volume, N is the number of particles, and a and b are new state parameters—renormalized van der Waals interaction parameters. The net heat that flows in/out of the system to change the internal energy of the van der Waals gas is given by:
đ Q = T d S N 2 V d a + N k B T V N b 1 d b .
The a and b state parameters appear in the net heat differential because, in a non-uniform system, the change in the density profile leads to the local absorption or release of heat. In general, under the mapping construction of global thermodynamics, all material parameters in the equilibrium equations of states become state parameters of the non-equilibrium state.
In a system kept at a constant temperature at equilibrium, energy is exchanged with the environment as heat only. The second law of equilibrium thermodynamics states [1,39] that the Helmholtz free energy, F ( T , V , N , x ) , is minimized with respect to x (the variable describing the internal constraint) at constant T , V , and N (temperature, volume, and number of particles, respectively). The minimization defines the equilibrium value of x. The change in F as a function of x is negative if the initial x does not correspond to the equilibrium state. Thus, the change in F when we move the system from the initial to the final x is given by Δ F 0 . Rewriting this equation using internal energy and entropy yields Δ U T Δ S 0 . In the infinitesimal form, we obtain d U T d S 0 . Finally, the net heat (the heat that enters or leaves the system and changes the internal energy) is đ Q = T d S . In general, the second law states that in spontaneous processes,
d U đ Q 0 .
Thus, the system minimizes part of its internal energy. This part does not account for the energy, which, at equilibrium, is continuously exchanged with the environment (here, the net heat).
In this paper, we elucidate non-equilibrium systems in a continuous heat flow. Out of equilibrium, net heat becomes the amount of energy that enters or leaves the system in the form of heat and changes the internal energy. We show the consequences of applying d U đ Q 0 in the systems’ stationary (steady) states. Section 2 discusses the ideal gas, van der Walls gas, and binary mixture of ideal gases enclosed by two fixed walls at different temperatures. We cite the functional forms of d U and đ Q as functions of the state parameters for a single compartment [35,36,37] and point out that đ Q differs from what we know from equilibrium. Next, as is done at equilibrium, we divide the enclosed gas with an internal diathermic wall and verify that (2) determines the spontaneous motion of the wall toward the stationary position. Section 3 extends the problem to account for gravity and heat flow. We present the energy differential and the fundamental relation for a single compartment in the case of the ideal gas [36] and expand them to the case of the van der Waals gas. Then, we introduce the internal diathermic wall between two compartments and verify that (2) holds during the spontaneous motion of the wall. Finally, we assume that the wall has a finite thickness and is permeable (porous). This leads to the derivation of Archimedes’ principle from the second law of thermodynamics. Section 4 demonstrates that the relative stability of multiple stable states can be compared within the proposed framework. We cite a system where the ideal gas is heated volumetrically [16] and calculate the integral of external work necessary to switch between the states. In Section 5, using the example of the ideal gas in shear flow [38], we introduce a further generalization of the proposed second law to the case where, in addition to heat, the system also exchanges energy in the form of work. The work originates from the shear forces imposed on the gas flow. We conclude this paper by generalizing the second law in the discussion section.

2. The Ideal Gas, Gas Mixtures, and van der Waals Gas in Heat Flow

Ideal gas. The model geometry is defined by two parallel walls at z 1 = 0 and z 2 = L , as shown in Figure 1. For now, let us focus on the arrangement that does not include the inner wall.
The area of each wall, A, is large enough so that the system is translationally invariant in the x and y directions. The walls are in contact with thermostats, which are maintained at different temperatures, T 1 < T 2 . In the presence of a heat flux, quantities specifying equilibrium thermodynamic states, such as number density n = N / V , pressure p, or temperature T, become space-dependent. The profiles n ( z ) , p ( z ) , and T ( z ) of these quantities can be determined using, for example, the irreversible thermo-hydrodynamics approach, which relies on the local equilibrium assumption and represents the conservation of mass, momentum, and energy, supplemented with the equations of states and relations between fluxes and thermodynamic forces. For a stationary state with a vanishing velocity field, the thermo-hydrodynamics problem reduces to the constant pressure condition p z = p .
The construction of global thermodynamics describing the non-equilibrium steady state is achieved by mapping it to a uniform system characterized by a finite number of state parameters. Mapping involves averaging the local pressure and energy over the system’s volume so that, when averaged, the global equations of state have the same form as in equilibrium. For an ideal gas with f degrees of freedom, this procedure gives, for a steady state with no mass flow,
p = 1 L 0 L n z k B T z d z = n k B T , U = f 2 A 0 L n z k B T z d z = f 2 V n k B T ,
where n = N / V = ( 1 / L ) 0 L n ( z ) d z is the average number density and T is given by
T = A 0 L d z n z T z A 0 L d z n z .
The new state parameter T can be interpreted as the average temperature of the uniform system onto which we have mapped the original, non-equilibrium one. The explicit expression for T depends on the specific temperature and number density profile. In the case of constant pressure [35], the equations of irreversible thermo-hydrodynamics yield T ( z ) = T 1 + ( T 2 T 1 ) ( z / L ) , n ( z ) = n T 2 T 1 / log ( T 2 / T 1 ) T ( z ) , and as a result,
T = ( T 2 T 1 ) / ln ( T 2 / T 1 ) .
Because mapping gives us the same formal structure as we know from equilibrium, the variable S conjugate to T is given by
S U , V , N / k B = N log U ( f / 2 ) k B N f / 2 V N + N s 0 / k B ,
where the constant s 0 is chosen such that S for T 2 = T 1 gives the equilibrium expression for entropy [1]. The internal energy of a non-equilibrium steady state is thus a function of three state parameters U S , V , N , with the thermodynamic relations
S U V , N = 1 T , S V U , N = p T .
Binary mixture. Mixing two ideal gases, a and b, with f a and f b degrees of freedom, respectively, adds to the temperature and density profile the concentration profile induced by a heat flux. As a result, the mapping procedure leads to a uniform system with two new state parameters. In addition to T , we have two effective degrees of freedom f a and f b , which, however, are not independent. These parameters allow us to write the internal energy and pressure in the same form as in equilibrium:
p = n a k B T + n b k B T , U = f a 2 V n a k B T + f b 2 V n b k B T ,
where n a = N a / V and n b = N b / V are the average number densities of the first and second components of the mixture, respectively. The effective degrees of freedom are obtained as
f i = f i x i 1 L 0 L x i ( z ) d z , i = a , b ,
where x i = n i / n is the average number fraction of the i component and x i ( z ) is its profile; x a + x b = 1 in the absence of chemical reactions. f a and f b are related via f a / f a x a + f b / f b x b = 1 . The non-equilibrium entropy S has the same form as in equilibrium but with T replaced with T and f i with f i . It is the sum of the entropy of the two components of the mixture considered separately and the entropy of mixing:
S U , V , N a , N b , f a / k B = N a f a 2 + 1 + ln V N a 2 U k B ( N a f a + N b f b ) f a / 2 + N b f b 2 + 1 + ln V N b 2 U k B ( N a f a + N b f b ) f b / 2 N a ln N a N N b ln N b N + S 0 / k B .
Here, S 0 may depend on N a and N b . Again, we choose S 0 such that S reduces to the equilibrium entropy of a binary mixture in the absence of heat flux. Thus, the internal energy of a binary mixture of ideal gases in the non-equilibrium steady state is a function of five state parameters: U S , V , N a , N b , f a .
Van der Waals gas. A gas of interacting particles obeying the van der Waals equations of state is described with two additional interaction parameters, a and b. As a result of the mapping procedure, we obtain
p = n k B T 1 n b a n 2 , U = f 2 V n k B T a V n 2 .
Three new state parameters describe the van der Waals gas in the steady state. Besides the effective temperature T , defined by the same expression (4) as for the ideal gas, we have the effective interaction parameter a given by
a = 1 L 0 L a n z 2 n 2 d z ,
and b defined by the formula
n k B T 1 n b = 1 L 0 L n z k B T ( z ) 1 b n ( z ) d z .
Because Equation (11) has the same structure as in equilibrium, the non-equilibrium entropy S has the same form as in equilibrium but with T replaced with T , a with a , and b with b . Specifically,
S U , V , N , a , b / k B = N log U + a N 2 V ( f / 2 ) k B N f / 2 V N b N + N s 0 / k B .
Thus, the internal energy of the van der Waals gas in the non-equilibrium steady state is a function of five state parameters: U S , V , N , a , b .
Note that for all systems discussed above, the non-equilibrium entropy S is part of the total entropy S t o t , defined as the integral of the volume entropy density s ( z ) over the system volume. It contains information about the heat absorbed/released in the system in addition to the dissipative background (temperature profile).
Net heat. In the case of a very slow transition between stationary states, e.g., by a slight change in temperature T 2 or a change in the distance between the confining walls L, the energy changes only through mechanical work and heat flow:
d U = đ Q + đ W , with đ W = p d V ,
where đ Q is the net heat transferred to the system during a small change between two non-equilibrium stationary states. Thus, the above equation can be considered as the first law of non-equilibrium thermodynamics. Using the fundamental relations for each of the above examples (Equations (6), (10) and (14)), we can determine the net heat from the energy balance d U = đ Q p d V . For an ideal gas, we have
d U = T d S p d V , and đ Q = T d S ,
which has the same formal structure as in equilibrium.
This is different for a binary mixture of ideal gases, where the net heat acquires additional terms. In general,
d U = T d S p d V + μ a d N a + μ b d N b + F a d f a + F b d f b ,
and for a fixed number of components N a and N b , we have
đ Q = T d S + F a d f a + F b d f b ,
where d f b = x b f a x a f b d f a and F i T = S f i U , V , N a , N b , f j i , with S given by Equation (10).
Additional terms in the net heat also occur for the van der Waals gas, where
d U = T d S p d V N 2 V d a + n k B T V N b 1 d b ,
so that
đ Q = T d S N 2 V d a + n k B T V N b 1 d b .
Movable internal wall. We introduce into the considered systems an internal constraint in the form of a wall parallel to the bounding walls (Figure 1). We assume the wall is thin, freely movable, diathermal, and impenetrable. The internal constraint divides the system into two subsystems, 1 and 2, each with a fixed number of particles, N 1 and N 2 , and volumes, V 1 and V 2 . We choose the volume of one subsystem, V 1 , as the parameter representing the constraint. The total energy is the sum of the energies of the two subsystems:
U ( S 1 , V 1 , N 1 , S 2 , V V 1 , N 2 , ) = U 1 ( S 1 , V 1 , N 1 , ) + U 2 ( S 2 , V V 1 , N 2 , ) ,
where … denotes possible other new state variables, such as f a for the binary mixture of ideal gases or a and b for the van der Waals gas. The steady position of the wall, which moves without friction, is determined by the equality of the pressures exerted by each subsystem: p 1 = p 2 .
In the case of a monoatomic ideal gas, the effective temperatures of both subsystems are [35]
T 1 = V 1 V T 2 T 1 log T 1 + V 1 V T 2 T 1 T 1 , T 2 = 1 V 1 V T 2 T 1 log T 2 T 1 + V 1 V T 2 T 1 .
Using the explicit forms of d U 1 , d U 2 and đ Q 1 , đ Q 2 , we can see that this condition can be obtained from the following minimum principle:
d U 1 + d U 2 đ Q 1 đ Q 1 = p 1 p 2 d V 1 0 ,
which we call the second law of non-equilibrium thermodynamics. Thus, the difference between the total energy of the system and the heat exchanged with the environment is minimized while the internal wall moves to the new position. The equality defines the condition of the stationary position of the wall, given by the equality of pressures p 1 = p 2 in the two subsystems. For the specific case of an ideal gas, we have
p 1 = k B N 1 V T 2 T 1 ln T 1 + ( T 2 T 1 ) V 1 V T 1 , p 2 = k B N 2 V T 2 T 1 ln T 2 T 1 + ( T 2 T 1 ) V 1 V .
Demanding p 1 = p 2 and solving for V 1 gives us the position of the internal wall at the stationary state for given temperatures T 2 and T 1 :
V 1 V = T 2 T 1 N 2 / N 1 T 2 T 1 1 .
We note that for the diathermal internal wall, there is only a single solution for V 1 for all possible values of T 2 > T 1 . In Section 4, we consider an adiabatic internal wall with volumetric heating and show that, in this case, there might be more solutions for V 1 , among which we can compare stability.

3. Gravitational Field: Ideal Gas, van der Waals Gas, and Archimedes’ Principle

Ideal gas in gravity. We insert the gas column into a constant, external gravitational field e ^ z g that acts in the direction of the heat flux (see Figure 2). Here, the situation of interest is without macroscopic fluid motion, which means there is no Rayleigh–Benard convection [40,41,42]. The meanings of symbols from previous sections remain valid. The non-equilibrium entropy of a single gas column of an ideal gas is [36]
S = 3 2 N k B ln U N M g L 2 N M g z 1 U 0 N M 0 g L 0 2 N M g z 1 , 0 + N k B ln V V 0 + S 0 ,
where L = z 2 z 1 is the height of the gas column; z 1 is the coordinate of the column base; z 2 is the coordinate of the column top; V = L A is the volume of the column, with A being the surface area of the column base; and M denotes the molecular mass of the gas. The lower index 0 indicates reference values. In the above, the state parameter M is the renormalized mass defined with the help of the potential gravitational energy of the gas [36]
E pot = N M g L 2 M = 2 E pot N g L ,
which is contained inside the column
E pot = g A L 2 z 1 / L z 2 / L ρ z d z .
The state parameter M is coupled to the gravitational field and informs about the system’s potential energy with respect to the position of the column base. To account for all effects of gravity, we also need to include the potential energy of the base of the column
E pot , 0 = N M g z 1 .
The total internal energy of the gas column
U = 3 2 N k B T + N M g L 2 + N M g z 1 ,
consists of thermal and gravitational contributions. The thermal contribution (through the mapping procedure) defines the renormalized temperature T
E T = 3 2 N k B T = A L z 1 / L z 2 / L 3 k B 2 M ρ T d z T = V N M z 1 / L z 2 / L ρ T d z .
The average pressure inside the column is the same as in the case without gravity (3)
p = N k B T V .
This is because changing the surface area of the column base does not affect the gravitational energy, and all work is done against the thermal motion of gas particles. However, the presence of an external field introduces anisotropy in the pressure profile. The pressures at the top
p ( z 2 ) = N k B T V N M g L 2 V ,
and bottom
p ( z 1 ) = p ( z 2 ) + N M g L V
of the column are different. The presented formulas are consistent with the functional form of S (26) and the internal energy differential
d U = T d S + N g L 2 d M N k B T V N M g L 2 V V L d L N k B T V d V + N M g d z 1 = T d S + N g L 2 d M + đ W = đ Q + đ W .
In the above, we distinguish two separate ways of changing the internal energy. The first is through heat,
đ Q = T d S + N g L 2 d M ,
which includes entropic components and gravitational interactions between the external field and the mass inside the column. The second is through mechanical interactions exerted on the system boundaries through mechanical work
đ W = N k B T V N M g L 2 V V L d L N k B T V d V + N M g d z 1 .
We deliberately wrote the work differential to distinguish between the column elongation, changes in the column volume, and changes in the column base position. The column elongation is coupled to the pressure at its top, while volume changes affect the average pressure.
In the two-compartment system (Figure 2a), the bottom part contains N 1 particles, and the top part contains N 2 particles of an ideal gas. We separate the compartments with a thin, diathermal wall at z = z w . We assume there are no changes in the column cross-sectional area ( d A = 0 ) and set this contribution to the energy differential to 0. We rewrite the energy balances for the bottom ( U 1 ) and top ( U 2 ) segments explicitly to account for the displacement of the separating wall
d U 1 = T 1 d S 1 + N 1 g ( z w z 1 ) 2 d M 1 N 1 k B T 1 V 1 N 1 M 1 g ( z w z 1 ) 2 V 1 V 1 ( z w z 1 ) d z w = T 1 d S 1 + N 1 g ( z w z 1 ) 2 d M 1 p 1 ( z w ) A d z w = đ Q 1 + đ W 1 ,
d U 2 = T 2 d S 2 N 2 g ( z 2 z w ) 2 d M 2 + N 2 k B T 2 V 2 N 2 M 2 g ( z 2 z w ) 2 v 2 V 2 ( z 2 z w ) d z w + N 2 M g d z w = T 2 d S 2 + N 2 g ( z 2 z w ) 2 d M 2 + p 2 ( z 2 ) A d z w + N 2 M g d z w = đ Q 2 + đ W 2 .
The sign in front of p 2 is consistent with the direction of changes in d z w . Following the second law of thermodynamics, we check whether the change in internal energies, when subtracted from all energies exchanged with the environment during the process, leads to the proper spontaneous state of the system, even beyond equilibrium. For the column subject to gravity with a fixed base and top, this is represented by the heat balance
d U 1 + d U 2 đ Q 1 đ Q 2 = A p 2 ( z 2 ) p 1 ( z w ) + N 2 M g d z w 0 .
We find that the wall will position itself where the force necessary to expand the bottom part of the column balances the sum of the forces required to compress the top part of the column and, additionally, to raise its weight.
Van der Waals gas in gravity. We combine the results from the previous sections and construct an analogous description of the van der Waals gas in a gravitational field. For the single compartment, the stationary entropy is
S = 3 2 N k B ln U N M g L 2 N M g z 1 + a N 2 V U 0 N 0 M 0 g L 0 2 N 0 M g z 1 , 0 + a 0 N 0 2 V 0 + N k B ln V N b V 0 N 0 b 0 + S 0 .
The internal energy of the van der Waals gas in a gravitational field is given by
U = 3 2 N k B T + N M g L 2 + N M g z 1 a N 2 V .
The pressures are as follows:
p av = N k B T V N b a N 2 V 2
is the average pressure,
p ( z 2 ) = N k B T V N b a N 2 V 2 N M g L 2 V
is the pressure at the top of the column, and
p ( z 1 ) = p ( z 2 ) + N M g L V
is the pressure at the bottom of the column. The energy differential has the form
d U = T d S + N g L 2 d M N k B T V N b a N 2 V 2 N M g L 2 V V L d L N k B T A V N b a N 2 V 2 d V N 2 V d a + N k B T V N b 1 d b + N M g d z 1 = T d S + N g L 2 d M N 2 V d a + N k B T V N b 1 d b + đ W = đ Q + đ W ,
where we distinguish the heat differential, including the gravitational effects,
đ Q = T d S + N g L 2 d M N 2 V d a + N k B T V N b 1 d b ,
and the mechanical work
đ W = N k B T V N b a N 2 V 2 N M g L 2 V V L d L N k B T A V N b a N 2 V 2 d V + N M g d z 1 .
Likewise, for the ideal gas, we write the second law of stationary thermodynamics as
d U 1 + d U 2 đ Q 1 đ Q 2 = A p 2 ( z 2 ) p 1 ( z w ) + N 2 M g d z w 0 .
Again, it reads that the wall will rest where the force necessary to extend the bottom compartment matches the force necessary to compress the top compartment plus the force to lift its weight.
Archimedes’ principle. Novel phenomena appear when the wall separating compartments is permeable. Assume that the wall has a finite thickness L w with its own mass M b and is pierced with small channels (pores) of negligible volume that allow for the passage of gas molecules (Figure 2b). We parametrize the position of the bottom of the wall at z w . The porous wall separating the two compartments constitutes a separate thermodynamic system. We make additional assumptions regarding the wall: the gas that fills the pores inside the wall has a volume V 3 V 2 , V 1 , but each pore is large enough so that the gas is in the thermal (not Knudsen) regime, and the wall is an excellent thermal conductor (temperature is constant inside the wall).
Each compartment containing a perfect gas is described by the fundamental relation (35). Since we allow for particle exchange, the energy balance (37) has to account for the resulting changes in energy U N d N = μ d N . We do not elucidate this term and treat it formally because it will cancel out during further derivation. The amended energy differential is
d U = T d S + N g L 2 d M N k B T V N M g L 2 V V L d L + N M g d z 1 + μ d N = T d S + N g L 2 d M + đ W + μ d N = đ Q + đ W + μ d N .
We rewrite the energy differential explicitly for both compartments. Compartment 1 is at the bottom, so its z 1 stays fixed while z 1 = z w + L w of the upper compartment 2 is subject to change:
d U 1 = T 1 d S 1 + N 1 g L 1 2 d M 1 N 1 k B T 1 V 1 N 1 M 1 g L 1 2 V 1 V 1 z w z 1 d z w + μ 1 d N ,
d U 2 = T 2 d S 2 + N 2 g L 2 2 d M 2 + N 2 k B T 2 V 2 N 2 M 2 g L 2 2 V 2 V 1 z 2 z w L w d z w + N 2 M g d z w μ 2 d N .
Here, we set d L = d z w and d N = d N 1 = d N 2 , which results from the direction of motion of the platform and the flux of the gas.
This case is different from the previous situations because, apart from the motion of the wall, we also allow for the exchange of particles, which introduces an additional degree of freedom. As previously guided by the second law of equilibrium thermodynamics, we sought to subtract from the total energy contributions exchanged with the environment. We perform the following heuristic consideration to elucidate the structure of the emerging second law of steady-state thermodynamics.
To write the total energy balance during the process of wall motion, we need to calculate the energetic cost of a passage of a particle batch d N from compartment 1 to compartment 2 through the porous wall. One way to achieve this is to use the metaphor of an external agent, similar to Maxwell’s demon. We calculate the energy the demon has to spend to displace d N particles from the bottom to the top compartment. The passage starts when d N particles leave compartment 1, which generates an energy gain in the demon’s account:
đ E b 1 w = μ 1 d N .
Inserting d N particles into the pore end at temperature T ( z w ) requires, in addition to the volumetric work, the thermal energy
đ E w 1 p = 3 2 k B T ( z w ) d N ,
which the demon has to subtract from the account. Transition through the region border is accompanied by the volumetric work performed on both sides of the border. When a batch of d N leaves compartment 1, the remaining gas has to fill the space. Inside the pore, space for the incoming batch has to be accommodated. Both processes—expansion and compression—happen under the same temperature and pressure, which are continuous over the border. Therefore, having the same magnitude with opposite signs, they cancel out in the demon’s account. Next, the batch is pushed through the pore, and during that process, it changes the potential energy,
đ E w , pot = M g L w d N ,
and undergoes decompression, performing work
đ W w = z w z w + L w p ( z ) V N ( z ) z d z d N .
Both quantities need to be supplied by the demon and thus subtracted from the account. Once the end of the pore is reached, thermal energy
đ E w 2 p = 3 2 k B T ( z w + L w ) d N
is released and added to the demon’s account. Simultaneously, the particle needs to equilibrate to the bulk conditions by taking energy
đ E b 2 w = μ 2 d N ,
which has to be paid by the demon. During passage through the region border, again, space has to be emptied and filled on both sides of the pore end at T ( z w + L w ) , which does not influence the demon’s account. The final state of the account and, thus, the energetic cost of the passage results from summing all the listed effects:
đ E pass = μ 2 μ 1 + M g L w d N + 3 2 k B T ( z w + L w ) T ( z w ) d N + W w d N .
To complete the energy balance, we also need to account for the potential energy change due to the motion of the mass of the wall M w in the gravitational field:
đ E w , pot = M w g d z w .
With the use of đ E pass and đ E w , pot , we can write the second law of stationary thermodynamics:
d U 1 + d U 2 + đ E w , pot đ Q 1 đ Q 2 + đ E pass = A p 2 ( z w + L w ) p 1 ( z w ) + M w g d z w + 3 2 k B T ( z w + L w ) T ( z w ) M g L w + W w d N 0 .
The first three terms represent the changes in the system. The following two terms represent the fluxes of energies that can be exchanged with the environment, and the last term accounts for the additional energetic costs to pass through the obstacle. We stress that all terms involving the chemical potential μ are canceled out.
To calculate W w , we use the assumption that the wall is an excellent heat conductor so that T ( z w + L w ) = T ( z w ) . This means that the demon had to support the cost of the isothermal expansion
W w = k B T ( z w ) ln p 2 ( z w + L w ) p 1 ( z w ) ,
which we can explicitly insert into the second law of stationary thermodynamics:
d U 1 + d U 2 + đ E w , pot đ Q 1 đ Q 2 + đ E pass = A p 2 ( z w + L w ) p 1 ( z w ) + M w g d z w k B T ( z w ) ln p 2 ( z w + L w ) p 1 ( z w ) + M g L w k B T ( z w ) d N 0 .
From the first bracket, one can see that the wall will move until the pressure at the top of the bottom compartment (1) matches the pressure at the bottom of the top compartment (2) plus the weight of the wall divided by its surface area:
p 1 ( z w ) p 2 ( z w + L w ) = M w g A .
This is the same as saying that the wall will move until the force necessary to expand the bottom compartment matches the force necessary to compress the top compartment plus the weight of the top compartment ( p 2 ( z w + L w ) = p 2 ( z 2 ) + N 2 M g A ; see (36d)) plus the weight of the wall. The second bracket informs about the conditions for the gas flow through the pores. The stationary state is reached when
p 2 ( z w + L w ) = p 1 ( z w ) e M g L w k B T ( z w ) .
This means that the pressures on both sides of the wall satisfy the jump given by the hydrostatic pressure drop inside the column of gas in the wall. When the expressions in the first and second brackets vanish, we find the condition for the position of the cylinder inside the whole column:
p 1 ( z w ) = M w g A 1 e M g L w k B T ( z w ) .
The above equation, for the stationary state, binds the mass of the wall, its thickness, temperature, and height above the ground through the pressure at p 1 ( z w ) . We see that the pressure needed to squeeze an appropriate amount of gas into the pore depends on the wall’s temperature. A warm wall requires more gas.
Finally, we expand the fraction in (55) assuming that M g L w k B T ( z w ) 1 and substitute p 1 ( z w ) = ρ 1 ( z w ) k B T ( z w ) M to find
p 1 ( z w ) = ρ 1 ( z w ) k B T ( z w ) M = M w g A 1 e M g L w k B T ( z w ) = M w g A k B T ( z w ) M g L w +
and Archimedes’ principle emerges:
ρ 1 ( z w ) L w A = M w .
It states that a body submerged in water will rest exactly where the weight of the displaced fluid is equal to that of the submerged body. Here, equivalently, punching another pore through the wall will not make it move.

4. Volumetric Heating of an Ideal Gas Separated by an Adiabatic Wall

Volumetric heating. We consider a system similar to the one in Section 2 but with volumetric heating, [16] and start by describing a single compartment. The additional energy supplied uniformly throughout the volume is λ per unit of time and unit of volume. It extends the set of parameters controlling the non-equilibrium steady state of this system to ( T 1 , T 2 , λ , V , N ) . In the steady state, the pressure p and, consequently, the energy density ϵ = ( f / 2 ) p are constant, similar to the systems discussed in Section 2. For volumetric heating, the outflow of the heat is balanced by the absorption of heat throughout the system. Therefore, the temperature profile is obtained from the following local energy continuity equation:
κ 2 T ( r ) = λ ,
which assumes Fourier’s law of heat conduction. The coefficient κ is the thermal conductivity, which we assume to be temperature-independent. We apply the mapping of a non-equilibrium steady-state system to a homogeneous equilibrium system. Solving Equation (58) for the temperature profile with the boundary conditions T ( z = 0 ) = T 1 and T ( z = L ) = T 2 , we find
T ( z ) = λ 2 κ z 2 + T 2 T 1 L + λ L 2 κ z + T 1 .
After the mapping procedure given by Equation (3), we obtain the effective temperature
T = T 1 F T 2 T 1 , λ L 2 κ T 1 ,
where F ( u , w ) is a dimensionless function given by
F ( u , w ) = ( u 1 + w / 2 ) 2 + 2 w 2 tanh 1 u 1 + w / 2 ( u 1 + w / 2 ) 2 + 2 w tanh 1 u 1 w / 2 ( u 1 + w / 2 ) 2 + 2 w .
In the limit λ 0 , the effective temperature T reduces to the value given by Equation (5) for an ideal gas in the heat flux induced by the temperature gradient T 2 T 1 (This can be seen by using the definition of inverse hyperbolic tangent in terms of logarithms). On the other hand, if the temperatures on both outer walls are equal (i.e., T 2 = T 1 ), the system reduces to that described in our previous article [16]. The entropy S ( U , V , N ) is given by Equation (14) with U ( S , V , N ) = ( f / 2 ) N k B T = ( f / 2 ) p V . For a fixed number of molecules N, the energy change is determined by the control parameters ( T 1 , T 2 , λ , V ) through Equation (60). For an incremental and slow change between steady states that does not disturb the pressure uniformity in the system, we have [35]
d U = đ Q p d V and đ Q = T d S .
Movable wall. We introduce a constraint into the system in the form of a movable adiabatic wall parallel to the bounding walls located at z = z w (see Figure 3). The adiabatic wall does not conduct heat, yet it allows for momentum transport [1]. In the microscopic picture, both are tied to molecular motion, and the mechanism of the position equilibration of the wall depends on the statistics of the thermal collisions on both sides of the wall. In the case of infinite compartments with equal pressure p but different temperatures, the wall surprisingly moves toward the hotter region [43]. The molecular collision mechanism is crucial in the case of an array of compartments separated by moving adiabatic walls that fill the space between two fixed walls kept at different temperatures. As a result of its existence, we obtain Fourier’s law [44]. Here, although it moves, we treat the adiabatic wall according to the macroscopic picture. The stable position is set by pressure equality on both sides, and the temperature profile has an additional boundary condition, i.e., T ( z ) / z = 0 , at the internal wall. Such a wall does not allow for energy transport between compartments as a result of molecular collisions.
The temperature profiles T ( 1 ) ( z ) and T ( 2 ) ( z ) for subsystems (1) and (2) with N 1 and N 2 molecules, respectively, are
T ( 1 ) ( z ) = λ 2 κ z z w 2 + λ 2 κ z w 2 + T 1 , T ( 2 ) ( z ) = λ 2 κ z z w 2 + λ 2 κ ( L z w ) 2 + T 2 .
The mapping procedure yields
T 1 = T 1 G λ z w 2 κ T 1 , T 2 = T 2 G λ ( L z w ) 2 κ T 2 ,
where
G ( w ) = w ( w + 2 ) 2 tanh 1 w w + 2 .
Thus, the energy of the subsystems of N 1 and N 2 particles under the constraint is
U = U 1 + U 2 = f 2 N 1 k B T 1 G λ z w 2 κ T 1 + f 2 N 2 k B T 2 G λ ( L z w ) 2 κ T 2 ,
with N 1 + N 2 = N .
Our second law of non-equilibrium thermodynamics states that
d U 1 + d U 2 đ Q 1 đ Q 1 = p 1 p 2 d V 1 0 ,
and, therefore, the condition for the whole system to reach a steady state is that the pressure exerted by each subsystem is equal. This condition is equivalent to
1 z w G λ z ˜ w 2 = r N 2 N 1 1 1 z ˜ w G λ r ( 1 z ˜ w ) 2 ,
where z ˜ w = z w / L , r is the temperature ratio r = T 2 / T 1 , and λ = λ / ( κ T 1 ) . At equilibrium, i.e., for r = 1 and λ = 0 , the internal wall is located precisely in the middle of the system at z w = 0 . For λ = 0 but r > 1 , the position of the internal wall is given by Equation (25).
An analysis of Equation (68) shows that the number of solutions for z w varies from one to three depending on the values of λ and r. A typical course of the variability of the function p 1 ( z w ) p 2 ( z w ) is shown in Figure 4, compared to the case of no temperature gradient ( r = 1 ). For r = 1 and small values of λ , the curve p 1 ( z w ) p 2 ( z w ) is a monotonically decreasing function with a single zero-crossing point at the midpoint of the system. A local stability analysis shows that this is a stable position of the internal wall. Upon increasing the volumetric heating, the function develops a minimum and a maximum but remains symmetrical with respect to z w = 1 / 2 with three zero-crossing points at z w ( 1 ) = 0 , z w ( 2 ) = z w , and z w ( 3 ) = L z w , as illustrated in Figure 4a. Among these three solutions, further analysis shows that locally, both z w ( 2 ) and z w ( 3 ) are stable, whereas z w ( 1 ) is unstable. Moreover, the work required to change the position of the wall from z w to L z w
W z w ( 2 ) > z w ( 3 ) = z w ( 2 ) z w ( 3 ) p 1 ( z ) p 2 ( z ) d z
is equal to zero ( W z w ( 2 ) > z w ( 3 ) = 0 ). This, in turn, means that these two stable steady states coexist.
In the case of a non-zero temperature gradient, symmetry is broken, and the situation changes qualitatively. The single solution that exists for smaller values of λ shifts away from the midpoint toward the lower temperature wall, as might be expected. The difference p 1 ( z ) p 2 ( z ) becomes non-monotonic as λ increases, similar to the behavior observed in the absence of a temperature gradient. Initially, the extremes of the function develop below the y = 0 axis. Above the critical value λ c ( r ) , at which the maximum of the function p 1 ( z ) p 2 ( z ) touches the y = 0 axis, three zero crossing points z w ( 1 ) < z w ( 2 ) < z w ( 3 ) appear, as shown in Figure 4b. These correspond to three non-equilibrium steady states. To evaluate their stability, we assume that the internal wall is displaced from a specific steady state in both directions. The pressure difference will push the inner wall toward z w ( 1 ) or z w ( 3 ) but push it away from z w ( 2 ) . Thus, the middle position is locally unstable, whereas both positions close to the external walls are locally stable. Considering the work needed to move the wall from position z w ( 1 ) to position z w ( 3 ) allows us to determine which steady state is globally stable. If the work W z w ( 1 ) > z w ( 3 ) done on the system during this process is negative, the final position z w ( 3 ) corresponds to a globally stable steady state. If this work is positive, the initial position z w ( 1 ) is globally stable. We calculate this work numerically up to machine accuracy. We find that for r > 1 , due to the asymmetry, the work is positive for all values of λ . Therefore, z w ( 1 ) is always the globally stable steady-state position of the internal wall. Thus, the internal wall is globally stable when moved closer to the colder external wall.

5. Couette Flow of an Ideal Gas

The last example is the Couette flow of an ideal gas [38]. The gas flows between two walls, as shown in Figure 5. The wall at z 2 = L moves at a constant speed v 2 in the direction of the imposed shear force. The wall at z 1 = 0 is fixed. The velocity profile is linear: v ( z ) = v 2 z / L . Inside the system there is an adiabatic massless wall at location z w , which moves at the local speed of the fluid v w = v 2 z w / L in the same direction as the upper wall. The wall divides the system into two parts: 1 and 2. The location of the wall is an internal constraint in the system. We ask the following question: if we allow the wall to move along the z direction, what will be the final stationary location of this internal wall? This system differs from the previous ones in two ways. Firstly, the total energy includes both internal energy and kinetic energy. The latter was absent in the cases discussed in the previous paragraphs. Secondly, the system exchanges energy with the external world in the form of heat and work. Continuous input of work maintains a constant flow and constant kinetic energy despite dissipation due to shear. In the steady state, with the internal wall fixed in place, the work dissipates as heat within the system. When we release the internal wall, it starts moving in the z direction. During its motion, the external machine must perform extra work to overcome extra shear (the constant velocity of the upper wall) and dissipation.
The internal wall divides the system into two subsystems: 1 and 2. According to the equations of irreversible thermodynamics [38], the change in internal energy d U i for subsystem i ( i = 1 , 2 ) is given by
d U i = đ Q i + đ W i + đ D i ,
where đ Q i is the excess heat in subsystem i ( i = 1 , 2 ), đ D i is the excess dissipation in subsystem i, and đ W i = p i d V i is the volumetric work performed during the transition in subsystem i. In this non-equilibrium case, we have T d S = đ Q + đ D . Thus, the entropy change is due to both excess heat exchange and dissipation due to shear within the system. d V i is the change in volume in subsystem i. From the equation of irreversible thermodynamics [38], it follows that the change in kinetic energy, E k , i in subsystem i ( i = 1 , 2 ) has the following form:
d E k , i = đ W w , i đ D i .
The excess work performed by the surface forces, đ W w , i , during the motion of the internal wall in the z direction maintains a constant velocity profile in the system despite changes in the shear forces at the walls. Again, motivated by the second law of equilibrium thermodynamics, we generalize the second law as follows: From the sum of the total energy of the two subsystems, we subtract the heat and work exchanged with the environment during the processes and require that the difference be non-positive:
d U 1 + d U 2 + d E k , 1 + d E k , 2 đ Q 1 đ Q 2 đ W w , 1 đ W w , 2 0 .
From this equation, we obtain the expected form of the second law:
( p 1 p 2 ) d V 0 .
The equality in the above expresion defines the condition for the stationary position of the wall, given by the equality of pressures p 1 = p 2 in the two subsystems.

6. Discussion

The system at equilibrium, which exchanges heat with the environment, satisfies the inequality d U đ Q 0 at a constant temperature. In this contribution, we showed that the same inequality holds for non-equilibrium states with constant boundary conditions, including boundary temperature. This inequality is the second law of non-equilibrium thermodynamics for systems coupled to the environment via heat flux.
In equilibrium thermodynamics, at constant temperature T and pressure p, the second law states that d U T d S p d V 0 . We can write it in a more general form as d U đ Q đ W 0 . In the example of the Couette flow, a similar inequality sets the direction of spontaneous processes in the system. In general, in a system that exchanges energy via heat đ Q and different forms of work đ W j , the following inequality should govern spontaneous processes:
d E đ Q j đ W j 0 .
E is the system’s total energy, including internal, potential, and kinetic energy. The most challenging part of our study was identifying the various terms in this equation and calculating the net heat and work performed during the process.
We are now in a position to discuss the second law in more detail, including all its subtleties. The second law of equilibrium thermodynamics states that the entropy of an isolated system reaches its maximum value at equilibrium. However, S ( U , V , N ) has a well-defined value for fixed U , V , N . Therefore, we need to compare the entropy of the system in its equilibrium state to the entropy of the system in states that are not reached simply by changing the state parameters U , V , N . We introduce an extra parameter x, usually in the form of a constraint in the system (like a movable wall), and state that S ( U , V , N , x ) is maximized as a function of x for fixed state parameters U , V , N (internal energy, volume, and number of particles, respectively). Similarly, when the system has a constant temperature, the Helmholtz free energy F ( T , V , N , x ) is minimized as a function of x (the variable describing the internal constraint) at constant T , V , N (temperature, volume, and number of particles, respectively). We treat this extra parameter as arising from an external device, i.e., an external device that performs work on the system by coupling to this internal parameter x and moving the system in a reversible manner (i.e., via a sequence of equilibrium states) between states that are not accessible by simply changing T , V , N . Let us define the work performed on the system by this external device as đ W z . In the process of transitioning between states, the first law of thermodynamics must be obeyed (conservation of energy). Thus, we have d U = đ Q + đ W z . the second law states that if we move a system from a stable state to a less stable state, the external device will perform work on our system, i.e., đ W z > 0 . On the other hand, if we move from a less stable state to a more stable state, the system will perform work on the device, and đ W z 0 . This inequality is yet another statement of the second law. For example, in the case of the movable wall, we can apply an external force to the wall and move the system between states, where the pressures on both sides of the internal wall are different and satisfy the equality p 1 p 2 = F / A , where A is the area of the wall and F is the external force. Now, we can rewrite the first law of thermodynamics in the following form: d U đ Q = đ W z 0 , and thus obtain d U đ Q 0 . This form is more convenient than đ W z 0 because we do not need to create a new device for each case or introduce new parameters. Everything can be calculated from the system’s state if we know d U and đ Q . However, the form d U đ Q compares neighboring states and is therefore local. In non-equilibrium states, this form is insufficient for predicting the direction of spontaneous processes. In many non-equilibrium situations, there may be many stable local states. For example, in Section IV, we described the volumetric heating of a gas. There, we identified three different states that satisfied d U đ Q = 0 . The question is how to compare these states and determine which is the most stable. We propose in this paper to calculate the total work performed along the path, i.e., đ W z . If the total work is negative, it means that we have moved from a less stable state to a more stable state. Thus, the second law would be ( d U đ Q ) 0 . Because work depends on the chosen path, we additionally require that this work be maximal, i.e., max ( d U d Q ) 0 . The calculations must be performed under constant boundary conditions. In simple terms, the second law states that the direction of a spontaneous process involves removing some energy from the system by performing work on an external device that keeps the system in a less stable state.
This contribution, together with our previous works on the first law of non-equilibrium thermodynamics, constitutes a good starting point for applying the presented second law of non-equilibrium thermodynamics to systems undergoing chemical reactions in photo-reactors and flow reactors, Rayleigh-Benard convection, thermoosmosis, and, finally, lift force in hydrodynamic flows.

Author Contributions

All authors have contributed substantially to the reported work including the methodology, conceptualization, and writing; P.J.Ż. prepared the visualization. All authors have read and agreed to the published version of the manuscript.

Funding

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 Sklodowska-Curie Grant Agreement No. 847413. This project was also part of an international co-financed initiative funded by the program of the Minister of Science and Higher Education, entitled “PMW”, for the years 2020–2024, Agreement No. 5005/H2020-MSCA-COFUND/2019/2. This research was funded in part by the Polish National Science Center (Opus Grant No. 2022/45/B/ST3/00936).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Callen, H.B. Thermodynamics and an Introduction to Thermostatistics; Wiley: New York, NY, USA, 1991. [Google Scholar]
  2. Prigogine, I. Introduction to Thermodynamics of Irreversible Processes; Interscience Publishers: New York, NY, USA, 1967. [Google Scholar]
  3. Oono, Y.; Paniconi, M. Steady state thermodynamics. Prog. Theor. Phys. Suppl. 1998, 130, 29–44. [Google Scholar] [CrossRef]
  4. Sekimoto, K. Langevin equation and thermodynamics. Prog. Theor. Phys. Suppl. 1998, 130, 17–27. [Google Scholar] [CrossRef]
  5. Hatano, T.; Sasa, S.i. Steady-state thermodynamics of Langevin systems. Phys. Rev. Lett. 2001, 86, 3463. [Google Scholar] [CrossRef]
  6. Sasa, S.i.; Tasaki, H. Steady state thermodynamics. J. Stat. Phys. 2006, 125, 125–224. [Google Scholar] [CrossRef]
  7. Guarnieri, G.; Morrone, D.; Çakmak, B.; Plastina, F.; Campbell, S. Non-equilibrium steady-states of memoryless quantum collision models. Phys. Lett. A 2020, 384, 126576. [Google Scholar] [CrossRef]
  8. Boksenbojm, E.; Maes, C.; Netočnỳ, K.; Pešek, J. Heat capacity in nonequilibrium steady states. EPL (Europhys. Lett.) 2011, 96, 40001. [Google Scholar] [CrossRef]
  9. Netz, R. Approach to equilibrium and nonequilibrium stationary distributions of interacting many-particle systems that are coupled to different heat baths. Phys. Rev. E 2020, 101, 022120. [Google Scholar] [CrossRef]
  10. Mandal, D. Nonequilibrium heat capacity. Phys. Rev. E 2013, 88, 062135. [Google Scholar] [CrossRef] [PubMed]
  11. 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]
  12. Jona-Lasinio, G. Thermodynamics of stationary states. J. Stat. Mech. Theory Exp. 2014, 2014, P02004. [Google Scholar] [CrossRef]
  13. Speck, T.; Seifert, U. Integral fluctuation theorem for the housekeeping heat. J. Phys. A Math. Gen. 2005, 38, L581. [Google Scholar] [CrossRef]
  14. Mandal, D.; Jarzynski, C. Analysis of slow transitions between nonequilibrium steady states. J. Stat. Mech. Theory Exp. 2016, 2016, 063204. [Google Scholar] [CrossRef]
  15. Maes, C.; Netočnỳ, K. Nonequilibrium calorimetry. J. Stat. Mech. Theory Exp. 2019, 2019, 114004. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Litniewski, M.; Makuch, K.; Żuk, P.J.; Maciołek, A.; Hołyst, R. Continuous nonequilibrium transition driven by heat flow. Phys. Rev. E 2021, 104, 024102. [Google Scholar] [CrossRef]
  17. Glansdorff, P.; Prigogine, I. On a general evolution criterion in macroscopic physics. Physica 1964, 30, 351–374. [Google Scholar] [CrossRef]
  18. Maes, C.; Netočnỳ, K. A nonequilibrium extension of the Clausius heat theorem. J. Stat. Phys. 2014, 154, 188–203. [Google Scholar] [CrossRef]
  19. Ruelle, D.P. Extending the definition of entropy to nonequilibrium steady states. Proc. Natl. Acad. Sci. USA 2003, 100, 3054–3058. [Google Scholar] [CrossRef]
  20. Sasa, S.i. Possible extended forms of thermodynamic entropy. J. Stat. Mech. Theory Exp. 2014, 2014, P01004. [Google Scholar] [CrossRef]
  21. Komatsu, T.S.; Nakagawa, N. Expression for the stationary distribution in nonequilibrium steady states. Phys. Rev. Lett. 2008, 100, 030601. [Google Scholar] [CrossRef]
  22. Komatsu, T.S.; Nakagawa, N.; Sasa, S.i.; Tasaki, H. Steady-state thermodynamics for heat conduction: Microscopic derivation. Phys. Rev. Lett. 2008, 100, 230602. [Google Scholar] [CrossRef]
  23. Komatsu, T.S.; Nakagawa, N.; Sasa, S.i.; Tasaki, H. Entropy and nonlinear nonequilibrium thermodynamic relation for heat conducting steady states. J. Stat. Phys. 2011, 142, 127–153. [Google Scholar] [CrossRef]
  24. 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]
  25. Nakagawa, N.; Sasa, S.i. Liquid-gas transitions in steady heat conduction. Phys. Rev. Lett. 2017, 119, 260602. [Google Scholar] [CrossRef] [PubMed]
  26. Nakagawa, N.; Sasa, S.i. Global Thermodynamics for Heat Conduction Systems. J. Stat. Phys. 2019, 177, 825–888. [Google Scholar] [CrossRef]
  27. Sasa, S.i.; Nakagawa, N.; Itami, M.; Nakayama, Y. Stochastic order parameter dynamics for phase coexistence in heat conduction. Phys. Rev. E 2021, 103, 062129. [Google Scholar] [CrossRef] [PubMed]
  28. Nakagawa, N.; Sasa, S.i. Unique extension of the maximum entropy principle to phase coexistence in heat conduction. Phys. Rev. Res. 2022, 4, 033155. [Google Scholar] [CrossRef]
  29. Van den Broeck, C.; Esposito, M. Ensemble and trajectory thermodynamics: A brief introduction. Phys. A Stat. Mech. Its Appl. 2015, 418, 6–16. [Google Scholar] [CrossRef]
  30. Rao, R.; Esposito, M. Conservation laws and work fluctuation relations in chemical reaction networks. J. Chem. Phys. 2018, 149, 245101. [Google Scholar] [CrossRef]
  31. Falasco, G.; Esposito, M. Macroscopic Stochastic Thermodynamics. arXiv 2024, arXiv:2307.12406. [Google Scholar]
  32. Glansdorff, P.; Prigogine, I. Thermodynamic Theory of Structure, Stability and Fluctuations; Wiley-Interscience: New York, NY, USA, 1971. [Google Scholar]
  33. De Groot, S.R.; Mazur, P. Non-Equilibrium Thermodynamics; Dover Publications: New York, NY, USA, 1984. [Google Scholar]
  34. Jou, D.; Casas-Vázquez, J.; Lebon, G. Extended Irreversible Thermodynamics; Springer: Berlin, Germany, 1996. [Google Scholar]
  35. 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]
  36. 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]
  37. 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. [Google Scholar] [CrossRef]
  38. Makuch, K.; Hołyst, R.; Giżyński, K.; Maciołek, A.; Żuk, P. Steady-state thermodynamics of a system with heat and mass flow coupling. J. Chem. Phys. 2023, 159, 194113. [Google Scholar] [CrossRef] [PubMed]
  39. Hołyst, R.; Poniewierski, A. Thermodynamics for Chemists, Physicists and Engineers; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  40. Chandrasekhar, S. Hydrodynamic and hydromagnetic stability; Courier Corporation: North Chelmsford, MA, USA, 2013. [Google Scholar]
  41. Getling, A.V. Rayleigh-Bénard Convection: Structures and Dynamics; World Scientific: Singapore, 1998; Volume 11. [Google Scholar]
  42. Mizerski, K.A. Foundations of Convection with Density Stratification; Springer: Cham, Switzerland, 2021. [Google Scholar]
  43. Gruber, C.; Piasecki, J. Stationary motion of the adiabatic piston. Phys. A Stat. Mech. Its Appl. 1999, 268, 412–423. [Google Scholar] [CrossRef]
  44. Caprini, L.; Cerino, L.; Sarracino, A.; Vulpiani, A. Fourier’s law in a generalized piston model. Entropy 2017, 19, 350. [Google Scholar] [CrossRef]
Figure 1. Schematic of a gas confined between two parallel walls with area A . The system is assumed to be translationally invariant in the x and y directions. The thin interior wall at z w is treated as an internal constraint. We assume that the internal wall is diathermic, impenetrable, and can move freely. The inner wall divides the system into two compartments containing N 1 and N 2 gas particles ( N 1 + N 2 = N ). The outer wall placed at z 1 = 0 has a fixed temperature T 1 , and the outer wall placed at z 2 = L is kept at a fixed temperature T 2 > T 1 . The resulting linear temperature profile is indicated by the dash-dot line.
Figure 1. Schematic of a gas confined between two parallel walls with area A . The system is assumed to be translationally invariant in the x and y directions. The thin interior wall at z w is treated as an internal constraint. We assume that the internal wall is diathermic, impenetrable, and can move freely. The inner wall divides the system into two compartments containing N 1 and N 2 gas particles ( N 1 + N 2 = N ). The outer wall placed at z 1 = 0 has a fixed temperature T 1 , and the outer wall placed at z 2 = L is kept at a fixed temperature T 2 > T 1 . The resulting linear temperature profile is indicated by the dash-dot line.
Entropy 26 00713 g001
Figure 2. Gas column in gravity. A column of height L contains N moles of gas particles. The base ( z 1 ) and top ( z 2 ) are in contact with temperature reservoirs kept at two different temperatures: T 1 at the base, and T 2 at the top. The column’s base and top surface area equal A, and the system is translationally invariant in the x and y directions. (a) Impermeable and diathermal wall at z w separates the bottom and top sections with N 1 and N 2 particles of gas ( N 1 + N 2 = N ). (b) Permeable, thick ( L w ) diathermal wall at z w separates sections containing N 1 and N 2 particles of gas ( N 1 + N 2 = N ).
Figure 2. Gas column in gravity. A column of height L contains N moles of gas particles. The base ( z 1 ) and top ( z 2 ) are in contact with temperature reservoirs kept at two different temperatures: T 1 at the base, and T 2 at the top. The column’s base and top surface area equal A, and the system is translationally invariant in the x and y directions. (a) Impermeable and diathermal wall at z w separates the bottom and top sections with N 1 and N 2 particles of gas ( N 1 + N 2 = N ). (b) Permeable, thick ( L w ) diathermal wall at z w separates sections containing N 1 and N 2 particles of gas ( N 1 + N 2 = N ).
Entropy 26 00713 g002
Figure 3. Schematic of a gas confined between two parallel walls of large area A subjected to a homogeneous external energy input of density λ . The external wall located at z 1 = 0 is maintained at a temperature T 1 , while the external wall located at z 2 = L is maintained at a temperature T 2 > T 1 . The mobile, thin internal wall is adiabatic, resulting in a jump in the temperature profile.
Figure 3. Schematic of a gas confined between two parallel walls of large area A subjected to a homogeneous external energy input of density λ . The external wall located at z 1 = 0 is maintained at a temperature T 1 , while the external wall located at z 2 = L is maintained at a temperature T 2 > T 1 . The mobile, thin internal wall is adiabatic, resulting in a jump in the temperature profile.
Entropy 26 00713 g003
Figure 4. The difference between the pressures in the two subsystems, normalized by the equilibrium pressure p 0 = N k b T 1 / V , where T 1 is the temperature of the external wall located at z 1 = 0 , as a function of z w for three values of λ . In (a), the temperatures at both external walls are equal ( r = T 2 / T 1 = 1 ). In (b), the temperatures at both external walls are different ( r = T 2 / T 1 = 1.5 ). The vertical lines mark the positions of the steady states. In (a), these are z w ( 1 ) for λ = 1 , and z w ( 1 ) , z w ( 2 ) , and z w ( 3 ) for λ = 40 and 300. In (b), they are z w ( 1 ) for λ = 1 and 40, and z w ( 1 ) , z w ( 2 ) , and z w ( 3 ) for λ = 300 . Here, λ = λ L 2 / ( κ T 1 ) .
Figure 4. The difference between the pressures in the two subsystems, normalized by the equilibrium pressure p 0 = N k b T 1 / V , where T 1 is the temperature of the external wall located at z 1 = 0 , as a function of z w for three values of λ . In (a), the temperatures at both external walls are equal ( r = T 2 / T 1 = 1 ). In (b), the temperatures at both external walls are different ( r = T 2 / T 1 = 1.5 ). The vertical lines mark the positions of the steady states. In (a), these are z w ( 1 ) for λ = 1 , and z w ( 1 ) , z w ( 2 ) , and z w ( 3 ) for λ = 40 and 300. In (b), they are z w ( 1 ) for λ = 1 and 40, and z w ( 1 ) , z w ( 2 ) , and z w ( 3 ) for λ = 300 . Here, λ = λ L 2 / ( κ T 1 ) .
Entropy 26 00713 g004
Figure 5. Schematic of a gas confined between two parallel moving walls of large area A . The external wall located at z 1 = 0 is fixed and is kept at a temperature T 1 , while the external wall located at z 2 = L is moving with velocity v 2 and kept at a temperature T 2 > T 1 . The mobile, thin internal wall is adiabatic, resulting in the temperature profile jump.
Figure 5. Schematic of a gas confined between two parallel moving walls of large area A . The external wall located at z 1 = 0 is fixed and is kept at a temperature T 1 , while the external wall located at z 2 = L is moving with velocity v 2 and kept at a temperature T 2 > T 1 . The mobile, thin internal wall is adiabatic, resulting in the temperature profile jump.
Entropy 26 00713 g005
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

Hołyst, R.; Żuk, P.J.; Maciołek, A.; Makuch, K.; Giżyński, K. Direction of Spontaneous Processes in Non-Equilibrium Systems with Movable/Permeable Internal Walls. Entropy 2024, 26, 713. https://doi.org/10.3390/e26080713

AMA Style

Hołyst R, Żuk PJ, Maciołek A, Makuch K, Giżyński K. Direction of Spontaneous Processes in Non-Equilibrium Systems with Movable/Permeable Internal Walls. Entropy. 2024; 26(8):713. https://doi.org/10.3390/e26080713

Chicago/Turabian Style

Hołyst, Robert, Paweł J. Żuk, Anna Maciołek, Karol Makuch, and Konrad Giżyński. 2024. "Direction of Spontaneous Processes in Non-Equilibrium Systems with Movable/Permeable Internal Walls" Entropy 26, no. 8: 713. https://doi.org/10.3390/e26080713

APA Style

Hołyst, R., Żuk, P. J., Maciołek, A., Makuch, K., & Giżyński, K. (2024). Direction of Spontaneous Processes in Non-Equilibrium Systems with Movable/Permeable Internal Walls. Entropy, 26(8), 713. https://doi.org/10.3390/e26080713

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