A parallel algorithm and code MVFT (multi-viscous-fluid and turbulence) of large-eddy simulation (LES) is developed from our MVPPM (multi-viscous-fluid piecewise parabolic method), and performed to solve the multi compressible Navier-Stokes (N-S) equations. The effect of the unresolved subgrid-scale (SGS) motions on the large scales is represented by different SGS stress models in LES. A Richtmyer-Meshkov instability experiment of the evolution of a rectangular block of SF 6 , which occupies half of the height of the shock tube test section, following the interaction with a planar shock wave, is numerically and exhaustively simulated by this code. The comparison between experimental and simulated images of the evolving SF 6 block shows that they are consistent. The numerical simulations reproduce the complex developing process of SF 6 block, which grows overturningly. The geometric quantities that characterize the extents of SF 6 block are also compared in detail between numerical simulations and experiment with good agreements between them, a quantitative demonstration of the developing law of SF 6 block. There is an evident discrepancy between the three numerical simulations for the maximum position of the right edge of block at the late stage, because the right interface grows complicated and the dissipation is different for different SGS models. The SGS turbulent dissipation, molecular viscosity dissipation and SGS turbulent kinetic energy have been studied and analyzed. They have a similar distribution to the large eddy structures. The SGS turbulent dissipation is much greater than the molecular viscosity dissipation; the SGS turbulent dissipation of Vreman model is smaller than the Smagorinsky model. In general, the simulated results of Vreman SGS model are better compared with the dynamic viscosity and Smagorinsky SGS model. The vorticity and circulation deposition on the block interface have also been investigated.large-eddy simulation, Navier-Stokes equations, subgrid-scale, Richtmyer-Meshkov instability, turbulent dissipation, turbulent kinetic energy, circulation deposition PACS: 47.10.ad, 47.20.Ma, 47.27.ep The instability of an accelerated interface separating two fluids with different properties is a fundamental hydrodynamic problem. The Richtmyer-Meshkov instability (RMI) is driven by impulsive acceleration such as a shock wave [1]. This instability is typically investigated in shock tube experiments. Any perturbations initially presented on the interface will, in most cases, be amplified following the passage of the shock. The physical mechanism for the amplification of perturbations on the interface is the baroclinic vorticity production as a result of the misalignment of the pressure gradient (∇p) of the shock and the local density gradient (∇ρ) at the interface. Nonlinearities come into play when the interface between two fluids becomes more distorted, and secondary baroclinic or shear-driven instabilities [2], such as the Kelvin-Helmholtz instability (KHI), also start developing. This combi...