Numerical analysis and implicit time stepping for high-order, fluid flow simulations on GPU architectures