Abstract
The freezing phenomena in supercooled liquid droplets are important for many engineering applications. For instance, a theoretical model of this phenomenon can offer insights for tailoring surface coatings and for achieving icephobicity to reduce ice adhesion and accretion. In this work, a mathematical model and hybrid numerical-analytical solutions are developed for the freezing of a supercooled droplet immersed in a cold air stream, subjected to the three main transport phenomena at the interface between the droplet and the surroundings: convective heat transfer, convective mass transfer and thermal radiation. Error-controlled hybrid solutions are obtained through the extension of the generalized integral transform technique to the transient partial differential formulation of this moving boundary heat transfer problem. The nonlinear boundary condition for the interface temperature is directly accounted for by the choice of a nonlinear eigenfunction expansion base. Also, the nonlinear equation of motion for the freezing front is solved together with the ordinary differential system for the integral transformed temperatures. After comparisons of the solution with previously reported numerical and experimental results, the influence of the related physical parameters on the droplet temperatures and freezing time is critically analysed.