For smooth accretion discs in SPH, the default quadratic artificial viscosity beta = 2 adds too much numerical viscosity, and linking beta to the switched alpha (beta = 2 alpha) gives the least numerical spreading.
Smoothed Particle Hydrodynamics
1 Pith paper cite this work. Polarity classification is still indexing.
abstract
I present a review of Smoothed Particle Hydrodynamics (SPH), with the aim of providing a mathematically rigorous, clear derivation of the algorithms from first principles. The method of discretising a continuous field into particles using a smoothing kernel is considered, and also the errors associated with this approach. A fully conservative form of SPH is then derived from the Lagrangian, demonstrating the explicit conservation of mass, linear and angular momenta and energy/entropy. The method is then extended to self-consistently include spatially varying smoothing lengths, (self) gravity and various forms of artificial viscosity, required for the correct treatment of shocks. Finally two common methods of time integration are discussed, the Runge-Kutta-Fehlberg and leapfrog integrators, along with an overview of time-stepping criteria.
fields
astro-ph.SR 1years
2025 1verdicts
CONDITIONAL 1representative citing papers
citing papers explorer
-
Minimising the numerical viscosity in Smoothed Particle Hydrodynamics simulations of discs
For smooth accretion discs in SPH, the default quadratic artificial viscosity beta = 2 adds too much numerical viscosity, and linking beta to the switched alpha (beta = 2 alpha) gives the least numerical spreading.