A heat polynomials method for the two-phase inverse Stefan problem

Abstract

In this paper, we extend the heat polynomials method (HPM) proposed by the authors for one-dimensional one-phase inverse Stefan problem to the two-phase case. The solution of the problem is presented in the form of linear combination of heat polynomials. The coefficients of this combination can be determined by satisfying the initial and boundary conditions or by the least square method for the boundary of a domain. The inverse problem is ill-posed, therefore, the regularization will be taken into account. Our numerical results are compared with results obtained by another method and show good enough accuracy.