A two-stage method for reconstruction of parameters in diffusion equations

Parameter reconstruction for diffusion equations has a wide range of applications. In this paper, we proposed a two-stage scheme to efficiently solve conductivity reconstruction problems for steady-state diffusion equations with solution data measured inside the domain. The first stage is based on total variation regularization of the log diffusivity and the split Bregman iteration method. In the second stage, we apply the K-means clustering for the reconstruction of ``blocky'' conductivity functions. The convergence of the scheme is theoretically proved and extensive numerical examples are shown to demonstrate the performance of the scheme.