Variational time discretization for mixed finite element approximations of nonstationary diffusion problems

We develop and study numerically two families of variational time discretization schemes for mixed finite element approximations applied to nonstationary diffusion problems. Continuous and discontinuous approximations of the time variable are encountered. The solution of the arising algebraic block system of equations by a Schur complement technique is described and an efficient preconditioner for the iterative solution process is constructed. The expected higher order rates of convergence are demonstrated in numerical experiments. Moreover, superconvergence properties are verified. Further, the efficiency and stability of the approaches are illustrated for a more sophisticated three-dimensional application of practical interest with discontinuous and anisotropic material properties.