2D Relativistic Euler Equations

• Two-Dimensional Relativistic Euler Equations are hyperbolic PDEs modeling perfect fluids in Minkowski spacetime with significant relativistic effects.
• They incorporate rich geometric, entropy, and Hamiltonian formulations to capture shock formation, ultra-relativistic regimes, and complex fluid dynamics.
• Recent advances demonstrate local well-posedness and effective entropy-stable numerical discretizations, critical for high-energy astrophysical and plasma simulations.

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 $(t,x)\in\mathbb{R}\times\mathbb{R}^2$0 for some $(t,x)\in\mathbb{R}\times\mathbb{R}^2$1 (Thein et al., 29 Aug 2025, Zhang, 18 Dec 2025). Specific forms incorporate barotropic, polytropic, or isothermal state functions, as in $(t,x)\in\mathbb{R}\times\mathbb{R}^2$2 ($(t,x)\in\mathbb{R}\times\mathbb{R}^2$3) (Shao et al., 2024).

For numerical and geometric analysis, conserved variables such as $(t,x)\in\mathbb{R}\times\mathbb{R}^2$4, $(t,x)\in\mathbb{R}\times\mathbb{R}^2$5, $(t,x)\in\mathbb{R}\times\mathbb{R}^2$6, and enthalpy $(t,x)\in\mathbb{R}\times\mathbb{R}^2$7 are convenient (Abbrescia et al., 2023).

2. Geometric and Analytic Structures

The system is strictly hyperbolic in the regime $(t,x)\in\mathbb{R}\times\mathbb{R}^2$8, $(t,x)\in\mathbb{R}\times\mathbb{R}^2$9, $\varrho(t,x)$0, with eigenvalues given by characteristic speeds: $\varrho(t,x)$1

$\varrho(t,x)$2

(Abbrescia et al., 2023). For generic equations of state, a second-order geometric formulation is available, involving wave operators associated with the acoustical metric

$\varrho(t,x)$3

and its inverse. Fluid variables (logarithmic enthalpy, four-velocity) satisfy covariant wave equations with quadratic null forms as nonlinearities (Disconzi et al., 2018). Associated transport and div–curl equations for entropy gradients and vorticity show elliptic regularity gain.

In ultra-relativistic regimes, primitive conserved variables $\varrho(t,x)$4 lead to entropy functionals $\varrho(t,x)$5, main (entropy) variables, and entropy fluxes, facilitating entropy-stable discretization (Thein et al., 29 Aug 2025). The system admits stream-function reductions and a noncanonical Poisson bracket structure, generalizing the Lie–Poisson formalism (Takeda et al., 2024).

3. Entropy and Hamiltonian Structures

Convex entropy functionals enable control of weak solutions and stability. For ultra-relativistic cases, the physical entropy $\varrho(t,x)$6 and fluxes $\varrho(t,x)$7 satisfy compatible relations: $\varrho(t,x)$8 with explicit entropy variables and potentials (Thein et al., 29 Aug 2025). Discrete entropy inequalities can be imposed by entropy-stable numerical fluxes.

The 2D relativistic Euler equations for $\varrho(t,x)$9-barotropic fluids possess a noncanonical Hamiltonian structure with the Poisson bracket

$u^\alpha$0

where $u^\alpha$1, $u^\alpha$2, and $u^\alpha$3 is a Casimir (Takeda et al., 2024). The Hamiltonian is

$u^\alpha$4

where $u^\alpha$5.

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 $u^\alpha$6 of their inverse foliation vanishes where shocks appear, which signals formation of a Cauchy horizon beyond the classical solution (Abbrescia et al., 2023).

Self-similar imploding solutions of the isothermal relativistic Euler equations were constructed for $u^\alpha$7 in $u^\alpha$8, revealing smooth profiles across the sonic point, finite-time blow-up, and explicit asymptotics: $u^\alpha$9 illustrating imploding singularity rather than classical shock (Shao et al., 2024).

In domains with vacuum boundary, solutions to polytropic equations of state ($u^\alpha u_\alpha = -1$0 or $u^\alpha u_\alpha = -1$1) maintain finite acceleration near vacuum, with the sound speed vanishing monotonically as $u^\alpha u_\alpha = -1$2 but $u^\alpha u_\alpha = -1$3, ensuring bounded acceleration (Oliynyk, 2011).

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 (Thein et al., 29 Aug 2025). 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

$u^\alpha u_\alpha = -1$4

with Strichartz and energy bounds (Zhang, 18 Dec 2025).

For the stiff fluid case ($u^\alpha u_\alpha = -1$5), the acoustic metric is flat; local well-posedness requires only $u^\alpha u_\alpha = -1$6 for $u^\alpha u_\alpha = -1$7 and $u^\alpha u_\alpha = -1$8 for $u^\alpha u_\alpha = -1$9. If the flow is irrotational ($m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$0), local solutions exist for $m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$1 initial data, and global well-posedness holds for small initial data in $m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$2 (Zhang, 18 Dec 2025).

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 (Disconzi et al., 2018).

7. Physical Implications and Comparative Structures

The two-dimensional relativistic Euler equations extend and deform classical incompressible Euler theory: Lorentz weighting by $m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$3 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 $m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$4 (Takeda et al., 2024).

Tables below summarize key formulations:

Formulation	Unknowns	Key Structure	Reference
Symmetric-Hyperbolic	$m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$5	Conservation laws; strict hyperbolicity	(Zhang, 18 Dec 2025, Abbrescia et al., 2023)
Geometric (Wave–Transport)	$m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$6	Acoustical metric; null forms; div–curl	(Disconzi et al., 2018, Abbrescia et al., 2023)
Ultra-Relativistic	$m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$7	Convex entropy, entropy potentials	(Thein et al., 29 Aug 2025)
Hamiltonian	Stream function $m_{\alpha\beta} = \operatorname{diag}(-1,1,1)$8	Noncanonical Poisson, Casimir enstrophy	(Takeda et al., 2024)

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.