In glass container manufacturing (e.g., production of glass bottles and jars) an important process step is the blowing of the final product. This process is fast and is characterized by large deformations and the interaction of a hot glass fluid that gets into contact with a colder metal, the mould. The objective of this paper is to create a robust finite-element model to be used for industrial purposes that accurately captures the blowing step of glass containers. The model should be able to correctly represent the flow of glass and the energy exchange during the process. For tracking the geometry of the deforming inner and outer interface of glass, level set technique is applied on structured and unstructured fixed mesh. At each time step the coupled problem of flow and energy exchange is solved by the model. Here the flow problem is only solved for the domain enclosed by the mould, whereas in the energy calculations, the mould domain is also taken into account in the computations. For all the calculations the material parameters (like viscosity) are based on the glass position, i.e., the position of the level sets. The velocity distribution, as found from this solution procedure, is then used to update the two level sets by means of solving a convection equation. A reinitialization algorithm is applied after each time step in order to let the level sets reattain the property of being a signed distance function. The model is validated by several examples focusing on both the overall behavior (such as conservation of mass and energy) and the local behavior of the flow (such as glass-mould contact) and temperature distributions for different mesh size, time step, level set settings and material parameters.