Journal Article

SPHS: smoothed particle hydrodynamics with a higher order dissipation switch

J. I. Read and T. Hayfield

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 422, issue 4, pages 3037-3055
Published in print June 2012 | ISSN: 0035-8711
Published online May 2012 | e-ISSN: 1365-2966 | DOI: http://dx.doi.org/10.1111/j.1365-2966.2012.20819.x
SPHS: smoothed particle hydrodynamics with a higher order dissipation switch

Show Summary Details

Preview

We present a novel implementation of smoothed particle hydrodynamics that uses the spatial derivative of the velocity divergence as a higher order dissipation switch. Our switch – which is second order accurate – detects flow convergence before it occurs. If particle trajectories are going to cross, we switch on the usual SPH artificial viscosity, as well as conservative dissipation in all advected fluid quantities (e.g. the entropy). The viscosity and dissipation terms (that are numerical errors) are designed to ensure that all fluid quantities remain single valued as particles approach one another, to respect conservation laws, and to vanish on a given physical scale as the resolution is increased. SPHS alleviates a number of known problems with ‘classic’ SPH, successfully resolving mixing, and recovering numerical convergence with increasing resolution. An additional key advantage is that – treating the particle mass similarly to the entropy – we are able to use multimass particles, giving significantly improved control over the refinement strategy. We present a wide range of code tests including the Sod shock tube, Sedov–Taylor blast wave, Kelvin–Helmholtz Instability, the ‘blob test’ and some convergence tests. Our method performs well on all tests, giving good agreement with analytic expectations.

Keywords: hydrodynamics; instabilities; methods: numerical

Journal Article.  12976 words.  Illustrated.

Subjects: Astronomy and Astrophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.