Mimetic discretization of two-dimensional magnetic diffusion equations

In case of non-constant resistivity, cylindrical coordinates, and highly distorted polygonal meshes, a consistent discretization of the magnetic diffusion equations requires new discretization tools based on a discrete vector and tensor calculus. We developed a new discretization method using the mimetic finite difference framework. It is second-order accurate on arbitrary polygonal meshes and a consistent calculation of the Joule heating is intrinsic within it. The second-order convergence rates in L 2 and L 1 norms were verified with numerical experiments.