A front-fixing method for the nonlinear and coupled phase change model in a moving domain

This study presents a mathematical model for a sublimation problem in a two-phase region which involves temperature and concentration-dependent thermophysical properties to investigate the temperature distribution and mass density in both regions. The numerical solution to the problem is acquired by employing a front-fixing explicit finite difference method. The consistency and stability of the numerical scheme are theoretically analyzed. The accuracy of the presented method is validated through a comparison with the exact solution achieved in a particular case and it is identified that current results are sufficiently near to them. The consequences of different dimensionless parameters on temperature and sublimation curve are presented graphically. The findings from this study provide some ways to accelerate the speed of phase transition with minimal energy absorption of particles during the sublimation process.