Conservative, special-relativistic smoothed particle hydrodynamics

Abstract

We present and test a new, special-relativistic formulation of Smoothed Particle Hydrodynamics (SPH). Our approach benefits from several improvements with respect to earlier relativistic SPH formulations. It is self-consistently derived from the Lagrangian of an ideal fluid and accounts for special-relativistic "grad-h terms". In our approach, we evolve the canonical momentum and the canonical energy per baryon and thus circumvent some of the problems that have plagued earlier formulations of relativistic SPH. We further use a much improved artificial viscosity prescription which uses the extreme local eigenvalues of the Euler equations and triggers selectively on a) shocks and b) velocity noise. The shock trigger accurately monitors the relative density slope and uses it to fine-tune the amount of artificial viscosity that is applied. This procedure substantially sharpens shock fronts while still avoiding post-shock noise. If not triggered, the viscosity parameter of each particle decays to zero. None of these viscosity triggers is specific to special relativity, both could also be applied in Newtonian SPH. The performance of the new scheme is explored in a large variety of benchmark tests where it delivers excellent results. Generally, the grad-h terms deliver minor, though worthwhile, improvements. The scheme performs close to perfect in supersonic advection tests, but also in strong relativistic shocks, usually considered a particular challenge for SPH, the method yields convincing results. For example, due to its perfect conservation properties, it is able to handle Lorentz-factors as large as $\gamma =50000$ in the so-called wall shock test. Moreover, we find convincing results in a rarely shown, but challenging test that involves so-called relativistic simple waves and also in multi-dimensional shock tube tests.

39 pages, 19 figures, Journal of Computational Physics in press, reference upated

References in corpus (9)

• The Physics of Gamma-Ray Bursts
• Turbulent gas motions in galaxy cluster simulations: The role of SPH viscosity
• Astrophysical Smooth Particle Hydrodynamics
• Resolving mixing in Smoothed Particle Hydrodynamics
• MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger applications
• Kelvin-Helmholtz instabilities in Smoothed Particle Hydrodynamics
• Smoothed Particle Hydrodynamics
• 3D Relativistic SPH
• Special-relativistic Smoothed Particle Hydrodynamics: a benchmark suite

Cited by in corpus (15)

• Binary neutron-star mergers: a review of Einstein's richest laboratory
• A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
• THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics
• SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch
• TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh
• Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
• Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests
• Hydrodynamics in full general relativity with conservative AMR
• The Lagrangian hydrodynamics code MAGMA2
• General relativistic smoothed particle hydrodynamics
• SPHINCS_BSSN: A general relativistic Smooth Particle Hydrodynamics code for dynamical spacetimes
• SPH Methods in the Modelling of Compact Objects
• Tensor classification of structure in smoothed particle hydrodynamics density fields
• Initial data for general relativistic SPH with Centroidal Voronoi Tessellations
• Special Relativistic Smoothed Particle Hydrodynamics Based on Riemann Solver