A posteriori discontinuous finite element error estimation for two-dimensional hyperbolic problems Original Research Article

We analyze the discontinuous finite element errors associated with p-degree solutions for two-dimensional first-order hyperbolic problems. We show that the error on each element can be split into a dominant and less dominant component and that the leading part is O(hp+1) and is spanned by two (p+1)-degree Radau polynomials in the x and y directions, respectively. We show that the p-degree discontinuous finite element solution is superconvergent at Radau points obtained as a tensor product of the roots of (p+1)-degree Radau polynomial. For a linear model problem, the p-degree discontinuous Galerkin solution flux exhibits a strong O(h2p+2) local superconvergence on average at the element outflow boundary. We further establish an O(h2p+1) global superconvergence for the solution flux at the outflow boundary of the domain. These results are used to construct simple, efficient and asymptotically correct a posteriori finite element error estimates for multi-dimensional first-order hyperbolic problems in regions where solutions are smooth.