Some properties of a newly developed polynomial preserving gradient recovery technique are discussed. Both practical and theoretical issues are addressed. Bounded-ness property is considered especially under anisotropic grids. For even-order finite element space, an ultra-convergence property is established under translation invariant meshes; for linear element, a superconvergence result is proven for unstructured grids generated by the Delaunay triangulation.