Your browser has javascript turned off or blocked. This will lead to some parts of our website to not work properly or at all. Turn on javascript for best performance.

The browser you are using is not supported by this website. All versions of Internet Explorer are no longer supported, either by us or Microsoft (read more here: https://www.microsoft.com/en-us/microsoft-365/windows/end-of-ie-support).

Please use a modern browser to fully experience our website, such as the newest versions of Edge, Chrome, Firefox or Safari etc.

Oscar Agertz. Profile photo.

Oscar Agertz

Assistant Professor / Associate senior university lecturer / Wallenberg Academy Fellow

Oscar Agertz. Profile photo.

Resolving mixing in smoothed particle hydrodynamics

Author

  • J. I. Read
  • T. Hayfield
  • O. Agertz

Summary, in English

Standard formulations of smoothed particle hydrodynamics (SPH) are unable to resolve mixing at fluid boundaries. We use an error and stability analysis of the generalized SPH equations of motion to prove that this is due to two distinct problems. The first is a leading order error in the momentum equation. This should decrease with an increasing neighbour number, but does not because numerical instabilities cause the kernel to be irregularly sampled. We identify two important instabilities: the clumping instability and the banding instability, and we show that both are cured by a suitable choice of kernel. The second problem is the local mixing instability (LMI). This occurs as particles attempt to mix on the kernel scale, but are unable to due to entropy conservation. The result is a pressure discontinuity at boundaries that pushes fluids of different entropies apart. We cure the LMI by using a weighted density estimate that ensures that pressures are single-valued throughout the flow. This also gives a better volume estimate for the particles, reducing errors in the continuity and momentum equations. We demonstrate mixing in our new optimized smoothed particle hydrodynamics (OSPH) scheme using a Kelvin-Helmholtz instability (KHI) test with a density contrast of 1:2, and the 'blob test'- a 1:10 density ratio gas sphere in a wind tunnel - finding excellent agreement between OSPH and Eulerian codes.

Publishing year

2010-07-01

Language

English

Pages

1513-1530

Publication/Series

Monthly Notices of the Royal Astronomical Society

Volume

405

Issue

3

Document type

Journal article

Publisher

Oxford University Press

Keywords

  • Hydrodynamics
  • Instabilities
  • Methods: numerical

Status

Published

ISBN/ISSN/Other

  • ISSN: 0035-8711