Global Thermodynamics for Heat Conduction Systems
Naoko Nakagawa, Shin-ichi Sasa
TL;DR
This paper introduces a new concept of global temperature for non-uniform heat conduction systems, extending thermodynamics to such systems and predicting interface temperature deviations and super-cooled gas formation.
Contribution
It develops a global thermodynamics framework and a variational principle for phase coexistence in heat conduction systems, extending classical equilibrium concepts.
Findings
Global temperature concept successfully extends thermodynamics.
Predicts deviation of interface temperature from equilibrium.
Identifies stable super-cooled gas near the interface.
Abstract
We propose the concept of global temperature for spatially non-uniform heat conduction systems. With this novel quantity, we present an extended framework of thermodynamics for the whole system such that the fundamental relation of thermodynamics holds, which we call "global thermodynamics" for heat conduction systems. Associated with this global thermodynamics, we formulate a variational principle for determining thermodynamic properties of the liquid-gas phase coexistence in heat conduction, which corresponds to the natural extension of the Maxwell construction for equilibrium systems. We quantitatively predict that the temperature of the liquid-gas interface deviates from the equilibrium transition temperature. This result indicates that a super-cooled gas stably appears near the interface.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
∎
11institutetext: N. Nakagawa 22institutetext: Department of Physics, Ibaraki University, Mito 310-8512, Japan
22email: [email protected] 33institutetext: S. Sasa 44institutetext: Department of Physics, Kyoto University, Kyoto 606-8502, Japan
44email: [email protected]
Global Thermodynamics for Heat Conduction Systems
Naoko Nakagawa
Shin-ichi Sasa
Abstract
We propose the concept of global temperature for spatially non-uniform heat conduction systems. With this novel quantity, we present an extended framework of thermodynamics for the whole system such that the fundamental relation of thermodynamics holds, which we call “global thermodynamics” for heat conduction systems. Associated with this global thermodynamics, we formulate a variational principle for determining thermodynamic properties of the liquid-gas phase coexistence in heat conduction, which corresponds to the natural extension of the Maxwell construction for equilibrium systems. We quantitatively predict that the temperature of the liquid-gas interface deviates from the equilibrium transition temperature. This result indicates that a super-cooled gas stably appears near the interface.
Keywords:
Thermodynamics Heat conduction Liquid-gas transition Super-cooled gas
††journal: Journal of Statistical Physics
Contents
-
3 Global thermodynamics for single-phase systems in the linear response regime
-
3.4 Correspondence of global thermodynamic quantities to equilibrium quantities
-
5 Variational principle for determining the liquid-gas interface position
-
6 Properties of the steady states with a liquid-gas interface
-
7 Global thermodynamics for systems with a liquid-gas interface
-
9 Other assumptions leading to the results of the variational principle
-
10.1 Global thermodynamic quantities in an arbitrarily shaped container
-
10.2 Equivalence of non-equilibrium global quantities with equilibrium quantities
1 Introduction
The behavior of liquids and gases close to equilibrium have been extensively studied for a long time. As a macroscopic universal theory describing these phases of matter, hydrodynamics is believed to be well-established Landau-Lifshitz-Fluid , and the connection of the hydrodynamics with the classical and quantum mechanics of atoms have been discussed for over a century Zubarev ; Mclennan . Nevertheless, in this paper, we construct a new universal theory for thermodynamic properties in the linear response regime. There are two main messages. First, a new concept, global temperature, is found, with which a novel framework of global thermodynamics is constructed to describe the whole of non-uniform non-equilibrium systems with local equilibrium thermodynamics. This outcome provides a fresh viewpoint for the description of systems out of equilibrium. Second, this formulation provides non-trivial quantitative predictions. As an example, let us consider pure water under a pressure of Pa, where two heat baths of temperature K and K are in contact with the sides of the system. See Fig. 1 for an illustration. Recall that the liquid-gas transition temperature is K. Based on global thermodynamics, in this paper, we predict that the interface temperature of the liquid-gas coexistence is K. This result means that super-cooled gas stably appears in heat conduction, which may be tested in experiments.
In the remaining part of this introduction, we first present a brief summary of development in non-equilibrium statistical mechanics, confirming that the phenomenon described above has never been discussed by established theories. We then provide a review of extended frameworks of thermodynamics so that readers can understand how the global thermodynamics proposed in this paper is different from previous frameworks. At the end of the introduction, we summarize the achievements of this paper.
1.1 Non-equilibrium statistical mechanics
In the universal theory of macroscopic irreversible dynamics built by Onsager, a steady-state current is assumed to be linear in the thermodynamic forces Onsager1931 . The linear coefficients form a non-negative symmetric matrix, which is referred to as the Onsager matrix. This theory with a variational principle for determining a current and/or a thermodynamic force is called irreversible thermodynamics, which is valid near equilibrium Groot-Mazur . Onsager theory can be interpreted as a universal theory for the dynamical fluctuations of thermodynamic variables, where the properties of the Onsager matrix are connected to the stability of the equilibrium state and the time-reversal symmetry in microscopic systems. In short, this fluctuation theory is represented by a simple stochastic process with a symmetry property Onsager1931 ; Schmitz . In accordance with this theory, statistical mechanics of trajectories has been studied based on microscopic dynamics, which leads to the expression of probability densities, linear and non-linear response formulas, and macroscopic deterministic dynamics Zubarev ; Mclennan ; Kubo ; Nakano ; Zwanzig ; Mori ; Kawasaki-Gunton .
After these developments, in the last two decades, our understanding of phenomena related to thermodynamics has progressed greatly because of the following two reasons. First, the development of experimental techniques for the measurement and manipulation of biological molecular machines Ashkin ; Bustamante ; Svoboda ; Chu ; Noji naturally leads to the extension of thermodynamics to mesoscopic scales, which is now called stochastic thermodynamics SeifertRPP ; Sekimoto-book ; Kalges-Just-Jarzynski . Second, the newly discovered relations, such as the fluctuation theorem Evans-Cohen-Morriss and the Jarzynski equality JarzynskiPRL , have emerged as universal Gallavotti ; Kurchan ; LS ; Maes ; Crooks ; JarzynskiDFT , providing a new starting point for the re-organization of previous theories with greatly simplified derivations of the formulas CrooksNRL ; KN ; KNST-rep ; Maes-rep ; SasaFluid . These two developed directions, stochastic thermodynamics and new universal relations, are related to each other and have given rise to a new connection with large deviation theory Bodineau-Derrida ; BertiniETAL ; Maes-LD ; Nemoto ; Bertini-rev , and information theory Information .
These results are quite useful when we specify a phenomenon under study. Indeed, the phase coexistence in heat conduction can be studied through a tough calculation based on the non-equilibrium statistical mechanics SNIN . However, we could not conjecture the non-trivial nature of this phenomenon immediately from principles of non-equilibrium statistical mechanics. For example, although one expects a variational formulation for the phase coexistence in heat conduction by recalling the minimum entropy production principle min-ent , we cannot obtain thermodynamic properties directly from the principle because the variational principle basically determines the statistical ensemble in the linear response regime as the minimizer of the entropy production Klein . In order to develop a variational principle for thermodynamic quantities, we have to start with an extended framework of thermodynamics.
1.2 Extended frameworks of thermodynamics
Equilibrium thermodynamics provides a unified description of thermodynamic properties of materials at equilibrium. It also formalizes the second law, which leads to a variational principle for determining equilibrium states Callen ; Prigogine-Kondepudi . The variational principle naturally suggests the law of fluctuation of thermodynamic variables Einstein , which is formulated as a large deviation theory. From this viewpoint, a framework of statistical mechanics may be constructed in a consistent manner with the fluctuation theory of thermodynamic variables Oono .
Therefore, in order to establish a universal theory for thermodynamic properties out of equilibrium, it is natural to consider an extended framework of thermodynamics. A naive attempt is to extend the equilibrium fundamental relation of thermodynamics
[TABLE]
for the case of simple fluids, where is the Helmholtz free energy and is the entropy. Examples of such attempts can be seen in Refs. Keizer ; Eu ; Jou-book . The heart of the problem for the extension is to confirm the two conditions: First, the theory is self-consistent and self-contained; second, new predictions specific to the extension are presented. However, the two conditions are not confirmed in many studies. The derivation of the extended thermodynamics from the microscopic theory, regardless of its importance, does not make sense unless the extended framework satisfies the above two conditions. Here, we do not give a complete review of previous studies on the extended framework of thermodynamics, but study heat conduction systems from the viewpoint of the extended framework of thermodynamics.
Heat conduction is described as a spatially extended system in which thermodynamic variables slowly vary in space. The local subsystems, which are small but macroscopic, are regarded as a local equilibrium state Landau-Lifshitz-Fluid ; Groot-Mazur ; Bedeauz86 . If we followed this standard description, we could not have a strong motivation to seek a thermodynamic framework because all thermodynamic properties can be calculated by the heat conduction equation and the local equilibrium thermodynamics. Nevertheless, we still have two possibilities for considering an extended framework of thermodynamics.
The first attempt is to go beyond the local equilibrium thermodynamics. In this approach, which may be adapted in Refs. Jou ; Sasa-Tasaki , the equation of state for the local subsystems is modified to contain the influence of the heat flux. Since such a contribution is quite small in the linear response regime, the theory is not useful even if it is correct. Although this approach may still be effective for spatially homogeneous driven systems, in which intensive parameters may describe the balance of extensive variables Sasa-Tasaki ; Hayashi-Sasa ; Bertin ; Seifert-contact ; Dickman , we do not deal with such systems in this paper.
In the second approach, we retain the standard description for a spatially extended system with local equilibrium thermodynamics. On this basis, we then seek a thermodynamic framework for the whole system. As an example, we first consider the extension of the second law, by which a state variable “entropy” is defined, and we then derive the thermodynamic relation which corresponds to an extension of the fundamental relation of thermodynamics. A concrete procedure of the first step was proposed by introducing the concept of excess heat Jou-book ; Landauer ; Oono-Paniconi . This idea was also studied from semi-macroscopic and microscopic theories Hatano-Sasa ; Ruelle ; KNST ; NN ; Jona-thermo ; Maes-thermo ; Spinney-Ford . However, the fundamental relation of thermodynamics has not been derived from the extended entropy. One reason may be that an interesting phenomenon associated with the extended entropy was not addressed. Nevertheless, one can continue to seek a possible framework without considering certain phenomena. Indeed, for a specific model of heat conduction, the extended entropy was numerically estimated Chiba-Nakagawa , which suggests that the extended entropy is close to the spatial integration of the local equilibrium entropy density. Since the extended entropy can be obtained in experiments, a natural question is “What is temperature?” If the fundamental relation of thermodynamics is formulated, there should be a temperature satisfying it. We consider this question seriously, putting aside specific phenomena.
1.3 Summary of results
We first focus on single-phase systems (either liquid or gas) in heat conduction. By assuming local equilibrium thermodynamics, we define the entropy and the Helmholtz free energy in heat conduction by the spatial integration of the local density fields corresponding to these variables. The pressure field is homogeneous in space. Then, the problem is to find a temperature satisfying the fundamental relation
[TABLE]
in the linear response regime. In §3, we solve this problem by defining as the kinetic temperature averaged over particles in the system. We call this temperature global temperature. It should be noted that the local temperature , which depends on the position , satisfies the thermodynamic relations for each point . In contrast to the local relations, (1.2) is a global relation applied for the whole system as if the system is at equilibrium. We call such a thermodynamic framework global thermodynamics. This formulation is also interpreted as a mapping of each heat conduction system to an equilibrium system through the novel quantity of the global temperature, .
The formulation for a single-phase system is indeed derived from fluid dynamics with local equilibrium thermodynamics. That is, this formulation is interpreted as a different formulation for describing thermodynamic quantities in the linear response regime. The prediction by using global thermodynamics can also be predicted by hydrodynamics with local equilibrium thermodynamics, in principle. Here, we go one step further. We consider a thermodynamic phenomenon that cannot be described by the standard hydrodynamics with local equilibrium thermodynamics. This is the phenomenon of phase coexistence in heat conduction. In §4, we study this phenomenon and we show how existing theories are not appropriate for determining the thermodynamic properties. The essential point is that the condition for the connection of the two phases is outside of the local equilibrium thermodynamics.
We study the phase coexistence in heat conduction with the framework of global thermodynamics. We first formulate a variational principle for determining the local temperature of the liquid-gas interface. The idea is quite simple. Fixing the global temperature of the whole system, we naturally extend the variational principle for equilibrium systems to that for heat conduction systems. This idea was proposed in our previous paper NS . Remarkably, by using the solution of the variational equation, in §5, we derive a universal relation among the interface temperature , the equilibrium transition temperature , the global temperature , and the mean temperature of the two heat baths , which is
[TABLE]
We call this relation the temperature relation. From the temperature relation, we find that the temperature of the liquid-gas interface deviates from the equilibrium transition temperature. That is, super-cooled gas stably appears near the interface in heat conduction. This is a qualitatively new phenomenon that has never been considered in previous studies.
Since the steady state as determined by the variational principle is expressed as a function of the global temperature, the prediction of measurable quantities for given conditions is indirect. Thus, in §6, we re-express all quantities in the steady state in terms of the two temperatures of the heat baths. We present several formulas of the interface temperature by directly using measurable quantities. Furthermore, we illustrate examples of quantitative results for a van der Waals fluid and pure water.
For the steady state determined by the variational principle, we further develop thermodynamics for heat conduction systems with phase coexistence. First, in §7, we derive the fundamental relation associated with the Gibbs free energy . At first sight, the result does not seem to contain non-equilibrium extensions. This is because is not differentiable at the transition point at equilibrium. By performing a careful analysis near equilibrium, we find that the fundamental relation holds in an appropriate equilibrium limit. The heat capacity and the compressibility, which are singular at equilibrium, are also obtained as a regularized form while breaking the additivity. In §8, we derive the fundamental relation associated with the Helmholtz free energy . Since the free energy is defined for the coexistence phase in equilibrium cases, its non-equilibrium extension can be written as a perturbation from the equilibrium form. We derive this expression explicitly.
It should be noted that the quantitative prediction is made based on a fundamental assumption of the variational principle. As is often observed in universal theories, one may replace the fundamental assumption by another one. In §9, we formulate the theory starting from assumptions other than the variational principle. For example, when we assume the fundamental relation of thermodynamics for the whole system by using the global temperature, we can derive the results of the variational principle. As another example, one may focus on how the volume change near the transition temperature at constant pressure exhibits the singularity in the equilibrium limit. Supposing the simplest form of a singularity, we can derive the results of the variational principle. These findings indicate that the theory itself possesses an elegant structure.
In §10, we extend the theory in §3 to cover steady states with an arbitrarily shaped container beyond the linear response regime, while restricting our attention to single-phase systems. Our theory quantitatively predicts a new relation among the global quantities in this setup.
2 Preliminaries
2.1 Equilibrium thermodynamics
We consider a macroscopic material at equilibrium. As the simplest example, we focus on a simple fluid whose thermodynamic state is characterized by temperature , the volume of a container, and the amount of material . For a system in contact with a heat bath of temperature which may be controlled externally, there exists a state variable , called entropy, which satisfies the Clausius equality
[TABLE]
for infinitesimal quasi-static heat from the heat bath. The infinitesimal change of internal energy of the material is determined by
[TABLE]
which is referred to as the first law of thermodynamics. is the infinitesimal quasi-static work required in the infinitesimal change of the volume , which is given by
[TABLE]
The substitution of (2.1) and (2.3) into (2.2) leads to the fundamental relation of thermodynamics:
[TABLE]
From this expression, we find that it is useful to consider as a function of . Indeed, a single function leads to all thermodynamic properties such as equation of state
[TABLE]
and heat capacity
[TABLE]
We then define chemical potential as
[TABLE]
Various thermodynamic functions equivalent to can be defined by the Legendre transformation of :
[TABLE]
The fundamental relations associated with these functions are
[TABLE]
Substituting into (2.13), we find
[TABLE]
This form together with (2.13) leads to the Gibbs-Duhem relation
[TABLE]
The extensivity of thermodynamic quantity leads to the concept of density defined as the quantity per unit volume, similarly to particle density . For instance, entropy density is defined as
[TABLE]
By using , which is obtained from (2.5), one may consider the entropy density as a function of , which is expressed as
[TABLE]
following the convention in thermodynamics. Similarly to the entropy density, for any extensive quantity , its density is defined as
[TABLE]
We here consider free energy density . Substituting into (2.14) and (2.15), we obtain
[TABLE]
where . From these, we have
[TABLE]
where , and are not independent, because , and are connected by the equation of state (2.5) such that
[TABLE]
Furthermore, we define extensive quantities per one particle, as
[TABLE]
Note that is equivalent to the chemical potential . We then have
[TABLE]
where is specific volume .
2.2 Setup of heat conduction systems
Throughout this paper, we consider a system of -particles, which are packed in a rectangular container with lengths of side , and , as shown in Fig. 2(a). We ignore the effect of gravity. and are fixed throughout this paper. is fixed for a constant volume system, or not fixed at constant pressure. We study heat conduction states driven by the temperature difference between two heat baths. As schematically described in Fig. 2(b), a heat bath of temperature is attached to the left end () and another heat bath of temperature to the right end (). Other four boundaries are thermally insulating. We take without loss of generality, and is assumed to be so small that the system reaches a unique nonequilibrium steady state in the linear response regime. In such an idealized steady state without any convection, the system is regarded as a one-dimensional system. Local states are homogeneous inside any section perpendicular to the axis, and therefore, local thermodynamic quantities are considered as functions of .
We introduce a dimensionless parameter that indicates the degree of non-equilibrium as
[TABLE]
where is the temperature difference and is the mean temperature of the heat baths, which are defined by
[TABLE]
When we focus on the linear response regime around an equilibrium state, we ignore the contribution of .
2.3 Local equilibrium thermodynamics
For such macroscopic non-equilibrium systems, the hypothesis of local equilibrium thermodynamics works well. That is, local thermodynamic quantities are assumed to satisfy local thermodynamic relations at each space and time Landau-Lifshitz-Fluid ; Groot-Mazur . Now, suppose that the temperature profile , the density profile and the pressure are determined by experimental observation. Any local thermodynamic quantities, such as and , are expressed as
[TABLE]
where functions and are determined in thermodynamics, as described in the previous subsection. According to (2.19), may be also written as via . Then, local equilibrium thermodynamics means
[TABLE]
which correspond to the relations (2.20), (2.21) and (2.22), respectively. For the quantities per one particle, we also have
[TABLE]
which correspond to the local version of the thermodynamic relations (2.25) and (2.26).
2.4 Global conditions for steady states
For steady state heat conduction, the local pressure satisfies
[TABLE]
for any and . Especially, for the system at the constant pressure ,
[TABLE]
holds for any . These equalities for the local pressure may be regarded as global relations because they are not obtained in the local thermodynamics. Furthermore, heat flux is uniform in as expressed in the form
[TABLE]
for any and , which may work as another global relation for the local quantities. In order to connect the heat flux with the local thermodynamic quantities, we assume a heat conduction equation
[TABLE]
in which is heat conductivity as a function of . We assume that there is no temperature gap at the boundaries, i.e.,
[TABLE]
Finally, the conservation of particle number is written as
[TABLE]
These global relations, (2.37), (2.38), (2.39), (2.42), together with the equation of state (2.5), the heat conduction equation (2.40), and the boundary condition (2.41), are sufficient to determine the profiles of the local temperature and the local density provided that the systems consist of a single phase.
As we will see in §4, there is a case where liquid and gas coexist in the container. For this special situation, the above global relations are not sufficient to determine the local states. This problem will be seriously studied in later sections.
2.5 Global thermodynamic quantities
Since thermodynamic quantities in heat conduction are not uniform in the space, they are basically described as local fields. Nevertheless, from the fact that the local states are governed by the global relations as explained in §2.4, we expect that some properties of the heat conduction systems are explained from a global point of view. Toward the characterization of global properties, we here define global thermodynamic quantities for heat conduction systems.
For extensive variables originally defined for equilibrium states, we define the following global quantities as an extension to those for heat conduction states:
[TABLE]
where we have used the same notation as that for equilibrium states. When we interpret as a state function of heat conduction states, we explicitly write which is in contrast to for equilibrium states. For instance, global entropy and global Gibbs free energy are defined as
[TABLE]
Here, we investigate the reference state dependence of these global thermodynamic quantities. In equilibrium thermodynamics, entropy density and internal energy density are defined up to an additive constant which depends on the choice of their reference state. Thus, the entropy density and the internal energy density possess this property. Concretely, let and be the shift of the additive arbitrary constants of entropy and energy per one particle, respectively, for the change of the reference state. We express this transformation by
[TABLE]
The shift of other thermodynamic quantities are induced as
[TABLE]
Note that local thermodynamic relations are invariant under the transformation.
Now, we consider the transformation of the global thermodynamic quantities, which are defined by the spatial integral of the local quantities. It is obvious that
[TABLE]
while for the global Gibbs free energy , we have
[TABLE]
The dependence of is far from trivial. Here, we introduce the global temperature such that the transformation is written as
[TABLE]
Explicitly, is given as
[TABLE]
which means that global temperature corresponds to the kinetic temperature averaged over particles. The transformation (2.53) suggests the consistency among , and . Indeed, we will show the global thermodynamic relations for these quantities. Similarly, we also define the global chemical potential as
[TABLE]
Since , the global chemical potential corresponds to the Gibbs free energy per one particle.
3 Global thermodynamics for single-phase systems
in the linear response regime
In this section, we restrict ourselves to a single phase where the heat conduction system is occupied by either liquid or gas at the constant pressure . When the environmental parameters , , and are fixed, the steady state is uniquely determined from the equation of state (2.5), the heat conduction equation (2.40), and the conservation law (2.42). That is, the values of are determined by
[TABLE]
Since the local thermodynamic quantities obey equilibrium thermodynamics, the profile leads to any local thermodynamic quantities such as and . Then, the global quantities such as , , and are calculated for the steady state.
Now, we show that the global thermodynamic functions satisfy
[TABLE]
This means that the global free energy and entropy, which are functions of , are expressed as the equilibrium free energy and the equilibrium entropy when we use the global temperature . This result leads to relations among the global quantities as
[TABLE]
which represents the fundamental relation of thermodynamics extended to heat conduction states in the linear response regime. We call such a framework the global thermodynamics for heat conduction systems. These results indicate that, as schematically shown in Fig. 3, the global nature of the heat conduction system is equivalent to that of the equilibrium system. Any global quantity defined by (2.43) is connected to a corresponding equilibrium quantity as
[TABLE]
3.1 Proof of (3.4) and (3.5)
When the environmental parameters are slightly changed to , the solution of the equations (3.1), (3.2), and (3.3) is modified slightly. We express the corresponding change as
[TABLE]
which leads to the change of any local thermodynamic quantity. For instance, the change of the local Gibbs free energy density is given by
[TABLE]
The change of the local quantities brings the change of the global quantities as
[TABLE]
where
[TABLE]
Here, since (2.34) holds for the local densities, the variation of each density satisfies a thermodynamic relation such as
[TABLE]
By substituting the local relations (2.32) and (3.15) into (3.14), we have
[TABLE]
where and we have defined
[TABLE]
Note that is the particle number in and that is estimated as . is the deviation of the local temperature from the global temperature , and is estimated as .
Since any local thermodynamic quantity is regarded as , we may estimate . Then, the integrant in (3.16) is estimated as
[TABLE]
where for and fixed. We find that
[TABLE]
holds from (2.54). Obviously, the conservation of the particle number leads to
[TABLE]
Thus, the integral in (3.16) is estimated as
[TABLE]
and then (3.16) becomes
[TABLE]
We emphasize that the definition of in (2.54) is essential for this relation.
Here, let us recall that the value of is determined for a given . Since there is one-to-one correspondence between and , is given as a function of . The relation (3.22) implies that is independent of in the linear response regime, and thus given as a function of . Therefore, we conclude (3.4). Next, we show (3.5). From the relation (3.22), the left-hand side of (3.5) is expressed as
[TABLE]
while the equilibrium fundamental relation results in
[TABLE]
Combining these two equalities, we obtain (3.5).
The definition (2.55) for the global chemical potential indicates . Then, (3.4) leads to
[TABLE]
We thus find
[TABLE]
where
[TABLE]
are entropy per one particle and specific volume of the heat conduction system.
3.2 Various global thermodynamic functions
From a local relation , we have
[TABLE]
By using (3.4), we obtain
[TABLE]
We thus have the fundamental relation for
[TABLE]
in the linear response regime.
Next, we consider the global internal energy and the global enthalpy . From local thermodynamic relations and , we write
[TABLE]
Remembering that , the integral is estimated as
[TABLE]
where we have used (3.19). Thus, we obtain
[TABLE]
Substituting (3.4),(3.5),(3.29) into (3.34) and (3.35), we have
[TABLE]
3.3 Clausius equality
Since the thermodynamic relations are extended with keeping the same forms as those in equilibrium thermodynamics, we may define quasi-static heat in an infinitely small quasi-static process as
[TABLE]
This corresponds to the absorbed heat by the system, which is the sum of the heat from the left and right heat baths during the infinitesimal change. There, by summing the two heats, the net heat flow is canceled. If the heat absorbed from the right heat bath is exactly the same as the heat released from the left heat bath at every moment during the process, the net heat vanishes for the heat conduction. Such systems are thought to be adiabatic in the sense of global thermodynamics, and
[TABLE]
in such a quasi-static adiabatic process in the linear response regime.
Then, constant-pressure heat capacity is defined as
[TABLE]
Applying the thermodynamic relation, we obtain
[TABLE]
Note that is defined as the response to the change of the global temperature . By changing the temperatures of heat baths, may change in accordance with absorbing heat corresponding to .
3.4 Correspondence of global thermodynamic quantities to
equilibrium quantities
In the previous subsections, we have obtained thermodynamic relations in the linear response regime and found that they are equivalent to the equilibrium ones. The key concept for extending thermodynamics is the global temperature and the global chemical potential , by which the heat conduction system of is mapped to the equilibrium system of . The connection between the two thermodynamic frameworks becomes clearer in the argument below by considering the estimation method of global thermodynamic quantities.
In a single phase system, local temperature is a continuous monotonic function of . Any extensive quantity , which is defined as a spatial integral of local density , can be transformed into an integral over temperature:
[TABLE]
where
[TABLE]
We expand around in the form
[TABLE]
The second term is canceled when the integral (3.43) is performed. We thus have an estimation
[TABLE]
For , we have for all . By combining this with (3.47), we find
[TABLE]
We then have arrived at a universal estimation of the global quantity as
[TABLE]
where . The result (3.49) directly connects the global quantity in the heat conduction state to the corresponding equilibrium quantity at .
We may apply the similar method to the global temperature. The global temperature is also written as
[TABLE]
in which
[TABLE]
By substituting it into (3.51), we obtain
[TABLE]
Thus, in the linear response regime, the global temperature in the heat conduction system is equal to the mean temperature of the heat baths. The relation (3.49), together with (3.54), concludes (3.8). We also note that (3.53), which is rewritten as with , corresponds to the global version of the equation of state
[TABLE]
We here remark that the relation in the linear response regime is a specific feature of systems in a rectangular container. We will show in §10.2 that the global temperature deviates from the mean temperature when the shape of the container is not a rectangle, whereas the relation (3.8) still holds.
4 Liquid-gas coexistence in heat conduction
From now on, we study liquid-gas coexistence in heat conduction. In §4.1, we review thermodynamics under equilibrium conditions while paying attention to the description of metastable states. In §4.2, we describe heat conduction states based on local equilibrium thermodynamics. We explain how the temperature profile and the density profile are determined for a given position of the liquid-gas interface. In §4.3, we present global thermodynamic quantities as a function of the interface position. In §4.4, we address the main problem of thermodynamics for the heat conduction systems with the liquid-gas coexistence. Hereafter, superscripts and are attached to quantities related to liquid and gas, respectively.
4.1 Thermodynamics under equilibrium conditions
We first consider equilibrium thermodynamics. For a given , the pressure is uniquely determined even when phase coexistence is observed. This equation of state is expressed as
[TABLE]
so as to emphasize that the pressure is the equilibrium value. See Fig. 4(a). A remarkable fact is that there is a plateau region for a given temperature below the critical temperature. From this equation of state, we find that the density exhibits the discontinuous change at when the temperature is changed with the pressure fixed, as shown in Fig. 4(a). That is, the whole system is occupied by liquid when , while by gas when .
In experiments, metastable states, which also may be called quasi-equilibrium states, are often observed. The most typical approach to the description of such meta-stable states is to employ the van der Waals equation of state:
[TABLE]
where is the Boltzmann constant. The constants and are van der Waals parameters. When is fixed as that below the critical temperature, the equation of state contains the unstable states . See Fig. 4(b). Following the standard procedure of the Maxwell construction, we obtain the equilibrium equation of state from (4.2). Here, we extract the two stable branches satisfying from the curve defined by (4.2), each of which is connected to the equilibrium liquid phase or to the equilibrium gas phase. By solving (4.2) in for each region, we obtain
[TABLE]
When , the spatially homogeneous phase of corresponds to a meta-stable gas. Similarly, when , the homogeneous phase of corresponds to a meta-stable liquid. By plotting as a function of with fixed as shown in Fig. 4(b), one finds that the meta-stable gas and the meta-stable liquid are super-cooled gas and super-heated liquid, respectively.
More generally, without using the van der Waals equation of state, we may obtain
[TABLE]
by experimental measurements. See Fig. 4(c). Below we assume (4.5) and (4.6), whereas one may interpret them as (4.3) and (4.4). We then define densities for the liquid phase and the gas phase as
[TABLE]
and define quantities per one particle for the liquid phase and the gas phase as
[TABLE]
4.2 Heat conduction states
When , we may observe liquid-gas coexistence at the constant pressure , which is in contrast with equilibrium cases. Concretely, the liquid occupies a lower temperature region, while the gas does a higher temperature region, as shown in Fig. 5(a). We assume that both the interface width and the temperature gap at the interface are negligible. See Fig. 5(b). Then, let denote the position of the interface. Although the value of should be uniquely determined under given conditions, we shall formally treat as if it were an independent variable in what follows. The density profile for a given is then expressed as
[TABLE]
See Fig. 5(b). The temperature profile is determined by the heat conduction equation
[TABLE]
where is the steady heat current which is constant in due to the conservation of energy. The boundary conditions for the equation (4.14) are
[TABLE]
Concretely, is obtained as follows. We first express the interface temperature as :
[TABLE]
By setting
[TABLE]
we find that the local heat conductivity is estimated as
[TABLE]
Then, (4.14) leads to
[TABLE]
Solving the second equality of (4.19) in , we obtain
[TABLE]
Defining the scaled position of the interface and an effective heat conductivity as
[TABLE]
we have
[TABLE]
which is obtained by substituting (4.20) into (4.19).
By using given by (4.20), the temperature profile is written as
[TABLE]
Substituting this profile into the equation of state (4.11), we have explicitly as exemplified in Fig. 5(c).
4.3 Global quantities as a function of
The system with the phase coexistence is interpreted as a composite system of a liquid region and a gas region. For each subsystem, we may define global thermodynamic quantities introduced in §3. In this section, we express global thermodynamic quantities as a function of .
First, the global temperatures in the liquid region and in the gas region are defined by
[TABLE]
In each region, the temperature profile and the density profile are smooth. We describe each subsystem as a single phase system discussed in §3. The boundary temperatures of the liquid are and , and those of the gas are and . Thus, applying the trapezoidal rule explained in §3.4 to (4.26), the global temperature in each region is expressed as
[TABLE]
where we explicitly write the dependence on the interface position . See Fig. 5(b).
Next, the volume of the whole system is expressed as
[TABLE]
and thus the volume of the liquid region and the volume of the gas region are written as
[TABLE]
Then, the particle number in the liquid region, , is obtained as
[TABLE]
The particle number in the gas region, , is immediately given by
[TABLE]
By using , , and , we can express all global thermodynamic quantities in the liquid region and the gas region. Explicitly, the extensive quantities are defined as
[TABLE]
Noting (3.8) in §3 with (2.24), we conclude that the global extensive quantities are written as
[TABLE]
Finally, we define the global temperature for the whole system by the formula (2.54), which is rewritten as
[TABLE]
4.4 Problem
We have shown that all thermodynamic quantities are determined for a given position of the liquid-gas interface. Since the liquid-gas interface in heat conduction may be at rest at constant pressure, the position is uniquely determined for .
For equilibrium cases, the phase coexistence is observed in a container with and fixed. When we start with the van der Waals equation of state, we can determine the interface position of the two phases by the Maxwell construction, which is equivalent to the continuity of the chemical potential at the interface. Recalling this theory, one may impose the continuity of the local chemical potential at the interface even in heat conduction Bedeauz86 . Since , it leads to . That is, under this assumption, the interface temperature is equal to the equilibrium transition temperature.
Here, the important observation is that there is no justification of this condition for heat conduction systems. To our best knowledge, there are no experimental measurements on this issue, no numerical simulations of sufficiently large systems, and no reliable theory for supporting this condition. To be more precise, although the continuity of the chemical potential was reported in numerical simulations Bedeaux00 ; Ogushi , the system size was too small to draw a definite conclusion. One may also recall from a viewpoint of non-equilibrium statistical mechanics that the local equilibrium distribution may be the leading contribution Bertini-rev . However, it should be noted that in the standard approach, fluctuations are described as those around the most probable profile, while the most probable profile itself is undetermined in the current problem. Therefore, the previous studies starting with the local equilibrium distribution do not result in the continuity of the chemical potential without imposing an additional assumption.
Furthermore, as another theoretical approach, one may study the stationary solution of the deterministic equation for the generalized Navier-Stokes equation with the interface thermodynamics Bedeaux03 ; Onuki . We can estimate the discontinuous jump of the chemical potential at the interface as based on reasonable assumptions SNIN , where represents the interface width. This result indicates the validity of the continuity of the chemical potential at the interface, because we may assume . However, as carefully argued in the paper SNIN , fluctuation effects should be taken into account for the generalized Navier-Stokes equation, so as to quantitatively describe the thermodynamic behavior.
Now, based on these facts, we reconsider the continuity of the chemical potential at the liquid-gas interface within the framework of equilibrium thermodynamics. We then find that the continuity is equivalent to the thermodynamic variational principle. If the principle were applied to a local sub-system including the interface, the continuity of the chemical potential would be concluded. However, the constraint of a variational principle is generally applied to a whole system but not to the local subsystem, and therefore no variational principle is expected for the local subsystem near the interface. Thus we cannot justify the continuity of the chemical potential. This suggests an alternative approach to determine the interface position . Namely, we formulate a variational principle for the whole heat conduction system as a natural extension of the equilibrium variational principle.
5 Variational principle for determining the
liquid-gas interface position
All our discussion on global thermodynamics so far has been firmly based on well-established local equilibrium thermodynamics. Now, we formulate the variational principle for determining the interface position and solve the variational problem. In §5.1, we start with the simplest example of the variational principle for equilibrium systems and naturally extend it to that for heat conduction systems. In §5.2 we derive the steady-state value by solving the minimization of the variational function. The solution is then reformulated as the temperature relation in §5.3. By using these results, in §5.4, we express the global quantities in terms of .
5.1 Variational principle
5.1.1 Equilibrium systems
We start with an example of the variational principle for equilibrium systems. We consider a system at the constant pressure in contact with a heat bath of the temperature . Since the volume of the system is not fixed, we determine for a given . When the equation of state is assumed, it is determined from the pressure balance equation
[TABLE]
We attempt to derive this condition from a variational principle. We take a variational function as
[TABLE]
This corresponds to the Helmholtz free energy minimum principle for the composite system of the system and the volume reservoir whose Helmholtz free energy is . We then find that the minimization of with respect to for fixed leads to (5.1). The volume in the equilibrium state satisfies
[TABLE]
with which the Gibbs free energy is determined as
[TABLE]
It should be noted that this principle may be applied to the case that is given by the van der Waals free energy:
[TABLE]
The variational principle for this case is equivalent to the Maxwell construction.
5.1.2 Extension to heat conduction systems
We study heat conduction systems at the constant pressure . First, we focus on single-phase systems studied in §3. From the equivalence between a heat conduction system and the corresponding equilibrium system, the volume of this system is determined by the minimization of
[TABLE]
where is the global temperature of the system. Since , we rewrite (5.7) as
[TABLE]
where depends on and should be fixed in the variation of . Then, the variational principle (5.3) holds for single-phase systems. Since this formula relies on the equivalence to the equilibrium system, we do not need to fix the degree of nonequilibrium in the variation.
Next, we consider a system with a liquid-gas interface. In this case, is not determined when we use the equation of state for each region, because the interface position is not determined yet as described in the previous section. We then assume that the minimization of (5.8) is also valid for such cases, where is also fixed in the variation . The last property is necessary to determine the value uniquely. That is, we define a variational function
[TABLE]
and propose that the volume of the system is determined as
[TABLE]
with
[TABLE]
This variational principle was first proposed in Ref. NS .
5.1.3 Remarks
Since the volume is uniquely determined by the interface position with fixed, as explicitly shown in (4.28), the variational principle (5.10) with (5.11) is expressed as
[TABLE]
Similarly, since the particle number in the liquid region is uniquely determined by , we also have another variational principle
[TABLE]
As the second remark, we provide a physical interpretation of (5.12). Let us consider a fluctuation of the interface position, . We assume that the motion of the interface is slowest and that the interface motion can be observed in a hypothetical system where is controlled such that and are fixed. See Fig. 6. It is then expected that fluctuations in this hypothetical system may be described by equilibrium statistical mechanics for the system with . Thus, the most probable position is given by the left equation of (5.12) and the stability of the interface position is expressed as the right inequality in (5.12).
Moreover, we interpret from an operational viewpoint. We virtually attach another heat bath of the temperature to the rigid top plate as shown in Fig. 7. Using the idealized assumption that the motion of the top plate is sufficiently slower than the other dynamical degrees of freedom, we control and with fixed such that the heat flux to the heat bath of is zero. Here, noting that the total kinetic energy in the bulk is given by with the total degrees of freedom in the system, is proportional to in the linear response regime. That is, is equivalent to . When we focus on the position of the top plate, the motion would be described as if it were in the equilibrium state at the temperature . From this picture for the specific setting, it is expected that (5.12) is equivalent to the equilibrium variational principle with .
As the last remark, we comment on the condition that is fixed. Since there are two temperatures and , two variables associated with should be fixed in the variation of . Although we adopt as a fixed variable in addition to the global temperature , one may conjecture that would be a more plausible fixed-variable than . As is formulated in (4.22), is proportional to but the proportional constant depends on the interface position. Thus, the variational principle with fixed results in a different steady state from the solution of (5.10) with (5.11). There are two reasons for our choice. The first is that the final result becomes simplest among our trials. Second, when is fixed in the variation, the enthalpy is conserved in the variation. In this case, the corresponding variational principle may be an extension of the equilibrium variational principle for adiabatic (thermally isolated) systems. From these two aspects, we assume that is fixed in the variation.
5.2 Steady state determined from the variational principle
We solve the variational equation (5.13). Hereafter, we consider the variation of with fixed. We abbreviate as or for the notational simplicity.
We express the variational function as
[TABLE]
where the free energy of the liquid , the free energy of the gas , and the volume of the system are given as functions of with fixed. Here, for any global quantity in each region, such as and , we define
[TABLE]
We then have
[TABLE]
From the argument in §3, we obtain
[TABLE]
Since every global thermodynamic relation obtained in §3 holds in each region, we have
[TABLE]
with
[TABLE]
Since the variation of the free energy in each region obeys the fundamental relation of thermodynamics for , we have
[TABLE]
where
[TABLE]
By using (5.20), we rewrite the variation of in (5.14) as
[TABLE]
where we used and .
The global temperature in the liquid region and in the gas region satisfy the relation
[TABLE]
which are derived from the relations
[TABLE]
and
[TABLE]
Furthermore, (5.23) leads to
[TABLE]
We have used in the first line. Especially, the formulas (5.26) and (5.27), with , , and fixed, bring
[TABLE]
which simplifies (5.22) as
[TABLE]
Here, we confirm
[TABLE]
and
[TABLE]
Let us define a specific temperature as
[TABLE]
By using (5.23), we rewrite (5.32) as
[TABLE]
Thus, the variation (5.29) is further simplified as
[TABLE]
This is equivalent to
[TABLE]
Since the functional form of is different from that of due to the crucial difference between the liquid and the gas, is not identically equal to . The equality
[TABLE]
holds only when
[TABLE]
Now, (5.32) and (5.37) yield the unique value as
[TABLE]
Thus we conclude that the variational principle (5.13) results in the unique steady value formulated by (5.38).
Next, we consider the second derivative of . By using (5.35) with (5.32), the second derivative of is obtained as
[TABLE]
At the steady state satisfying (5.37), the above entropy difference is connected to the latent heat per one particle
[TABLE]
Therefore, we estimate the second derivative as
[TABLE]
Since and , we conclude that
[TABLE]
in the linear response regime.
5.2.1 Careful analysis of
The solution (5.38) of the variational principle includes an undetermined term of . Nevertheless, (5.38) provides additional information to equilibrium behavior. To clarify the situation, we consider the limit with keeping the phase coexistence. This may be formalized by fixing a parameter defined as
[TABLE]
whereas (5.38) indicates
[TABLE]
As , satisfies
[TABLE]
In Fig. 8, we draw straight lines connecting from to for several values of in the parameter region where the liquid-gas coexistence is observed. The convergence of these lines at indicates that the equilibrium state behaves as a singular point. By specifying the value of , we identify the corresponding line which undoes the degeneration at equilibrium. is determined by considering the equilibrium limit , while it is not uniquely determined at equilibrium.
5.3 Temperature relation
We obtained the steady state value of as (5.38) by solving the variational principle (5.13). Substituting this solution into (5.33) with (5.37), we have the global temperature for each region as
[TABLE]
We sum up the two relations in (5.46) and substitute (4.27) into . We then obtain
[TABLE]
which we call the temperature relation. This non-trivial relation, which results from the variational principle, is a simple condition that the steady states satisfy. Generally, is not equal to , because the particle density of liquid is larger than that of gas, . Therefore, the relation (5.47) implies that the liquid-gas coexistence temperature in a heat conduction system must be strictly lower than the equilibrium coexistence temperature . This means that metastable states stably appear near the liquid-gas interface as schematically shown in Fig. 9.
Now, suppose that (5.47) holds. Summing up the two equalities in (5.23), we obtain
[TABLE]
Here, it should be noted that we do not use the variational principle for the derivation of (5.23). By substituting the temperature relation (5.47) into this result, we derive (5.38) without the variational principle. In this sense, we may say that the temperature relation (5.47) is equivalent to the variational principle.
5.4 Global quantities as functions of
We determine the steady state values of several quantities. Let denote the equilibrium saturated volume defined as
[TABLE]
with
[TABLE]
The volume of the liquid (gas) region is given by
[TABLE]
Since is determined with an error of in (5.38), (5.51) is the best estimate of in . Substituting (5.38) into (5.51), we obtain
[TABLE]
which corresponds to the solution of the variational principle (5.10). Then, the steady position of the liquid-gas interface,
[TABLE]
is determined as
[TABLE]
which is the solution of the variational principle (5.12).
Then, the interface temperature is obtained from the elimination of from (5.47) and (4.20) as
[TABLE]
The deviation of from and the relation between and will be demonstrated in §6 for explicit examples.
6 Properties of the steady states with a
liquid-gas interface
In §5.4, we derived the steady state values for given such as in (5.38), in (5.52) and in (5.54). However, the global temperature is not easily controlled in experiments, differently from the heat bath temperatures and . In this section, we express any quantity in the steady state as a function of . In §6.1, we derive as a function of . Since we already have
[TABLE]
the relations
[TABLE]
lead to the expression
[TABLE]
Particularly, we consider the interface temperature and the jump of the chemical potential at the interface as a function of in §6.2 as (6.12) and §6.3 as (6.27), respectively. In §6.4, we discuss how the super-cooled gas appears near the interface quantitatively. The last subsection §6.5 is devoted to the demonstration of examples. In what follows, we use the notations
[TABLE]
for simplicity.
6.1 and as functions of
Substituting (5.55) into (5.47), we obtain
[TABLE]
where the dimensionless position of the steady interface is
[TABLE]
as is derived in (5.54). Eliminating from (6.6) and (6.7), we have
[TABLE]
Solving this equation in , we derive
[TABLE]
with
[TABLE]
(6.9) provides a concrete functional form of (6.3).
Substituting (6.8) into (6.6), we obtain a concrete form of (6.2),
[TABLE]
with the form of in (6.9).
6.2 Interface temperature as a function of
Applying the temperature relation (5.47) to (6.11), we find that the interface temperature is written as
[TABLE]
where is given by (6.9). This formula clarifies that the interface temperature generally deviates from when the temperature gradient is imposed. Since the particle density satisfies , i.e., , and the heat conductivity is expected to be , the interface temperature is deviated from the equilibrium transition temperature in the order of , which is not negligible.
We have another formula for the interface temperature , which is expressed solely by experimentally accessible quantities:
[TABLE]
where . The derivation of this formula is shown in §6.2.1.
Up to here, we have studied the systems with , where and . More generally, from the left-right symmetry of the system in Fig. 2, we notice that the interface temperature is invariant for the transformation to . That is, (6.13) is expressed as
[TABLE]
for any .
One might guess that the coexistence temperature is uniquely determined by local quantities that characterize the phase boundary, namely, the pressure and the heat current per unit area. Rather interestingly, this is not the case. Reflecting the global nature of our variational principle, the coexistence temperature explicitly depends on global conditions of the system, namely, the temperatures of the two heat baths.
6.2.1 Derivetion of (6.13)
We first recall (4.19) and (4.20). From these,
[TABLE]
Combining these two relations, we have another form of as
[TABLE]
with an error of .
Second, the result of the variational principle (5.38) is written as
[TABLE]
As and with and , the relation (6.18) is further expressed as
[TABLE]
with an error of . By subtracting (6.17) from (6.19) and using the temperature relation (5.47), we have two forms of the interface temperature:
[TABLE]
By substituting (6.20) and (6.21) into the first term and the second term of the right-hand side of the trivial identity
[TABLE]
we obtain (6.13).
6.3 as a function of
Next, we discuss the chemical potential jump
[TABLE]
at the interface. Note that
[TABLE]
which gives the definition of the transition temperature . Thus, means the imbalance of the chemical potential at the interface, which is quantitatively expressed as
[TABLE]
Since the entropy difference is connected to the latent heat per one particle as
[TABLE]
the jump of the chemical potential is proportional to . Now, substituting (6.12) into (6.25), the amount of the jump is estimated as
[TABLE]
Another expression of is also obtained by substituting (6.13) into (6.25). The result is
[TABLE]
6.4
Super-cooled gas in the liquid-gas coexistence
When , we find the position at which the local temperature satisfies . Suppose . Then, the position of the interface satisfies
[TABLE]
We observe the liquid in the region and the gas in . Then the local temperature of the gas in is less than . This means that a super-cooled gas is observed in , which is not stable in equilibrium.
We plot (6.9) and (6.12) in Fig. 10 for three sets of . It is clearly seen that the interface temperature deviates from the equilibrium transition temperature . The deviation of the global temperature from shows the same figure with Fig. 10(b) due to the temperature relation (5.47). The jump of the chemical potential at the interface also exhibits the same dependence on since it is proportional to as shown in (6.27). In the presented three sets of , a super-cooled gas region appears near the right side of the interface, as shown schematically in Fig. 9.
According to Fig. 10(a), the stable position of the interface is shifted continuously as we increase the temperature of both heat baths. This feature is different from the numerical solution of the variational equation reported in Ref. NS , where the interface discontinuously appears. A more accurate numerical calculation was necessary for the cases with high .
6.5 Examples
In this subsection, we illustrate some examples of the liquid-gas coexistence in heat conduction quantitatively, where the particle density is measured in mol.
First, we consider the van der Waals equation of state
[TABLE]
of Pa m6/mol2 and m3/mol for at a high pressure Pa CO2 . The gas constant is J/Kmol. Substituting these parameters into (6.30), we obtain the transition temperature as K, and the mol density at as mol/m3 and mol/m3, which results in . From the database NIST , we set the heat conductivity of the liquid as W/m K and that of the gas as W/m K, which yields . The result is close to the result which was already shown in Fig. 10. The temperature profile becomes linear both in the liquid and gas regions in this example due to the constant heat conductivity. Then the volume fraction of the super-cooled gas is simply obtained as
[TABLE]
where is the volume of the super-cooled gas, i.e., . The interface temperature is given by (6.12) and results from (6.8) as
[TABLE]
where . In Fig. 11, we show the interface temperature and the volume fraction as a function of the position of the interface when K. We notice that the interface temperature may deviate about K from , and that the volume fraction of the super-cooled gas may exceed .
Second, more importantly, we do not need the equation of state, but have only to know the value of , , , when we obtain the phase diagram for a given material. For example, consider pure water at Pa, where K. From the database NIST , we have and . From these, we can predict the interface temperature and the volume fraction as a function of the position of the interface as shown in Fig. 12. As an example, when K and K, we obtain K and .
7 Global thermodynamics for systems with a liquid-gas interface
Hereafter, all the quantities are defined in the steady state. When a global quantity is determined for , we express this relation as
[TABLE]
Recall that for single-phase systems, where is the equilibrium function, as discussed in §3. For the systems with a liquid-gas interface, dependence appears even in the linear response regime. In this section, we extend the global thermodynamics introduced in §3 so as to describe the liquid-gas coexistence determined in the previous sections. In §7.1, we derive formulas as preliminaries for later sections. In §7.2, we formally derive a fundamental relation of thermodynamics for systems with a liquid-gas interface. However, since the derivative of is not defined at , the formal expression is not properly defined. Then, in §7.3, we perform a careful analysis of the limit . In §7.4, we show the form of the entropy and the volume in the appropriate limit. In §7.5, we present formulas for constant pressure heat capacity and compressibility.
7.1 Preliminaries
We first study global extensive quantities
[TABLE]
It is obvious that they possess the additivity
[TABLE]
where and are defined as global thermodynamic quantities in the liquid region and the gas region, respectively. By using and , we write
[TABLE]
Now, we consider an infinitely small change , which results in
[TABLE]
Accordingly, any quantity for the steady state changes as
[TABLE]
Then, (7.4) and (7.5) bring the following relations among the infinitely small changes:
[TABLE]
We sum up the above two relations (7.8) and (7.9) and substitute (5.26) and (5.27) into them. Then, the sum becomes
[TABLE]
The fourth line is transformed into
[TABLE]
By summarizing the terms proportional to in (7.10) and (7.11), the coefficient of becomes
[TABLE]
Here, we note that the sum of the two relations for and in (5.23) results in
[TABLE]
where the transformation from the first line to the second line is found by noting and from (3.54). Now, using the temperature relation (5.47), we find that the right-hand side of (7.13) is . We thus obtain
[TABLE]
Here, the last term should be expressed as a linear combination of , , and . Indeed. from (5.38) which gives as a function of , we derive
[TABLE]
We leave the formula (7.14) and (7.15) as they are for later convenience.
7.2 Formal derivation of fundamental relation
We set
[TABLE]
in the formulas in the previous subsection. Then, we recall the thermodynamic relations
[TABLE]
where is the specific volume. We also define
[TABLE]
where is the latent heat defined in (5.40). By using them, we rewrite (7.14) as
[TABLE]
Here, is the global chemical potential defined by (2.55). It should be noted that
[TABLE]
where is the global chemical potential in each region such that and .
Now, noting the relations
[TABLE]
and (7.15), we obtain
[TABLE]
By formally considering the infinitely small change , we obtain
[TABLE]
which corresponds to the fundamental relation of thermodynamics.
At first sight, one may be afraid that the relation (7.23) does not provide a non-equilibrium extension because the relation involves the error of . However, it includes non-trivial information. The point is that some equilibrium thermodynamic quantities are singular at . For example, and are not uniquely defined at in equilibrium. Nevertheless, (7.23) provides a definition of at as the limit of the derivative of with respect to while fixing . We study (7.23) more carefully by going back to (7.22).
7.3 Careful analysis of
We consider the limit with keeping the phase coexistence. The parameter defined in (5.43) identifies the line terminating to the equilibrium state as exemplified in Fig. 8, and undoes the degeneration of the equilibrium state. Since (5.43) is rewritten as
[TABLE]
a global quantity is considered as a function of through
[TABLE]
Then, let us define the equilibrium limit of entropy and volume by
[TABLE]
where we explicitly write only the dependence for and .
Substituting (7.24) into (7.22), we obtain
[TABLE]
where we have used (5.44). This relation leads to
[TABLE]
Since we find from (7.24) that the left-hand side is equal to
[TABLE]
we obtain
[TABLE]
Similarly, (7.28) leads to
[TABLE]
Here, from
[TABLE]
with (7.24), we have the identity
[TABLE]
Substituting this into (7.32), we obtain
[TABLE]
Combining (7.31) and (7.35), we may write
[TABLE]
which is the more precise expression of (7.23) and convey non-trivial information of thermodynamics.
7.4 Explicit forms of and
By using (5.44), we express any extensive quantity as
[TABLE]
Letting , we have
[TABLE]
This leads to
[TABLE]
where we have used the formula (5.40) for the latent heat . Similarly, we obtain
[TABLE]
where we have used the Clausius-Clapeyron relation
[TABLE]
Note that in (7.40) is consistent with in (5.52) at the limit .
7.5 Heat capacity and compressibility
Let be the enthalpy for the heat conduction system with the phase coexistence. Since is given by the spatial integral of the local enthalpy, it satisfies
[TABLE]
This leads to
[TABLE]
in an infinitely small quasi-static process at constant pressure. We interpret this as the quasi-static heat in this process. We thus define heat capacity at constant pressure as
[TABLE]
Substituting into (7.14) with fixed, we have
[TABLE]
where the specific heat for the liquid (gas) is defined as
[TABLE]
By using (5.38) and (7.44), we rewrite (7.45) as
[TABLE]
with the latent heat . Note that the third term mainly contributes to in (7.47) because it diverges as . This is consistent with the singularity at for equilibrium cases, where is given by the derivative of a discontinuous function (enthalpy) at . The first and the second terms correspond to the heat capacity of the liquid region and the gas region, respectively. Thus, the third term may be interpreted as the interface contribution. The expression (7.47) also clarifies the violation of additivity for .
Similarly to the heat capacity, we can also study the additivity of other response functions. As one example, we show the singular nature of the compressibility from the viewpoint of the violation of the additivity. The compressibility in heat conduction is defined as
[TABLE]
By taking in the formula (7.14), we obtain
[TABLE]
Similarly to , the main contribution in is the non-additive term, which diverges in the equilibrium limit.
8 Thermodynamic relations in
The fundamental relation (7.23) derived in the previous section is expressed in terms of the state variable . As is familiar with thermodynamics, one may consider other fundamental relations. The most familiar one to physicists may be the form using the state variable , because this form is directly related to experimental configurations with fixed, instead of fixed in the previous sections. Since the same steady state can be realized when the fixing condition is replaced from to , there is one-to-one correspondence between and . See Fig. 13. In this section, we show thermodynamic relations with volume . In §8.1, we derive the fundamental relation in . In §8.2, we express the free energy by using the saturated pressure for equilibrium systems. Since the value of is not easily controlled in experiments, in §8.3, we formulate and as functions of and . Below, we fix without loss of generality and sometimes abbreviate as .
8.1 Fundamental relation
Since and are given as the spatial integral of the local free energies, the uniformity of the pressure leads to
[TABLE]
Substituting this into (7.22), we have
[TABLE]
where the error of in (7.22) is regarded as which results in the error of besides and so on. The most important thing here is the order estimate of , which has been estimated as when is a function of . See (7.15). As a function of , becomes of because the liquid-gas coexistence stably appears at equilibrium () in a certain range of and, therefore, the term is not singular in the limit . Indeed, the second line is estimated as
[TABLE]
where we have used
[TABLE]
Thus, the last term in (8.2) remains to be negligible as term, and we conclude
[TABLE]
It should be noted that there are the errors of in contrast to (7.23).
8.2 Free energy for the coexistence phase
In the equilibrium coexistence phase, as shown in Fig. 14, the saturated volumes, and , and the saturated pressure are defined for a given temperature . Even for the coexistence phase in heat conduction, we define the saturated volumes, and , at which the liquid and the gas start to coexist. Note that the pressure is not kept constant in , differently from the equilibrium case.
We first express , and in terms of , and . Since the steady state satisfies (5.38), we have
[TABLE]
Solving this in , we obtain
[TABLE]
For equilibrium cases, the relations and lead to
[TABLE]
By using these expressions, we rewrite (8.7) as
[TABLE]
where . Figure 14 shows a - curve described by (8.9), in which is linearly decreasing with in the range . Configurations in the phase coexistence are exemplified in Fig. 15. Then, the saturated volumes under heat conduction are expressed as
[TABLE]
where and are the equilibrium equation of state for the liquid and the gas, respectively.
Next, we study the free energy . We first consider the free energy difference
[TABLE]
Since the saturated state is in a single phase, we have
[TABLE]
We then transform
[TABLE]
where we have used the Taylor expansion in to obtain the second line, and substituted the equilibrium relation
[TABLE]
to obtain the last line.
Second, from (8.5), we have
[TABLE]
Performing the integration with (8.9), we obtain
[TABLE]
with
[TABLE]
Here, by using (8.8), we can confirm
[TABLE]
Combining it with the Clausius-Clapeyron relation
[TABLE]
we derive
[TABLE]
where was introduced in (7.18). Summarizing these, we schematically draw the graph in Fig. 16. For equilibrium cases, the free energy of the liquid-gas coexistence phase is expressed as a common tangent at and as shown in Fig. 16(a). For heat conduction cases, the free energy of the coexistence phase is expressed as a common tangential quadratic curve at and as shown in Fig 16(b), and therefore keeps the convexity on .
Since the convexity of is concluded, the Gibbs free energy is expressed as the Legendre transform of :
[TABLE]
The graph in Fig. 16(b) yields the graph in Fig. 17(b), while it contains an error of as well as formulated explicitly in (8.16). The liquid-gas coexistence states are expressed as a common tangential quadratic curve to the graph at and .
8.3 , and as functions of
We have studied the dependence of the steady states on for fixed. However, when we fix values of and , the change of yields the change of . Here we propose a protocol to control and , by which is fixed when changing . In other words, we derive and as functions of . We also give the interface temperature and its deviation from as functions of .
In (6.8), is expressed in terms of , , , , and , where the value of and are determined for each . For the coexistence phase in heat conduction, the pressure is determined as (8.9) and can be replaced by for the argument of and . We thus represent and , and simply write and . Then, we rewrite (6.8) as
[TABLE]
with
[TABLE]
Furthermore, from
[TABLE]
and at equilibrium, we obtain
[TABLE]
These expressions of , and lead to
[TABLE]
with which (8.22) is rewritten as
[TABLE]
Here, let us recall that defined in (5.43) and that
[TABLE]
as used in (8.9). Substituting this formula into (8.29) and noting
[TABLE]
we finally arrive at
[TABLE]
Since and , we obtain
[TABLE]
These formulas inform how to operate and in order to fix the global temperature when changing and .
Substituting the temperature relation (5.47) into (8.29), we obtain the interface temperature as
[TABLE]
Similarly, (8.32) and (5.47) lead to the deviation from the transition temperature,
[TABLE]
9 Other assumptions
leading to the results of the variational principle
To this point, it has been assumed that the steady states are determined by the variational principle proposed in §5. In this section, we introduce other assumptions from which the result of the variational principle is obtained. First, in §9.1, we formulate the variational principle for the constant volume systems and show that its solution is equivalent to the previous one (5.38). In §9.2, we start with the fundamental relation of thermodynamics without assuming the variational principle. We then derive the result of the variational principle. In §9.3, we notice a singularity relation which is a simple assumption for the singularity in the limit . We confirm that this is equivalent to the result of the variational principle. Finally, in §9.4 we argue some scaling behavior of the system which also leads to the result of the variational principle.
9.1 Variational principle for constant volume systems
In this subsection, we study constant volume systems, where the steady state is characterized by . As shown in Fig. 13, this steady state is equivalent to the steady state of the system at the constant pressure when the value of is chosen as the pressure . We assume the mechanical balance everywhere. That is,
[TABLE]
for any and in , where . This condition determines and when and are determined. The problem is then to determine . After briefly reviewing the corresponding variational principle for equilibrium cases, we propose the variational principle for heat conduction systems and determine the value of .
9.1.1 Equilibrium systems
We define the variational function as
[TABLE]
in which is the Helmholtz free energy of the liquid (gas). and are determined by
[TABLE]
and . Then, the variational principle for determining is
[TABLE]
which is equivalent to the equality of the chemical potential.
9.1.2 Heat conduction systems
Let us consider a variational principle for the heat conduction systems characterized by . We assume that the steady state with a liquid-gas interface is determined by the variational function
[TABLE]
where we have used a similar formula as (7.3) with (7.4) and (7.5) to obtain the second line. As we formulated in §4, thermodynamic quantities are expressed as a function of such that . The variation of is then written as
[TABLE]
Remembering that the relation (5.20) holds in each region, we have
[TABLE]
Substituting , and (5.28) into (9.7), we rewrite (9.6) as
[TABLE]
where is given in (5.32). Then, the variational principle
[TABLE]
determines as that satisfying
[TABLE]
This is equivalent to (5.36) so that (9.10) leads to (5.38). Therefore, the variational principle for the constant volume system is equivalent to that for the system at constant pressure.
9.2 Thermodynamic relation
In this subsection, we do not assume any variational principles, but assume the fundamental relation of thermodynamics
[TABLE]
for the liquid-gas coexistence phase in the linear response regime.
First of all, it should be noted that (7.10), (7.11), (7.12) and (7.13) hold regardless of the variational principle. Then, we have
[TABLE]
instead of (7.14). Letting and , and using , we obtain
[TABLE]
where
[TABLE]
Since , the assumption (9.11) leads to
[TABLE]
This means the temperature relation (5.47). As shown in §5.3, (5.47) leads to (5.38) so that the thermodynamic relation is equivalent to the variational principle.
9.3 Singularity relation
We here consider a protocol to shift the liquid-gas interface from to with fixing and . This protocol is obtained by varying as
[TABLE]
as shown in Figs. 18(b) and (c), where is changed as
[TABLE]
On the phase diagram of for a given , the change occurs along the line displayed by the arrow in Fig. 18(a) when is sufficiently small. Therefore, it is natural to assume that
[TABLE]
for sufficiently small . This differential equation (9.18) is written as
[TABLE]
which we call a singularity relation. Solving (9.18) with the boundary conditions
[TABLE]
in Fig. 18(b) and
[TABLE]
in Fig. 18(c), we obtain
[TABLE]
which is the same form as the result of the variational principle (5.38).
The relation (9.19) indicates
[TABLE]
which corresponds to an expression of the singularity associated with the first-order transition for equilibrium cases. This singularity is consistent with the discontinuous change of at . It should be noted that (9.23) is connected to the singularity of constant pressure heat capacity and compressibility as shown in §7.5.
9.4 Scaling relation
We characterize the coexistence phase by . As examples, we write
[TABLE]
From the homogeneity in the direction perpendicular to , and are invariant for . We thus write
[TABLE]
where
[TABLE]
By noting , we express and as
[TABLE]
Now, we attempt to extend the gas region while keeping the thermodynamic state in the liquid region. See Fig. 19. First, we increase as with fixed, where we take such that the system is still in the linear response regime. We estimate the dependence of quantities in this asymptotic regime. First, the heat flux is proportional to as a common value to the liquid and the gas region. Since is fixed, it is plausible that the temperature difference in the liquid region remains to be , and then this leads that the horizontal length of the liquid region is proportional to . On the other hand, in the gas region, the temperature difference is proportional to and the horizontal length is in the leading order estimate. Thus, the total volume saturates at the finite value with which the whole system is occupied by the gas with the temperature less than . From these, we find that the volume of the liquid region is proportional to . Although its proportional coefficient is not determined in the asymptotic region, we assume the scaling relation that the volume of the liquid is as it is consistent with the case . Since and is fixed in this operation, the density of the liquid is also fixed. Thus, the particle number of the liquid becomes . Next, we increase as . Then, we have and . That is, for a series of processes
[TABLE]
we have
[TABLE]
Since this condition cannot be concluded from the local equilibrium thermodynamics, we impose (9.32) as a requirement for the steady state.
By using (9.29) for (9.32), we have
[TABLE]
for . Expanding
[TABLE]
and substituting it to (9.33), we obtain
[TABLE]
for . This gives for . Thus,
[TABLE]
which restricts the dependence of .
Next, we extend the liquid region while keeping the gas region. In this case, we consider the transformation , while fixing . From (9.30), we have
[TABLE]
Expanding
[TABLE]
in the limit , we obtain and for . This leads to
[TABLE]
By combining (9.36) and (9.39), we obtain
[TABLE]
Since this holds for any , we obtain . By substituting (9.36) with into (9.29), we have arrived at
[TABLE]
Substituting (5.23) into this relation, we obtain the solution (5.38) of the variational principle.
10 Generalization for single phase systems
In previous sections, global thermodynamics has been developed for systems in a rectangular container as shown in Fig. 2. In this section, we consider an arbitrarily shaped container in contact with two heat baths of and . Here, we assume that the response of a local quantity is not singular for the change of the system parameters , where is the volume of the container. We also assume that local steady states inside the container are well described by the local equilibrium thermodynamics of local quantities , , , and so on. In §10.1, we define global thermodynamic quantities as they are consistent with the definitions given in §3. We construct global thermodynamics in the linear response regime in §10.2. Then, we extend the thermodynamic formulas beyond the linear response regime in §10.3. In what follows, we restrict to single-phase systems without liquid-gas interfaces.
10.1 Global thermodynamic quantities in an arbitrarily
shaped container
We define all global thermodynamic quantities similarly to the case of rectangular containers. Explicitly, the global temperature and the global chemical potential are defined as
[TABLE]
and
[TABLE]
More generally, global extensive quantities are defined by the spatial integration of local thermodynamic quantities per unit volume or per one particle as
[TABLE]
It should be noted that all global quantities transform consistently for the change of the reference state on entropy and internal energy as discussed in §2.5.
Due to the mechanical balance, the local pressure is homogeneous. That is,
[TABLE]
We here note that the pressure does not satisfy the local equilibrium equation of state beyond the linear response regime. Concretely, the contribution out of the local equilibrium pressure was explicitly calculated in the kinetic regime by analyzing the Boltzmann equation Kim-Hayakawa , and also the long-range correlation yields the additional pressure depending on the system size Aoki ; Wada ; Sengers2 . Thus, we express
[TABLE]
where is the pressure determined from the local equilibrium equation of state such that
[TABLE]
with . Let us introduce a global version of the local equilibrium pressure as
[TABLE]
This may deviate from the global pressure as
[TABLE]
10.2 Equivalence of non-equilibrium global quantities with
equilibrium quantities
We show that the results in §3 hold in the container with an arbitrary shape in the linear response regime. As we did in §3.1, we introduce a function as
[TABLE]
From the definition of , we have . Any local quantity per one particle is expanded in as
[TABLE]
We then derive
[TABLE]
This indicates that all global thermodynamic quantities are equivalent to those in equilibrium by adopting the global temperature in (10.1), and that global thermodynamics for heat conduction systems are generally mapped to equilibrium thermodynamics, regardless of the shape of the container. It should be noted that we cannot apply the trapezoidal rule argued in §3.4 to such general configurations. Thus, in contrast to the case of rectangular containers, we have
[TABLE]
even in the linear response regime.
10.3 Fundamental relation beyond the linear response regime
Since global thermodynamics relies on local equilibrium thermodynamics, its validity depends on the extent how the local states are close to the local equilibrium. For instance, the local equilibrium thermodynamic relations may be valid up to and a non-equilibrium driving may bring a slight deviation of from the local equilibrium. We here expect that the deviation from the local equilibrium is small even beyond the linear response regime, and attempt to extend the global thermodynamics to the second order of . Indeed, in §10.3.1 and §10.3.2, we prove
[TABLE]
where is either , or . is the conjugate variable to , which is identified as
[TABLE]
This is a quantity of with constant pressure specific heat . is a geometrical factor independent of and as examined in §10.3.3. That is, we can extend the thermodynamic framework with an error of beyond the linear response regime.
10.3.1 Derivation of
(10.13)
Any local density per one particle is expanded as
[TABLE]
Then, we have
[TABLE]
from the definition of and . Since and , the integral of the third term in the right-hand side of (10.17) is expressed as
[TABLE]
As the integral of (10.18) vanishes, the third term is estimated as . Next, by applying to the second term of (10.17), we obtain
[TABLE]
where
[TABLE]
is a geometrical factor depending on the configuration of the container as will be examined in §10.3.3. Because , we conclude that global quantities are extended to those including contribution:
[TABLE]
With this formula, we write the global thermodynamic quantities as
[TABLE]
Performing the spatial integration of the local relation , we have the global relation
[TABLE]
For the later purpose, we define a new quantity as
[TABLE]
Using (10.24), we obtain
[TABLE]
which is the form given in (10.15).
Now, consider a parameter change , and . Accordingly, the local quantities are changed as
[TABLE]
for , and the global quantities as
[TABLE]
This leads to
[TABLE]
where we have used (10.25). From (10.24), we derive
[TABLE]
Substituting this into (10.30), we obtain
[TABLE]
where is rewritten by the linear combination of , and . Moreover, for the change with fixed, (10.25) leads to
[TABLE]
Substituting (10.31) into this, we obtain
[TABLE]
By combining (10.32) and (10.34), we have arrived at the formula (10.13).
10.3.2 Derivation of (10.14)
Since , we have
[TABLE]
Substituting the expansion
[TABLE]
and into (10.35), we obtain
[TABLE]
where the first term in the right-hand side of the first line turns out to be zero from the definition of . Thus, (10.35) is rewritten as
[TABLE]
which corresponds to (10.14).
10.3.3 Factor
determined by the geometry of containers
Hereafter, we show that is a function of the geometry of the container and independent of and . We consider a container as shown in Fig. 20, which is in contact with the first heat bath of in a region of the surface and with the second bath of in . The rest of the surface of the container is thermally insulating. The steady profiles of temperature, density, and heat flux are , , and , respectively.
The profile is determined by the condition for the steady state, which is written as
[TABLE]
with an error of . By introducing a scale transformation
[TABLE]
as shown in Fig. 20, we express (10.41) as
[TABLE]
where . Letting the solution of (10.45) be , the global temperature is written as
[TABLE]
We here expand in :
[TABLE]
With these expansions, (10.45) is written as
[TABLE]
It is obvious that the solution of (10.51) depends only on the geometry of the surface but is independent of the scale of the container and all physical parameters. More explicitly, we confirm that is independent of and . (10.46) yields
[TABLE]
which indicates is independent of and .
At last, we examine the parameter which is rewritten as
[TABLE]
Using (10.47) and (10.52), we have
[TABLE]
which is determined by the function given by (10.51). For instance, calculating for rectangular containers, we obtain
[TABLE]
11 Concluding Remarks
In this long paper, we have developed a thermodynamic framework for heat conduction states, which we call global thermodynamics. The key concept of the framework is the global temperature defined by (2.54). We describe spatially inhomogeneous systems by a set of global thermodynamic variables. Although the method is rather different from a standard continuum description, global thermodynamics is built on local equilibrium thermodynamics and the transportation equation. Thus, it is not contradictory with established theories. Nevertheless, our framework provides a new prediction on the thermodynamic properties of a liquid-gas interface in heat conduction. As we explained in §4.4, the rule of connecting two phases at the interface is out of the standard hydrodynamics with local equilibrium thermodynamics. Even for this case, global thermodynamics can determine the way how to connect the two phases by the variational principle that is a natural extension of the Maxwell construction for equilibrium cases.
Although global thermodynamics is self-consistent, self-contained, and natural, it does not ensure the validity of the framework. Similarly to other universal theories such as thermodynamics, hydrodynamics, and statistical mechanics, global thermodynamics involves a fundamental assumption on the determination of the steady state. In other words, if our quantitative prediction on the phase coexistence is denied by experiments, global thermodynamics could not be valid. In this case, we will attempt to understand the experimental results carefully so that we have a correct description of phenomena out of equilibrium. This is of course quite stimulating, and we sincerely wait for experiments. Numerical simulations of microscopic dynamics such as molecular dynamics also provide useful materials for further studies. To the present, no deviation of the interface temperature from the transition temperature was observed Bedeaux00 ; Ogushi . We conjecture that this is due to the insufficiency of the separation of scales, because the deviation should be observed in systems with enough separation of scales, as discussed in Ref. SNIN .
Putting aside the ultimate validity of global thermodynamics, we emphasize that the framework itself is quite fruitful. This paper presents only the minimum backbone, but exhibits non-trivial simple relations among several quantities. These somewhat miraculous statements may suggest a deeper structure behind global thermodynamics. Moreover, importantly, the story does not end. We did not discuss the Clausius equality for the cases with a liquid-gas interface, while we confirmed it for single-phase systems in §3.3. The equality connects “heat” to the entropy change, and the “heat” should be defined as the renormalized one, called “excess heat”, after subtracting the persistent heat as we reviewed in Introduction. It is quite natural to study the liquid-gas coexistence by considering an extended form of the Clausius relation. Since we find that this provides another source of further development of the framework, we have decided to argue it in a separate paper.
From a theoretical viewpoint in a broad context, we expect that we may derive global thermodynamics based on a more microscopic description. Since we know the basic framework of non-equilibrium statistical mechanics, e.g. the linear response theory, we obtain a linear response formula of the temperature of the interface. It is quite formal, and it is not evident that the result satisfies the temperature relation (5.47). Related to this problem, recently, we have calculated the temperature of the interface for a stochastic order parameter equation describing the phase coexistence in heat conduction SNIN . The result shows that the interface temperature deviates from the equilibrium transition temperature. It should be noted that global thermodynamics can be formulated for systems that exhibits an order-disorder transition in heat conduction. By finding a connection between the two theories, we will understand a microscopic mechanism of the interesting macroscopic phenomenon.
Our theory may be naturally generalized to that for heat conduction systems where a liquid flows in and a gas flows out. When a steady state is observed, the system exhibits perpetual evaporation at the interface. The latent heat is persistently generated at the interface and as a result additional heat flux occurs from the interface to the liquid region and gas region. This phenomenon was experimentally studied for pure water Fang-Ward , and a macroscopic temperature gap as much as about K was found at the liquid-gas interface. As a related phenomenon, the temperature gradient in the inverse direction to the imposed temperature difference was observed for the liquid-gas coexistence of a heat conducting Helium far from equilibrium Urban . These phenomena may be closely related to our prediction on the phase coexistence in heat conduction. We will attempt to explain such phenomena in global thermodynamics. Furthermore, we also expect that the description of global quantities may be useful for more complicated systems. Even for phenomena far from equilibrium such as motility-induced phase separation mips , there may be an appropriate set of global variables on the basis of local variables whose space-time evolution may be described by another evolution rule. We will accumulate such examples.
We hope that this paper is a good starting point of many studies.
Acknowledgment
The authors thank Masato Itami, Michikazu Kobayashi, Christian Maes, Yohei Nakayama, Kazuya Saito, Hal Tasaki, and Satoshi Yukawa for their useful comments. One of the authors (N.N.) also thanks Yoshiyuki Chiba, Daiki Takei, Yuya Kai, Kyosuke Tachi and Akira Yoshida for their efforts on numerical experiments of subjects related to §3. The present study was supported by KAKENHI (Nos. JP17H01148, JP19H05496, JP19H05795 and JP19K03647).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) Landau, L. D., Lifshitz, E. M. : Fluid Mechanics. Pergamon Press, Oxford (1959).
- 2(2) Zubarev, D. N. : Nonequilibrium Statistical Thermodynamics. Consultants Bureau, New York (1974).
- 3(3) Mclennan, J. A. : Phys. Fluids 3 , 493-502 (1960); Introduction to Non-equilibrium Statistical Mechanics. Prentice-Hall (1988).
- 4(4) Onsager, L. : Reciprocal relations in irreversible processes. I. Phys. Rev. 37 , 405-426 (1931); Reciprocal relations in irreversible processes. II. Phys. Rev. 38 , 2265-2279 (1931).
- 5(5) de Groot, S. and Mazur, P. : Non-equilibrium Thermodynamics. North-Holland, Amsterdam (1962)
- 6(6) Schmitz, R. : Fluctuations in nonequilibrium fluids, Physics Reports 171 , 1-58 (1988).
- 7(7) Kubo, R., Toda, M., Hashitsume, N. : Statistical Physics II. Springer, (1985).
- 8(8) Nakano, H. : Linear response theory - historical perspective. Int. J. Mod. Phys. B 7:2397-2467 (1993).