Numerical solution of moving boundary problem in a finite domain

The temperature distribution in two regions and the location of moving inte rface during freezing in a finite domain is studied numerically. The differ ential equations governing the process of heat transfer in two regions are converted to initial value problem of vector matrix form. The solution of t his initial value problem is utilized iteratively in the interface heat flu x equation to determine interface location as well as the temperature distr ibution in two regions. The whole analysis is presented in a nondimensional form and the results thus obtained are discussed in detail.