Implicit large eddy simulation (ILES) of compressible turbulence with shock capturing schemes requires wide investigations and numerical experiments. In this study, a newly proposed P n T m− BVD (polynomial of n-degree and THINC function of m-level reconstruction based on BVD algorithm) shock capturing scheme is introduced to simulate compressible turbulence flow with ILES. The new scheme is designed by employing high-order linear-weight polynomials and THINC (Tangent of Hyperbola for INterface Capturing) functions with adaptive steepness as the reconstruction candidates. The final reconstruction function in each cell is determined with a multi-stage BVD (Boundary Variation Diminishing) algorithm so as to effectively control numerical oscillation and dissipation. Numerical tests involving shock waves and broadband turbulence are conducted in comparison with WENO (Weighted Essentially Non-oscillatory) schemes which are widely used in ILES. The results demonstrate performing ILES with P n T m− BVD scheme is able to obtain higher resolution and more faithful results than WENO does. Importantly, the superiority of P n T m− BVD becomes more notable in high wave-number region. Thus this paper provides and verifies a new scheme which is promising in providing high-resolution results for real-case ILES of compressible turbulence flow.