Thermodynamics and stability of equilibrium/non-equilibrium steady states in thermodynamically isolated/open systems -- case study for compressible heat conducting fluid

This paper reviews the necessary calculations for constructing a Lyapunov-like functional to analyze the nonlinear stability of both equilibrium and non-equilibrium steady states in thermodynamically isolated and open systems composed of compressible, heat-conducting fluids.

The Gist

Imagine you have a giant, sealed jar filled with a hot, swirling gas. You shake it up, stir it with a spoon, or heat one side. What happens next?

According to our everyday experience, the gas will eventually stop swirling, the temperature will even out, and the gas will settle into a calm, uniform state. It will reach "equilibrium."

This paper is a mathematical proof that explains why this happens and how we can predict it, even when the system is messy, complex, and far from calm. The author, Vít Průša, uses the tools of thermodynamics and fluid mechanics to build a special "stability meter" that proves the system must settle down.

Here is the breakdown of the paper's ideas using simple analogies.

1. The Two Scenarios: The Sealed Jar vs. The Leaky Jar

The paper looks at two different situations:

• The Isolated System (The Sealed Jar): Imagine a perfectly insulated box. No heat gets in or out, and the walls don't move. If you stir the gas inside, it will eventually stop moving and become uniform.
• The Open System (The Leaky Jar): Imagine a box where the walls are kept at a specific temperature (like a radiator), but no gas can escape. The gas inside might settle into a pattern where one side is hot and the other is cool, but it stops changing over time. This is a "steady state."

The big question is: If we start with a chaotic mess, will the system always find its way to this calm state?

2. The Problem with "Entropy" (The Messy Room Analogy)

In physics, we often use a concept called Entropy to measure disorder. The Second Law of Thermodynamics says entropy always increases (the room gets messier). You might think, "If entropy is always going up, we can just use it as a ruler to measure how close we are to the calm state."

The Catch: The paper explains that total entropy is a bad ruler for this job.

• Analogy: Imagine a messy bedroom. You could have a pile of clothes on the floor (high entropy) or clothes scattered everywhere on the bed (also high entropy). Just knowing the "total messiness" is high doesn't tell you how messy it is or if it's getting better. You need a specific way to measure the distance between "chaos" and "order."

3. The Solution: The "Lyapunov" Stability Meter

To solve this, the author builds a special mathematical tool called a Lyapunov functional. Think of this as a specialized "Distance-to-Home" meter.

• How it works: Imagine you are lost in a foggy forest (the chaotic state). You need a device that always points toward home (the calm state) and always shows a number that gets smaller as you get closer.
• The Magic Formula: The author constructs this meter by combining three things:
  1. Kinetic Energy: How fast the gas is moving (the swirling).
  2. Internal Energy: How hot the gas is.
  3. Entropy: How disordered it is.

By mixing these ingredients with specific "weights" (mathematical multipliers), he creates a single number.

• Rule 1: This number is always positive (you are never "negative distance" from home).
• Rule 2: This number is zero only when the gas is perfectly calm and uniform.
• Rule 3: As time passes, this number always decreases. It never goes up.

Because the number is always dropping and can't go below zero, the gas must eventually reach zero. This proves mathematically that the chaotic gas will inevitably settle down.

4. The "Affine Correction" Trick (The Moving Target)

The first part of the paper deals with the sealed jar where the final state is simple (uniform temperature everywhere).

The second part deals with the "leaky jar" where the final state is complex (hot on one side, cool on the other). This is harder because the "home" is moving or shaped differently.

• The Trick: The author uses a clever mathematical shortcut called the Affine Correction Trick.
• Analogy: Imagine you are trying to measure how far a runner is from the finish line.
  • First, you measure their distance from the start line.
  • Then, you realize the finish line isn't at the start.
  • Instead of building a whole new map, you just take your "distance from start" map and subtract the "distance from start to finish" map.
  • The result is a new map that shows the distance from the runner to the finish line.

The author does this with the math. He takes the "calm state" formula he already built and subtracts the "steady state" formula. This instantly gives him the correct "Distance-to-Home" meter for the more complex, open system.

5. Why This Matters

This isn't just about gas in a jar. This logic applies to:

• Weather patterns: Why storms eventually die down or settle into predictable flows.
• Engineering: Designing engines or heat exchangers that won't explode or become unstable.
• Biology: Understanding how cells maintain stable internal temperatures despite external changes.

The Bottom Line

The paper proves that nature has a built-in "brake." Even if you push a fluid system into a wild, chaotic state, the laws of thermodynamics (specifically the way heat and pressure interact) act like a giant, invisible hand that slowly pushes the system back to a stable, calm state.

The author didn't just guess this; he built a mathematical "speedometer" that proves the system is always slowing down and heading toward peace, whether it's in a sealed box or an open one. It's a victory for the idea that order eventually wins over chaos, provided you have the right tools to measure it.

Technical Summary

Technical Summary: Thermodynamics and Stability of Equilibrium/Non-Equilibrium Steady States in Thermodynamically Isolated/Open Systems

1. Problem Statement

The paper addresses the fundamental problem of nonlinear stability analysis for compressible heat-conducting fluids governed by the Navier-Stokes-Fourier (NSF) equations. The author investigates whether solutions to these equations, starting from arbitrary initial conditions, converge to specific steady states as time approaches infinity.

The study focuses on two distinct scenarios:

1. Thermodynamically Isolated Systems: A fluid in a container with no mechanical or thermal exchange with the surroundings (zero velocity and zero heat flux at boundaries). The target state is a spatially homogeneous rest state (constant density $\hat{\rho}$, zero velocity, constant temperature $\hat{\theta}$).
2. Thermodynamically Open Systems: A fluid in a mechanically isolated container but with fixed boundary temperatures (Dirichlet conditions). The target state is a spatially inhomogeneous steady state (constant pressure, spatially varying temperature and density, zero velocity).

The core challenge is to construct a Lyapunov functional (a scalar quantity that decreases monotonically over time and vanishes only at the target state) to prove that the system is globally stable and converges to these steady states, regardless of the magnitude of the initial perturbations (nonlinear regime).

2. Methodology

The author employs a rigorous thermodynamic approach combined with continuum mechanics, moving beyond linearized stability analysis. The methodology involves:

• Linearized Analysis: First, the NSF equations are linearized around the steady state to identify necessary conditions for stability (positivity of specific heat and pressure derivative).
• Construction of Lyapunov Functionals:
  • Isolated Systems: The author constructs a functional based on the maximization of entropy subject to constraints of constant total energy and mass. This involves identifying Lagrange multipliers ($\lambda_1, \lambda_2$) corresponding to the inverse temperature and a chemical potential-like term.
  • Open Systems: The author utilizes an "Affine Correction Trick" (related to Bregman divergence). This technique constructs a functional for a non-equilibrium steady state by taking the functional for the homogeneous equilibrium state and subtracting its value and linear approximation at the steady state.
• Thermodynamic Potentials: The analysis relies heavily on the Helmholtz free energy ($\psi$) and the modified internal energy ($e(\eta, \rho)$), where $\eta$ is the specific entropy and $\rho$ is density. The convexity of these potentials is central to proving the non-negativity of the Lyapunov functionals.
• Bregman Divergence: The paper explicitly links the constructed Lyapunov functionals to the concept of Bregman divergence, a measure of distance generated by a strictly convex function. This provides a geometric interpretation of the stability measure.

3. Key Contributions

A. Derivation of the Lyapunov Functional for Isolated Systems

The paper derives a specific functional, denoted $V_{meq}$, for isolated systems:
$$V_{meq} = \int_\Omega \left[ \frac{1}{2}\rho|v|^2 - \hat{\theta}(\rho\eta - \hat{\rho}\hat{\eta}) + (\rho e - \hat{\rho}\hat{e}) - \left(\frac{\hat{p}_{th}}{\hat{\rho}} + \hat{\psi}\right)(\rho - \hat{\rho}) \right] dv$$

• Identification of Multipliers: The author rigorously identifies the Lagrange multipliers required to make the functional vanish at the equilibrium state and have a zero first variation.
• Non-negativity Proof: It is proven that $V_{meq} \geq 0$ and $V_{meq} = 0$ if and only if the system is at the equilibrium state, provided the standard thermodynamic stability conditions hold ($c_V > 0$ and $\partial p_{th}/\partial \rho > 0$).
• Time Derivative: The time derivative is shown to be strictly non-positive:
  $$\frac{dV_{meq}}{dt} = -\hat{\theta} \int_\Omega \left( \frac{\tilde{\lambda}(\text{div } v)^2 + 2\nu D^\delta:D^\delta + \kappa \nabla \theta \cdot \nabla \theta}{\theta} \right) dv \leq 0$$
  This confirms that the system dissipates energy and converges to the rest state.

B. Extension to Open Systems via Affine Correction

For open systems with spatially inhomogeneous steady states ($\hat{\theta}_{std}(x)$), the paper demonstrates that the Lyapunov functional can be constructed using the Affine Correction Trick:
$$V_{non-eq}(W \| \hat{W}_{std}) = V_{meq}(W \| \hat{W}_{eq}) - V_{meq}(\hat{W}_{std} \| \hat{W}_{eq}) - D_W V_{meq}|_{\hat{W}_{std}} [W - \hat{W}_{std}]$$

• This construction effectively creates a Bregman divergence between the current state and the non-equilibrium steady state.
• The paper proves that for Dirichlet boundary conditions (constant or spatially varying temperature), the time derivative of this new functional remains non-positive, ensuring stability even in the presence of heat flux through the boundaries.

C. Connection to Classical Thermodynamics and Convexity

• The paper bridges the gap between continuum mechanics stability and classical thermodynamics. It shows that the stability conditions ($c_V > 0$, $\partial p/\partial \rho > 0$) are equivalent to the convexity of the modified internal energy $e(\eta, \rho)$ with respect to its natural variables.
• It clarifies that the "relative entropy" or "ballistic free energy" used in previous mathematical literature (e.g., Feireisl's work) is mathematically equivalent to the thermodynamically derived Lyapunov functional presented here.

4. Key Results

1. Global Nonlinear Stability: The paper provides a constructive proof that spatially homogeneous rest states in isolated compressible fluids are globally asymptotically stable. Any initial perturbation (provided it respects mass and energy conservation) decays to the equilibrium state.
2. Stability of Non-Equilibrium Steady States: The methodology successfully extends to open systems, proving the stability of spatially inhomogeneous steady states where temperature varies spatially but velocity is zero.
3. Universality of the Functional: The derived functional is shown to be identical to the "relative energy" functionals used in rigorous mathematical analysis of weak-strong uniqueness, validating the physical intuition behind those mathematical tools.
4. Role of Convexity: The stability of the system is fundamentally tied to the convexity of the thermodynamic potentials. If the material satisfies the standard thermodynamic stability inequalities, the Lyapunov functional is guaranteed to be non-negative and decreasing.

5. Significance

• Theoretical Unification: The paper unifies the physical concept of entropy maximization with the mathematical requirement for Lyapunov stability in infinite-dimensional systems (PDEs). It resolves the question of how to construct a functional that controls the "size" of perturbations in a way that simple entropy does not.
• Rigorous Foundation for Continuum Thermodynamics: It provides a detailed, step-by-step derivation of stability conditions for compressible fluids, moving beyond linearized approximations to full nonlinear stability.
• Methodological Advancement: The use of the "Affine Correction Trick" and the explicit connection to Bregman divergence offers a powerful, generalizable framework for analyzing stability in other complex thermodynamic systems (e.g., mixtures, viscoelastic fluids).
• Validation of Mathematical Models: By showing that the physically derived functional matches the "ballistic free energy" used in advanced mathematical analysis (Feireisl et al.), the paper validates the mathematical techniques used to prove existence and uniqueness of solutions for the Navier-Stokes-Fourier equations.

In summary, this paper provides a comprehensive thermodynamic framework for proving the stability of both equilibrium and non-equilibrium steady states in compressible heat-conducting fluids, demonstrating that the laws of thermodynamics (specifically entropy maximization and convexity) naturally yield the Lyapunov functionals required for rigorous stability analysis.