Numerical Study on the Effects of Wall Wettability on Heat Transfer and Flow Characteristics during Flow Boiling in Microchannels

Resumen Técnico: Estudio Numérico sobre los Efectos de la Humectabilidad de la Pared en la Transferencia de Calor y las Características del Flujo durante la Ebullición en Flujo en Microcanales

Planteamiento del Problema
La miniaturización de chips de alto rendimiento ha provocado aumentos significativos en la densidad de flujo de calor local, que a menudo supera los 1000 W/cm², algo que los métodos de enfriamiento tradicionales luchan por disipar. La ebullición en flujo en microcanales es una solución prometedora debido a sus altos coeficientes de transferencia de calor y su capacidad de absorción de calor latente. Sin embargo, el rendimiento de estos sistemas está fuertemente influenciado por la dinámica de las burbujas y la evolución de los patrones de flujo. Un factor crítico, aunque complejo, en este proceso es la humectabilidad de la pared. Si bien la literatura existente ha explorado superficies hidrofílicas, hidrofóbicas y de humectabilidad mixta, los mecanismos acoplados que gobiernan cómo las condiciones extremas de humectabilidad afectan la interacción entre el comportamiento de las burbujas, la estabilidad del flujo, la eficiencia de la transferencia de calor y la irreversibilidad termodinámica (generación de entropía) aún no se comprenden suficientemente. Específicamente, existe la necesidad de cuantificar las compensaciones entre la mejora de la transferencia de calor y la estabilidad del flujo bajo diferentes regímenes de humectabilidad.

Metodología
Este estudio emplea un modelo numérico tridimensional basado en el método de Volumen de Fluido (VOF) para simular la ebullición en flujo en un microcanal ranurado. La simulación utiliza el modelo de Fuerza de Superficie Continua (CSF) para capturar los efectos de la tensión superficial y un modelo de evaporación-condensación (basado en Lee et al.) para gestionar el cambio de fase.

• Modelo Físico: La geometría consiste en un microcanal con ranuras triangulares en la pared lateral, diseñado para estudiar la nucleación y el crecimiento de burbujas. El fluido de trabajo es agua desionizada y el dominio sólido es silicio.
• Condiciones de Humectabilidad: Se comparan dos condiciones extremas de superficie: una pared hidrofílica (ángulo de contacto $\theta = 30^\circ$) y una pared hidrofóbica ($\theta = 150^\circ$).
• Enfoque Numérico: Las ecuaciones de gobierno (conservación de masa, momento y energía) se resuelven mediante el método de volúmenes finitos utilizando el algoritmo PISO para el acoplamiento presión-velocidad. Un estudio de independencia de malla confirmó que una malla de aproximadamente 275,000 elementos proporciona suficiente precisión. El modelo fue validado con datos experimentales, mostrando una discrepancia máxima de menos del 8% en el aumento de la temperatura de la base del chip.
• Métricas de Análisis: El estudio evalúa la dinámica de las burbujas, los patrones de flujo, las fluctuaciones de temperatura de salida, el número de Nusselt (Nu), la caída de presión y la generación de entropía (tanto componentes térmicas como friccionales).

Contribuciones Clave y Resultados

1. Dinámica de Burbujas y Evolución del Patrón de Flujo:

   • Superficies Hidrofílicas ($\theta = 30^\circ$): Las fuertes interacciones líquido-sólido causan que las burbujas se adhieran a la pared y se fusionen, formando estructuras de película líquida. Esto promueve una transición hacia el flujo anular pero dificulta la desprendimiento de las burbujas.
   • Superficies Hidrofóbicas ($\theta = 150^\circ$): Las burbujas se desprenden rápidamente y son transportadas hacia la corriente principal. Esto acelera la acumulación de vapor y conduce a una transición más temprana del flujo de burbujas a flujo de tapón (slug flow) y flujo anular.

2. Rendimiento Térmico:

   • Eficiencia de Transferencia de Calor: La superficie hidrofóbica demuestra un rendimiento de transferencia de calor superior. El rápido desprendimiento de las burbujas facilita la eliminación del calor latente, manteniendo temperaturas de pared más bajas y números de Nusselt más altos. Por el contrario, la superficie hidrofílica sufre una acumulación de calor debido a las burbujas adheridas, lo que resulta en temperaturas de pared más altas y un declive continuo en el número de Nusselt.
   • Estabilidad de Temperatura: Aunque las superficies hidrofóbicas mejoran la transferencia de calor promedio, exhiben una menor inestabilidad de la temperatura de salida en comparación con las superficies hidrofílicas. La superficie hidrofílica muestra mayores fluctuaciones de temperatura debido al impedimento de la disipación de calor por las persistentes películas líquidas.

3. Características del Flujo e Inestabilidad:

   • Caída de Presión: El canal hidrofóbico experimenta caídas de presión significativamente mayores y oscilaciones de presión más pronunciadas. Esto se atribuye a la mayor fracción de volumen de vapor y a la descarga periódica de burbujas que alteran la estructura del flujo local.
   • Campos de Velocidad: El canal hidrofóbico exhibe perturbaciones de flujo local más fuertes y gradientes de velocidad más altos debido a la dinámica activa de las burbujas, mientras que el canal hidrofílico muestra variaciones de velocidad más suaves y bajas.

4. Análisis de Generación de Entropía:

   • Térmica vs. Friccional: El estudio distingue entre la generación de entropía térmica (impulsada por gradientes de temperatura) y la generación de entropía friccional (impulsada por la disipación viscosa).
   • Mecanismo Dominante: La generación de entropía térmica es el contribuyente dominante a la generación total de entropía. La superficie hidrofílica genera una entropía térmica significativamente mayor debido a los grandes gradientes de temperatura causados por el deterioro de la transferencia de calor.
   • Compensación (Trade-off): Aunque la superficie hidrofóbica genera más entropía friccional debido a los intensos gradientes de velocidad, su generación total de entropía permanece más baja que la de la superficie hidrofílica debido a la reducción sustancial de la irreversibilidad térmica.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona una guía teórica valiosa para el diseño de la humectabilidad y la optimización de los sistemas térmico-hidráulicos en disipadores de calor de microcanales bajo condiciones de alto flujo de calor. Los hallazgos resaltan una compensación fundamental: mientras que las superficies hidrofóbicas mejoran significativamente la transferencia de calor por ebullición y reducen la irreversibilidad térmica, introducen una mayor resistencia al flujo y fluctuaciones de presión. En contraste, las superficies hidrofílicas, aunque ofrecen caídas de presión más estables, sufren de un deterioro severo de la transferencia de calor y mayores pérdidas termodinámicas.

El artículo concluye que la optimización de la humectabilidad requiere equilibrar la mejora de la transferencia de calor por cambio de fase frente a la penalización del aumento de la inestabilidad del flujo. Los resultados sugieren que para aplicaciones que priorizan el rendimiento térmico y la minimización de la entropía, las superficies hidrofóbicas son ventajosas, siempre que el sistema pueda acomodar la caída de presión y las oscilaciones de flujo asociadas. El estudio no propone aplicaciones futuras específicas o validaciones experimentales más allá del alcance de la investigación numérica actual, sino que enfatiza la necesidad de comprender estos mecanismos acoplados para un diseño eficaz de la gestión térmica.