Reconstruction of Voronoi diagrams in inverse potential problems

In this paper we propose and analyze a numerical method for the recovery of a piecewise constant parameter with multiple phases in the inverse potential problem. The potential is assumed to be constant in each phase, and the phases are modeled by a Voronoi diagram generated by a set of sites, which are used as control parameters. We first reformulate the inverse problem as an optimization problem with respect to the position of the sites. Combining techniques of non-smooth shape calculus and sensitivity of Voronoi diagrams, we are able to compute the gradient of the cost function, under standard non-degeneracy conditions of the diagram. We provide two different formulas for the gradient, a volumetric and an interface one, which are compared in numerical experiments. We provide several numerical experiments to investigate the dependence of the reconstruction on the problem parameters, such as noise, number of sites and initialization.