Do you have a minus sign to much in your RHSheatspatial? When I remove one of them, i.e., du_dt = c**2 * dd_u, I get the following solution:
I am not sure if this is the correct solution to your heat equation, though, but it looks more plausible, I believe.