ON THE STEFAN PROBLEM WITH INTERNAL HEAT GENERATION AND PRESCRIBED HEAT FLUX CONDITIONS AT THE BOUNDARY

We study the evolution of the solid-liquid interface during melting and solidification of a material with constant internal heat generation and prescribed heat flux at the boundary for a plane wall and a cylinder. The equations are solved by splitting them into transient and steady-state components and then using separation of variables. This results in an ordinary differential equation for the interface that involves infinite series. The initial value problem is solved numerically, and solutions are compared to the previously published quasi-static solutions. We show that when the internal heat generation and the heatflux at the boundary are close in value to each other, the motion of the phase change front takes longer to reach steady-state than when the values are farther apart. As the difference between the internal heat generation and the heat flux increases, the transient solutions become more dominant and the numerical solution of the phase change front does not reach steady-state before the outer boundary or centerline is reached. The difference between the internal heat generation and the heatflux at the boundary can be used to control the motion and speed of the interface. The problem has applications for a nuclear fuel rod during meltdown.