A new first-order hyperbolic system (FOHS) is formulated for the compressible Navier-Stokes equations. The resulting hyperbolic Navier-Stokes system (HNS), termed HNS20G in this paper, introduces the gradients of density, velocity, and temperature as auxiliary variables. Efficient, accurate, compact and robust reconstructed discontinuous Galerkin (rDG) methods are developed for solving this new HNS system. The newly introduced variables are recycled to obtain the gradients of the primary variables. The gradients of these gradient variables are reconstructed based on a newly developed variational formulation in order to obtain a higher order polynomial solution for these primary variables without increasing the number of degrees of freedom. The implicit backward Euler method is used to integrate solution in time for steady flow problems, while the third-order explicit first stage singly diagonally Runge-Kutta (ESDIRK) time marching method is implemented for advancing solutions in time for unsteady flows. The flux Jacobian matrices are obtained with an automatic differentiation toolkit TAPENADE. The approximate system of linear equations is solved with either symmetric Gauss-Seidel (SGS) method or general minimum residual (GMRES) algorithm with a lower-upper symmetric Gauss-Seidel (LU-SGS) preconditioner. A number of test cases are presented to assess accuracy and performance of the newly developed HNS+rDG methods for both steady and unsteady compressible viscous flows. Numerical experiments demonstrate that the developed HNS+rDG methods are able to achieve the designed order of accuracy for both primary variables and the their gradients, and provide an attractive and viable alternative for solving the compressible Navier-Stokes equations.