On the accuracy of high-order discretizations for underresolved turbulence simulations