On high-order conservative finite element methods

We describe and analyze a volumetric and residual-based Lagrange multipliers saddle point reformulation of the standard high-order finite method, to impose conservation of mass constraints for simulating the pressure equation on two dimensional convex polygons, with sufficiently smooth solution and mobility phase. We establish high-order a priori error estimates with locally conservative fluxes and numerical results are presented that confirm the theoretical results.