A mass invariant in a compressible turbulent medium

Predicting the measurable statistical properties of density fluctuations in a supersonic compressible turbulent flow is a major challenge in physics. In 1951, Chandrasekhar derived an invariant under the assumption of the statistical homogeneity and isotropy of the turbulent density field and stationarity of the background density. Recently, Jaupart & Chabrier (2021) extended this invariant to non-isotropic flows in a time-evolving background and showed that it has the dimension of a mass. This invariant $M_{\rm inv}$ is defined by $M_{\rm inv} = \mathbb{E}(\rho)\text{Var}\left(\frac{\rho}{\mathbb{E}(\rho)}\right)(l_{\rm c}^\rho)^3$ where $\rho$ is the density field and $l_{\rm c}^\rho$ is the correlation length. In this article, we perform numerical simulations of homogeneous and isotropic compressible turbulence to test the validity of this invariant in a medium subject to isotropic decaying turbulence. We study several input configurations, namely different Mach numbers, injection lengths of turbulence and equations of state. We confirm that $M_{\rm inv}$ remains constant during the decaying phase of turbulence. Furthermore, we develop a theoretical model of the density field statistics which predicts without any free parameter the evolution of the correlation length with the variance of the logdensity field beyond the assumption of the gaussian field for the logdensity. Noting that $M_{\rm inv}$ is independent of the Mach number, we show that this invariant can be used to relate the non-gaussian evolution of the logdensity probability distribution function to its variance with no free parameters.