In this paper, an inverse problem of glass forming process is studied to determine a number of unknown heat transfer coefficients which are imposed as boundary conditions. An analysis program for transient heat conduction of axi-symmetric dimension is developed to simulate the forming and cooling process. The analysis is repeated until it attains periodic state, which requires at least 30 cycles of iteration. Measurements are made for the temperatures at several available time and positions of glass and moulds in operation. Heat removal by the cooling water from the plunger is also recorded. An optimization problem is formulated to determine heat transfer coefficients which minimize the difference between the measured data and analysis results. Significant time savings are achieved in finite difference based sensitivity computation during the optimization by employing distributed computing technique. The analysis results by the optimum heat transfer coefficients are found to agree well with the measured data.