John Crepeau

We analyze the evolution of the solid-liquid front during melting and solidification in materials with constant internal heat generation, under prescribed temperature and heat flux conditions at the boundary of an infinite cylinder. We employ a sharp interface approach and assume that the motion of the front is slow relative to the temperature changes in both phases of the material. We derive infinite series solutions for the temperature in each phase and a nonlinear first-order differential equation for the evolution of the interface. Additionally, we solve the problem using the catching of the front into a node method and the Ansys Fluent enthalpy-porosity method. The latter incorporates a mushy zone that is a mixed solid-liquid transition zone. All three methods provide consistent results, especially when the mushy zone is taken into account. The series and front catching solutions develop a finite time overheated zone during melting, whereas the enthalpy solutions do not exhibit this phenomenon. We show that the evolution of the overheated and mushy zones is very similar in shape and time for both boundary conditions.

A weakly time-dependent equation for the evolution of the solid–liquid interface in spherical coordinates, driven by internal heat generation, is derived for constant surface temperature boundary conditions. The derivation comes by making an assumption that the interface moves slowly compared to the changes in the temperature so that the technique of separation of variables may be applied for Stefan numbers less than one. Under this approximation, we can separate the nonhomogeneous heat diffusion equation into transient and steady-state terms, and then integrate to get the temperature relations. With the temperature equations in hand, the derivatives are inserted into the interface equation giving a first-order differential equation for the location of the solid–liquid interface as a function of time. The results are compared to a previously derived quasi-static solution and a numerical simulation generated using the method of catching of the front. This method allows for direct tracking of a moving boundary via the calculation of the time it takes to move from node to node in a discretized grid characteristic of classical finite difference methods. Plots of the interface evolution show excellent agreement between the three methods, especially for lower Stefan numbers. The quality of the approximation decreases as the Stefan number increases, but the model is more accurate than the previously studied quasi-static model. For the Stefan numbers St = 1.0 and 10.0, the weakly time-dependent solutions are in better agreement with the numerical results than the quasi-static solutions.