SUMMARYComputation of non-isothermal ow problems involving moving interfaces is presented. A Lagrangian interface technique, deÿned in the context of a ÿxed-mesh ÿnite element formulation for incompressible ows, is employed to update the interface position. A global mass-corrector algorithm is used to accurately enforce the global mass conservation. The Navier-Stokes equations are solved with an improved sub-element integration technique to more accurately account for sudden changes in the uid properties across the interface. The method described is applied to two mould ÿlling problems.