“…Because it has many advantages, such as high order accuracy, flexibility for mesh refinement, localizability, stability, parallelizability and less numerical diffusion/dipersion, the DG method has been widely extended and used to solve many partial differential equations. For solving the nonlinear systems of the general conservation laws, the total variation bounded (TVB) RungeKutta local projection discontinuous Galerkin (RKDG) method with high order of accuracy is developed [4−9] Another popular DG method is the local discontinuous Galerkin method, see [10][11][12][13][14][15][16][17] and references therein. Additionally, many other DG methods, such as generalized DG (GDG) method [18−19] , coupled DG method [20] , DG finite volume element method [21] , EulerianLagrangian DG method [22] , characteristic DG method [23] and compact DG method [24] , have been developed.…”