This paper summarises the findings of a study on mathematical modeling of heat transfer in glass melts using the commonly employed diffusion or Rosseland approximation (R) for optically thick media and the more general, spectral discrete ordinates (DO) approach. The systems considered were the melting of an ultra-clear glass in a furnace (optical thickness ~ 10) and heat transfer in 2 glass melts with different levels of iron oxide flowing through a delivery channel (optical thicknesses of ~17 and ~80). The results showed that the use of the R approach for modeling heat transfer in the ultra-clear will lead to serious errors in the prediction of important practical parameters such as the refractory temperature, energy efficiency, and the quality of glass melt.