Richardson extrapolations of Galerkin finite element methods  for parabolic partial differential equations