A DENSITY-INDEPENDENT FORMULATION OF SMOOTHED PARTICLE HYDRODYNAMICS

• Published in: The Astrophysical journal Vol. 768; no. 1; pp. 1 - 24
• Main Authors: Saitoh, Takayuki R, Makino, Junichiro
• Format: Journal Article
• Language: English
• Published: United States 01-05-2013
• Subjects: ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; Computational fluid dynamics; Contact; DENSITY; DIFFUSION; Discontinuity; DISTRIBUTION; ENERGY CONSERVATION; ENERGY DENSITY; EVOLUTION; EXPLOSIONS; Fluid flow; GALAXIES; High density; HYDRODYNAMICS; MASS; Mathematical models; PARTICLES; RAYLEIGH-TAYLOR INSTABILITY; SHOCK TUBES; SURFACE TENSION
• ISSN: 0004-637X; 1538-4357
• DOI: 10.1088/0004-637X/768/1/44

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-density (high-density) side is overestimated (underestimated). 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 the 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 as a special case. Our formulation can be extended to handle a non-ideal gas easily.