Process parameters of Precision Glass Molding (PGM) are often sought by Finite Element (FE) simulation. Mechanical as well as thermal boundary conditions (BCs) are necessary for FE simulation in which mechanical BCs are usually known or easily determinable. However, most of the thermal BCs are generally assumed in the FE simulation as they cannot be measured directly. The focus of this article is to propose a novel method for evaluating the thermal BC of glass-N2 gas. FE simulations as well as thermal cycling experiments are carried out for a glass disk specimen for three different cooling rates. CFD analysis of N2 flow in the PGM machine is performed to understand the heat extraction mechanism. Based on this, adhoc values for equivalent heat transfer coefficient (heqv) are obtained by lumped system analysis. A novel methodology is then proposed for obtaining accurate heqv values by measurement of integrated residual birefringence in glass using digital photoelasticity. FE simulation is repeated for different values of heqv until the integrated birefringence based on simulation matches with that of the experiment. For the same cooling rates, two aspherical glass lenses are molded and their residual birefringence is measured and compared with the glass disk specimen. © 2015 The American Ceramic Society.