SUMMARYCompact ÿnite di erence methods feature high-order accuracy with smaller stencils and easier application of boundary conditions, and have been employed as an alternative to spectral methods in direct numerical simulation and large eddy simulation of turbulence. The underpinning idea of the method is to cancel lower-order errors by treating spatial Taylor expansions implicitly. Recently, some attention has been paid to conservative compact ÿnite volume methods on staggered grid, but there is a concern about the order of accuracy after replacing cell surface integrals by average values calculated at centres of cell surfaces. Here we introduce a high-order compact ÿnite di erence method on staggered grid, without taking integration by parts. The method is implemented and assessed for an incompressible shear-driven cavity ow at Re = 10 3 , a temporally periodic ow at Re = 10 4 , and a spatially periodic ow at Re = 104 . The results demonstrate the success of the method.