Liquid-solid phase transition is accompanied by a latent heat release, both in isothermal and non-isothermal phase transformations. The latent heat and the discontinuous phase change function increase the difficulty of obtaining a solution for the Fourier heat conduction equation. Celentano et al. [Int. J. Numer. Meth. Eng. 37(20), 1994] proposed a temperature-based finite element model for solving multidimensional transient heat conduction involving phase change. The present work addresses the computational aspects of the Celentano et al. model. The importance of a line search algorithm for improving the convergence of the Newton-Raphson iterations are explained in detail. While performing the iterations in this kind of fixed domain methods, the phase front moves back and forth fictitiously. The introduced phase change matrix handles the latent effect efficiently. The phase fractions are evaluated at the integration points instead of the nodal points. Several numerical examples are presented and the benefits and difficulties of the solution technique are elaborately discussed.
- Heat Transfer Division
Finite Element Technique for Phase-Change Heat Conduction Problem
- Views Icon Views
- Share Icon Share
- Search Site
Nallathambi, AK, Specht, E, & Bertram, A. "Finite Element Technique for Phase-Change Heat Conduction Problem." Proceedings of the ASME 2009 Heat Transfer Summer Conference collocated with the InterPACK09 and 3rd Energy Sustainability Conferences. Volume 2: Theory and Fundamental Research; Aerospace Heat Transfer; Gas Turbine Heat Transfer; Computational Heat Transfer. San Francisco, California, USA. July 19–23, 2009. pp. 667-674. ASME. https://doi.org/10.1115/HT2009-88106
Download citation file: