⚡ ingeniería eléctrica

Improvement of Two-Phase Flow Simulation in a Thermosyphon Using a New Saturation Temperature Model

Este estudio introduce y valida un nuevo modelo de temperatura de saturación media dentro de un solver de CFD basado en OpenFOAM que mejora significativamente la precisión de las simulaciones de flujo bifásico transitorio en termosifones al contabilizar dinámicamente las variaciones internas de presión y temperatura, superando los enfoques convencionales de saturación constante con errores inferiores al 6%.

Autores originales: Quang-Hat Cao, Shehnaz Akhtar, Minseong Choi, Kyoungsik Chang, Sang-Wook Lee

Publicado 2026-11-04

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Quang-Hat Cao, Shehnaz Akhtar, Minseong Choi, Kyoungsik Chang, Sang-Wook Lee

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

El calor se desplaza por el mundo de formas predecibles, pero cuando un líquido se convierte en gas y viceversa dentro de un tubo sellado, la física se convierte en una compleja danza de presión, temperatura y gravedad. Estos tubos sellados, conocidos como termosifones, son dispositivos simples pero potentes que se utilizan para mover el calor de un lugar a otro sin partes móviles como bombas o ventiladores. Funcionan porque un líquido en la parte inferior absorbe calor, hierve convirtiéndose en vapor, sube a la parte superior, libera ese calor y vuelve a convertirse en líquido para caer de nuevo. Este ciclo es esencial para enfriar desde naves espaciales hasta sistemas de energía solar. Sin embargo, predecir exactamente cómo se comportan estos tubos, especialmente cuando apenas están arrancando, ha sido un desafío persistente para los ingenieros. Durante décadas, los modelos informáticos utilizados para diseñar estos sistemas dependían de un supuesto simplificador: que la temperatura a la que el líquido hierve permanece constante durante todo el proceso. Aunque esto facilitaba las matemáticas, ignoraba la realidad de que la presión y la temperatura dentro del tubo cambian constantemente a medida que el sistema se calienta desde un arranque en frío.

Un equipo de investigadores de la Universidad de Ulsan en Corea del Sur ha desarrollado una nueva forma de simular este comportamiento que captura esas condiciones cambiantes. En lugar de asumir un punto de ebullición fijo, su nuevo modelo informático calcula la temperatura de saturación —el punto específico donde el líquido se convierte en gas— como un valor dinámico que cambia según las temperaturas promedio de la sección caliente inferior y la sección fría superior. Integraron este modelo en una herramienta de software gratuita y de código abierto llamada OpenFOAM, que permite a los científicos visualizar el flujo invisible de fluidos dentro del tubo. Al rastrear cómo se mezclan y se mueven el líquido y el vapor, y al contabilizar cómo el calor viaja a través de las paredes metálicas del propio tubo, los investigadores crearon una imagen más realista de lo que sucede dentro de un termosifón. Su trabajo sugiere que, al permitir que la temperatura de ebullición cambie conforme el sistema opera, pueden predecir el rendimiento del tubo con una precisión mucho mayor que los métodos anteriores.

Los investigadores probaron este nuevo enfoque simulando el arranque de un termosifón y comparando los resultados con experimentos del mundo real y modelos informáticos antiguos. En el pasado, las simulaciones estándar a menudo tenían dificultades para predecir la temperatura de las paredes del tubo durante la fase crítica de arranque, fallando por un margen significativo en ocasiones. El nuevo modelo, sin embargo, mostró una alineación mucho más estrecha con la realidad. Cuando el equipo comparó sus predicciones con datos experimentales, el error en las lecturas de la temperatura de la pared cayó a menos del seis por ciento en todos los escenarios que probaron. Esta mejora es particularmente importante en la sección del evaporador, la parte inferior del tubo donde primero se absorbe el calor. Aquí, los modelos antiguos tendían a sobreestimar la temperatura porque no podían contabilizar los cambios rápidos en el punto de ebullición a medida que el sistema se calentaba. El nuevo modelo, por el contrario, rastreó estos cambios locales, proporcionando una visión más clara de cómo el calor se mueve desde la pared hacia el líquido.

Más allá de simplemente acertar con los números, las simulaciones revelaron cómo la cantidad de líquido dentro del tubo afecta la rapidez con la que el sistema se establece en un ritmo constante. Los investigadores descubrieron que la "relación de llenado", o cuánto del tubo está lleno de líquido, desempeña un papel importante en la velocidad de arranque. Cuando el tubo contenía solo una pequeña cantidad de líquido, tardaba mucho más en alcanzar un estado operativo estable. En estos escenarios de bajo llenado, el líquido luchaba por formar una película continua a lo largo de las paredes, lo que provocaba parches de superficies secas que debían humedecerse repetidamente antes de que el flujo se estabilizara. Esta inestabilidad retrasaba el sistema, tardando hasta 55 segundos en estabilizarse. A medida que el equipo aumentaba la cantidad de líquido, el sistema se estabilizaba mucho más rápido. Cuando el tubo se llenaba a su capacidad, el estado estacionario se alcanzaba en tan solo 18 segundos. El exceso de líquido aseguraba un suministro constante a la superficie calentada, evitando los parches secos y permitiendo que la circulación se suavizara casi de inmediato.

El estudio también analizó cómo el tubo maneja diferentes cantidades de entrada de calor. Independientemente de si el calor aplicado era moderado o alto, el nuevo modelo superó consistentemente a los enfoques estáticos anteriores. Predijo correctamente que, a medida que se añade más calor, la resistencia general al flujo de calor disminuye, lo que significa que el tubo se vuelve más eficiente al mover la energía. Esto sucede porque los niveles de calor más altos hacen que el líquido hierva con más vigor, creando más vapor que sube rápidamente y transporta la energía hacia la parte superior. Los investigadores señalaron que, si bien su modelo es un paso significativo hacia adelante, todavía tiene límites. Supone que el flujo es suave y bidimensional, y aún no tiene en cuenta los gases no condensables que podrían quedar atrapados dentro de los tubos reales. No obstante, el trabajo demuestra que tratar la temperatura de ebullición como una variable en lugar de una constante es la clave para comprender estos sistemas. Al refinar la forma en que simulamos la danza invisible del cambio de fase, esta nueva herramienta ofrece una forma fiable para que los ingenieros diseñen mejores sistemas de refrigeración para el futuro, asegurando que el calor se gestione de manera eficiente en todo, desde la electrónica hasta la exploración espacial.

Resumen Técnico: Mejora de la Simulación de Flujo Bifásico en un Termosifón mediante un Nuevo Modelo de Temperatura de Saturación

Planteamiento del Problema
Los termosifones (tubos de calor sin mecha) son fundamentales para la gestión térmica en aplicaciones aeroespaciales, electrónicas y solares. Si bien la Dinámica de Fluidos Computacional (CFD) ofrece un método rentable para investigar su flujo interno y transferencia de calor, las simulaciones numéricas existentes suelen basarse en una suposición simplificadora: una temperatura de saturación ( $T_{sat}$ ) constante. Este artículo identifica que tal suposición es inadecuada para modelar operaciones transitorias, particularmente durante el arranque. Durante el arranque, la presión y la temperatura interna del termosifón evolucionan dinámicamente hasta que se alcanza un estado estacionario; una $T_{sat}$ fija no logra capturar este comportamiento transitorio, lo que conduce a imprecisiones en la predicción de las temperaturas de las paredes y la resistencia térmica. Además, los autores señalan que los marcos de trabajo de CFD de código abierto como OpenFOAM poseen actualmente capacidades limitadas para simular la transferencia de calor conjugada que involucra cambios de fase.

Metodología
Para abordar estas limitaciones, los autores desarrollaron un nuevo solver basado en el solver interFoam de OpenFOF (versión 8). La metodología integra los siguientes componentes:

Ecuaciones de Gobierno: El solver emplea el método de Volumen de Fluidos (VOF) para el seguimiento de la interfaz y asume un flujo laminar, incompresible y bidimensional. Resuelve las ecuaciones de continuidad, cantidad de movimiento y energía de la mezcla, incluyendo términos de calor latente para el cambio de fase y la transferencia de calor conjugada entre el fluido y las paredes sólidas.
Modelo de Cambio de Fase: Se utiliza el modelo de evaporación-condensación de Lee para calcular las tasas de transferencia de masa.
Nuevo Modelo de Temperatura de Saturación: La innovación central es un modelo de temperatura de saturación media derivado de las condiciones de equilibrio del modelo de Lee. En lugar de un valor constante, la $T_{sat}$ se calcula como una función ponderada de las temperaturas medias en las secciones del evaporador ( $\bar{T}_e$ ) y del condensador ( $\bar{T}_c$ ):
$T_{sat} = \frac{A\bar{T}_e + \bar{T}_c}{A + 1}$
donde $A$ es un coeficiente adimensional dependiente de los coeficientes de evaporación/condensación y las densidades de fase.
Esquema Numérico: El algoritmo PIMPLE gestiona el acoplamiento presión-velocidad, y el Limitador Universal Multidimensional para Soluciones Explícitas (MULES) realiza el seguimiento de la fracción de volumen. Las simulaciones utilizan un paso de tiempo adaptativo con un número de Courant de 0.8.
Validación: El modelo fue validado contra datos experimentales publicados y resultados numéricos previos a través de múltiples condiciones de flujo de calor y geometrías (Casos 1–4).

Contribuciones Clave

Desarrollo de un Modelo de Saturación Transitorio: El artículo introduce una formulación donde la $T_{sat}$ evoluciona dinámicamente basándose en el estado térmico del evaporador y el condensador, en lugar de permanecer estática.
Implementación en OpenFOAM: El estudio proporciona un solver funcional de código abierto capaz de manejar la transferencia de calor conjugada y el cambio de fase en termosifones, abordando una brecha en las capacidades actuales de OpenFOAM.
Mejora de la Precisión: El nuevo modelo se enfoca específicamente en la mejora de las predicciones de la temperatura de la pared en la sección del evaporador, la cual se identifica como la región más crítica para el diseño térmico.

Resultados

Comportamiento de Arranque: En el caso de validación (Caso 1), el modelo capturó con éxito la evolución transitoria de la temperatura desde el arranque hasta el estado estacionario (aproximadamente 120 segundos). Los perfiles de temperatura predichos en las secciones del evaporador, adiabática y del condensador mostraron una estrecha concordancia con los datos experimentales, con una desviación máxima de menos del 2%.
Precisión de la Temperatura de la Pared: Comparado con las simulaciones que utilizan una temperatura de saturación constante, el modelo de saturación media propuesto redujo significativamente los errores en las predicciones de la temperatura de la pared del evaporador. Los errores se mantuvieron por debajo del 6% en todos los casos probados. El modelo constante tendió a sobrepredecir las temperaturas de la pared en el evaporador debido a su incapacidad para capturar el sobrecalentamiento local que impulsa la ebullición nucleada.
Resistencia Térmica: El nuevo modelo predijo la resistencia térmica con un error máximo del 34%, una mejora notable frente al error del 54.7% observado en las simulaciones previas de saturación constante. Se observó que la resistencia térmica disminuye a medida que aumenta la potencia de calentamiento, de acuerdo con la literatura establecida.
Efectos de la Relación de Llenado: El estudio investigó la relación entre la relación de llenado y el tiempo requerido para alcanzar el estado estacionario. Los resultados indicaron que las relaciones de llenado bajas (por ejemplo, 25%) requieren tiempos significativamente más largos para estabilizarse (55 s) debido a ciclos intermitentes de secado y rehumectación. A medida que la relación de llenado aumentó al 100%, el tiempo de estado estacionario disminuyó a 18 s, logrando el sistema un perfil de temperatura más uniforme e isotérmico.

Significancia y Reivindicaciones
Los autores afirman que el solver propuesto basado en OpenFOAM ofrece una "herramienta fiable y accesible" para simular la transferencia de calor por cambio de fase en sistemas de termosifón. La principal significancia radica en la capacidad demostrada del nuevo modelo de temperatura de saturación para mejorar la precisión de las simulaciones transitorias, particularmente durante la fase crítica de arranque y en la sección del evaporador. Al alejarse de la suposición de saturación constante, el modelo proporciona una descripción físicamente más representativa de la dinámica interna del termosifón.

El artículo mantiene un tono modesto respecto a su alcance, reconociendo limitaciones como la suposición de flujo laminar y bidimensional, así como la omisión de gases no condensables, que pueden contribuir a las desviaciones en la sección del condensador. Los autores sugieren que trabajos futuros extenderán el modelo a tres dimensiones y refinarán el cierre de la temperatura de saturación para incluir la ponderación de la fracción de fase local, pero no proponen nuevas aplicaciones específicas más allá de la simulación de sistemas de termosifón.

Más como este