Skip to content

Conservation Laws and Hydrodynamic Fields

Relativistic hydrodynamics evolves expectation values of the stress tensor and conserved currents. A velocity field is introduced to decompose those observables into local energy, momentum, stress, charge, and spatial flux; it is not itself a uniquely defined microscopic observable away from equilibrium.

Required background. Current Sources and Generating Functionals derives currents by varying background sources. Spacetime Currents, Stress Tensors, and Charge Algebras fixes the stress-tensor and charge conventions used here.

Helpful background. Hydrodynamic Effective-Theory Architecture explains the derivative expansion that later closes these equations.

For a nonanomalous U(1)U(1) current coupled to a background gauge field AμA_\mu and a stress tensor coupled to gμνg_{\mu\nu},

μJμ=0,μTμν=FνλJλ,Fμν=μAννAμ.\nabla_\mu J^\mu=0, \qquad \nabla_\mu T^{\mu\nu}=F^{\nu\lambda}J_\lambda, \qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu.

The second Ward identity includes the work and force supplied by the external field. In a local inertial frame its time component is the power density and its spatial components are the Lorentz-force density, once the displayed definition of FμνF_{\mu\nu} is used consistently. If the current is anomalous, if translations are explicitly broken, or if dynamical electromagnetism is included, the right-hand sides and field content must be changed rather than hidden inside a constitutive coefficient.

Choose a future-directed unit timelike field uμu^\mu, u2=1u^2=1, and the positive rest-space metric Pμν=uμuνgμνP^{\mu\nu}=u^\mu u^\nu-g^{\mu\nu}. The most general symmetric stress tensor and current decompose as

Tμν=Euμuν+PPμν+uμqν+uνqμ+πμν,Jμ=Nuμ+jμ,\begin{aligned} T^{\mu\nu} &= \mathcal E u^\mu u^\nu +\mathcal P P^{\mu\nu} +u^\mu q^\nu+u^\nu q^\mu +\pi^{\mu\nu},\\ J^\mu &=\mathcal N u^\mu+j^\mu, \end{aligned}

with

uμqμ=uμjμ=0,uμπμν=0,πμμ=0.u_\mu q^\mu=u_\mu j^\mu=0, \qquad u_\mu\pi^{\mu\nu}=0, \qquad \pi^\mu{}_{\mu}=0.

This irreducible decomposition and its relativistic-fluid interpretation are developed in Kovtun 2012, §2.1, pp. 15–21, Open PDF and Rezzolla and Zanotti 2013, ch. 2.

The scalars and fluxes are obtained directly from the physical tensors:

E=uμuνTμν,qμ=TμνuνEuμ,N=uμJμ,jμ=JμNuμ.\mathcal E=u_\mu u_\nu T^{\mu\nu}, \quad q^\mu=T^{\mu\nu}u_\nu-\mathcal E u^\mu, \quad \mathcal N=u_\mu J^\mu, \quad j^\mu=J^\mu-\mathcal N u^\mu.

In the local rest frame, T00=ET^{00}=\mathcal E, T0i=qiT^{0i}=q^i, and J0=NJ^0=\mathcal N. The spatial stress splits into its average P\mathcal P and traceless part πij\pi^{ij}. These relations are identities for any chosen uμu^\mu; hydrodynamics enters only when the components are expressed locally in terms of slow fields and sources.

The table is the canonical convention and translation reference for Chapters 11–12. Every equality is in the site’s (+)(+---) convention; dd denotes spacetime dimension.

ObjectDefinition in this volumeLocal-rest-frame or invariant checkStatus under a first-order frame redefinition
Metric and velocitygμν=diag(1,1,)g_{\mu\nu}=\mathrm{diag}(1,-1,\ldots), u2=1u^2=1uμ=(1,0)u^\mu=(1,\mathbf0)Convention fixed; uμu^\mu itself is frame dependent
Rest-space metricPμν=uμuνgμνP^{\mu\nu}=u^\mu u^\nu-g^{\mu\nu}Pij=δijP^{ij}=\delta^{ij}, Pμνuν=0P^{\mu\nu}u_\nu=0Changes when uμu^\mu changes
DerivativesD=uμμD=u^\mu\nabla_\mu, μ=Pμνν\nabla_\perp^\mu=P^{\mu\nu}\nabla_\nuD=tD=\partial_t, i=i\nabla_\perp^i=\partial_i at restBasis elements mix only at the next retained order
Expansion and accelerationθ=μuμ\theta=\nabla_\mu u^\mu, aμ=Duμa^\mu=Du^\muua=0u\cdot a=0Frame dependent
Shearσμν=PμαPνβ(αuβ+βuα+2d1Pαβθ)\sigma^{\mu\nu}=P^{\mu\alpha}P^{\nu\beta}(\nabla_\alpha u_\beta+\nabla_\beta u_\alpha+\tfrac{2}{d-1}P_{\alpha\beta}\theta)transverse and traceless; ησij\eta\sigma^{ij} is the usual negative viscous shear stressTensor coefficient η\eta is invariant at first order
Vorticityωμν=12PμαPνβ(αuββuα)\omega^{\mu\nu}=\tfrac12P^{\mu\alpha}P^{\nu\beta}(\nabla_\alpha u_\beta-\nabla_\beta u_\alpha)transverse and antisymmetricNot an independent parity-even first-order stress term
Stress and currentTμν=Euu+PP+2u(μqν)+πμνT^{\mu\nu}=\mathcal E uu+\mathcal P P+2u^{(\mu}q^{\nu)}+\pi^{\mu\nu}; Jμ=Nuμ+jμJ^\mu=\mathcal N u^\mu+j^\mureconstructing T,JT,J must give the original tensorsTμνT^{\mu\nu} and JμJ^\mu are invariant; components are not
Ideal dataT(0)μν=ϵuu+pPT_{(0)}^{\mu\nu}=\epsilon uu+pP, J(0)μ=nuμJ_{(0)}^\mu=n\,u^\muw=ϵ+p=Ts+μnw=\epsilon+p=Ts+\mu nEquation of state is invariant; its field coordinates change
Landau frameqμ=0q^\mu=0 through the chosen orderuμu^\mu is the timelike energy-flow eigenvector when it existsA matching convention, not an observable condition
Eckart framejμ=0j^\mu=0 through the chosen order, requiring n0n\ne0uμu^\mu follows charge flowSingular at n=0n=0 and not automatically stable or causal
Invariant diffusion currentJμ=jμ(n/w)qμ\mathcal J^\mu=j^\mu-(n/w)q^\muunchanged when qqwδuq\mapsto q-w\delta u, jjnδuj\mapsto j-n\delta uFrame invariant through first order
Invariant scalar stressB=δPpϵδEpnδN\mathcal B=\delta\mathcal P-p_\epsilon\delta\mathcal E-p_n\delta\mathcal Nremoves shifts along the equilibrium equation-of-state surfaceFrame invariant through first order
Physical linear dataretarded poles, residues of fixed operators, Kubo coefficientscompute from T,JT,J correlatorsInvariant through the consistently transformed truncation order

This is the chapter’s canonical translation table. Its scientific check is reconstruction: after any field redefinition, rebuild TμνT^{\mu\nu} and JμJ^\mu and verify equality through the retained derivative order.

Local equilibrium supplies the zeroth-order functions

E=ϵ(T,μ),P=p(T,μ),N=n(T,μ),\mathcal E=\epsilon(T,\mu), \qquad \mathcal P=p(T,\mu), \qquad \mathcal N=n(T,\mu),

and the thermodynamic identities

dϵ=Tds+μdn,dp=sdT+ndμ,w=Ts+μn.\mathrm d\epsilon=T\,\mathrm ds+\mu\,\mathrm dn, \qquad \mathrm dp=s\,\mathrm dT+n\,\mathrm d\mu, \qquad w=Ts+\mu n.

The entropy current is not an additional exactly conserved microscopic current in a generic interacting theory. At ideal order sμ=suμs^\mu=su^\mu follows from energy and charge conservation and obeys μsμ=0\nabla_\mu s^\mu=0 for smooth solutions. Dissipative corrections produce entropy and allow improvement ambiguities; the inequalities derived from them are treated on Onsager Reciprocity and Entropy Production.

In flat spacetime about rest, take T00=ϵT^{00}=\epsilon, T0i=wviT^{0i}=wv^i, and Jμ=(n,J)J^\mu=(n,\mathbf J). The Ward identities reduce to

tϵ+ ⁣(wv)=E ⁣ ⁣J,\partial_t\epsilon+\boldsymbol\nabla\!\cdot(w\mathbf v) =\mathbf E\!\cdot\!\mathbf J, t(wv)+p=nE+J×B,\partial_t(w\mathbf v)+\boldsymbol\nabla p =n\mathbf E+\mathbf J\times\mathbf B,

up to the order at which products of perturbations are retained. The first equation checks the sign of external work; the second checks that momentum of the matter sector alone is not conserved in a fixed electromagnetic background. If the electromagnetic field is dynamical, its stress tensor must be included so that the total stress is conserved.

Treating a decomposition as a constitutive law. Writing T=Euu+T=\mathcal E uu+\cdots is kinematics. Predictive content starts only after the components are specified as functionals of the slow fields and sources.

Calling entropy an independently conserved charge. Entropy advection is an ideal, smooth-flow consequence. Dissipation, shocks, and fluctuations change the statement.

Forgetting improvements and magnetization. Local representatives of TT and JJ can contain identically conserved pieces. Integrated charges and properly defined transport response remain the comparison targets.

Local Equilibrium and Hydrostatic Constraints determines which parts of these constitutive tensors follow from a stationary generating functional. Hydrodynamic Frames and Constitutive Data then derives the transformations summarized in the reference table.

The schematic places the Ward identities at the entrance to the constitutive construction. Inspect especially the dashed frame-change branch: it changes the fields used to decompose TμνT^{\mu\nu} and JμJ^\mu, not those physical tensors themselves.

Conserved stress and charge currents feed a local-equilibrium and hydrostatic construction, then a constitutive tensor and derivative basis; a dashed frame-change branch leaves physical transport data unchanged, while the final box lists ideal sound, shear, and charge modes.

The physical inputs are the source-covariant conservation laws and the tensors TμνT^{\mu\nu} and JμJ^\mu. Local equilibrium and hydrostatics constrain their constitutive form; a hydrodynamic frame only chooses coordinates on that form. The final mode box refers to the ideal-order truncation. The diagram is schematic and does not encode coefficient values or a causal initial-value theorem.

In text: decompose the fixed tensors relative to a chosen uμu^\mu, match ϵ\epsilon, nn, pp, energy flux, charge flux, and viscous stress, and translate those components order by order under field redefinitions. Ward identities and observable correlators survive that translation even though individual constitutive coefficients need not.

  • Kovtun, Pavel. 2012. “Lectures on Hydrodynamic Fluctuations in Relativistic Theories.” Journal of Physics A: Mathematical and Theoretical 45: 473001. DOI. Open PDF.

  • Rezzolla, Luciano, and Olindo Zanotti. 2013. Relativistic Hydrodynamics. Oxford University Press. DOI.