I tested on 2D and 3D models and the response was the same: temperature field not symmetric. This sounds a bit strange since Hgen applied in a node should provide spherical (3D) or cylindrical (2D) distribution. This happens only for transient hgen load.