“…The rDG method is a general framework for constructing efficient high-order schemes with reconstruction techniques, having the finite-volume (FV) and discontinuous Galerkin (DG) methods as special cases [3,4]. As we have shown in our previous developments [5,6,7,8,9], the two approaches can be combined systematically to simplify the discretization of diffusion/viscous terms, improve gradient accuracy, accelerate iterative convergence, and achieve higher-order accuracy with fewer numbers of degrees of freedom than DG methods. These advantages have been demonstrated for diffusion with scalar and tensor diffusion coefficients [5,6], nonlinear diffusion [7] advection-diffusion equations [8], and the NS equations [9].…”