Discretization errors in finite element models of hyperbolic cooling