Action based approach to dissipative relativistic fluid systems

Abstract: We develop an action principle for a relativistic two-fluid system with dissipation. The specific constituents of the model - which serves as a proof of principle - are particles and entropy. The linchpin of the action is the assertion that a given flux is dissipative if its covariant divergence is non-zero. For our model, the particle flux is taken to be conservative while the entropy flux is dissipative. This allows for a "top-down" approach where the general question is geometric. Previous work has shown that new terms (the proper time derivative of matter space "metrics") must be included in the Lagrangian in order to produce equations of motion with terms representing bulk and shear viscosity. In addition to including these terms we show that further terms - interpreted as velocities - can be included. The new action-based model recovers known relativistic formulations of the Cattaneo equation, which results in causal heat propagation. We further advance our understanding by exploring the single-fluid limit by locking the entropy four-velocity to that of the matter component. This reduces the system to a single field equation along with a constraint equation. We show that this constraint leads to a dynamical extension of the standard Tolman red-shift condition. Finally, we provide three example actions (of increasing complexity) which demonstrate that the model is able to reproduce (in the single-fluid limit) the anticipated terms from the relativistic Navier-Stokes equations. In the general case, the action based approach allows for a much richer structure, which may be relevant for realistic models of non-linear dissipative relativistic fluid systems.

• The paper presents a variational formulation that embeds dissipation in relativistic two-fluid systems by modeling separate fluxes for particles and entropy.
• The methodology involves varying the action with respect to dual fluxes and the metric, yielding evolution equations that generalize the Navier-Stokes and Cattaneo heat conduction models.
• The approach compares multiple model hierarchies to capture shear, bulk viscosities, and heat flow while ensuring positive-definite entropy production in near-equilibrium regimes.

Action-Based Modeling of Dissipative Relativistic Two-Fluid Systems

Introduction and Motivation

This paper provides a detailed formulation of a variational (action-based) approach to modeling dissipative relativistic fluid systems with two distinct fluxes: particles and entropy. Traditional treatments of relativistic hydrodynamics encounter deep conceptual and technical obstacles when including dissipation, largely because classical thermodynamics and energy-momentum conservation do not transfer straightforwardly into general relativistic contexts, where global symmetries are absent except in trivial cases. The approach undertaken in this work reframes the problem in geometric and variational terms, distinguishing itself from more conventional, phenomenological treatments of dissipative effects and from linear perturbative expansions around equilibrium.

A core assertion underlying the construction is that dissipative behavior arises if a flux vector field's covariant divergence is nonzero. The model treats the particle flux as conservative and the entropy flux as allowed to possess non-zero divergence, representing dissipation. This setup enables the embedding of dissipative fluxes into solutions of the Einstein field equations in a manner dictated by an underlying action principle.

Theoretical and Mathematical Framework

Two-Fluid Action Principle and Matter Spaces

The model extends the multi-fluid paradigm, representing fluids as fluxes $n^a$ (particles) and $s^a$ (entropy), both treated covariantly. The action is constructed from all scalar combinations of the metric and these fluxes; for particles and entropy, this gives three independent scalars: $n^2$, $s^2$, $x^2 = -g_{ab} n^a s^b$. Essential for general relativity, the state variables and Lagrangian density are formulated in a frame-invariant, coordinates-free fashion.

To facilitate the variation of fluxes while respecting topological constraints (incompressibility, conservation, or non-conservation), the framework associates each fluid with a three-dimensional matter space, tracking the assignment of fluid elements to spacetime worldlines via scalar mappings $\{X^A, X^{\bar{A}}\}$ for particle and entropy spaces, respectively.

Dissipation enters the formalism through the functional dependence of the entropy three-form $s_{abc}$ on fields beyond mere spacetime position—specifically, on variables encoding velocities, metric derivatives (i.e., time derivatives of the matter-space metrics), and geometric deformation. Prior work established the necessity of including terms proportional to these time derivatives to retrieve the structure of bulk and shear viscosities; this work extends the set of allowed variables, introducing terms interpretable as relative velocities and generalizing the action's possible structure.

Equations of Motion and Stress-Energy Tensor

Variation of the action with respect to the independent degrees of freedom (including the particle and entropy fluxes and the metric) yields evolution equations for each flux and the energy-momentum tensor. For the non-dissipative sector, these equations reduce to familiar multi-fluid conservation laws and the generalized Euler equation. When dissipative terms are included, the entropy flux equation contains nonvanishing divergence, and extra evolution/constraint equations emerge, consistent with the presence of relaxation timescales and finite-speed propagation of thermal signals.

A key analytic result is that, for appropriate choices of functional dependence, the variational equations reduce in the appropriate limits to well-known causal heat conduction models, in particular those featuring the Cattaneo equation—a manifestly hyperbolic heat-flux evolution law, unlike the pathological, acausal parabolic Fourier law. In the single-fluid limit, where the entropy and particle four-velocities "lock," the model recovers the expected relativistic generalizations of Navier-Stokes including shear and bulk viscosity, with a natural derivation of the Tolman temperature red-shift relation found as a dynamical constraint.

Thermodynamics, Entropy Production, and the Second Law

An important observation is that the second law of thermodynamics is not automatically enforced by the variational structure; it must be imposed by selecting actions that guarantee positive-definite entropy production. At the fully nonlinear level, this constraint can only be rigorously achieved in the small entropy drift regime, where the action-based formalism can be shown to reduce to a positive-definite entropy evolution analogous to the Navier-Stokes system. This is a significant, non-trivial result, given the notorious difficulty of constructing manifestly causal and stable relativistic dissipative fluid models.

Model Hierarchies and Results

Three classes of actions with increasing complexity are analyzed:

• Model (i): Entropy depends purely on matter space coordinates. Dissipation enters through the entropy flux four-divergence alone. The resulting equations reproduce the causal Cattaneo formulation for heat flow in the relativistic regime. In the single-fluid limit, this yields an extended Tolman law for temperature in curved spacetime.
• Model (ii): Additional dependence on relative velocities and metric degrees of freedom. The model generates new resistive and dissipative vector and tensor structures, modifying the constitutive relation for heat flow but lacking time derivatives necessary for positive-definite entropy creation when far from equilibrium.
• Model (iii): Inclusion of matter-space metric time derivatives. This construction produces the complete set of dissipative channels (shear and bulk viscosities, in addition to heat flow), and, importantly, in the entropy drift/small deviation regime it guarantees positive-definite entropy production, automatically reducing to the relativistic Navier-Stokes system with explicit coefficients for all dissipative mechanisms.

The explicit equations allow for separation of the full nonlinear evolution (with independent particle and entropy flows) from the single-fluid limit, elucidating how relaxation to equilibrium and causal propagation of dissipative signals arise naturally. The constraint equations combine to yield generalizations of well-known results (e.g., the Tolman temperature relation) for nontrivial, even dynamical, spacetimes.

Implications and Future Directions

The variational, action-based formulation outlined in this work provides a principled framework for constructing and analyzing dissipative relativistic fluids, ameliorating issues of causality, stability, and coordinate invariance endemic to previous approaches. In the single-fluid (locked) limit, all expected structure of relativistic hydrodynamics is recovered, but the "top-down" approach explored here uncovers the necessary and sufficient structure of the action to admit physically reasonable (causal, thermodynamically consistent) dissipative dynamics even in the strongly nonlinear regime.

From a practical standpoint, this formulation is directly relevant to the modeling of highly dynamical, relativistically gravitating fluids—e.g., in neutron star mergers, core-collapse supernovae, accretion disks, and cosmological applications—where dissipative effects can be both nonlinear and strongly coupled to the gravitational field, and where traditional near-equilibrium, phenomenological models are likely inadequate or inconsistent.

The theoretical implications are substantial: the construction demonstrates how action-based methods can provide unique insight into the global structure and possible generalizations of relativistic non-equilibrium thermodynamics, especially in regimes where microscopic and statistical underpinnings may be ambiguous or ill-defined.

Avenues for future research involve:

• Exploring the full nonlinear regime of the extended model, particularly situations far from equilibrium
• Application to multi-constituent systems beyond particles and entropy (e.g., superfluids, charged fluids interacting with electromagnetic fields)
• Detailed comparison with kinetic-theory-based derivations of relativistic dissipative hydrodynamics to validate and possibly calibrate the phenomenological coefficients emerging from the action-based model
• Investigation of the stability and causality properties of the derived equations, especially in highly dynamical and inhomogeneous spacetimes

Conclusion

This paper advances the program of reconciling dissipation with general covariance by formulating a broad class of relativistic two-fluid models within a rigorous variational framework. The analysis systematically characterizes the minimal and sufficient structure of the action required to produce causal, stable, and thermodynamically consistent dissipative dynamics for fluids embedded in curved spacetime. The approach yields both the standard phenomenology in familiar limits and a framework capable of probing new regimes inaccessible to traditional methods, providing a robust foundation for theoretical and computational studies of relativistic hydrodynamics under realistic, non-ideal conditions.