High-order Discontinuous Galerkin (DG) hydrodynamics can efficiently capture supersonic turbulence, provided the method is augmented with an appropriate artificial viscosity treatment and a projection of primitive variables onto the polynomial basis.
