A Density Independent Formulation of Smoothed Particle Hydrodynamics

The standard formulation of the smoothed particle hydrodynamics (SPH) assumes that the local density distribution is differentiable. This assumption is used to derive the spatial derivatives of other quantities. However, this assumption breaks down at the contact discontinuity. At the contact discontinuity, the density of the low-density side is overestimated while that of the high-density side is underestimated. As a result, the pressure of the low (high) density side is over (under) estimated. Thus, unphysical repulsive force appears at the contact discontinuity, resulting in the effective surface tension. This tension suppresses fluid instabilities. In this paper, we present a new formulation of SPH, which does not require the differentiability of density. Instead of the mass density, we adopt the internal energy density (pressure), and its arbitrary function, which are smoothed quantities at the contact discontinuity, as the volume element used for the kernel integration. We call this new formulation density independent SPH (DISPH). It handles the contact discontinuity without numerical problems. The results of standard tests such as the shock tube, Kelvin-Helmholtz and Rayleigh-Taylor instabilities, point like explosion, and blob tests are all very favorable to DISPH. We conclude that DISPH solved most of known difficulties of the standard SPH, without introducing additional numerical diffusion or breaking the exact force symmetry or energy conservation. Our new SPH includes the formulation proposed by Ritchie & Thomas (2001) as a special case. Our formulation can be extended to handle a non-ideal gas easily.Comment: 24 pages, 21 figures. Movies and high resolution figures are available at http://v1.jmlab.jp/~saitoh/sph/index.htm

rpSPH: a novel Smoothed Particle Hydrodynamics Algorithm

We suggest a novel discretisation of the momentum equation for Smoothed Particle Hydrodynamics (SPH) and show that it significantly improves the accuracy of the obtained solutions. Our new formulation which we refer to as relative pressure SPH, rpSPH, evaluates the pressure force in respect to the local pressure. It respects Newtons first law of motion and applies forces to particles only when there is a net force acting upon them. This is in contrast to standard SPH which explicitly uses Newtons third law of motion continuously applying equal but opposite forces between particles. rpSPH does not show the unphysical particle noise, the clumping or banding instability, unphysical surface tension, and unphysical scattering of different mass particles found for standard SPH. At the same time it uses fewer computational operations. and only changes a single line in existing SPH codes. We demonstrate its performance on isobaric uniform density distributions, uniform density shearing flows, the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, the Sod shock tube, the Sedov-Taylor blast wave and a cosmological integration of the Santa Barbara galaxy cluster formation test. rpSPH is an improvement these cases. The improvements come at the cost of giving up exact momentum conservation of the scheme. Consequently one can also obtain unphysical solutions particularly at low resolutions.Comment: 17 pages, 13 figures. Final version. Including section of how to break i

Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH

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

A general class of Lagrangian smoothed particle hydrodynamics methods and implications for fluid mixing problems

Various formulations of smoothed particle hydrodynamics (SPH) have been proposed, intended to resolve certain difficulties in the treatment of fluid mixing instabilities. Most have involved changes to the algorithm which either introduces artificial correction terms or violates what is arguably the greatest advantage of SPH over other methods: manifest conservation of energy, entropy, momentum and angular momentum. Here, we show how a class of alternative SPH equations of motion (EOM) can be derived self-consistently from a discrete particle Lagrangian – guaranteeing manifest conservation – in a manner which tremendously improves treatment of these instabilities and contact discontinuities. Saitoh & Makino recently noted that the volume element used to discretize the EOM does not need to explicitly invoke the mass density (as in the ‘standard’ approach); we show how this insight, and the resulting degree of freedom, can be incorporated into the rigorous Lagrangian formulation that retains ideal conservation properties and includes the ‘∇h’ terms that account for variable smoothing lengths. We derive a general EOM for any choice of volume element (particle ‘weights’) and method of determining smoothing lengths. We then specify this to a ‘pressure–entropy formulation’ which resolves problems in the traditional treatment of fluid interfaces. Implementing this in a new version of the GADGET code, we show it leads to good performance in mixing experiments (e.g. Kelvin–Helmholtz and ‘blob’ tests). And conservation is maintained even in strong shock/blastwave tests, where formulations without manifest conservation produce large errors. This also improves the treatment of subsonic turbulence and lessens the need for large kernel particle numbers. The code changes are trivial and entail no additional numerical expense. This provides a general framework for self-consistent derivation of different ‘flavours’ of SPH

Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests

We study the impact of different discretization choices on the accuracy of SPH and we explore them in a large number of Newtonian and special-relativistic benchmark tests. As a first improvement, we explore a gradient prescription that requires the (analytical) inversion of a small matrix. For a regular particle distribution this improves gradient accuracies by approximately ten orders of magnitude and the SPH formulations with this gradient outperform the standard approach in all benchmark tests. Second, we demonstrate that a simple change of the kernel function can substantially increase the accuracy of an SPH scheme. While the "standard" cubic spline kernel generally performs poorly, the best overall performance is found for a high-order Wendland kernel which allows for only very little velocity noise and enforces a very regular particle distribution, even in highly dynamical tests. Third, we explore new SPH volume elements that enhance the treatment of fluid instabilities and, last, but not least, we design new dissipation triggers. They switch on near shocks and in regions where the flow --without dissipation-- starts to become noisy. The resulting new SPH formulation yields excellent results even in challenging tests where standard techniques fail completely.Comment: accepted for publication in MNRA

ISPH modeling of Rayleigh–Taylor instability

This paper presents a Smoothed Particle Hydrodynamics (SPH) solution to a Rayleigh-Taylor Instability (RTI) problem in an incompressible viscous two-phase immiscible fluid with an interfacial tension. The evolution of the fluid-fluid interface is numerically investigated for four different density ratios. The simulation outcomes are compared with existing results in literature. Three stages of instability, namely the exponential growth rate, the formation of circular form at the crest of spike and the appearance of the final shape of instability, are discussed for different density ratios. It is shown that the numerical algorithm used in this work is capable of capturing the complete physics behind the RTI, such as interface evolution, growth rate and secondary instability accurately, and successfully

An improved SPH scheme for cosmological simulations

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.Comment: 21 figures, 2 tables, accepted to MNRA