Vortex-p: a Helmholtz-Hodge and Reynolds decomposition algorithm for particle-based simulations

Abstract

Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relatively simple to perform numerically for uniformly-sampled data, such as the one emerging from Eulerian, fix-grid simulations; but their computation is remarkably more complex in the case of non-uniformly sampled data, such as the one stemming from particle-based or meshless simulations. In this paper, we describe, implement and test vortex-p, a publicly available tool evolved from the vortex code, to perform both these decompositions upon the velocity fields of particle-based simulations, either from smoothed particle hydrodynamics (SPH), moving-mesh or meshless codes. The algorithm relies on the creation of an ad-hoc adaptive mesh refinement (AMR) set of grids, on which the input velocity field is represented. HHD is then addressed by means of elliptic solvers, while for the RD we adapt an iterative, multi-scale filter. We perform a series of idealised tests to assess the accuracy, convergence and scaling of the code. Finally, we present some applications of the code to various SPH and meshless finite-mass (MFM) simulations of galaxy clusters performed with OpenGadget3, with different resolutions and physics, to showcase the capabilities of the code.

Accepted for publication in Computer Physics Communications. 20 pages, 13 figures + appendices. The code is available in https://github.com/dvallesp/vortex-p and the documentation can be accessed in https://vortex-particles.github.io

References in corpus (16)

• The Statistics of Supersonic Isothermal Turbulence
• On the Star Formation Efficiency of Turbulent Magnetized Clouds
• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• The Fractal Density Structure in Supersonic Isothermal Turbulence: Solenoidal versus Compressive Energy Injection
• Resolved magnetic dynamo action in the simulated intracluster medium
• Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence
• Turbulence and Vorticity in Galaxy Clusters Generated by Structure Formation
• Hydrodynamical adaptive mesh refinement simulations of turbulent flows - II. Cosmological simulations of galaxy clusters
• The Matryoshka Run (II): Time Dependent Turbulence Statistics, Stochastic Particle Acceleration and Microphysics Impact in a Massive Galaxy Cluster
• Physical viscosity in smoothed particle hydrodynamics simulations of galaxy clusters
• Dynamical evolution of magnetic fields in the intracluster medium
• Shaping the X-ray spectrum of galaxy clusters with AGN feedback and turbulence
• Hydrodynamical adaptive mesh refinement simulations of turbulent flows - I. Substructure in a wind
• A multifiltering study of turbulence in a large sample of simulated galaxy clusters
• Troubled cosmic flows: turbulence, enstrophy and helicity from the assembly history of the intracluster medium
• Unravelling cosmic velocity flows: a Helmholtz-Hodge decomposition algorithm for cosmological simulations