Approximation of diffusion operators with discontinuous tensor coefficients on distorted meshes Original Research Article

A new finite volume method is presented for discretizing diffusion operators with discontinuous tensor coefficients. The main advantage of this method is that arbitrary distorted meshes can be used without the numerical results being altered. The matrices that need to be inverted are positive definite, so the most powerful linear solvers can be applied. Moreover the numerical experiments show that the method is second-order accurate. The method has been tested on a few elliptic and parabolic equations.