THE SOLUTION OF A HEAT TRANSFER PROBLEM WITH CHANGE OF STATE INSIDE A SPHERICAL SHELL

The heat transfer problem is solved for the case of cooling, below the freezing temperature, an initially liquid material inside a spherical shell. The shell is limited by a fixed inner surface and by an outer surface, free to radially expand or contract. As boundary conditions it is imposed that the inner surface is kept constant below the freezing temperature of the liquid and the outer surface is maintained constant above it. The solution represents the formation of a solidified layer that expands outward, separated from the liquid by a spherical surface kept constant at the freezing temperature. The problem is solved in the form of two closed form solutions, written in non-dimensional variables: one for the heat conduction equation in the solid layer and the second for the heat conduction - advection equation in the liquid layer. The solutions formally depend of and are linked by the time dependent radius of the spherical solid liquid interface and its time derivate, which are, at first, unknown. A differential equation describing the solid liquid interface radius as function of time is obtained imposing the conservation of energy through the interface during the heat exchange process. This equation is non-linear and has to be numerically solved. Substitution of the interface radius and its time derivate for a particular instant in the heat transfer equations solutions furnishes the temperature distribution inside the spherical shell at that moment. The solution is illustrated with numerical examples.