Arbitrary-order mixed methods for heterogeneous anisotropic diffusion on general meshes

We devise mixed methods for heterogeneous anisotropic diffusion problems supporting general polyhedral meshes. For a polynomial degree k≥0, we use as potential degrees of freedom the polynomials of degree at most k inside each mesh cell, whereas for the flux we use polynomials of degree at most k for the normal component on each face and fluxes of polynomials of degree at most k inside each cell. The method relies on three ideas: a flux reconstruction obtained by solving independent local problems inside each mesh cell, a discrete divergence operator with a suitable commuting property and a stabilization enjoying the same approximation properties as the flux reconstruction. Two static condensation strategies are proposed to reduce the size of the global problem, and links to existing methods are discussed. We carry out a full convergence analysis yielding flux-error estimates of order (k+1) and L²-potential estimates of order (k+2) if elliptic regularity holds. Numerical examples confirm the theoretical results.