Numerical approximation and analysis of the inverse Robin problem using the Kohn-Vogelius method

In this work, we numerically investigate the inverse Robin problem of recovering a piecewise constant Robin coefficient in an elliptic or parabolic problem from the Cauchy data on a part of the boundary, which arises commonly in several applications, e.g., non-destructive corrosion detection. We employ a Kohn-Vogelius type variational functional for the regularized reconstruction, and discretize the resulting optimization problem using the Galerkin finite element method on a graded mesh. We establish rigorous error estimates on the recovered Robin coefficient in terms of the mesh size, temporal step size and noise level. This is achieved by combining the approximation error of the direct problem, a priori estimates on the functional, and suitable conditional stability estimates of the continuous inverse problem. We present several numerical experiments to illustrate the approach and to complement the theoretical findings.