Numerical simulations of freezing behaviors of water droplets impacting cold hydrophobic surfaces

Freezing, as a prominent phase transition phenomenon, occurs frequently in nature and industry. This study employs a combination of the Volume of Fluid(VOF) model, solidification model, and a dynamic contact angle model to simulate the dynamics and solidification process of droplets impacting supercooled surfaces. Two distinct freezing morphologies, central-concave and central-convex, are observed. These different morphologies are attributed to the competition between the time scales associated with droplet impact dynamics and solidification (retraction time t_r and solidification time t_f,p). Thus, we propose a scaling law to predict the occurrence of different morphologies: when t_r < t_f,p, the final freezing morphology exhibits a central-convex shape, while when t_r > t_f,p, it exhibits a central-concave shape. In order to explore effects of nucleation on solidification, the classical nucleation theory is adopted to determine the nucleation time, and thus the freezing behaviors can be well captured. Furthermore, we examine the effect of volume expansion during phase transition on the freezing behaviors, and find that the increase in freezing time due to volumetric expansion can be well predicted by the scaling law t_f ∼ Δ(H²)/ΔT.