“…The discontinuous Galerkin (DG) methods are locally conservative, stable, and high-order accurate methods which can easily handle complex geometries, irregular meshes with hanging nodes, and approximations that have polynomials of different degrees in different elements. These properties, which render them ideal to be used with ¢ ¡ -adaptive strategies, not only have brought these methods into the main stream of computational fluid dynamics, for example, in gas dynamics [11,39,14], compressible [10,64,66,65] and incompressible [13,33,32]flows, turbomachinery [12], magneto-hydrodynamics [77], granular flows [53,52], semiconductor device simulation [25,24], viscoplastic crack growth and chemical transport [21], viscoelasticity [51,5,8] and transport of contaminant in porous media [41,1,28,29], but have also prompted their application to a wide variety of problems for which they were not originally intended like, for example, Hamilton-Jacobi equations [60,59,62], second-order elliptic problems [38,67,6,70,20,22,31,4], elasticity [54,46], and Korteweg-deVries equations [73,72].…”