Flux recovery scheme for elliptic interface problems