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

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.