Global Linear-Irreversible Principle (GLIP)

• GLIP is a foundational variational principle in linear irreversible thermodynamics that maximizes entropy production under conservation constraints near equilibrium.
• It unifies Onsager’s least-dissipation and maximum entropy production principles using linear flux-force relations derived from fundamental conservation laws.
• The principle provides analytic models for finite and cyclic systems, informing applications from thermoelectric devices to climate models.

The Global Linear-Irreversible Principle (GLIP) is a foundational variational principle in linear irreversible thermodynamics. It asserts that, for systems sufficiently close to equilibrium and subject only to conservation-law constraints (such as the first law of thermodynamics and boundary conditions), the physically realized state is the one that maximizes the total (global) entropy production. GLIP unifies several classical principles—including Onsager's least-dissipation principle, the maximum entropy production principle, and global flux-force relations—across continuous media, cyclic machines, and finite-sized devices, creating a rigorous link between microscopic transport laws and macroscopic observable behavior (Zupanovic et al., 2010, Johal, 2017, Kaur et al., 2020, Izumida, 2021).

1. Formal Statement and Mathematical Structure

GLIP applies to continuous thermodynamic systems near equilibrium, defined by a set of local intensive forces $\{X_i(\mathbf{r})\}$ and conjugate fluxes $\{J_i(\mathbf{r})\}$. The total entropy production is

$$\Sigma = \int_V \sigma(\mathbf{r})\,dV = \int_V \sum_i J_i(\mathbf{r}) X_i(\mathbf{r})\,dV,$$

where the local entropy-production density is $$\sigma(\mathbf{r}) = \sum_i J_i(\mathbf{r}) X_i(\mathbf{r})$$. The system is constrained by global conservation laws, such as total energy conservation, and physically warranted boundary conditions.

The fluxes are assumed to obey linear constitutive relations (Onsager relations): $J_i(\mathbf{r}) = \sum_j L_{ij} X_j(\mathbf{r}),$ with a positive-definite, symmetric matrix $L_{ij}$ of Onsager coefficients. The only constraints considered are those dictated by fundamental conservation laws and imposed boundary conditions—no additional "ad hoc" constraints (e.g., fixed fluxes) are allowed within the GLIP framework.

2. Variational Derivation and Onsager–Maximum Entropy Production Equivalence

GLIP establishes that Onsager's principle of least dissipation and the maximum entropy production (MEP) principle are mathematically equivalent under appropriate constraints. The maximization of $\Sigma$ under energy conservation yields the physically realized state.

The Lagrangian functional is constructed as

$$\mathcal{F}[J] = \Sigma[J] - \lambda \left\{ \int_V \left( \sum_i X_i J_i - \sum_{ij} R_{ij} J_i J_j \right) dV \right\},$$

where $R_{ij} = (L^{-1})_{ij}$ and $\lambda$ is a Lagrange multiplier enforcing energy balance. Variation with respect to $J_k(\mathbf{r})$ yields, after imposing the constraint, the linear flux–force relations: $J_i(\mathbf{r}) = \sum_j L_{ij} X_j(\mathbf{r}).$ Second variation analysis demonstrates that this stationary point is a maximum of $\Sigma$, confirming the MEP character of GLIP (Zupanovic et al., 2010).

3. Global Flux–Force Relations in Finite and Cyclic Systems

GLIP generalizes the local, infinitesimal linear-response theory to global, finite-size, and/or cyclic systems:

• Finite Thermoelectric Devices: For a thermoelectric generator of finite length and area, global fluxes (total electric current $I$, hot-side heat flow $\dot{Q}_h$) and global forces ($X_1 = -\Delta V/T_c$, $X_2 = 1/T_c - 1/T_h$) are defined via the global entropy production, $\dot S = I X_1 + \dot Q_{h} X_2$. The driven flux (current) remains exactly linear, $I = \mathcal L_{11} X_1 + \mathcal L_{12} X_2$, while the driver flux ($\dot Q_h$) acquires an additional quadratic term reflecting Joule dissipation. The global kinetic coefficients maintain Onsager reciprocity, $\mathcal L_{12} = \mathcal L_{21}$, inherited from the local coefficients (Kaur et al., 2020).
• Cyclic Heat Engines: For cyclic thermal machines, GLIP prescribes that the total cycle duration is determined by

$$\tau = \frac{\bar Q^2}{\lambda\,\Delta_{\rm tot} S},$$

where $\bar Q$ is an effective heat (an algebraic mean of heats supplied and rejected), $\lambda$ is an effective conductance, and $\Delta_{\rm tot} S$ is the total entropy generated (Johal, 2017). This allows for explicit analytic expressions for efficiency, power, and other figures of merit at optimal performance, reducing many optimization problems to one-parameter maximizations.

4. Hierarchical Structure and Protocol Dependence

For cyclic and adiabatically driven systems, GLIP elucidates a hierarchical structure connecting local and global Onsager coefficients:

$$L_{ij}^{\rm(g)} = \frac{1}{\tau_{\rm cyc}} \int_0^{\tau_{\rm cyc}} \alpha_i(t) A_{ij}(t) \alpha_j(t) dt,$$

where $A_{ij}(t)$ are the instantaneous local Onsager coefficients and $\alpha_i(t)$ are modulation functions encoding protocol shape (Izumida, 2021). Onsager reciprocity at the local level (detailed balance and microreversibility) is inherited by the global coefficients, establishing a structural link between microscopic and macroscopic irreversibility.

Efficiencies in such engines are tightly bounded: $$\eta \leq \eta_C \left(1 - \frac{T_c [\Delta S_{\rm ad}]^2}{L_{ww}^{(g)}} \right),$$ with the gap to Carnot controlled by the diagonal global Onsager coefficient and the protocol-dependent adiabatic entropy change.

5. Applications and Physical Examples

GLIP encompasses and rigorously explains a wide array of physical phenomena and engineering systems:

• Kinetic Theory: The stationary solution to the linearized Boltzmann equation maximizes entropy-production rate under the condition of local balance of collisional and flux-induced entropy production (Zupanovic et al., 2010).
• Thermomechanics: Classical Ziegler thermomechanics is equivalent to maximizing global entropy production under energy constraints, reproducing Onsager's constitutive equations.
• Climate Models: Paltridge's planetary heat transport model employs GLIP at the planetary scale, with the observed distribution maximizing global entropy production (Zupanovic et al., 2010).
• Resistor Networks: For planar electrical networks, GLIP recovers Kirchhoff's current-division law when only energy conservation is imposed; fixing boundary currents (an ad hoc constraint) instead corresponds to a minimization principle and lies outside the GLIP regime.
• Anisotropic Heat Conduction: Maximizing entropy production with respect to temperature gradients yields the anisotropic Fourier law, equivalent to energy-conservation-constrained maximization.

6. Universality, Limitations, and Broader Implications

GLIP provides a unifying framework for diverse linear irreversible systems, both continuous and discrete, steady and cyclic. Near equilibrium, universal predictions for optimal efficiency and power (such as Curzon–Ahlborn efficiency and first- and second-order universality coefficients) emerge solely from the linear structure and symmetry of Onsager coefficients (Johal, 2017, Izumida, 2021). This minimal embedding clarifies that details of finite-time protocols influence only higher-order deviations.

The limitations of GLIP are explicit:

• Nonlinear flux–force relations, finite reservoir size effects, heat leaks, and multi-reservoir interactions fall outside the scope of the current principle.
• The selection of effective means (for heat flows in cyclic engines) remains phenomenological.
• Breakdown of local reciprocity (e.g., under magnetic fields or time-reversal symmetry breaking) destroys the global reciprocity structure.

A plausible implication is that extending GLIP to fully nonlinear, far-from-equilibrium, or strongly correlated settings is an open direction.

7. Relationship to Other Extremum Principles and Theorems

Unlike Prigogine's minimum entropy production theorem, which assumes both fixed forces and fluxes (creating an over-constrained scenario), GLIP is anchored solely in conservation-law constraints appropriate to physical energy balance. No contradiction arises between GLIP and the minimum production theorem as long as one respects the correct constraint structure. GLIP thus stands as the foundational extremum principle for near-equilibrium, linear irreversible systems under physically consistent global constraints, unifying phenomenological laws, kinetic theory, and optimization of real devices (Zupanovic et al., 2010).