
arXiv: 1206.5006
Various formulations of smooth-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 introduce artificial correction terms or violate 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 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 'Grad-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 & 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 sub-sonic 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 'flavors' of SPH.
16 pages, 16 figures, MNRAS. Matches accepted version (tests added, typo fixed). A version of the public GADGET-2 modified to use the 'pressure-entropy' formulation proposed here can be downloaded at http://ascl.net/1305.006
Cosmology and Nongalactic Astrophysics (astro-ph.CO), turbulence, Fluid Dynamics (physics.flu-dyn), FOS: Physical sciences, Physics - Fluid Dynamics, Computational Physics (physics.comp-ph), 530, methods: numerical, 510, instabilities, cosmology: theory, hydrodynamics, Astrophysics - Instrumentation and Methods for Astrophysics, Physics - Computational Physics, Instrumentation and Methods for Astrophysics (astro-ph.IM), Astrophysics - Cosmology and Nongalactic Astrophysics
Cosmology and Nongalactic Astrophysics (astro-ph.CO), turbulence, Fluid Dynamics (physics.flu-dyn), FOS: Physical sciences, Physics - Fluid Dynamics, Computational Physics (physics.comp-ph), 530, methods: numerical, 510, instabilities, cosmology: theory, hydrodynamics, Astrophysics - Instrumentation and Methods for Astrophysics, Physics - Computational Physics, Instrumentation and Methods for Astrophysics (astro-ph.IM), Astrophysics - Cosmology and Nongalactic Astrophysics
| selected citations These citations are derived from selected sources. This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | 387 | |
| popularity This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network. | Top 1% | |
| influence This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | Top 1% | |
| impulse This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network. | Top 1% |
