Polynomial preserving recovery for PHT-splines

We propose a polynomial preserving recovery method for PHT-splines within isogeometric analysis to obtain more accurate gradient approximations. The method fully exploits the local interpolation properties of PHT-splines and avoids the need for information on gradient superconvergent points. By leveraging the superconvergence argument of difference quotients and the interior error estimate, we establish the superconvergence property of the recovered gradient on translation invariant meshes. As a byproduct, a recovery-based a posteriori error estimator is developed for adaptive refinement. Numerical results confirm the theoretical findings and demonstrate the effectiveness of the proposed method.