When u is an order 1 H1 function, its gradient is piecewise constant, but not continuous, so you get different gradients from different elements. I think that explains the picture from your first post.
As for visualizing the rhs, did you just copy the vector into a gridfunction? That does not work, the right hand side is a dual vector, not a FEM function.