---
title: 2D Relativistic Euler Equations
url: https://www.emergentmind.com/topics/two-dimensional-relativistic-euler-equations
type: topic
---

# 2D Relativistic Euler Equations

The two-dimensional relativistic Euler equations describe the evolution of perfect fluids in $\mathbb{R}^{1+2}$ Minkowski spacetime, accounting for relativistic effects arising from large velocities and strong coupling between pressure, energy density, and momentum. These equations govern many phenomena in high-energy astrophysical and plasma settings, including those where thermal and kinematic energies are comparable to rest mass energy. In two dimensions, the system admits both classical and ultra-relativistic equations of state, rich geometric structures, entropy-stable discretizations, precise shock descriptions, well-posedness at low regularity, and a noncanonical Hamiltonian formulation.

## 1. Fundamental Equations and Variants

The two-dimensional relativistic Euler system for a perfect fluid consists of conservation laws for energy-momentum and (where relevant) particle number. In Eulerian coordinates $(t,x)\in\mathbb{R}\times\mathbb{R}^2$, canonical variables are:
- Proper energy density $\varrho(t,x)$,
- Four-velocity $u^\alpha$ normalized by $u^\alpha u_\alpha = -1$ under the Minkowski metric $m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$,
- Pressure $p=p(\varrho)$ specified by the equation of state.

The system is compactly written as:
\[
\partial_\alpha\left[(p+\varrho)u^\alpha u^\beta + p\,m^{\alpha\beta}\right] = 0,
\]
and admits a symmetric-hyperbolic formulation:
\[
\begin{cases}
u^\kappa\partial_\kappa\varrho +(p+\varrho)\,\partial_\kappa u^\kappa =0, \\
(p+\varrho)\,u^\kappa\partial_\kappa u^\alpha +(m^{\alpha\kappa}+u^\alpha u^\kappa)\,\partial_\kappa p =0.
\end{cases}
\]
Ultra-relativistic and isentropic reductions yield closed systems without explicit particle conservation, often with $e=3p$ or $p=K\varrho^A$ for some $A\geq1$ [2508.21427][2512.16090]. Specific forms incorporate barotropic, polytropic, or isothermal state functions, as in $p=\frac{1}{\ell}\varrho$ ($\ell>1$) [2403.11471].

For numerical and geometric analysis, conserved variables such as $D = \rho W$, $m_i = \rho h W^2 v_i$, $W = (1-|v|^2)^{-1/2}$, and enthalpy $h(\rho)$ are convenient [2308.07289].

## 2. Geometric and Analytic Structures

The system is strictly hyperbolic in the regime $\rho>0$, $0< c_s^2 < 1$, $|v|<1$, with eigenvalues given by characteristic speeds:
\[
\lambda_{\pm} = \frac{v\cdot n\pm c_s}{1\pm(v\cdot n)c_s},
\]
\[
c_s^2 = \frac{dp}{d\rho} \Big/ \Big(h+\frac{p}{\rho}\Big)
\]
[2308.07289]. For generic equations of state, a second-order geometric formulation is available, involving wave operators associated with the acoustical metric
\[
g_{\alpha\beta} = c^{-2} m_{\alpha\beta} + (c^{-2}-1) u_\alpha u_\beta,
\]
and its inverse. Fluid variables (logarithmic enthalpy, four-velocity) satisfy covariant wave equations with quadratic null forms as nonlinearities [1809.06204]. Associated transport and div–curl equations for entropy gradients and vorticity show elliptic regularity gain.

In ultra-relativistic regimes, primitive conserved variables $(p,u_x,u_y)$ lead to entropy functionals $\eta = p^{3/4}\sqrt{1+u_x^2+u_y^2}$, main (entropy) variables, and entropy fluxes, facilitating entropy-stable discretization [2508.21427]. The system admits stream-function reductions and a noncanonical Poisson bracket structure, generalizing the Lie–Poisson formalism [2410.22612].

## 3. Entropy and Hamiltonian Structures

Convex entropy functionals enable control of weak solutions and stability. For ultra-relativistic cases, the physical entropy $\eta(\mathbf{w})$ and fluxes $q^k$ satisfy compatible relations:
\[
\nabla_{\mathbf{w}} q^k(\mathbf{w})^T = \nabla_{\mathbf{w}}\eta(\mathbf{w})^T D_{\mathbf{w}}\mathbf{F}^k(\mathbf{w}),\quad k=1,2,
\]
with explicit entropy variables and potentials [2508.21427]. Discrete entropy inequalities can be imposed by entropy-stable numerical fluxes.

The 2D relativistic Euler equations for $\gamma$-barotropic fluids possess a noncanonical Hamiltonian structure with the Poisson bracket
\[
\{F,G\}_{2D} = \int_\Omega \left[\Delta^{-1}\frac{\delta F}{\delta\Psi}, \Psi\right]\Delta^{-1}\frac{\delta G}{\delta\Psi} - (F\leftrightarrow G),
\]
where $\Delta = \nabla\cdot(\gamma\rho\nabla)$, $\gamma = (1 - |u|^2/c^2)^{-1/2}$, and $E[\Psi]=\tfrac12\int_\Omega \gamma\rho\,\omega^2$ is a Casimir [2410.22612]. The Hamiltonian is
\[
H[\Psi] = \int_\Omega \gamma\rho c^2 + S(\rho) \,dx,
\]
where $S'(\rho) = p'(\rho)/\rho$.

## 4. Shock Formation, Blow-Up, and Free Boundaries

Shock waves form by breakdown of strict hyperbolicity and gradient blow-up. Level sets of a suitably defined eikonal function define acoustic characteristics; the density $\mu$ of their inverse foliation vanishes where shocks appear, which signals formation of a Cauchy horizon beyond the classical solution [2308.07289].

Self-similar imploding solutions of the isothermal relativistic Euler equations were constructed for $p=\frac{1}{\ell}\varrho$ in $d=2$, revealing smooth profiles across the sonic point, finite-time blow-up, and explicit asymptotics:
\[
\varrho(t,x) = \frac{\hat{\varrho}(r/(T_*-t))}{(T_*-t)^\beta} \sim \frac{\varrho_*}{|x|^\beta},\quad u(t,x)\sim\frac{x}{|x|}u_*,
\]
illustrating imploding singularity rather than classical shock [2403.11471].

In domains with vacuum boundary, solutions to polytropic equations of state ($p=K\rho^2$ or $p=K\rho^{(\gamma+1)/\gamma}$) maintain finite acceleration near vacuum, with the sound speed vanishing monotonically as $\rho\to0$ but $\nabla s^2\neq0$, ensuring bounded acceleration [1102.4276].

## 5. Discretization and Numerical Methods

Entropy-stable Discontinuous Galerkin (DG) methods provide high-fidelity finite element frameworks for the 2D (and 3D) Euler equations. DG discretization on triangulated domains employs basis polynomials at quadrature nodes, numerical flux differencing for Summation-By-Parts properties, and entropy-stable interface fluxes [2508.21427]. Shock-capturing is managed by switching to finite-volume subcell schemes in troubled elements.

Benchmarking with radially symmetric data in 2D (self-similar shocks, rarefactions, pressure bubbles) demonstrates that DG solutions sharply resolve shock waves, pressure blow-ups, and smooth expansions, closely matching one-dimensional reference codes (Trixi.jl vs. RadSymS).

## 6. Well-Posedness and Regularity

Recent advances established local and global existence, uniqueness, and regularity at lower thresholds than previously believed possible. By introducing log-enthalpy, rescaled velocity, and vorticity as good variables, the Euler equations reduce to a coupled wave–transport system. Local well-posedness holds when
\[
(h_0, v_0) \in H^{7/4+}(\mathbb{R}^2),\quad w_0 \in H^{3/2+},\quad \nabla w_0 \in L^8(\mathbb{R}^2),
\]
with Strichartz and energy bounds [2512.16090].

For the stiff fluid case ($p(\varrho)=\varrho$), the acoustic metric is flat; local well-posedness requires only $H^{7/4+}$ for $h_0, v_0$ and $H^{1+}$ for $w_0$. If the flow is irrotational ($w=0$), local solutions exist for $H^{1+}$ initial data, and global well-posedness holds for small initial data in $\dot{B}^1_{2,1}$ [2512.16090].

Coupled covariant wave, transport, and elliptic equations show enhanced regularity for vorticity and entropy: elliptic div–curl structure upgrades their regularity by one Sobolev derivative, a property crucial to the analysis of shock formation [1809.06204].

## 7. Physical Implications and Comparative Structures

The two-dimensional relativistic Euler equations extend and deform classical incompressible Euler theory: Lorentz weighting by $\gamma$ modulates all derived quantities, including the Laplacian and enstrophy. Relativistically, localized vorticity transports nontrivial energy-momentum density profiles. Conservation of relativistic enstrophy remains a Casimir invariant, controlling turbulence and inverse cascade, with heightened effects as $|u|\to c$ [2410.22612].

Tables below summarize key formulations:

| Formulation        | Unknowns          | Key Structure                | Reference           |
|--------------------|-------------------|------------------------------|---------------------|
| Symmetric-Hyperbolic| $u^\alpha, \varrho$ | Conservation laws; strict hyperbolicity | [2512.16090], [2308.07289] |
| Geometric (Wave–Transport) | $(h, u^\alpha)$ | Acoustical metric; null forms; div–curl | [1809.06204], [2308.07289] |
| Ultra-Relativistic | $(p, u_x, u_y)$  | Convex entropy, entropy potentials | [2508.21427]       |
| Hamiltonian        | Stream function $\Psi$ | Noncanonical Poisson, Casimir enstrophy | [2410.22612]       |

In summary, the two-dimensional relativistic Euler equations admit highly structured PDE, geometric, entropy, and Hamiltonian formulations, with well-understood regularity, shock mechanisms, and numerical methods. This facilitates rigorous analysis and simulation of relativistic fluids in astrophysical, plasma, and high-energy contexts, underpinning both local well-posedness and complex nonlinear phenomena such as shock formation and turbulence.

Source: https://www.emergentmind.com/topics/two-dimensional-relativistic-euler-equations