The affluent literature of finite element methods applied to linear parabolic problems, generally, produces numerical procedures with satisfactory properties from a numerical approximation point of view. However, some initial, boundary conditions or abrupt sink/source terms may cause large gradients at some points and consequently jumps in the solution that usually needs a certain period of tim...