Thermomechanical modeling of glass cooling including thermal radiation using P 1 method