Nota Bene: The present document constitutes a 'postprint' version of the published paper, in which a few errors in the proofs (on the present pages 2, 5, 15 and 17) have been corrected and an incomplete reference on p. 22 has been completed. These minor corrections are marked out in red. Results are unchanged.
AbstractA modified Green operator is proposed as an improvement of Fourier-based numerical schemes commonly used for computing the electrical or thermal response of heterogeneous media. Contrary to other methods, the number of iterations necessary to achieve convergence tends to a finite value when the contrast of properties between the phases becomes infinite. Furthermore, it is shown that the method produces much more accurate local fields inside highly conducting and quasi-insulating phases, as well as in the vicinity of phase boundaries. These good properties stem from the discretization of Green's function, which is consistent with the pixel grid while retaining the local nature of the operator that acts on the polarization field. Finally, a fast implementation of the 'direct scheme' of Moulinec et al. (1994) that allows for parsimonious memory use is proposed.