Numerical Analysis of an Adaptive FEM for Distributed Flux Reconstruction