Abstract. In recent years the increasing attention to high-order Finite Volume (FV), Finite Element (FE) and spectral methods and the growth of computing power promote the development of high-order temporal schemes to perform robust, accurate and efficient unsteady long-time simulations. In this context, some features of the Discontinuous Galerkin finite element (DG) methods, e.g. compactness and flexibility, can be advantageous both for explicit and implicit time integration approaches. Explicit schemes can achieve very high accuracy, but are limited by time-step restrictions, while implicit schemes, even if memory consuming, can overcome time-step limitations, thus improving the time integration efficiency. During last decades several high order implicit temporal schemes have been developed, and some of them have been successfully coupled with DG methods. However these schemes can show the order reduction if applied to very stiff problems or problems with time-dependent boundary conditions. To overcome these limitations, high-order linearly implicit two-step peer methods have been proposed and successfully applied to the numerical solution of differential-algebraic equations. The aim of the present work is to implement high-order two-step peer methods in a DG code and assess their performance for the unsteady solution of the incompressible Navier-Stokes (INS) and Reynolds Averaged Navier-Stokes equations (RANS) equations.