CONVECCIÓN LIBRE

ANDRÉS SANTIAGO BARRETO

WILLIAM FELIPE MARTÍNEZ VERGARA
ALEJANDRO FIGUEROA JARAMILLO

UNIVERSIDAD NACIONAL DE COLOMBIA

FACULTADA DE INGENIERÍA
TRANSFERENCIA DE CALOR Y MASA
BOGOTÁ
2019

1
INTRODUCCIÓN

La convección es una de las tres formas de transferencia de calor, esta se da mediante la

transferencia de energía entre una superficie y el líquido o gas contiguos que están en
movimiento y comprende los efectos combinados de la conducción y el movimiento de
fluidos. Mientras el movimiento del fluido sea más rápido, mayor es la transferencia de calor
por convección.

Existen dos tipos de convección, la libre o natural y la forzada. En la convección libre, el

movimiento del fluido se debe a las fuerzas de “empuje” dentro de éste, mientras que la
convección forzada se impone de forma externa. El “empuje” de la convección libre se debe a
la presencia combinada de un gradiente de densidad del fluido y de una fuerza de cuerpo que
es proporcional a la densidad. Normalmente la fuerza de cuerpo es gravitacional en la
práctica, no obstante, también puede ser una fuerza centrífuga en una maquinaria de fluido
giratoria.

OBJETIVOS

1) Realizar el balance de energía y determinar sus pérdidas

2) Para los datos obtenidos determinar el tipo de convección dada en la práctica,

utilizando para ello los criterios de la literatura.

3) Utilizando las ecuaciones apropiadas determine los coeficientes de película teóricos

para el agua y el vapor

4) Comparar los coeficientes obtenidos en el literal anterior con los datos dados por la
literatura

5) Determine los coeficientes de película prácticos para el agua y el vapor

6) Representar gráficamente el log (Nu) vs log [(Gr)(Pr)] y desarrollar una ecuación

empírica aproximada mediante una técnica de ajuste

MATERIALES Y MÉTODOS

En la práctica de convección libre tiene como fin hallar los diferentes tipos de coeficientes de
película, así como realizar el análisis de cada uno de los datos obtenidos en la presente
práctica. Esta práctica se dividió en tres partes, debido a que se realizó el mismo proceso en
todas las tres, pero para la escala de medida se usaron 3 niveles diferentes (-½”,0, +½”). Los
materiales utilizados en la práctica son:

2
Figure 1, Equipo Modelo 9054 para prácticas de
convección
Figure 2, Probetas

Figure 3, Termómetro
Figure 4, Cronómetro

3
Figure 5, Tuberías y Válvulas Figure 6, Tuberias y Válvulas

Figure 7, Escala de medida (Nivel)

Como primera medida se hizo una revisión el equipo, observando las válvulas que debían
estar abiertas y cerrada. A continuación, se coloca la escala de medida en el nivel cero para el
primer ensayo. Con el tanque de alimentación abierto a la atmósfera, para que el equipo
funcionara a convección libre, se abrió lentamente v-9 (válvula), hasta llegar a un nivel
constante en el tanque. Después se abre la válvula de purga (v-10) a la vez que se agrega
vapor, el cual su presión debía ser regulada (v-1), para no sobrepasar el nivel indicado de la
máquina ubicado en la parte superior de esta.

Una vez terminada la purga del sistema, asegurando que no hay presencia de gases dentro del
sistema, las válvulas v-1 y v-10 se cerraron simultánea e inmediatamente se siguió regulando
la presión con v-1. Bajo estas condiciones el sistema ya funcionaba a convección libre. Se dio

4
paso a conectar el termómetro digital al selector de termopares con el fin de observar el
momento en el que la temperatura de salida del agua era estable. En ese momento se toman
los datos de la temperatura de entrada del agua, la temperatura de la pared inferior del tubo, la
temperatura de la pared superior del tubo, temperatura de salida del agua, la temperatura del
condensado, la temperatura de la cámara y los caudales. Para los dos siguientes ensayos solo
se modificó la posición de la escala de medida (+½’’ y -½’’ del nivel de referencia del ensayo
anterior) y se ajustó la alimentación de agua con las válvulas.

CÁLCULOS Y RESULTADOS

Para obtener los resultados se ejecutó primero una toma de datos de temperaturas y flujos de
agua, con el fin de sacar las diferentes temperaturas medias para determinar sus propiedades
con el fin de hacer los cálculos correspondientes para determinar los coeficientes de película
exteriores e interiores.

Los datos tomados se evidencian en la siguiente tabla

Table 1, Datos tomados en el laboratorio

Datos tomados en el laboratorio

Datos 1(+1/2'') 2(0'') 3(-1/2'')

Temperatura Entrada de T1 °C 31.6 21.3 42.9

Agua

Temperatura Pared Inferior T2 °C 87.5 83.1 92.7

Temperatura Pared Superior T3 °C 88.7 85.3 92.1

Temperatura Salida de Agua T4 °C 69.9 53.6 86.1

Temperatura Vapor T5 °C 82.6 78.0 88.1

Temperatura de Condensado Tc °C 74.8 73.8 97.1

Flujo de Agua ṁa m3/s 8.3E-03 1.3E-02 8.4E-03

Flujo de Condensado ṁv m3/s 4.3E-04 4.4E-04 6.9E-04

5
Cada dato de la temperatura fue tomado dependiendo donde estaba el termopar instalado
previamente en el equipo. Para la determinación de tanto el flujo de agua como el flujo de
condensado, se ejecutó un proceso de medir en un tiempo específico la cantidad de agua que
caía y así determinar el flujo volumétrico de las dos condiciones tanto del agua como del
condensado.

a). Dadas las temperaturas se ejecutó un balance energético, en donde a partir de las
siguientes ecuaciones se encontraron los calores del vapor y del agua.

𝑄𝑄 = 𝑄𝑄 ∗ ℎ𝑄𝑄 (1)
Donde Qv es el flujo de calor cedido por el vapor, 𝑄𝑄 es el flujo de condensado y ℎ𝑄𝑄 es la
entalpía de condensación de la cámara de vapor.

𝑄𝑄 = 𝑄𝑄 ∗ 𝑄𝑄 ∗ (𝑄4 − 𝑄1 ) (2)
Donde Qa es el flujo de calor tomado por el agua, 𝑄𝑄 es el flujo de agua, 𝑄𝑄 es el calor
específico del agua a la temperatura media de película y T1 y T4 son temperaturas tomadas
de la anterior tabla.

𝑄𝑄 = 𝑄𝑄 + 𝑄𝑄 (3)
Donde Qp es el flujo de calor perdido en el sistema.

Estos cálculos se ven evidenciados en las siguientes tablas.

Table 2, Cálculo del flujo de calor cedido por el vapor

Cálculo del flujo de calor cedido por el vapor

Dato Tv(°C) ṁv(m3/s) ρTv(m3/kg) hfg(J/kg) Qv(W)

1 82.60 4.25E-04 0.620 2647350.00 1.81E+03

2 78.00 4.39E-04 0.628 2639660.00 1.84E+03

3 88.10 6.94E-04 0.610 2656430.00 3.02E+03

Cálculo del flujo de calor tomado por el agua

6
 Dato T1(°C) T4(°C) ΔT(°C) Tf(°C) ρTf(m3/ Cp(J/k ṁa(m3/s Qa(
kg) g*K) ) W)

1 31.60 69.90 38.30 50.75 987.655 4183.9 8.28E-03 1.34

0

2 21.30 53.60 32.30 37.45 993.118 4182.0 1.25E-02 1.70

0

3 42.90 86.10 43.20 64.5 980.770 4182.5 8.42E-03 1.55

1
Table 3,Cálculo del flujo de calor tomado por el agua

Cálculo del flujo de calor perdido

Dato Qv(W) Qa(W) Qp(W)

1 1.81E+03 1.34 1813.17

2 1.84E+03 1.70 1841.78

3 3.02E+03 1.55 3022.33

Table 4, Cálculo del flujo de calor perdido

b). Como segunda parte del experimento se calculó el tipo de convección que se generó, es
decir si es libre, o combinada y si es un régimen laminar o turbulento.

Para determinar si el flujo es laminar o turbulento, se calcula el número de Grashof a partir de

la siguiente ecuación.

𝑄𝑄𝑄3 (𝑄𝑄 −𝑄𝑄 )𝑄2

𝑄𝑄 = (4)
𝑄2
Donde g es la gravedad, ꞵ es el coeficiente de expansión volumétrica del fluido (calculado en
la siguiente ecuación), D que es el diámetro, Tw que es la temperatura media entre T2 y T3,
Ta que es la temperatura media entre T1 y T4, ⍴ que es la densidad calculada a la
temperatura media entre Tw y Ta y 𝜇que es la viscosidad calculada a la temperatura media
entre Tw y Ta.

𝑄21 −𝑄24
𝛽= (5)
2𝑄1 𝑄4 (𝑄4 −𝑄1 )

7
Para el cálculo de ꞵ, se calculó G1 y G4 que son las gravedades del agua a las temperaturas 1
y 4 respectivamente y las temperaturas T1 y T4.

Para determinar si la convección es libre o combinada, se determina el número de Reynolds a

partir de la siguiente ecuación.

𝑄𝑄
𝑄𝑄 = (6)
𝑄
Donde D es el diámetro del tubo, G es la división entre el flujo de masa (ṁa) y el area de
flujo.

Por último, se calcula el número de Prandtl y el número de Nusselt, dependiendo de si es un

régimen laminar o régimen turbulento.

Para calcular el número Prandtl se utilizó la siguiente ecuación.

𝑄𝑄 𝑄
𝑄𝑄 = (7)
𝑄
Donde Cp es el calor específico evaluado en la temperatura media de Tw y Ta, la viscosidad
igualmente evaluado en la misma temperatura media y la conductividad térmica del fluido
evaluada en la misma temperatura media.

Para calcular el número de Nusselt se usa una ecuación genérica que varía dependiendo del
régimen si es laminar o turbulento.

𝑄𝑄 = 𝑄(𝑄𝑄 ∗ 𝑄𝑄)𝑄 (8)

Si es un régimen laminar C es igual 0.59 y n es igual a ¼.
Si es un régimen Turbulento C es igual a 0.1 y n es igual a ⅓.

Los resultados los encontramos en la siguiente tabla. Los tres datos coincidencialmente nos
entregaron resultados de régimen turbulento, así que fueron calculados dadas esas
condiciones.

Resultados de determinación del tipo de convección

Dato Gr Re2 Gr/Re2 Pr Ral Nu

1 2.8E+09 2.17E-03 1.3E+12 2.55 7.3E+09 193.57

2 2.7E+09 3.78E-03 7.0E+11 2.95 7.8E+09 198.61

8
 3 3.9E+09 2.82E-03 1.4E+12 2.29 8.9E+09 207.23

Table 5,Resultados de determinación del tipo de convección

c). Para la tercera parte del experimento se calculó lo que se denomina como los coeficientes
de película teóricos tanto exteriores como interiores.

Para calcular el coeficiente de película teórico interior se hizo uso de los resultados
calculados anteriormente dado que para calcular el coeficiente interior se usa la siguiente
ecuación.

𝑄𝑄𝑄
ℎ𝑄 = (9)
𝑄𝑄
Donde Nu es el número de Nusselt previamente calculado, k es la conductividad térmica del
agua y Di es el diámetro interior de la tubería.

Esta ecuación nos llevó a los siguientes resultados:

Coeficiente de película teórico interior

Dato Nu k(W/m*k) Di(m) hi(W/m2*K)

1 193.57 6.7E-01 1.3E-02 1.02E+04

2 198.61 6.6E-01 1.3E-02 1.02E+04

3 207.23 6.6E-01 1.3E-02 1.08E+04

Table 6, Coeficiente de película teórico interior

Ahora, para calcular el coeficiente de película exterior primero se debe determinar el número
Reynolds dado que si Reynolds es menor a 1800 se usa una ecuación y si es mayor a 1800 se
aplica otra ecuación. Las ecuaciones son las siguientes.

4ṁ𝑄
𝑄𝑄 = (10)
𝑄𝑄𝑄𝑄
Donde 𝜇𝑄 es la viscosidad del condensado Tc, ṁ𝑄 que es el flujo de condensado y D que es el
diámetro exterior del tubo.

9
Dadas las condiciones entregadas por el fluido a partir de la ecuación, se determina cuál de
las siguientes ecuaciones debe ser usada, la primera es para valores donde Re es menor 1800
y la segunda para mayores a 1800.

1
𝑄𝑄2𝑄 ℎ𝑄𝑄 𝑄3𝑄
1. ℎ𝑄 = 1.13( )4 (11)
𝑄𝑄 (𝑄5 −𝑄𝑄 )𝑄
Donde g es la gravedad, 𝜌𝑄 es la densidad a temperatura del condensado, ℎ𝑄𝑄 la entalpía de
vaporización, 𝑄𝑄 es la conductividad térmica a temperatura de condensado, T5 es la
temperatura del vapor y Tw que es la temperatura media T2 y T3.
1
𝑄𝑄 𝑄𝑄3𝑄 (𝑄𝑄 −𝑄𝑄 ) 3
2. ℎ𝑄 = 0.0076𝑄𝑄0.4 [ ] (12)
𝑄2𝑄
Donde g es la gravedad, 𝜌𝑄 es la densidad a temperatura del condensado, ℎ𝑄𝑄 la entalpía de
vaporización, 𝑄𝑄 es la conductividad térmica a temperatura de condensado, Re es el número
de Reynolds y 𝜌𝑄 es la densidad a temperatura del vapor.

Los resultados encontrados fueron los siguientes.

Coeficiente de película teórico Exterior

Dato Re h0(W/m2*k)

1 9.23E-02 9.52E+03

2 9.39E-02 9.21E+03

3 1.99E-01 1.08E+04

Table 7, Coeficiente de película teórico Exterior

d). Para la cuarta parte del experimento, se buscaron los coeficientes de película prácticos que
son calculados dependiendo de si se aplican en agua o en vapor.

Para la determinación del coeficiente interior dada las características del agua, se presenta la
siguiente ecuación.

𝑄𝑄
ℎ𝑄 = (13)
𝑄𝑄 (𝑄𝑄 −𝑄𝑄 )

10
Donde Qa es el calor tomado por el agua, Ai es el area interna del tubo, Tw es la temperatura
media de T2 y T3, y Ta es la temperatura media de T1 y T4.

Para la determinación del coeficiente exterior dada las características del vapor, se presenta la
siguiente ecuación.

𝑄𝑄
ℎ𝑄 = (14)
𝑄0 (𝑄𝑄 −𝑄𝑄 )
Donde Qv es el calor cedido por el vapor, A0 es el area interna del tubo, Tv es la temperatura
del vapor y Tw es la temperatura media entre T2 y T3.

Los resultados se presentan a continuación.

Calculo de coeficiente de película práctico para el agua

Dato Qa(W) Ai(m2) Tw(°C) Ta(°C) hi(W/m2*K)

1 1.34 0.024322 88.1 50.75 1.478

2 1.70 0.024322 84.2 37.45 1.495

3 1.55 0.024322 92.4 64.5 2.285

Table 8, Calculo de coeficiente de película práctico para el agua

Calculo de coeficiente de película práctico para el vapor

Dato Qa(W) A0(m2) Tv(°C) Tw(°C) h0(W/m2*K)

1 1.81E+03 0.030402 82.6 88.1 10851.50

2 1.84E+03 0.030402 78 84.2 9779.99

3 3.02E+03 0.030402 88.1 92.4 23130.64

Table 9, Calculo de coeficiente de película práctico para el vapor

e). Por último, se realizó un análisis gráfico a partir de los datos encontrados en el punto b, en
donde se buscó la constante C y la constante n, sin embargo en los resultados dado que los
tres datos encontrados fueron calculados con la misma ecuación debido a que se encontraban

11
en el mismo régimen, la ecuación que encuentra el número de Nusselt es la misma que se usa
1
para calcularlo, es decir, 𝑄𝑄 = 0.1(𝑄𝑄 ∗ 𝑄𝑄)3 . (15)

Los resultados fueron:

Número de Nusselt vs Número de

Rayleigh

Dato Ral Nu

1 7.25E+09 193.57

2 7.83E+09 198.61

3 8.90E+09 207.23

Table 10, Número de Nusselt vs Número de Rayleigh

Graph 1, Número de Nusselt vs Número de Rayleigh

DISCUSIÓN DE RESULTADOS

A partir de toda la experimentación se encontraron diferentes datos importantes para la

investigación de la convección natural o libre.

12
En primer lugar, se evidenció que el experimento fue exitoso debido a que a la hora de
determinar si se estaba ejecutando una acción de convección libre o combinada, el cociente
entre Grashof y Reynolds este último al cuadrado, da como resultado mayor a uno,
confirmando que la práctica si se ejecutó en condiciones de ser una convección natural o
libre.

En segundo lugar, se evidenció que cuando se contrastan los coeficientes de película teóricos
y prácticos exteriores se encuentra que tanto en el dato 1, como en el dato 2 no existe una
diferencia mayor al 15%, es decir, que se asemejan los datos, sin embargo cuando se van a
contrastar los datos entre los coeficientes de película interiores, los datos si tienen diferencia
y esto se debe a que el calor absorbido del agua es muy bajo y eso hace que el coeficiente sea
mucho menor el práctico al teórico.

Por último, se evidencia que a la hora de generar la gráfica de los resultados del número de
Nusselt en función de Grashof y Prandtl, se encuentra que debido a que los tres resultados
encontrados están en el mismo régimen, es decir, son turbulentos, usan la misma ecuación
para determinar Nusselt, dado eso cuando se hace la tabulación y se busca la ecuación que
rige a los resultados, entrega que C= 0.1 y n= 0.33, qué es lo mismo que cuando se calcula
con la ecuación inicial.

CONCLUSIONES

1. El número de Nusselt tiende a aumentar con el incremento del número de Reynolds y

el número de Rayleigh.
2. El coeficiente de película aumenta con el aumento del número de Nusselt y Reynolds.
3. El flujo turbulento está relacionado con las perturbaciones del fluido, la velocidad con
la que se mueve el fluido y el medio que lo rodea.

Bibliography
DeWitt, D., Incropera, F., Bergman, T., & Lavine, A. (2007). Fundamentals of Heat and Mass
Transfer (6 ed.). Hoboken, NJ, United States of America: John Wiley & Sons, Inc.
GARAVITO, N. G. (2009). OPERACIONES UNITARIAS MANUAL DE PRÁCTICAS. Bogotá:
Universidad Nacional de Colombia.

13