SUMMARYWe investigate implicit large eddy simulation of the Taylor-Green vortex, Comte-Bellot-Corrsin experiment, turbulent channel flow and transitional and turbulent flow over an SD7003 airfoil using the high-order unstructured correction procedure via reconstruction (CPR) scheme, also known as the flux reconstruction scheme. We employ P1 (second-order) to P5 (sixth-order) spatial discretizations. Results show that the CPR scheme can accurately predict turbulent flows without the addition of a sub-grid scale model. Numerical dissipation, concentrated at the smallest resolved scales, is found to filter high-frequency content from the solution. In addition, the high-order schemes are found to be more accurate than the low-order schemes on a per degree of freedom basis for the canonical test cases we consider. These results motivate the further investigation and use of the CPR scheme for simulating turbulent flows.