The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work, by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows. The turbulent statistics agree well with those from the existing direct numerical simulation (DNS) results, indicating that the LES method established in the present work is reliable. The analysis of the turbulent fluctuations computed by the present LES reveals that the flows considered in this work follow the Morkovin's hypothesis. Thus, the compressibility effects are dominated by the mean field properties, and the relevant statistical ratios are invariant to the variation of Mach number. The near-wall streamwise streaks are more coherent and the spacing between streaks is wider as the Mach number increases. This can be regarded as a direct feature characterized by the compressibility effects. The restrained influences of compressibility effects on the production and dissipation of the turbulence kinetic energy are also identified based on the present LES results. compressibility, channel flow, large eddy simulation, subgrid-scale model PACS number(s): 98.70.Sa, 90.50.sb, 13.85.Tp Citation: Wang S Z, Gao Z X, Lee C H. Numerical investigation of compressibility effects in turbulent channel flows using large eddy simulation.Compressible turbulence has been one of the areas extracting much of the attention in turbulence research in these past decades, due to its crucial applications in hypersonic technology. Effects of compressibility can be categorized into two types [1], namely, those associated with variations of the mean properties, and those due to fluctuations of thermodynamic quantities (so-called acoustic effects). According to Morkovin's hypothesis, compressibility effects influencing the mean flow properties become dominant when pressure fluctuation is much less than its mean value [2]. Bradshaw [3] and Spina et al. [4] found that only the mean effects are significant for transonic and supersonic wall-bounded flows when Ma<5.Along with the continuous improvement on computational fluid dynamics and the rapid development of computer technology, numerical simulation has become an important tool for turbulence research. There are currently three approaches available for turbulence numerical simulations, i.e., the direct numerical simulation (DNS), the large eddy simulation (LES) and the Reynolds averaged NavierStokes (RANS). For LES, the turbulent instantaneous movement is decomposed into large scales and small scales through a filter. The quantities of large scales are directly computed from the equations of motion, while the effects of unresolved small scales of motion on large scales are imposed by subgrid-scale (SGS) terms [5]. Compared with DNS on the one hand, LES can save considerable computational efforts, and be adaptable to high Reynolds number flows. By comparing it with RANS, on the other hand, the SGS model can be universal and provide resolution of the fluctuations of large scale...