Improved and reproducible heterodyne mixing (noise temperatures of 950 K at 2.5 THz) has been realized with NbN based hot-electron superconducting devices with low contact resistances. A distributed temperature numerical model of the NbN bridge, based on a local electron and a phonon temperature, has been used to understand the physical conditions during the mixing process. We find that the mixing is predominantly due to the exponential rise of the local resistivity as a function of electron temperature.