Numerical solution of phase change heat transfer problems by effective heat capacity model and element differential method

Abstract

The element differential method (EDM) is extended to solve the transient heat transfer problems with phase change. The governing equation of the phase change problem is established in the whole domain by using the effective heat capacity method. Based on the analytical expressions of spatial derivatives of the shape function with respect to the coordinate, discrete equations for internal nodes are generated by employing the governing equation. The interface nodes and the outer boundary nodes meet the flux equilibrium condition. Several numerical examples are given to validate the effectiveness and accuracy of the present method for solving phase change problems.