Based on a previously proposed dimensionless phase-change-driven frosting model, this study numerically investigates frost formation on a horizontal cold plate under natural convection using a Eulerian multiphase framework coupled with species transport. The model is validated against experimental data, showing errors within 5–18%; the maximum deviation of 17.07% occurs at Tw = −25 °C, possibly due to increased experimental uncertainty at very low temperatures. Results demonstrate that lower cold plate temperatures lead to greater frost thickness and higher ice volume fraction. A key physical insight is that under natural convection, local convective circulation causes enhanced frosting at the plate edges, resulting in spatial non-uniformity in both thickness and density. The study covers cold plate temperatures from −10 °C to −25 °C at relative humidity of 60%. The frost growth rate and density at both ends of the cold plate exceed those in the central region, and this difference intensifies with decreasing temperature. Within the frost layer, humid air velocity is nearly zero, while maximum velocity occurs near the sides due to natural convection. The simulation results show good agreement with experimental data, confirming the model’s reliability for natural convection scenarios.
Yang et al. (Tue,) studied this question.