A posteriori error estimates for mixed finite element approximations of parabolic problems

We derive residual based a posteriori error estimates for parabolic problems on mixed form solved using Raviart–Thomas–Nedelec finite elements in space and backward Euler in time. The error norm considered is the flux part of the energy, i.e. weighted L2(Ω) norm integrated over time. In order to get an optimal order bound, an elementwise computable post-processed approximation of the scalar variable needs to be used. This is a common technique used for elliptic problems. The final bound consists of terms, capturing the spatial discretization error and the time discretization error and can be used to drive an adaptive algorithm.