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
is a function of five state parameters, where
is the non-equilibrium entropy,
V is the volume,
N is the number of particles, and
and
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:
The
and
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,
, is minimized with respect to
x (the variable describing the internal constraint) at constant
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
. Rewriting this equation using internal energy and entropy yields
. In the infinitesimal form, we obtain
. Finally, the net heat (the heat that enters or leaves the system and changes the internal energy) is
. In general, the second law states that in spontaneous processes,
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
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
and
as functions of the state parameters for a single compartment [
35,
36,
37] and point out that
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
and
, 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, . In the presence of a heat flux, quantities specifying equilibrium thermodynamic states, such as number density , pressure p, or temperature T, become space-dependent. The profiles , , and 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 .
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,
where
is the average number density and
is given by
The new state parameter
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
depends on the specific temperature and number density profile. In the case of constant pressure [
35], the equations of irreversible thermo-hydrodynamics yield
,
, and as a result,
Because mapping gives us the same formal structure as we know from equilibrium, the variable
conjugate to
is given by
where the constant
is chosen such that
for
gives the equilibrium expression for entropy [
1]. The internal energy of a non-equilibrium steady state is thus a function of three state parameters
, with the thermodynamic relations
Binary mixture. Mixing two ideal gases,
a and
b, with
and
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
, we have two effective degrees of freedom
and
, which, however, are not independent. These parameters allow us to write the internal energy and pressure in the same form as in equilibrium:
where
and
are the average number densities of the first and second components of the mixture, respectively. The effective degrees of freedom are obtained as
where
is the average number fraction of the
i component and
is its profile;
in the absence of chemical reactions.
and
are related via
. The non-equilibrium entropy
has the same form as in equilibrium but with
T replaced with
and
with
. It is the sum of the entropy of the two components of the mixture considered separately and the entropy of mixing:
Here,
may depend on
and
. Again, we choose
such that
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:
.
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
Three new state parameters describe the van der Waals gas in the steady state. Besides the effective temperature
, defined by the same expression (
4) as for the ideal gas, we have the effective interaction parameter
given by
and
defined by the formula
Because Equation (
11) has the same structure as in equilibrium, the non-equilibrium entropy
has the same form as in equilibrium but with
T replaced with
,
a with
, and
b with
. Specifically,
Thus, the internal energy of the van der Waals gas in the non-equilibrium steady state is a function of five state parameters:
.
Note that for all systems discussed above, the non-equilibrium entropy is part of the total entropy , defined as the integral of the volume entropy density 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
or a change in the distance between the confining walls
L, the energy changes only through mechanical work and heat flow:
where
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
. For an ideal gas, we have
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,
and for a fixed number of components
and
, we have
where
and
, with
given by Equation (
10).
Additional terms in the net heat also occur for the van der Waals gas, where
so that
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,
and
, and volumes,
and
. We choose the volume of one subsystem,
, as the parameter representing the constraint. The total energy is the sum of the energies of the two subsystems:
where … denotes possible other new state variables, such as
for the binary mixture of ideal gases or
and
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:
.
In the case of a monoatomic ideal gas, the effective temperatures of both subsystems are [
35]
Using the explicit forms of
and
, we can see that this condition can be obtained from the following minimum principle:
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
in the two subsystems. For the specific case of an ideal gas, we have
Demanding
and solving for
gives us the position of the internal wall at the stationary state for given temperatures
and
:
We note that for the diathermal internal wall, there is only a single solution for
for all possible values of
. In
Section 4, we consider an adiabatic internal wall with volumetric heating and show that, in this case, there might be more solutions for
, 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
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]
where
is the height of the gas column;
is the coordinate of the column base;
is the coordinate of the column top;
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
is the renormalized mass defined with the help of the potential gravitational energy of the gas [
36]
which is contained inside the column
The state parameter
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
The total internal energy of the gas column
consists of thermal and gravitational contributions. The thermal contribution (through the mapping procedure) defines the renormalized temperature
The average pressure inside the column is the same as in the case without gravity (
3)
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
and bottom
of the column are different. The presented formulas are consistent with the functional form of
(
26) and the internal energy differential
In the above, we distinguish two separate ways of changing the internal energy. The first is through heat,
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
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
particles, and the top part contains
particles of an ideal gas. We separate the compartments with a thin, diathermal wall at
. We assume there are no changes in the column cross-sectional area (
) and set this contribution to the energy differential to 0. We rewrite the energy balances for the bottom (
) and top (
) segments explicitly to account for the displacement of the separating wall
The sign in front of
is consistent with the direction of changes in
. 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
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
The internal energy of the van der Waals gas in a gravitational field is given by
The pressures are as follows:
is the average pressure,
is the pressure at the top of the column, and
is the pressure at the bottom of the column. The energy differential has the form
where we distinguish the heat differential, including the gravitational effects,
and the mechanical work
Likewise, for the ideal gas, we write the second law of stationary thermodynamics as
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
with its own mass
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
. 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
, 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
. We do not elucidate this term and treat it formally because it will cancel out during further derivation. The amended energy differential is
We rewrite the energy differential explicitly for both compartments. Compartment 1 is at the bottom, so its
stays fixed while
of the upper compartment 2 is subject to change:
Here, we set
and
, 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
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
particles from the bottom to the top compartment. The passage starts when
particles leave compartment 1, which generates an energy gain in the demon’s account:
Inserting
particles into the pore end at temperature
requires, in addition to the volumetric work, the thermal energy
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
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,
and undergoes decompression, performing work
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
is released and added to the demon’s account. Simultaneously, the particle needs to equilibrate to the bulk conditions by taking energy
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
, 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:
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
in the gravitational field:
With the use of
and
, we can write the second law of stationary thermodynamics:
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
, we use the assumption that the wall is an excellent heat conductor so that
. This means that the demon had to support the cost of the isothermal expansion
which we can explicitly insert into the second law of stationary thermodynamics:
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:
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 (
; 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
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:
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
. 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
and substitute
to find
and Archimedes’ principle emerges:
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
. In the steady state, the pressure
p and, consequently, the energy density
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:
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
and
, we find
After the mapping procedure given by Equation (
3), we obtain the effective temperature
where
is a dimensionless function given by
In the limit
, the effective temperature
reduces to the value given by Equation (
5) for an ideal gas in the heat flux induced by the temperature gradient
(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.,
), the system reduces to that described in our previous article [
16]. The entropy
is given by Equation (
14) with
. For a fixed number of molecules
N, the energy change is determined by the control parameters
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]
Movable wall. We introduce a constraint into the system in the form of a movable adiabatic wall parallel to the bounding walls located at
(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.,
, at the internal wall. Such a wall does not allow for energy transport between compartments as a result of molecular collisions.
The temperature profiles
and
for subsystems (1) and (2) with
and
molecules, respectively, are
The mapping procedure yields
where
Thus, the energy of the subsystems of
and
particles under the constraint is
with
.
Our second law of non-equilibrium thermodynamics states that
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
where
,
r is the temperature ratio
, and
. At equilibrium, i.e., for
and
, the internal wall is located precisely in the middle of the system at
. For
but
, the position of the internal wall is given by Equation (
25).
An analysis of Equation (
68) shows that the number of solutions for
varies from one to three depending on the values of
and
r. A typical course of the variability of the function
is shown in
Figure 4, compared to the case of no temperature gradient (
). For
and small values of
, the curve
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
with three zero-crossing points at
,
, and
, as illustrated in
Figure 4a. Among these three solutions, further analysis shows that locally, both
and
are stable, whereas
is unstable. Moreover, the work required to change the position of the wall from
to
is equal to zero (
). 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
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
axis. Above the critical value
, at which the maximum of the function
touches the
axis, three zero crossing points
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
or
but push it away from
. 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
to position
allows us to determine which steady state is globally stable. If the work
done on the system during this process is negative, the final position
corresponds to a globally stable steady state. If this work is positive, the initial position
is globally stable. We calculate this work numerically up to machine accuracy. We find that for
, due to the asymmetry, the work is positive for all values of
. Therefore,
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
moves at a constant speed
in the direction of the imposed shear force. The wall at
is fixed. The velocity profile is linear:
. Inside the system there is an adiabatic massless wall at location
, which moves at the local speed of the fluid
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
for subsystem
is given by
where
is the excess heat in subsystem
i (
),
is the excess dissipation in subsystem
i, and
is the volumetric work performed during the transition in subsystem
i. In this non-equilibrium case, we have
. Thus, the entropy change is due to both excess heat exchange and dissipation due to shear within the system.
is the change in volume in subsystem
i. From the equation of irreversible thermodynamics [
38], it follows that the change in kinetic energy,
in subsystem
i (
) has the following form:
The excess work performed by the surface forces,
, 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:
From this equation, we obtain the expected form of the second law:
The equality in the above expresion defines the condition for the stationary position of the wall, given by the equality of pressures
in the two subsystems.
6. Discussion
The system at equilibrium, which exchanges heat with the environment, satisfies the inequality 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
. We can write it in a more general form as
. 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
and different forms of work
, the following inequality should govern spontaneous processes:
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, has a well-defined value for fixed . 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 . We introduce an extra parameter x, usually in the form of a constraint in the system (like a movable wall), and state that is maximized as a function of x for fixed state parameters (internal energy, volume, and number of particles, respectively). Similarly, when the system has a constant temperature, the Helmholtz free energy is minimized as a function of x (the variable describing the internal constraint) at constant (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 . Let us define the work performed on the system by this external device as . In the process of transitioning between states, the first law of thermodynamics must be obeyed (conservation of energy). Thus, we have . 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., . 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 . 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 , 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: , and thus obtain . This form is more convenient than 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 and . However, the form 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 . 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., . 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 . Because work depends on the chosen path, we additionally require that this work be maximal, i.e., . 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.