  In every turbulent flow with non-zero viscosity, heat is generated by viscous
friction. This heat is then mixed by the velocity field. We consider how heat
fluctuations generated this way are injected and distributed over length scales
in isotropic turbulence. A triadic closure is derived and numerically
integrated. It is shown how the heat fluctuation spectrum depends on the
Reynolds and Prandtl numbers.
