Improved numerical method for three dimensional diffusion equation with strongly discontinuous coefficients

The numerical heat barrier is induced when the traditional support operator method is used to solve the diffusion equation with strongly discontinuous coefficients. In this paper, a new support operator method is proposed. In this method, anisotropy conductivity tensors are introduced and defined on grid vertices. The method is applied to two diffusion cases with strongly discontinuous conductivity coefficients. The numerical results agree well with the analytical values, which proves the feasibility of the proposed scheme to solve the diffusion equation with strongly discontinuous coefficients.