Polynomial preserving recovery on boundary

In this paper, we propose two systematic strategies to recover the gradient on the boundary of a domain. The recovered gradient has comparable superconvergent property on the boundary as that in the interior of the domain. This superconvergence property has been validated by several numerical experiments.