High-order discontinuous finite element methods are becoming increasingly more popular for simulations of vortex dominated flows over complex geometries because they are inherently less dissipative than traditional second-order finite volume methods. In particular, the Flux Reconstruction (FR) approach has gained popularity because it not only offers less dissipation but is also eminently parallelizable on multi-core processors and accelerators. An extensive amount of research has been performed for accelerated explicit methods for FR but more realistic simulations often require high aspect ratio meshes where the maximum stable time step or CFL condition is determined by the smallest cell. This can lead to significant limitations in explicit time integration.