Higher-order method of lines for the numerical simulation of turbulence