Phase-field modeling and numerical simulation for ice melting

In this paper, we propose a mathematical model and present numerical simulations for ice melting phenomena. The model is based on the phase-field modeling for the crystal growth. To model ice melting, we ignore anisotropy in the crystal growth model and introduce a new melting term. The numerical solution algorithm is a hybrid method which uses both the analytic and numerical solutions. We perform various computational experiments. The computational results confirm the accuracy and efficiency of the proposed method for ice melting.