Removing the stiffness from interfacial flows with surface tension
Journal of Computational Physics
Journal of Computational Physics
A switch to reduce SPH viscosity
Journal of Computational Physics
Journal of Computational Physics
A multi-phase SPH method for macroscopic and mesoscopic flows
Journal of Computational Physics
An Hamiltonian interface SPH formulation for multi-fluid and free surface flows
Journal of Computational Physics
An unconditionally stable fully conservative semi-Lagrangian method
Journal of Computational Physics
Smoothed particle hydrodynamics and magnetohydrodynamics
Journal of Computational Physics
Constrained hyperbolic divergence cleaning for smoothed particle magnetohydrodynamics
Journal of Computational Physics
A comparison of SPH schemes for the compressible Euler equations
Journal of Computational Physics
Hi-index | 31.47 |
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 dissipative 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. [O. Agertz, B. Moore, J. Stadel, D. Potter, F. Miniati, J. Read, L. Mayer, A. Gawryszczak, A. Kravtsov, A. Nordlund, F. Pearce, V. Quilis, D. Rudd, V. Springel, J. Stone, E. Tasker, R. Teyssier, J. Wadsley, R. Walder, Fundamental differences between SPH and grid methods, MNRAS 380 (2007) 963-978] 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.