Global thermodynamics for heat-conducting fluids under weak gravity
Pith reviewed 2026-06-28 20:09 UTC · model grok-4.3
The pith
The free-energy decomposition shows effective gravity governs liquid-gas transitions under heat conduction while residual heat term ensures consistency.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We construct a variational free-energy function for the fixed-global-temperature description and decompose it into two parts. The first has the same configurational form as the equilibrium weak-gravity free energy with gravity replaced by the effective gravity, and it determines the first-order configurational transition between the two separated liquid-gas arrangements. The second is a residual excess-latent-heat contribution that vanishes without heat conduction. Although it does not decide which separated liquid-gas arrangement is thermodynamically favored, this residual part is needed to derive the fundamental relation in the laboratory variables and to recover thermodynamic observables
What carries the argument
The decomposition of the variational free-energy function into a configurational part using effective gravity and a residual excess-latent-heat part.
If this is right
- The first-order transition between liquid-gas arrangements is controlled by effective gravity.
- Thermodynamic observables including spatially averaged pressure require the residual contribution for correct recovery.
- The residual term modifies the geometry of barriers, ridges, valleys, and interfaces in the free-energy landscape.
- Van der Waals model calculations illustrate the altered landscape under heat conduction.
- Experimental setups can detect effective-gravity inversion at accessible scales.
Where Pith is reading between the lines
- The framework could apply to other gradient-driven phase separations beyond gravity.
- Measuring pressure averages might distinguish the residual contribution experimentally.
- The landscape changes could affect nucleation rates or interface dynamics in heat-conducting systems.
- Extensions to stronger gravity or different fluids may reveal further nonequilibrium effects.
Load-bearing premise
The variational free-energy function for the fixed-global-temperature description can be decomposed into the configurational part with effective gravity and the residual excess-latent-heat contribution that vanishes without heat conduction.
What would settle it
Observation of whether the spatially averaged pressure in a heat-conducting fluid under gravity matches the value predicted only when the residual excess-latent-heat term is included in the free-energy derivation.
read the original abstract
We study liquid-gas coexistence under gravity and heat conduction from the viewpoint of global thermodynamics. We construct a variational free-energy function for the fixed-global-temperature description and decompose it into two parts. The first has the same configurational form as the equilibrium weak-gravity free energy with gravity replaced by the effective gravity, and it determines the first-order configurational transition between the two separated liquid-gas arrangements. The second is a residual excess-latent-heat contribution that vanishes without heat conduction. Although it does not decide which separated liquid-gas arrangement is thermodynamically favored, this residual part is needed to derive the fundamental relation in the laboratory variables and to recover thermodynamic observables such as the spatially averaged pressure. The same residual contribution reshapes the barrier geometry, ridge/valley structure, and interfacial anomalies of the fixed-global-temperature free-energy landscape. Numerical examples based on the van der Waals model illustrate the resulting landscape structure, and estimates of experimental scales suggest a setup for detecting the effective-gravity inversion.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs a variational free-energy function for liquid-gas coexistence in heat-conducting fluids under weak gravity at fixed global temperature. It decomposes this free energy into a configurational part identical in form to the equilibrium weak-gravity free energy (with g replaced by an effective gravity g_eff) that controls the first-order configurational transition between separated liquid-gas arrangements, plus a residual excess-latent-heat contribution that vanishes without heat conduction. The residual term is required to recover the fundamental relation in laboratory variables and observables such as spatially averaged pressure; it also modifies the free-energy landscape. Numerical results from the van der Waals model illustrate the landscape, and experimental scales are estimated for detecting g_eff inversion.
Significance. If the asserted decomposition holds without additional assumptions, the work provides a systematic global-thermodynamic treatment of non-equilibrium phase coexistence under simultaneous gravity and heat flow, extending equilibrium weak-gravity results via an effective gravity while isolating the role of latent-heat transport. The separation into configuration-determining and observable-recovering parts, together with the van der Waals numerics and experimental estimates, offers a concrete route to testable predictions in steady-state thermodynamics.
major comments (3)
- [Abstract / §2–3] Abstract and the central construction (presumably §2–3): the decomposition of the fixed-global-temperature variational free energy into F_config(g_eff) + F_residual is asserted to be exact, with F_config identical in form to the known equilibrium functional and F_residual vanishing identically when the heat current is zero. No explicit derivation from the global-thermodynamic postulates or from the van der Waals equation of state is supplied; without this step the claim that F_config alone determines the first-order transition cannot be verified and remains the load-bearing assumption.
- [Abstract / §3] The definition and independence of the effective gravity g_eff: the manuscript introduces g_eff as the replacement that makes the configurational part match the equilibrium form, yet provides no independent determination (e.g., from local force balance or from the steady-state heat-current constraint) that would avoid circularity when locating the transition point.
- [§4] Recovery of laboratory observables: the claim that the residual term is required to obtain the fundamental relation in laboratory variables and the spatially averaged pressure is central, but the explicit mapping from the decomposed free energy to these quantities (including any integration over the interface or averaging procedure) is not shown in sufficient detail to confirm consistency with known limits (zero heat current, zero gravity).
minor comments (2)
- [Abstract] Notation for the residual excess-latent-heat contribution should be introduced with an explicit symbol and its functional dependence on the heat current stated at first use.
- [§5] The numerical examples would benefit from a brief statement of the discretization or minimization method used to obtain the landscape structure.
Simulated Author's Rebuttal
We thank the referee for the thorough reading and for identifying the points where additional detail would strengthen the manuscript. We address each major comment below and will incorporate the requested clarifications and derivations in a revised version.
read point-by-point responses
-
Referee: [Abstract / §2–3] Abstract and the central construction (presumably §2–3): the decomposition of the fixed-global-temperature variational free energy into F_config(g_eff) + F_residual is asserted to be exact, with F_config identical in form to the known equilibrium functional and F_residual vanishing identically when the heat current is zero. No explicit derivation from the global-thermodynamic postulates or from the van der Waals equation of state is supplied; without this step the claim that F_config alone determines the first-order transition cannot be verified and remains the load-bearing assumption.
Authors: We agree that the decomposition step requires an explicit derivation to be fully verifiable. In the revised manuscript we will add a dedicated subsection (new §2.3) that starts from the global-thermodynamic postulates, introduces the fixed-global-temperature variational functional, and carries out the decomposition algebraically, showing that the configurational term recovers the equilibrium weak-gravity form upon the replacement g → g_eff while the residual term is proportional to the heat current and vanishes identically when the current is zero. The same derivation will be repeated for the van der Waals equation of state to confirm that the transition location is indeed controlled solely by the configurational part. revision: yes
-
Referee: [Abstract / §3] The definition and independence of the effective gravity g_eff: the manuscript introduces g_eff as the replacement that makes the configurational part match the equilibrium form, yet provides no independent determination (e.g., from local force balance or from the steady-state heat-current constraint) that would avoid circularity when locating the transition point.
Authors: We will revise §3 to supply an independent expression for g_eff obtained directly from the steady-state force balance and the constraint of constant heat current, without reference to the location of the configurational transition. This expression will be derived from the local momentum and energy equations under the weak-gravity, steady-state assumptions and will be shown to reduce to the known equilibrium gravity when the heat current vanishes, thereby removing any appearance of circularity. revision: yes
-
Referee: [§4] Recovery of laboratory observables: the claim that the residual term is required to obtain the fundamental relation in laboratory variables and the spatially averaged pressure is central, but the explicit mapping from the decomposed free energy to these quantities (including any integration over the interface or averaging procedure) is not shown in sufficient detail to confirm consistency with known limits (zero heat current, zero gravity).
Authors: We accept that the mapping from the decomposed free energy to laboratory observables needs to be written out explicitly. In the revised §4 we will insert the intermediate steps that relate the total variational functional to the fundamental relation in laboratory variables, including the explicit surface integrals across the interface and the spatial averaging procedure for pressure. We will then verify that both the zero-heat-current and zero-gravity limits recover the standard equilibrium expressions, as required. revision: yes
Circularity Check
No significant circularity; construction of variational free energy is the core contribution
full rationale
The paper explicitly constructs the variational free-energy function for the fixed-global-temperature description and then decomposes it by definition into the configurational part (matching equilibrium form with g replaced by g_eff) plus residual term. This is presented as the modeling step rather than a derived claim that reduces to inputs by construction. Numerical illustrations use the independent van der Waals equation of state. No load-bearing self-citations, fitted inputs renamed as predictions, or uniqueness theorems imported from prior author work are indicated in the text. The framework is self-contained against external benchmarks such as the van der Waals model.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The system can be described using a fixed-global-temperature variational principle.
- ad hoc to paper The free-energy decomposes into configurational and residual parts as described.
invented entities (2)
-
effective gravity
no independent evidence
-
residual excess-latent-heat contribution
no independent evidence
Forward citations
Cited by 1 Pith paper
-
Pressure-drop localization and momentum insulation in liquid-gas coexistence Poiseuille flow
In liquid-gas coexistence Poiseuille flow, pressure difference concentrates across the interfacial region and ordinary Poiseuille particle current is strongly reduced due to small parameters A^L and A^G from squared m...
Reference graph
Works this paper leans on
-
[1]
J. R. Thome, Boiling in microchannels: A review of experiment and theory, Int. J. Heat Fluid Flow25, 128 (2004)
2004
-
[2]
Ahlers, L
G. Ahlers, L. I. Berge, and D. S. Cannell, Thermal convection in the presence of a first-order phase change, Phys. Rev. Lett.70, 2399 (1993)
1993
-
[3]
Zhong, D
J.-Q. Zhong, D. Funfschilling, and G. Ahlers, Enhanced heat transport by turbulent two-phase Rayleigh-Benard convection, Phys. Rev. Lett.102, 124501 (2009)
2009
-
[4]
Cristea, G
A. Cristea, G. Gonnella, A. Lamura, and V. Sofonea, A lattice Boltzmann study of phase separation in liquid-vapor systems with gravity, Commun. Comput. Phys.7, 350–361 (2010)
2010
-
[5]
L. E. Czelusniak, L. Cabezas-Gomez, and A. J. Wagner, Effect of gravity on phase transition for liquid–gas simula- tions, Phys. Fluids35, 043324 (2023)
2023
-
[6]
Davis and B
D. Davis and B. Sen Gupta, Kinetics of vapor–liquid and vapor–solid phase separation under gravity, Soft Matter 21, 1012–1023 (2025)
2025
-
[7]
Weiss and G
S. Weiss and G. Ahlers, Nematic-isotropic phase transition in turbulent thermal convection, J. Fluid Mech.737, 308–328 (2013)
2013
-
[8]
Urban, D
P. Urban, D. Schmoranzer, P. Hanzelka, K. R. Sreenivasan, and L. Skrbek, Anomalous heat transport and conden- sation in convection of cryogenic helium, Proc. Natl. Acad. Sci. U.S.A.110, 8036–8039 (2013)
2013
-
[9]
Muscatello, E
J. Muscatello, E. Chac´ on, P. Tarazona, and F. Bresme, Deconstructing temperature gradients across fluid interfaces: The structural origin of the thermal resistance of liquid-vapor interfaces, Phys. Rev. Lett.119, 045901 (2017)
2017
-
[10]
Fang and C
G. Fang and C. A. Ward, Temperature measured close to the interface of an evaporating liquid, Phys. Rev. E59, 417–428 (1999)
1999
-
[11]
Nakagawa and S.-i
N. Nakagawa and S.-i. Sasa, Liquid-gas transitions in steady heat conduction, Phys. Rev. Lett.119, 260602 (2017)
2017
-
[12]
Nakagawa and S.-i
N. Nakagawa and S.-i. Sasa, Global thermodynamics for heat conduction systems, J. Stat. Phys.177, 825–888 (2019). 56 N. Nakagawa and S.-i. Sasa
2019
-
[13]
Nakagawa and S.-i
N. Nakagawa and S.-i. Sasa, Unique extension of the maximum entropy principle to phase coexistence in heat conduction, Phys. Rev. Res.4, 033155 (2022)
2022
-
[14]
Kobayashi, N
M. Kobayashi, N. Nakagawa, and S.-i. Sasa, Control of metastable states by heat flux in the Hamiltonian Potts model, Phys. Rev. Lett.130, 247102 (2023)
2023
-
[15]
Sasa and N
S.-i. Sasa and N. Nakagawa, Non-equilibrium phase coexistence in boundary-driven diffusive systems, J. Stat. Phys. 192, 26 (2025)
2025
-
[16]
Nakagawa, S.-i
N. Nakagawa, S.-i. Sasa, T. Hirao, T. Shiina, K. Tachi, and A. Yoshida, Global thermodynamics for isothermal fluids under weak gravity, J. Stat. Phys.192, 91 (2025)
2025
-
[17]
N. Nakagawa and S.-i. Sasa, Thermodynamic variational principle unifying gravity and heat flow, Phys. Rev. Lett. 136, 057101 (2026), doi:10.1103/bbqy-hptc
-
[18]
Yoshida, N
A. Yoshida, N. Nakagawa, and S.-i. Sasa, Heat-induced liquid hovering in liquid-gas coexistence under gravity, Phys. Rev. Lett.133, 117101 (2024)
2024
-
[19]
Schweizer, H
M. Schweizer, H. C. ¨Ottinger, and T. Savin, Nonequilibrium thermodynamics of an interface, Phys. Rev. E93, 052803 (2016)
2016
-
[20]
Struchtrup and H
H. Struchtrup and H. C. ¨Ottinger, Nonequilibrium liquid-vapor interfaces: Linear and nonlinear descriptions, Phys. Rev. E108, 064801 (2023)
2023
-
[21]
S.-i. Sasa, N. Nakagawa, M. Itami, and Y. Nakayama, Stochastic order parameter dynamics for phase coexistence in heat conduction, Phys. Rev. E103, 062129 (2021)
2021
-
[22]
A. H. Persad and C. A. Ward, Expressions for the evaporation and condensation coefficients in the Hertz-Knudsen relation, Chem. Rev.116, 7727–7767 (2016)
2016
-
[23]
M. T. Rauter, A. Aasen, S. Kjelstrup, and Ø. Wilhelmsen, A comparative study of experiments and theories on steady-state evaporation of water, Chem. Thermodyn. Therm. Anal.8, 100091 (2022)
2022
-
[24]
E. W. Lemmon, I. H. Bell, M. L. Huber, and M. O. McLinden, Thermophysical Properties of Fluid Systems, inNIST Chemistry WebBook, NIST Standard Reference Database Number 69, eds. P. J. Linstrom and W. G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899,https://doi.org/10.18434/ T4D303, retrieved May 24, 2026
2026
-
[25]
Touchette, The large deviation approach to statistical mechanics, Phys
H. Touchette, The large deviation approach to statistical mechanics, Phys. Rep.478, 1–69 (2009)
2009
-
[26]
Bertini, A
L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, and C. Landim, Macroscopic fluctuation theory, Rev. Mod. Phys.87, 593–636 (2015)
2015
-
[27]
Uruno and N
S. Uruno and N. Nakagawa, in preparation (2026)
2026
-
[28]
J. W. Cahn and J. E. Hilliard, Free energy of a nonuniform system. I. Interfacial free energy, J. Chem. Phys.28, 258–267 (1958)
1958
-
[29]
J. S. Rowlinson and B. Widom,Molecular Theory of Capillarity(Clarendon Press, Oxford, 1982)
1982
-
[30]
Qu´ er´ e, Leidenfrost dynamics, Annu
D. Qu´ er´ e, Leidenfrost dynamics, Annu. Rev. Fluid Mech.45, 197–215 (2013)
2013
-
[31]
Bourrianne, C
P. Bourrianne, C. Lv, and D. Qu´ er´ e, The cold Leidenfrost regime, Sci. Adv.5, eaaw0304 (2019)
2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.