La formación de escarcha es un fenómeno común que se observa en muchos campos diferentes. La acumulación de escarcha en las superficies de los intercambiadores de calor perjudica significativamente la eficiencia de la transferencia de calor1, obstruye el flujo de fluidos e interrumpe el rendimiento general de los intercambiadores de calor2, lo que en última instancia dificulta su funcionamiento normal3. Por lo tanto, comprender los mecanismos y el comportamiento de la formación de heladas es fundamental para abordar este problema en los sistemas de refrigeración4. En las últimas décadas, se ha dedicado un cuerpo sustancial de investigación a investigar las causas y características de la formación de heladas en estos sistemas.
Los estudios experimentales han demostrado que la formación de escarcha está influenciada por varios factores, incluida la temperatura del aire, la humedad y la temperatura de la superficie fría 5,6,7,8,9,10. Numerosos hallazgos experimentales han indicado que las temperaturas más bajas del aire entrante tienden a dar lugar a capas de escarcha más gruesas6, mientras que los niveles de humedad más altos contribuyen a la formación de capas de escarcha más densas7. Song et al. estudiaron la formación de escarcha en superficies horizontales y encontraron que las variaciones cíclicas de temperatura de la superficie fría pueden causar derretimiento en la interfaz de la capa de escarcha, lo que influye significativamente en la tasa de formación de escarcha, el espesor de la capa de escarcha y la densidad dinámica de la escarcha8. Otras investigaciones han examinado tanto la morfología como la distribución de las heladas. Jeong et al. observaron en sus experimentos que la escarcha se forma inicialmente cerca de la ensenada, lo que lleva a la ocurrencia de un fenómeno conocido como colina helada9. Noorshams et al. investigaron la formación de escarcha en la superficie de un tubo circular horizontal y encontraron que las capas de escarcha en las superficies delantera y trasera del cilindro eran más gruesas que las de las superficies superiores10. Además, varios estudios 11,12,13,14,15,16 han desarrollado modelos para predecir el espesor unidimensional de la capa de escarcha, utilizando patrones experimentales de formación de escarcha junto con enfoques teóricos y empíricos. Jones y Parker desarrollaron un modelo predictivo para el espesor de las heladas basado en la teoría de la difusión molecular11. La discrepancia entre su modelo y los datos experimentales se mantuvo por debajo del 30% durante un período de 3 h. Con la progresión de la tecnología computacional, un número cada vez mayor de investigadores ha recurrido a la dinámica de fluidos computacional (CFD) para simular la formación de escarcha. A diferencia de los modelos unidimensionales tradicionales, las simulaciones CFD brindan beneficios significativos, particularmente en la visualización de la distribución del espesor de la escarcha y los perfiles de temperatura. Cui et al. realizaron simulaciones CFD de formación de escarcha basadas en la teoría de nucleación12. Sus predicciones de espesor de escarcha mostraron una desviación de menos del 13% de los datos experimentales proporcionados por Lenic et al.13. Paralelamente, los estudios CFD sobre condensación en tubos estructurados han demostrado que las características geométricas como los hoyuelos14 o los pasos helicoidales15 mejoran la transferencia local de calor y masa. Recientemente, You et al.16 desarrollaron un modelo CFD basado en malla dinámica que caracteriza la capa de escarcha como un medio poroso en crecimiento e incorpora la difusión de vapor directamente, logrando una desviación relativa de menos del 5% mientras se mantiene un bajo costo computacional. Estos hallazgos subrayan el potencial de CFD para resolver fenómenos complejos de cambio de fase, ofreciendo información valiosa para el modelado de formación de heladas.
En conclusión, un número significativo de estudios 5,6,7,8,9,10,11,12,13 han explorado la formación de escarcha en superficies frías, contribuyendo a una comprensión evolutiva de los patrones de formación de escarcha bajo varios parámetros. Aunque se han desarrollado múltiples modelos numéricos, que incorporan diversas dimensiones y mecanismos, con frecuencia carecen de una validación integral. La mayoría de los estudios 6,7,8,9,10,11,12,13 validan principalmente modelos que utilizan espesor de escarcha, lo que limita su aplicabilidad más amplia 16. Para superar estas limitaciones, este artículo presenta un modelo numérico que integra el modelo de cambio de fase de Lee y el modelo de flujo multifásico euleriano, con énfasis en los mecanismos centrales que subyacen a la formación de heladas. Además, se introduce un método novedoso para calcular el límite superior de la fracción de volumen de heladas, teniendo en cuenta los cambios dependientes del tiempo en la densidad de las heladas, abordando así las deficiencias de los modelos anteriores. La precisión y confiabilidad del modelo propuesto se evalúan desde varios ángulos, incluido el espesor de la escarcha, las fluctuaciones de densidad y los patrones de formación de escarcha en diferentes condiciones experimentales. Esta extensa validación ofrece un marco teórico sólido para predecir con mayor precisión el comportamiento de las heladas en aplicaciones del mundo real.