A numerical approach to a two-phase free boundary problem with MPC material in a finite domain

Purpose

The main purpose of this paper is given below: To present a mathematical model of a two-phase Stefan problem including a moving phase change material and variable thermophysical properties. To find a numerical solution of the problem to discuss the dependence of considered phase change problem on variable thermal conductivity, variable specific heat and Peclet number.

Design/methodology/approach

In this paper, a numerical solution of the problem is obtained using the front-fixing method in tandem with the explicit finite difference scheme. The authors have also discussed the consistency and stability of proposed numerical scheme.

Findings

In this study, it is observed that the considered scheme is an efficient tool that provides sufficiently accurate results for exploring the behaviors of moving interface (free boundary) and temperature profile for a nonclassical two-phase free boundary problem. In this study, the authors have observed that the parameters α1 and α2 influence the temperature profiles of the liquid region and the solid region. It is also found that the free boundary propagates faster when the authors increase the parameter α1 or decrease the parameter α2.

Originality/value

From the literature, it is seen that most of the two-phase problems with free boundary in an infinite domain are considered by the authors with constant thermophysical properties. Because it is possible to establish an analytical solution of two-phase problems with free boundary in case of an infinite domain. Moreover, a two-phase problem in a finite domain involving moving phase change material with the unidirectional speed is not considered. Therefore, the authors have considered a two-phase free boundary problem with variable thermal coefficients.