The bubble evolution and dissolution process in a glass furnace is an extremely complicated physico-chemical phenomenon. Most of the huge numbers of bubbles generated in molten glass during the decomposition of the glass batch are removed during the trajectory from the charging end to the forming section, but some remain in final products. On the basis of the two mechanisms, flotation and absorption, in the bubble removing process, a numerical simulator was developed to evaluate the influence of first, glass and batch composition, including refining agents, second, geometrical tank design and third, furnace operating conditions upon the bubble quality in products. Furthermore, an index was devised to synthetically assess the numbers of bubbles in final products.