In this paper we discuss the treatment of discontinuities in Smoothed
Particle Hydrodynamics (SPH) simulations. In particular we discuss the
difference between integral and differential representations of the fluid
equations in an SPH context and how this relates to the formulation of dissip
ative terms for the capture of shocks and other discontinuities.
This has important implications for many problems, in particular related to
recently highlighted problems in treating Kelvin-Helmholtz instabilities across
entropy gradients in SPH. The specific problems pointed out by Agertz et al.
(2007) are shown to be related in particular to the (lack of) treatment of
contact discontinuities in standard SPH formulations which can be cured by the
simple application of an artificial thermal conductivity term. We propose a new
formulation of artificial thermal conductivity in SPH which minimises
dissipation away from discontinuities and can therefore be applied quite
generally in SPH calculations.Comment: 31 pages, 10 figures, submitted to J. Comp. Phys. Movies + hires
version available at http://www.astro.ex.ac.uk/people/dprice/pubs/kh/ . v3:
modified as per referee's comments - comparison with Ritchie & Thomas
formulation added, quite a few typos fixed. No major change in metho