On p-Robust Saturation for hp-AFEM
  We consider the standard adaptive finite element loop SOLVE, ESTIMATE, MARK, REFINE, with ESTIMATE being implemented using the $p$-robust equilibrated flux estimator, and MARK being D\"orfler marking. As a refinement strategy we employ $p$-refinement. We investigate the question by which amount the local polynomial degree on any marked patch has to be increase in order to achieve a $p$-independent error reduction. The resulting adaptive method can be turned into an instance optimal $hp$-adaptive method by the addition of a coarsening routine. 
Claudio Canuto and Ricardo H. Nochetto and Rob Stevenson and Marco   Verani
math.NA