LABORATORIO #6 - PRÁCTICA DE EBULLICIÓN

FUNDAMENTOS TRANSFERENCIA DE CALOR 2024-1

UNIVERSIDAD NACIONAL DE COLOMBIA

Andrés Camilo Calderón Gómez, Juan Esteban Vizcaino Gomez

[email protected], [email protected]

RESUMEN El presente informe busca estudiar los Ebullición en convección natural: Se da con
fenómenos de transferencia de calor por ∆T_exceso muy bajos teniendo líquido ligeramente
ebullición. Es analizada la ebullición nucleada y la sobrecalentado. Se da una transferencia de calor
ebullición en película. Posteriormente, serán relativamente baja.
usadas relaciones empíricas para determinar los
coeficientes de transferencia de calor. Además, Ebullición nucleada: Se caracteriza por tener
será realizado un balance de energía en el burbujas sobre la superficie de calentamiento. En
condensador para hallar el coeficiente de este caso el ∆T_exceso es más alto que en
transferencia de calor de este y se analizarán los convección natural. Al tener una temperatura en la
coeficientes de transferencia de calor por superficie mayor que la temperatura de saturación
convección y el coeficiente global de transferencia del líquido se forman burbujas con mucha más
de calor del condensador. rapidez. Para este régimen se presenta la
correlación de Rohsenow [2].
Palabras Clave— Ebullición nucleada, Ebullición
en película, Condensador.

(2)

I. MARCO TEÓRICO Kutateladze y Zuber determinaron el flujo máximo

de calor de manera teórica con la siguiente
1A. Ebullición: Proceso de cambio de fase de expresión: [2]
líquido a vapor que ocurre en una interfase sólido
líquido, cuando un líquido se pone en contacto con
una superficie que se encuentra a una temperatura (3)
mayor a su temperatura de saturación. Al ser un
proceso de transferencia de calor por convección Donde C_cr dependen de la configuración
el calor transferido se calcula mediante la ley de geométrica. Los valores para C_cr se encuentran
enfriamiento de Newton: en la Tabla 1.

(4)
(1)
Ebullición de transición: en este régimen a
1A-1. Ebullición en estanque ocurre cuando no se medida que se aumenta el ∆T_exceso se van
tiene flujo másico del fluido, es decir, cuando el formando más burbujas que se quedan en la
fluido se encuentra en reposo y cualquier superficie de calentamiento, estas funcionan como
movimiento de este se debe solo a las corrientes aislante térmico debido a la baja conductividad del
de convección natural y al movimiento de las vapor. Una vez se alcanza el punto crítico o
burbujas formadas durante la ebullición. La máximo de transferencia de calor, se entra en la
ebullición en un estanque se caracteriza por medio región llamada transición, donde el ∆T_exceso
de la curva propuesta por Nukiyama presentada en aumenta, sin embargo, la transferencia de calor
la Figura 1. disminuye hasta el punto mínimo de flujo de calor
En esta curva se presentan cuatro zonas: (Punto de Leiderfrost).
 1D. Convección forzada interna en un serpentín
helicoidal.

(5) La convección se da en presencia de un fluido en

Así mismo en el VDI Heat Atlas [3] se tiene la movimiento, particularmente la convención
siguiente expresión: forzada interna se da cuando se forza a un fluido a
fluir por el interior de determinada geometría
usando una bomba, compresor, ventilador, entre
otros. Se puede determinar el calor transferido al
(6) fluido interno mediante la siguiente expresión.

Ebullición de película: Esta región se caracteriza

por que la superficie de calentamiento queda
completamente cubierta por la una película (9)
continua y estable de vapor, sin embargo, el En este tipo de convección se presenta en dos
∆T_exceso sigue aumentando y la radiación situaciones particulares. Cuando se tiene
empieza a tomar un papel importante en la temperatura superficial constante y cuando se
transferencia de calor. La correlación de Bromley tiene flujo de calor constante a través de la
[2] permite calcular el flujo de calor en este superficie de calentamiento o enfriamiento. Para el
régimen. caso de temperatura constante se tiene la
siguiente expresión.

(7)
Donde C_pel es 0.62 para cilindros horizontales y (10)
0.67 para esferas. Para el caso donde no se conoce T_s pero si T_∞,
se puede usar la ecuación 10 reemplazando T_s
1B. Coeficiente de transferencia de calor por por T_∞ y h por U.
ebullición: El coeficiente de transferencia se Existen distintas correlaciones para calcular el
determina a partir de la ley de enfriamiento de número de Nusselt, para ello hay que calcular el
Newton. número de Reynolds con la siguiente correlación.

(8)
(11)
1C. CONDENSACIÓN:
Como la ebullición, la condensación es un proceso
Una de las formas de calcular el número de Nusselt
de cambio de estado, pero esta vez se pasa de
es la propuesta por Churchill y Bernstein que se
vapor a líquido, y se presenta cuando el vapor
muestra a continuación.
tiene una temperatura menor que su temperatura
de saturación, por lo general cuando dicho vapor
entra en contacto con una superficie sólida que
tenga una temperatura superficial menor que la
temperatura de saturación del vapor.

La condensación se divide en: Condensación tipo

gota y tipo película. En la condensación tipo
película, el condensado forma una capa de líquido
sobre la superficie que fluye hacia abajo por acción (12)
de la gravedad, mientras que, en la condensación Donde RePr > 0.2. Determinando las propiedades a
tipo gota, el vapor forma gotas sobre la superficie T_prom.
en lugar de una película continua.
Otra correlación importante es la propuesta por
Zukauskas [2].
 Se debe cumplir que 20000 < Re < 150000 con
R/r < 5.
Para 1500 < Re < 20000 se recomienda usar la
correlación de Pratt [5]:
(13)

Donde C y n se determinan a partir de la Tabla 2. (20)

Otra correlación propuesta por Zukauskas [2] es: 2. INSTRUMENTOS Y MATERIALES [6]:

2.1 Equipo:
Para la realización de la práctica, se cuenta con un
módulo de ebullición y condensación H112S, este
(14) permite hacer diferentes variaciones de
Donde todas las propiedades se evalúan a T_∞ parámetros con el fin de analizar los fenómenos de
excepto P_rs la cual se evalúa a T_s. Si P_r <= 10, n transferencia de calor por ebullición y
es 0,37; mientras que Pr > 10, n es 0,36. condensación.

C y m dependen del número de Reynolds y se Modulo H112S: Presente en la Figura 2, cuenta con
determinan a partir de la Tabla 3. las siguientes partes:

Una vez se determine el valor del número de ● Cámara: Cilindro vertical de vidrio en la cual
Nusselt teórico, es posible determinar el flujo de se encuentra el fluido de transferencia de
calor usando la diferencia de temperatura media calor SES36 y el cilindro calefactor, el
logarítmica: condensador helicoidal y los termopares.

● Condensador helicoidal refrigerado por agua:

(15) Tubo helicoidal de diámetro interno de 7 mm
y de diámetro externo de 7.6 mm, así como
Donde la temperatura media logarítmica es: una longitud total de 18.0 cm , por el cual
fluye agua de enfriamiento.

● Cilindro de calentamiento: Elemento

encargado de disipar la potencia para
producir la ebullición. Tiene una longitud de
(16)
42 mm y un diámetro de 12.7 mm.
Teniendo en cuenta que:
● Rotámetro: El cual permite controlar cómo
medir el flujo del agua de refrigeración, antes
(17) de que entre en el condensador.
● Manómetro: Indica la presión en el interior
(18) de la cámara.
● Válvula de carga y drenaje: Utilizada para
Dichas correlaciones son usadas para tubos rectos. recargar la cámara con SES36.
Para el caso de intercambiadores de calor de tubo ● Válvula de ventilación: Usada para sacar el
helicoidal se utiliza la siguiente correlación de aire presente en la cámara.
Schmit [5].
A. Sistema de control:
Este dispositivo permite establecer el valor de la
potencia de trabajo, la cual será disipada por el
cilindro de calentamiento del módulo H112S, así
(19) como visualizar la temperatura de la superficie de
Donde Nu_c es el número de Nusselt para el tubo éste. Se puede observar en la Figura 3.
helicoidal, Nu_s para el tubo recto, r el radio del
tubo y R el radio de la hélice.
B. Módulo de servicio y adquisición de 4. RESULTADOS
datos:
Este módulo permite tanto monitorear cómo Se realizó un proceso de revisión de datos para
registrar los datos de temperatura, presión, caudal eliminar valores incongruentes obtenidos en la
y potencia del módulo H112S. El sistema cuenta
recolección de datos durante la prueba de
con 5 termopares:
● T1: cilindro calentador laboratorio.
● T2: Temperatura del líquido
● T3: temperatura del vapor A Partir de las Tablas del Refrigerante SES36 [7][8]
● T4: temperatura de la entrada del agua de se grafica la curva Presión vs T_sat del refrigerante,
refrigeración. se evidencia en la Figura 5. La ecuación fue usada
● T5: temperatura de la salida del agua de
para calcular T_sat del refrigerante a lo largo del
refrigeración.
experimento.

C. Sistema de refrigeración:
Este sistema cuenta con un depósito de agua 4.1 Comportamiento térmico del sistema durante
destilada, una bomba sumergible y las líneas de
la prueba.
conexión al módulo H112S. Se observa en la Figura
4.
Se grafican Δ T Exceso y Potencia a través del
3. METODOLOGÍA [6]: tiempo, este comportamiento se evidencia en las
1. Calibrar: Encender el equipo y calibrar Figura 6 y Figura 7.
temperatura, potencia, caudal y presión.
2. Abrir Software: Iniciar el software y seleccionar Se procede a graficar la curva de ebullición del
“Collect Data”. fluido Potencia vs Δ T Exceso para cada presión, se
3. Elegir Visualización: Seleccionar “Multiple Graph evidencian los regímenes de ebullición y puntos de
Display” para graficar datos. calor máximo y mínimo en las Figura 8, Figura 9,
4. Configurar Frecuencia: Revisar la frecuencia de Figura 12, Figura 13 como aproximación de la
toma de datos (5-30 segundos). curva de Nukiyama.
5. Seleccionar Variables: Elegir T1, T2, T3, T4, T5,
presión, potencia y flujo de agua. Se muestran los fenómenos observables que
6. Iniciar Registro: Oprimir “Begin recorder”. relacionan los momentos en el Anexo 2.
7. Configurar Caudal y Potencia: Encender la
bomba, ajustar el caudal y establecer potencia a 4.2 .Determinación de los coeficientes de
20W. transferencia de calor por ebullición.
8. Registrar Datos: Esperar estabilización y registrar
temperaturas, caudal, potencia y presión. Se calcula el C_sf despejando la ecuación 2, por
9. Aumentar Potencia: Incrementar potencia en otra parte, con la misma ecuación usando un n=1,7
pasos de 40W hasta 340W. se procede a calcular el flujo de calor. Los
10. Controlar Temperatura: Si T1 > 160°C, el resultados obtenidos se encuentran registrados en
controlador suspenderá la potencia. la Tabla 4.
11. Repetir Pruebas: Seleccionar 3 caudales
adicionales y repetir pasos 7-9. Una vez hallado el flujo de calor se calcula el valor
12. Finalizar Registro: Oprimir “End recorder” al teórico del coeficiente de transferencia de calor h
terminar. por ebullición. Se registran los datos en la Tabla 4.

Se comparan los coeficientes de calor h teórico y

experimental, se registran los resultados en la
Tabla 4 y Tabla 5 y su comportamiento se
evidencia en las Figura 10 y Figura 11.
4.3 Determinación del flujo de calor máximo y mientras que para 70kPa se presentaron errores
mínimo. más altos con valores de 29,33% y 29,99%
respectivamente.
A partir de los valores de potencia registrados se
determina el valor del flujo de calor máximo y En el cálculo de los coeficientes globales de
mínimo con el uso de las ecuaciones 3, ecuación 5. transferencia se compararon los coeficientes
ecuación 6. y las Tabla 1 , Tabla 2, Tabla 3 . Los teórico y experimental para cada una de las
resultados se registraron en las Tabla 8 - Tabla 10. presiones, presentando relaciones de 59% para
50kPa y 76% para 70kPa.
4.4 Determinación del coeficiente global de
transferencia de calor en el Se estimaron las pérdidas para el condensador de
condensador. 193W para una presión de 50kPa y 162W para 70
kPa.
Utilizando las ecuaciones para convección forzada
interna en un serpentín helicoidal (ecuaciones 6. CONCLUSIONES
10-20) se calcula el coeficiente global de
transferencia de calor para cada presión. Los La estabilización de la temperatura de exceso en el
resultados se muestran en las Tablas 4. Tabla 6. régimen de ebullición nucleada refleja una
Tabla 7. transferencia de calor eficiente, con un aumento
controlado de potencia que coincide con los
Adicionalmente en las tablas anteriores se regímenes esperados según la curva de Nukiyama.
calcularon las pérdidas de calor reales .
Las diferencias entre los coeficientes teóricos y
5. DISCUSIÓN DE RESULTADOS. experimentales, con variaciones entre 12% y 17%,
indican errores experimentales asociados a las
La temperatura de exceso aumenta gradualmente mediciones de temperatura y potencia,
al inicio por el calentamiento, estabilizándose especialmente en configuraciones de mayor
luego en el régimen nucleado. Este incremento presión.
refleja el aumento de potencia aplicada al sistema.
Los resultados muestran concordancia con los
Se evidencio dentro de los resultados cada uno de fenómenos observados en la práctica como la
los regímenes de el proceso de ebullición en las aparición de una película de vapor que ocasiona la
curvas de Nukiyama graficadas, estos resultados disminución del flujo de calor en la curva de
fueron acordes a los fenómenos observados Nukiyama.
dentro del experimento.
Durante la práctica se evidenció el fenómeno de
Se presenta una relación directa entre el ebullición para diferentes presiones.
coeficiente de transferencia h y la potencia
suministrada, evidenciado en los resultados de las Se evidencia que la fase de transición se presenta a
gráficas y tablas. una potencia cercana a 320W, donde hay un
aumento repentino de la temperatura.
La variación entre los coeficientes de transferencia
teórico y experimental presentan valores entre el REFERENCIAS:
12% y 17% para cada presión respectivamente
como se evidencia en la Tabla 5. mostrando
errores asociados en la práctica.

Para una presión de 50 kPa se obtuvieron errores

relativos para q_max y q_min de 24% y 15%
[6]Video - Guia de laboratorio - Práctica de
Ebullición.

COPIAR:
[7]P. Fröba, H. Kremer, A. Leipertz, F. Flohr and C.
Meurer. Thermophysical Properties of a Refrigerant
Mixture of R365mfc (1,1,1,3,3-Pentafluorobutane)
and Galden HT55 (Perfluoropolyether). Vol. 28, No.
2. International Journal of Thermophysics. 2007.

[8]Propiedades SES36. SOLKATHERM. 2009.

ANEXOS:

Figura 4: Sistema de refrigeración. [6]

Figura 1: Curva de Nukiyama. [1]

Tabla 1: Coeficientes C_cr.

Tabla 2: Valores de C y n para la ecuación 13.

Figura 2: Equipo H112S para la práctica de

ebullición. [6]

Tabla 3: Valores de C y m para la ecuación 14.

Figura 3: Módulo de control para la práctica de

ebullición. [6]
ANEXOS:

Figura 9: Potencia vs ΔT_exceso para una

presión de 70kPa.

Figura 5: Presión vs T_sat del refrigerante

SES36 mediante gráficas y tablas. [7][8]

Figura 10: Coeficiente de transferencia de

calor h vs ΔT_exceso para una presión de
50kPa.

Figura 6: Potencia vs Tiempo.

Figura 11: Coeficiente de transferencia de

calor h vs ΔT_exceso para una presión de
70kPa.
Figura 7: ΔT_exceso vs Tiempo.

Figura 8: Potencia vs ΔT_exceso para una Figura 12: Potencia vs ΔT_exceso para una
presión de 50kPa. presión de 50kPa - Aproximación curva de
Nukiyama.
 Tabla 8: Cálculo de Flujo de calor Máximo
teórico

Tabla 9: Cálculo de flujo de calor mínimo

Figura 13: Potencia vs ΔT_exceso para una teórico, q* máximo y mínimo experimental
presión de 70kPa - Aproximación curva de para 50 kPa.
Nukiyama.

Tabla 10: Cálculo de flujo de calor mínimo

teórico, q* máximo y mínimo experimental
para 70 kPa.

Tabla 4: Cálculos para flujo de calor,

coeficientes de transferencia.

Tabla 5: Comparación coeficiente de

transferencia teórico y experimental

Tabla 6: Resultados Globales a 50kPa

Tabla 7: Resultados Globales a 70kPa