A simple yet accurate iterative method for solving a one-dimensional phase change problem with convection boundary is described. The one-dimensional model takes into account the variation in the wall temperature along the direction of the flow as well as the sensible heat during preheating/precooling of the phase change material (PCM). The mathematical derivation of convective boundary conditions has been integrated into a phase change processor (PCP) algorithm that solves the liquid fraction and temperature of the nodes. The algorithm is based on the heat balance at each node as it undergoes heating or cooling which inevitably involves phase change. The paper presents the model and its experimental validation.
Halawa, E., Saman, W., & Bruno, F. (2010). A phase change processor method for solving a one-dimensional phase change problem with convection boundary. Renewable Energy, 35(8), 1688-1695. https://doi.org/10.1016/j.renene.2010.01.016