  Recent progress in PDE constrained optimization on shape manifolds is based
on the Hadamard form of shape derivatives, i.e., in the form of integrals at
the boundary of the shape under investigation, as well as on intrinsic shape
metrics. From a numerical point of view, domain integral forms of shape
derivatives seem promising, which rather require an outer metric on the domain
surrounding the shape boundary. This paper tries to harmonize both points of
view by employing a Steklov-Poincar\'e type intrinsic metric, which is derived
from an outer metric. Based on this metric, efficient shape optimization
algorithms are proposed, which also reduce the analytical labor, so far
involved in the derivation of shape derivatives.
