A hybride finite element/control volume method was employed to calculate the pressure field,to trace the moving boundaries,one for melt front,the other for gas melt interface.