“…Here, we consider approximation by means of finite elements, which has become the most popular choice for PDE-constrained shape optimization due to its flexibility for engineering applications. Nevertheless, it is worth mentioning that alternatives based on other discretizations have also been considered [6,9,18,43].…”