Validation of large eddy simulation method for study of flatness and skewness of decaying compressible magnetohydrodynamic turbulence

In the present work we study potential applicability of large eddy simulation (LES) method for prediction of flatness and skewness of compressible magnetohydrodynamic (MHD) turbulence. The knowledge of these quantities characterizes non-Gaussian properties of turbulence and can be used for verification of hypothesis on Gaussianity for the turbulent flow under consideration. Prediction accuracy of these quantities by means of LES method directly determines efficiency of reconstruction of probability density function (PDF) that depends on used subgrid-scale (SGS) parameterizations. Applicability of LES approach for studying of PDF properties of turbulent compressible magnetic fluid flow is investigated and potential feasibilities of five SGS parameterizations by means of comparison with direct numerical simulation results are explored. The skewness and the flatness of the velocity and the magnetic field components under various hydrodynamic Reynolds numbers, sonic Mach numbers, and magnetic Reynolds numbers are studied. It is shown that various SGS closures demonstrate the best results depending on change of similarity numbers of turbulent MHD flow. The case without any subgrid modeling yields sufficiently good results as well. This indicates that the energy pile-up at the small scales that is characteristic for the model without any subgrid closure, does not significantly influence on determination of PDF. It is shown that, among the subgrid models, the best results for studying of the flatness and the skewness of velocity and magnetic field components are demonstrated by the Smagorinsky model for MHD turbulence and the model based on cross-helicity for MHD case. It is visible from the numerical results that the influence of a choice subgrid parametrization for the flatness and the skewness of velocity is more essential than for the same characteristics of magnetic field.