In this paper, we consider the conductivity problem with piecewise‐constant conductivity and Robin‐type boundary condition on the interface of discontinuity. When the quantity of interest is the jump of the conductivity, we perform a local stability estimate for a parameterized non‐monotone family of domains. We give also a quantitative stability result of local optimal solution with respect to a perturbation of the Robin parameter. In order to find an optimal solution, we propose a Kohn–Vogelius‐type cost functional over a class of admissible domains subject to two boundary values problems. The analysis of the stability involves the computation of first‐order and second‐order shape derivative of the proposed cost functional, which is performed rigorously by means of shape‐Lagrangian formulation without using the shape sensitivity of the states variables.
