Approximate shape gradients for interface problems

Shape gradients of shape differentiable shape functionals constrained to an interface problem (IP) can be formulated in two equivalent ways. Both formulations rely on the solution of two IPs, and their equivalence breaks down when these IPs are solved approximatively. We establish which expression for the shape gradient offers better accuracy for approximations by means of finite elements. Great effort is devoted to provide numerical evidence of the theoretical considerations.