The phase change problem in materials with internal heat generation in a cylinder

• Published in: International communications in heat and mass transfer 2024-12, Vol.159, p.108150, Article 108150
• Main Authors: Barannyk, Lyudmyla L., Crepeau, John C., Paulus, Patrick, Sakhnov, Alexey Yu, Williams, Sidney D.V.
• Format: Article
• Language: English
• Subjects: Enthalpy-porosity method; Melting and solidification; Method of catching of the front into a space grid node; Mushy zone; Overheated zone; Separation of variables; Sharp interface; Solid-liquid phase change; Stefan problem
• ISSN: 0735-1933
• DOI: 10.1016/j.icheatmasstransfer.2024.108150

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. [Display omitted] •Solid-liquid phase change with internal heat generation is analyzed with two boundary conditions.•The problem is modeled using analytical, front catching and the enthalpy-porosity methods.•The overheated regions produced by the sharp interface closely resemble the mushy zone morphology.•The solutions show good agreement using both analytical and numerical techniques.