UZH-Logo

Maintenance Infos

Resolving mixing in smoothed particle hydrodynamics


Read, J I; Hayfield, T; Agertz, O (2010). Resolving mixing in smoothed particle hydrodynamics. Monthly Notices of the Royal Astronomical Society, 405(3):1513-1530.

Abstract

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.

Abstract

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.

Citations

79 citations in Web of Science®
69 citations in Scopus®
Google Scholar™

Altmetrics

Downloads

73 downloads since deposited on 02 Mar 2011
21 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Institute for Computational Science
Dewey Decimal Classification:530 Physics
Language:English
Date:July 2010
Deposited On:02 Mar 2011 15:03
Last Modified:05 Apr 2016 14:34
Publisher:Wiley-Blackwell
ISSN:0035-8711
Additional Information:The definitive version is available at www.blackwell-synergy.com
Publisher DOI:https://doi.org/10.1111/j.1365-2966.2010.16577.x
Related URLs:http://arxiv.org/abs/0906.0774

Download

[img]
Preview
Content: Accepted Version
Filetype: PDF (Accepted manuscript, Version 2)
Size: 937kB
View at publisher
[img]
Preview
Content: Accepted Version
Filetype: PDF (Accepted manuscript, Version 1)
Size: 879kB

TrendTerms

TrendTerms displays relevant terms of the abstract of this publication and related documents on a map. The terms and their relations were extracted from ZORA using word statistics. Their timelines are taken from ZORA as well. The bubble size of a term is proportional to the number of documents where the term occurs. Red, orange, yellow and green colors are used for terms that occur in the current document; red indicates high interlinkedness of a term with other terms, orange, yellow and green decreasing interlinkedness. Blue is used for terms that have a relation with the terms in this document, but occur in other documents.
You can navigate and zoom the map. Mouse-hovering a term displays its timeline, clicking it yields the associated documents.

Author Collaborations