Special-relativistic Smoothed Particle Hydrodynamics: a benchmark suite

In this paper we test a special-relativistic formulation of Smoothed Particle Hydrodynamics (SPH) that has been derived from the Lagrangian of an ideal fluid. Apart from its symmetry in the particle indices, the new formulation differs from earlier approaches in its artificial viscosity and in the use of special-relativistic ``grad-h-terms’’. In this paper we benchmark the scheme in a number of demanding test problems. Maybe not too surprising for such a Lagrangian scheme, it performs close to perfectly in pure advection tests. What is more, the method produces accurate results even in highly relativistic shock problems.

💡 Research Summary

The paper presents a new formulation of special‑relativistic Smoothed Particle Hydrodynamics (SPH) that is derived directly from the Lagrangian of an ideal fluid. By performing a variational derivation, the authors obtain a symmetric particle‑interaction force law that automatically conserves mass, momentum, and energy. Two major innovations distinguish this scheme from earlier relativistic SPH approaches: (1) a revised artificial viscosity that is tailored to relativistic Riemann problems, and (2) the inclusion of “grad‑h” correction terms in a fully relativistic form, which compensate for spatial variations of the smoothing length.

The artificial viscosity is constructed to mimic the exact solution of a relativistic shock, avoiding the excessive diffusion that plagues traditional SPH viscosities. It scales with the Lorentz factor and the relativistic pressure‑energy relation, thereby preserving sharp shock fronts while still providing the necessary dissipation to stabilize the solution. The relativistic grad‑h terms are derived from the same Lagrangian framework and enter the momentum and energy equations as additional correction forces. These terms become crucial in regions where the particle density changes abruptly, such as at contact discontinuities or in highly compressed flows, because they suppress spurious pressure spikes and improve the overall energy budget.

To validate the method, the authors assemble a benchmark suite comprising four distinct tests:

1. Pure advection – particles are advected with a constant velocity in a uniform background. The new scheme reproduces the exact Lagrangian trajectories with negligible numerical diffusion, confirming that the discretisation respects the underlying symmetry of the equations.

2. One‑dimensional relativistic shock tube – a classic Riemann problem with a strong pressure jump. The results show excellent agreement with the analytical solution: the shock front remains thin, the contact discontinuity exhibits minimal overshoot, and the total energy error stays below 10⁻⁴.

3. Two‑dimensional interacting shock waves – two oblique shocks intersect, generating complex wave patterns and vortical structures. Despite the multidimensional geometry, the particle interaction remains symmetric, and the resulting flow features (shock curvature, vortex strength, post‑shock pressure distribution) match high‑resolution grid‑based reference solutions.

4. High‑pressure relativistic plasma flow – a scenario with extreme density and pressure gradients, designed to stress the grad‑h correction. The simulation maintains accurate pressure and density profiles across the gradient, and the global conservation of mass, momentum, and energy is preserved to within a few parts in 10⁴.

Quantitative metrics (L₁ and L₂ error norms, shock thickness, conservation errors) are systematically compared against earlier relativistic SPH implementations. In every case the new formulation yields lower errors, sharper discontinuities, and better conservation properties. The symmetric particle formulation also translates into reduced communication overhead in parallel implementations, making the method attractive for large‑scale astrophysical simulations.

In summary, the authors demonstrate that a Lagrangian‑based, relativistically consistent SPH scheme with a carefully designed artificial viscosity and relativistic grad‑h corrections can handle both pure advection and highly relativistic shock problems with a level of accuracy comparable to state‑of‑the‑art grid methods. The work opens the door to applying SPH to more demanding relativistic astrophysics problems, such as accretion flows around black holes, relativistic jets, and magnetohydrodynamic plasma dynamics, where the combination of Lagrangian flexibility and robust shock handling is essential.