Fixed-grid sharp-interface numerical solutions to the three-phase spherical Stefan problem

TL;DR

This paper develops a fixed-grid sharp-interface numerical method to solve the complex three-phase Stefan problem in spherical coordinates, capturing simultaneous melting, boiling, and condensation in finite-sized particles.

Contribution

It extends previous work by addressing the three-phase Stefan problem in spherical coordinates with a novel numerical approach and validates the method against existing two-phase results.

Findings

Numerical solutions accurately predict interface positions and velocities.

Kinetic energy effects are significant in nano-sized particles.

The method demonstrates diminishing kinetic effects for larger particles.

Abstract

Many metal manufacturing processes involve phase change phenomena, which include melting, boiling, and vaporization. These phenomena often occur concurrently. A prototypical 1D model for understanding the phase change phenomena is the Stefan problem. There is a large body of literature discussing the analytical solution to the two-phase Stefan problem that describes only the melting or boiling of phase change materials (PCMs) with one moving interface. Density-change effects that induce additional fluid flow during phase change are generally neglected in the literature to simplify the math of the Stefan problem. In our recent work [1], we provide analytical and numerical solutions to the three-phase Stefan problem with simultaneous occurrences of melting, solidification, boiling, and condensation in Cartesian coordinates. Our current work builds on our previous work to solve a more challenging problem: the three-phase Stefan problem in spherical coordinates for finite-sized particles. There are three moving interfaces in this system: the melt front, the boiling front, and the outer boundary which is in contact with the atmosphere. Although an analytical solution could not be found for this problem, we solved the governing equations using a fixed-grid sharp-interface method with second-order spatio-temporal accuracy. Using a small-time analytical solution, we predict a reasonably accurate estimate of temperature (in the three phases) and interface positions and velocities at the start of the simulation. Our numerical method is validated by reproducing the two-phase nanoparticle melting results of Font et al. [2]. Lastly, we solve the three-phase Stefan problems numerically to demonstrate the importance of kinetic energy terms during phase change of smaller (nano) particles. In contrast, these effects diminish for large particles (microns and larger).