We consider PDE constrained shape optimization in the framework of finite element discretization of the underlying boundary value problem. Our approach employs (i) B-spline based representations of the deformation diffeomorphism, and (ii) superconvergent domain integral expressions for the shape gradient. We study the balance of the discretization errors of the finite element method and the B-s...