u[x_,y_,t_]=Exp[-((m^2)+(n^2)) k t] Sin[m x] Cos[n y]

FullSimplify[D[u[x,y,t],t]-k (D[D[u[x,y,t],x],x]+D[D[u[x,y,t],y],y])]