Abstract: Using the method of line and the variable domain finite element method, the numerical simulation of 2-D., two-phase Stefan problems has been carried out successfully. The advantage of the method is that two boundary conditions on the moving boundary are simultaneously satisfied. Honce, it has a higher precision even for the coarse meshes.
A melting problem of a square ice will be as a typical calculating example. For the different parameters (the thermal diffusivity, the latent heat of melting etc.and two types of the boundary condition, the numerical results of the temperature distribution and the position of the interface have been obtained. For the degenerated one-phase case, our results are compared with Meyer's results[7]. The agreement is quite good.