AbstractWe consider the corrosion detection problem in terms of the Laplace equation and study a simply connected bounded domain with Wentzell-type GIBC boundary condition. We derive the systems of integral equations and establish the equivalence to the inverse shape problem in a Sobolev space setting. For the direct problem, we use potential theory to simulate the Neumann data from Dirichlet data on Dirichlet boundary. Then we propose a Newton iterative approach based on the boundary integral equations derived from Green’s representation theorem. After describing the linearization and the iteration scheme for the inverse shape, we compute the Fréchet derivatives with respect to the unknowns. We conclude by presenting several numerical examples for shape reconstructions to show the validity of the proposed method.