PDE-constrained shape optimization with FEM-tailored discretization of diffeomorphisms