The authors study the effects on heat diffusion related to the presence of a quenched disorder of impurities at the boundaries. A simple lattice model summarizes the main characteristics of the surface through which the heat flux input is produced. The disorder induced temperature fluctuations are proportional to the imposed temperature gradient and depend on the distribution of impurities. The existence of a nonlinearity in the boundary conditions leads to the renormalization of the Nusselt number