We study large-eddy simulation (LES) of compressible decaying isotropic turbulence of dense gas at initial turbulent Mach numbers of 0.4 and 0.8. The unclosed subgrid scale (SGS) terms are approximated by the dynamic iterative approximate deconvolution (DIAD) model proposed by Yuan \emph{et al.} "Dynamic iterative approximate deconvolution models for large-eddy simulation of turbulence," Phys. Fluids 33, 085125 (2021), and compared with the dynamic Smagorinsky (DSM) model.In an \textit{a priori} test, the correlation coefficients of DIAD model for most SGS terms are larger than 0.98, and the relative errors are smaller than 0.2, except for the SGS internal energy flux. In an \textit{a posteriori} test, the DIAD model can well predict the probability density functions (PDFs) of SGS terms involving thermodynamic variables. Moreover, the DIAD model shows greater advantages than the DSM model in predicting various statistics and structures of compressible turbulence of dense gas, including spectra of velocity and thermodynamic variables, PDFs of SGS kinetic energy flux, deviatoric SGS stress and normalized strain-rate tensor, and the instantaneous spatial structures of vorticity.

The small-scale statistics and local flow topology of compressible homogeneous isotropic turbulence of dense gas are numerically investigated with the turbulent Mach number and Taylor Reynolds number, respectively, nearly equaling 1.0 and 153.0. The initial state of the flow field is in the inversion zone, where the fundamental derivative of gas dynamics is negative. After reaching the stationary state, the flow field includes three different gas regions: a Bethe–Zel'dovich–Thompson (BZT) region, a classical dense gas (CDG) region, and a usual gas region. In the present study, the effects of different gas regions on the statistical properties of the enstrophy production term are investigated. Based on Helmholtz decomposition, it is found that the enstrophy production mainly comes from its solenoidal component. The dense gas effect reduces the production of enstrophy in the compression region and weakens the loss of enstrophy in the expansion region. Furthermore, the properties of flow topology based on the three invariants of the velocity gradient tensor are studied. The expansion region is mainly occupied by the BZT and CDG regions. In the expansion region, the dense gas effect significantly reduces the expansive vortex structure and weakens the contribution of this structure to the enstrophy loss.

The effects of heat sources on kinetic energy transfer in compressible homogeneous shear turbulence are studied using numerical simulations at turbulent Mach numbers 0.1 and 0.4 for two levels of heat source. It is found that the strong heat source can significantly enhance both positive and negative components of subgrid-scale (SGS) kinetic energy flux and pressure–dilatation. After adding a strong heat source, compression motions enhance the positive SGS flux, and expansion motions enhance the negative SGS flux at a low turbulent Mach number. According to the Helmholtz decomposition, we found that the solenoidal and dilatational components of pressure–dilatation and SGS kinetic energy flux are increased greatly by a strong heat source at a low turbulent Mach number. The solenoidal mode plays a dominant role in the kinetic energy transfer process, but the contribution of the dilatational mode is not negligible. The dilatational component of the production term is increased by a strong heat source at a low turbulent Mach number, providing the main source of kinetic energy to the dilatational mode. The strong heat source also enhances the kinetic energy exchange between solenoidal mode and dilatational mode through nonlinear advection at a low turbulent Mach number. Moreover, the strong heat source enhances pressure anisotropy, redistribution of the kinetic energy of two transverse components, and energy transfer from internal energy to the kinetic energy through pressure–dilatation term. At a high turbulent Mach number, the strong heat source has little impact on the solenoidal and dilatational components of kinetic energy transfer terms.