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 $p=K\varrho^A$ for some $A\geq1$ (Thein et al., 29 Aug 2025, Zhang, 18 Dec 2025). Specific forms incorporate barotropic, polytropic, or isothermal state functions, as in $p=\frac{1}{\ell}\varrho$ ($\ell>1$) (Shao et al., 2024).

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 (Abbrescia et al., 2023).

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)$$

(Abbrescia et al., 2023). 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 (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 $(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 (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 $\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 (Thein et al., 29 Aug 2025). 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 (Takeda et al., 2024). 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 (Abbrescia et al., 2023).

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 (Shao et al., 2024).

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 (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

$$(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 (Zhang, 18 Dec 2025).

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}$ (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 $\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$ (Takeda et al., 2024).

Tables below summarize key formulations:

Formulation	Unknowns	Key Structure	Reference
Symmetric-Hyperbolic	$u^\alpha, \varrho$	Conservation laws; strict hyperbolicity	(Zhang, 18 Dec 2025, Abbrescia et al., 2023)
Geometric (Wave–Transport)	$(h, u^\alpha)$	Acoustical metric; null forms; div–curl	(Disconzi et al., 2018, Abbrescia et al., 2023)
Ultra-Relativistic	$(p, u_x, u_y)$	Convex entropy, entropy potentials	(Thein et al., 29 Aug 2025)
Hamiltonian	Stream function $\Psi$	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.