Melting glass furnace mathematical modelling currently includes couple
d simulations of electrical boosting, bubbling, flame radiation, etc..
. Increasing modelling complexity may also result from consideration o
f unsteady solutions. Different sources of thermoconvective unsteadine
ss can be considered: pull modification, time variation of bubbling an
d/or boosting, shift of the temperature profile or, finally, colour ch
ange [1], this paper will present the numerical solutions to the last
unsteady coupled mixing problem.