RESIDUAL-BASED A POSTERIORI ERROR ESTIMATES FOR hp-DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE BIHARMONIC PROBLEM

We introduce a residual-based a posteriori error estimator for a novel hp-version interior penalty discontinuous Galerkin method for the biharmonic problem in two and three dimensions. We prove that the error estimate provides an upper bound and a local lower bound on the error, and that the lower bound is robust to the local mesh size but not the local polynomial degree. The suboptimality in terms of the polynomial degree is fully explicit and grows at most algebraically. Our analysis does not require the existence of a C-1-conforming piecewise polynomial space and is instead based on an elliptic reconstruction of the discrete solution to the H-2 space and a generalized Helmholtz decomposition of the error. This is the first hp-version error estimator for the biharmonic problem in two and three dimensions. The practical behavior of the estimator is investigated through numerical examples in two and three dimensions.