EnglishAbstractWe consider PDE constrained shape optimization in the framework of finite element discretization of the underlying boundary value problem. We present an algorithm tailored to preserve and exploit the approximation properties of the finite element method, and that allows for arbitrarily high resolution of shapes. It employs (i) B-spline based representations of the deformation diffeomorphism, and (ii) superconvergent domain integral expressions for the shape gradient. We provide numerical evidence of the performance of this method both on prototypical well-posed and ill-posed shape optimization problems.