From particles to precision. Simulating subsonic turbulence with Smoothed Particle Hydrodynamics

Abstract

The numerical simulation of subsonic turbulence with smoothed particle hydrodynamics (SPH) has traditionally been hampered by zeroth-order (E0) errors, inaccurate gradient evaluations, and excessive numerical dissipation. We aim to investigate whether a modern SPH formulation can overcome these challenges by comparing its results to those obtained using moving-mesh and meshless finite-volume methods such as AREPO and GIZMO. For this purpose, we used SPH-EXA, a highly scalable, natively GPU-accelerated, state-of-the-art SPH code. Our results show that SPH-EXA accurately reproduces the Kolmogorov inertial range scaling in the subsonic regime with increasing resolution, closely matching the results of the reference methods. We also identify accurate grad-h terms as critical: a noisy standard implementation can imprint spurious granulation in the density field once dissipation is sufficiently reduced. These results demonstrate that, with appropriate methodological advances, SPH can achieve a level of fidelity in modeling subsonic turbulence comparable to the most advanced Eulerian and moving-mesh approaches.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…