A sharp numerical method for the simulation of Stefan problems with convective effects.

We present a numerical method for the solution of interfacial growth governed by the Stefan model coupled with incompressible fluid flow. An algorithm is presented which takes special care to enforce sharp interfacial conditions on the temperature, the flow velocity and pressure, and the interfacial velocity. The approach utilizes level-set methods for sharp and implicit interface tracking, hybrid finite-difference/finite-volume discretizations on adaptive quadtree grids, and a pressure-free projection method for the solution of the incompressible Navier-Stokes equations. The method is first verified with numerical convergence tests using a synthetic solution. Then, computational studies of ice formation on a cylinder in crossflow are performed and provide good quantitative agreement with existing experimental results, reproducing qualitative phenomena that have been observed in past experiments. Finally, we investigate the role of varying Reynolds and Stefan numbers on the emerging interface morphologies and provide new insights around the time evolution of local and average heat transfer at the interface. • Sharp numerical solver for the solution of Stefan problems with convective effects. • Finite-difference/finite-volume, level-set, and pressure-free projection methods. • Efficient computation using adaptive quadtree grids. • Apply sharp boundary conditions and calculate accurate gradients at the interface. • Validation/investigation of interface morphology and interfacial heat transfer. [ABSTRACT FROM AUTHOR]