A FIXED GRID METHOD FOR THE NUMERICAL-SOLUTION OF PHASE-CHANGE PROBLEMS

A fixed grid method using an updated iterative implicit scheme is developed to solve one-dimensional phase change problems. The temperature field is deduced from the resolution of the governing equations whose discretization takes into account the discontinuous variation of the temperature derivative at the phase change front. At each iteration an updated position of the moving front is found from the resolution of the energy conservation at the solid-liquid interface. The accuracy of the proposed numerical method has been checked on three test problems.