CONVECCIN LIBRE OBJETIVOS 1. GENERAL: Analizar las condiciones en las cuales se desarrolla la transferencia de calor por conveccin natural.

. Comprobar las ecuaciones desarrolladas para analizar este mecanismo. 2. ESPECIFICOS: Realizar los balances trmicos para cada ensayo. Calcular los coeficientes peliculares experimentales, para la condensacin del vapor de agua y para el agua. Calcular los coeficientes peliculares con base a ecuaciones propuestas en la bibliografa para el caso. Calcular los nmeros adimensionales de Un, Gr y Pr. Para el agua en cada ensayo. Graficar UN vs Gr y Un vs Re. Calcular las fuerzas de empuje o boyantes para cada ensayo.

PROCEDIMIENTO EXPERIMENTAL
Identificar las lneas de uso para la practica de conveccin libre dentro del equipo.

Graduar el nivel del tubo vertical que tiene el rebosadero, este tubo es el de en medio de los tubos que descienden de la cmara acumuladora, posteriormente quitar la tapa roscada que tiene esta cmara, para asegurar que el agua circulara a presin atmosfrica.

Purgar la lnea de vapor antes de conectarla al equipo de practica.

Abrir la vlvula de paso de agua fra, posteriormente abrir las vlvulas de paso primero del agua dentro del equipo y despus permitir el paso del vapor dentro del equipo

Esperar 15 minutos para que el equipo alcance el estado estable Tomar datos de temperaturas a la entrada salida del agua; as como tambin temperaturas de pared inferior y superior. Medir el peso de condensado y de agua cada 120 y 30 segundos respectivamente

Cuando los balances de calor para el agua y el condensado no presenten diferencias porcentuales superiores al 4% por lo menos para dos ensayos consecutivos, pasamos al segundo ensayo bajando o subiendo el nivel del tubo del rebosadero.

Repetir el procedimiento para dos ensayos mas

DEFINICIN DE VARIABLES T1= Temperatura de entrada del agua al equipo de conveccin C T2= Temperatura de pared inferior en el equipo de conveccin C T3= Temperatura de salida del agua del equipo de conveccin C T4= Temperatura de pared superior en el equipo de conveccin C

ALGORITMO DE CALCULO 1. Realizar los balances trmicos para cada ensayo. Los balances trmicos se realizan para el vapor y el agua respectivamente: Qvapor = mv' * vap Qagua = m L' * Cp * (T3-T1) Donde: Qagua y Q vapor = Calores ganado y cedido por el vapor (KJ/s) m v' y mL' =Flujo msico del vapor y del agua respectivamente. (Kg/s) Cp = Capacidad calorfica del agua (J/Kg s) El porcentaje de error debe ser menor o igual al 4% y se calcula de la siguiente manera: % de Error = (Qmayor - Qmenor) / Qmayor 2. Calculo de los coeficientes peliculares experimentales para el agua y el vapor.

El balance de calor para el agua viene dado por: Qagua = h * A * Tml

Donde: h= Coeficiente pelicular para el agua (J/s*m2*K) A= rea de transferencia de calor (m2) Tml = Delta medio logartmico para el agua y se define como sigue.

Tml = (T2 - T1) / Ln (T2/ T1) Donde: T1 = Tp - T1 T2 = Tp - T3 Tw= Temperatura promedio de la pared externa del tubo. =(T2 + T4)/2 DE y DI = Dimetro externo e interno del tubo (m) Tp=Tw-(2.3*Qagua/(L*2**K))*Log(DE/DI)

El balance de calor para el condensado viene dado por:

Qvapor = A * h * (Tsat - Tm) Donde: Tm = (T2 + T4)/2 Tsat = temperatura de saturacin del agua a la presin de Manizales (585 mmHg). Esta temperatura de saturacin la hallamos con las constantes de Antoine para el vapor del agua como se muestra a continuacin. Con la temperatura en C y la presin en Kpa Ln 78 = 16.262 - (3799.88/(T+226.34)) Tsat = 92.8C

3. Calculo de los coeficientes peliculares con base a ecuaciones propuestas para el caso.

Para el caculo de los coeficientes, primero calculamos los nmeros adimensionales Nu, Gr, Pr para cada ensayo. Los nmeros adimensionales se definen asi:

D3 * 2 * g * *T Gr= 2
donde: Tm = (T1 + T3) / 2 Tm1 = (T2 + T4) / 2 = viscosidad a Tm (N*s/m2) D = Dimetro interno (m) = densidad del fluido a Tm (Kg/m3) g = gravedad 9.8 m/s2 = Coeficiente de expansin trmica = 1/*(inf - ) / (T inf - T) = 1/T (hallado en tablas a una temperatura media de entrada y salida) T = Tm1 - Tm

Pr= C * / K
donde: C= Calor especifico a Tm (Kj/Kg K) = viscosidad aTm (N * S/m2) K= Conductividad trmica (N/ m*K)

Para conveccin libre dentro de tubos horizontales se utiliza la correlacin siguiente. Las correlaciones fueron encontradas en l Perry , en la pagina 10-12 del Tomo III. Para cilindros y placas verticales y con 1 < Pr < 40 Gr>109 Nu = 0.138 Gr0.36 (Pr0.175 - 0.55)

Gr<109 Nu = 0.683 Gr0.25 (Pr / (0.861 + Pr)) 0.25

Otras correlaciones como las halladas en el libro fenmenos de transporte del Bird.

Para 104 < Gr * Pr < 109

Nu = 0.59 (Gr * Pr)1/4

! Despus de haber definido el numero de Nusselt, conocemos que el numero se define as:

Nu = h * D / K Con esta ecuacin despejamos el coeficiente pelicular para el agua, siendo D el dimetro interior del tubo en metros.

4. Calculo de los coeficientes peliculares del condensado por ecuaciones empricas. Para la condensacin del vapor se encuentra gran variedad de correlaciones empricas, suponiendo que la condensacin que se presenta es del tipo de pelcula descendente de fluido. ! Con la ecuacin del Bird para tubos o paredes verticales.
hm = 4 K * *g *( ) (1) 3 3* *

donde: K = Conductividad trmica del fluido condensado a Tm1 (W/m k) = Densidad del fluido condensado a Tm1 (Kg/m3) = Viscosidad del fluido condensado a Tm1 (N/m*K) = Velocidad de flujo del condensado en el fondo de la superficie de condensacin por unidad de anchura de dicha superficie. Para un tubo vertical = W / *D siendo W la velocidad msica total de condensacin. Cuando se encuentran datos para el numero de Reynolds >350 se encuentra otro tipo de correlacin (tambin extradas del Bird pagina 13-31)

K3 *2 * g *(TdT0) *L hm = 0.003 ( * ) 3 *Hvap

Para un Reynolds < 2100 Re= 4 * / = W/( D * )

La ecuacin de Nusselt se puede escribir en las formas que siguen

Tipo Colburn:

h c * 5.35 = = c *G K 4*/
donde:
G= (3 * * / 2 * g)1/ 3

(2)

Tipo Nusselt :

L3 * 2 * g * L3 * 2 * g 1 / 3 h * L / K = 0.943 * ( ) = 0.925 ( ) K * * T *
(3)
5.

Calculo de la fuerza boyante o de empuje.

La fuerza boyante viene dad por: Fuerza boyante= g/ * (inf - ) = g * * (T3 - T1)

Donde: inf y = densidad del agua a T1 y T3 respectivamente.

TABLAS DE DATOS 1er Ensayo T1 T2 T3 T4 m vap' *103 0.3833 0.666 magua '*103 2.767 9.5 Qcond 0.8801 1.5154 Qagua 0.8745 1.4504 % Error 0.6406 4.3

15.3 93.1 90.8 92.4 15 92 85 93.3

2do Ensayo T1 T2 T3 T4 m vap' *103 magua'*103 7.675 7.483 Qcond 1.875 Qagua 1.8023 % Error 3.87 4.02

14.4 90.2 70.5 93.3 0.8166 14.7 90.2 71.2 93.4 0.8033 3er Ensayo T1 T2 T3 T4 94 93.8 m vap' *103 0.933 0.983 0.958

1.7697 1.7697

magua '*103 13.125 14.3 13.33

Qcond 2.1423 2.2571 2.2004

Qagua 2.0877 2.2567 2.1092

% Error 2.5 0.0019 4.14

14.2 87.4 52.2 14.3 87.3 52

14.4 87.1 52.2 93.8

INFORME 1. Calculo de los balances trmicos para cada ensayo. La forma de realizar los balances trmicos, las ecuaciones y el criterio para aceptar un conjunto de datos, se expusieron en el numeral 1 del preinforme. Los resultados de los balances trmicos se muestran en la tabla de datos para los tres ensayos, todos por debajo de 4% de error de desviacin entre el calor cedido y el calor recibido. 2. Calculo de los coeficientes peliculares experimentales para el agua. Se indica el calculo para el primer ensayo, el resto de clculos se resumen en la tabla subsiguiente. Tw= (93.1 + 92.9)/2 = 92.75C = 365.9K Tp=365.9-(2.3*0.8745/(0.6096*2**393))=365.899K T2 = 92.749 - 15.3 = 77.749C T1 = 92.749 - 90.8 = 1.949 C Tml = 20.56335K A= * DI * L = *0.01486 m * 0.6096m = 2.8456*10 -2 m2

h= Qagua / (A * Tml) = 0.8745 / (20.56335 * 2.8456*10-2) = 1.4943 KJ / s*m2*K

Tabla de resultados de los coeficientes experimentales para el agua.

h (KJ/K*s*m2)

Tw (K) 1er Ensayo

Tp (K)

Tml

1.4943

92.75

92.74

20.5633

1.6876

92.65 2do ensayo

92.64

30.2014

1.459 1.4535

91.75 91.8 3er ensayo

91.74 91.78

43.4101 42.786

1.3262 1.4353 1.3454

90.7 90.55 90.45

90.68 90.53 90.41

55.3217 55.2528 54.9604

3. Calculo de los coeficientes peliculares experimentales para el vapor. Necesitamos conocer la temperatura de saturacin del agua a la presin de 78 Kpa. Este calculo se realizo con la ecuacin de Antoine, para el vapor del agua con la temperatura en C. El resultado fue de 92.8 C. hvap = Qcond / A * (Tsat - Tm1) = 0.8801 / 0.3039 = =14.5231 KJ/m2*s*C

Para los otros ensayos los clculos se resumen en la siguiente tabla.

Qcond (Kj/s)

Tm1 (C) 1er ensayo

h (KJ/K*s*m2 )

0.8801 1.5154

93 92.65 2do ensayo

14.5231 33.2434

1.875 1.8445

91.75 91.8 3er ensayo

5.876 6.0694

2.1423 2.2571 2.2004

90.7 90.55 90.45

3.3568 3.3009 3.081

4. Calculo de los coeficientes peliculares con base a ecuaciones propuestas en la bibliografa para el caso.

Primero haremos el calculo de los nmeros adimensionales Nu, Gr, Pr. Para cada ensayo. ! Para el primer ensayo:

D3 * 2 * g * *T Gr= 2
0 .0149 3 * 987 .16 2 * 9 .8 * 0 .4712 * 10 2 * 39 .95 Gr = ( 0.528 * 10 3 ) 2

Gr = 2133114 .65742
! Para el primer ensayo:

Pr = C * / K
Pr = 4182 * 0 .528 * 10 3 / 0 .645

Pr= 3.4234
! Para el primer ensayo:

Para conveccin dentro de tubos se utiliza la correlacin encontrada en el Perry, referenciada en el preinforme.

Perry: Gr<109 Nu = 0.138 Gr0.36 (Pr0.175 - 0.55) Nu = 0.138 2133114.65740.36 (3.42340.175 - 0.55) Nu = 18.087

Bird: Gr*Pr =7302504.7182 Nu = 0.59 (Gr * Pr)1/4 Nu = 0.59 (7302504.7182) 1/4 Nu = 20.7935

Los resultados para el resto de ensayos, tanto para los nmeros adimensionales como el coeficiente pelicular se muestran en la siguiente tabla

Tm (C)

Tm1 (C)

* 103 (N*s/m2)

(Kg /m3)

C K Gr* 10-3 *103 (1/K) (J/Kg*K) (W/m*K) 1er ensayo

Pr

Nu Nu (Perry) (Bird)

53.05 50

93

39.95

0.5228 0.527

987.16 986.16

0.4712 0.4715

4182 4187

0.645 0.644

2133.114 3.4234 24.6783 2133.200 3.2021 23.6882

30.67 29.9

92.65 42.65

2do ensayo 42.45 42.95 91.8 91.8 49.35 48.85 0.631 0.577 991.08 989.1196 0.4004 0.4367 4179 4180 0.634 0.640 1580.246 4.1592 23.1033 29.8734 2028.999 3.7685 24.4849 31.025

3er ensayo 33.2 90.7 57.5 57.4 0.769 0.769 0.769 995.0248 995.0248 995.0248 0.3206 0.3206 0.3206 4178 4178 4178 0.620 0.620 0.620 1000.532 5.1821 20.7864 28.1537 1000.532 5.1821 20.7864 28.1537 1000.532 5.1821 20.7864 28.1537

33.15 90.55 33.3

90.45 57.15

5. Calculo de los coeficientes peliculares de acuerdo a ecuaciones de la bibliografa. Nu = h * D / K

Para el primer ensayo: por Perry h = 24.6783 * 0.645 / 0.0149 = 1068.2888 W/m*K Por Bird h= 1327.6611 W/m*K

Para el resto de los ensayos se resumen los clculos en la siguiente tabla

Nu Perry

K 1er ensayo

h Perry

H Bird

24.6783 23.6882

0.645 0.644

1068.2889 1023.8389 2do ensayo

1327.6611 1292.3221

23.1033 24.4849

0.634 0.640

983.053 1051.7004 3er ensayo

1271.1232 1332.6174

20.7864 20.786 20.7864

0.620 0.620 0.620

864.9374 864.9374 864.9374

1171.4962 1171.4962 1171.4962

6. Calculo de los coeficientes para el condensado a partir de ecuaciones empiricas

Tm1

K *103

*103

Re 1er ensayo

h (1)

h (2)

h (3)

92.75 92.65

677.4 677.4

962.83 962.83

0.3026 0.3026

0.02415 79.8083 0.0142 46.9266 2do ensayo

4210 4210

204.1130 6324.35 6344.1168 6063.17 126.368 8037.72 8062.8375 7237.29

91.75 91.8

677 677

963.39 963.39

0.306 0.306

0.05144 168.104 0.0506 168.04 3er ensayo

4209 4209

297.073 5213.0229 5232.4033 4696.66 297.083 5241.7115 5258.0919 4722.509

90.7 90.55 90.45

676.22 676.22 676.22

964.13 964.13 964.13

0.309 0.309 0.309

0.05877 190.1942 4207.8 323.7747 4967.227 0.0619 200.3236 4207.8 0.0604 195.4682 4207.8

4982.75

4403.59

335.17 4882.0517 4897.3081 4439.70 330.87 4939.139 4954.5737 4455.04

En la tabla anterior, las propiedades fueron tomadas a la temperatura Tm1, adems los coeficientes de pelcula estn dados en J/m2 *s C, habra que dividir los resultados por 1000 y as comparar con los datos calculados de forma experimental.

7. Calculo de la fuerza de empuje o boyante: Fuerza boyante= g/ * (inf - ) = g * * (T3 - T1) Los resultados para los diferentes ensayos se presentan en la siguiente tabla

inf a T1 (Kg/m3)

a T3 (Kg/m3)

*103 a tm

Fuerza boyante

1er ensayo 999 999 964.1714 964.1714 2do ensayo 999.49 999.49 977.335 977.335 3er ensayo 999.53 999.53 999.53 987.1668 987.1668 987.1668 0.3206 0.3206 0.3206 0.1194 0.1194 0.1194 0.4009 0.4009 0.2199 0.2199 0.471 0.71 0.3484 0.3484

CONCLUSIONES

La conveccin libre presenta 2 caractersticas que dominan su comportamiento fsico. Estos son las fuerzas viscosas y de diferencias

de densidad, estos ltimos son los que determinan las fuerzas boyantes, que permiten el flujo de un fluido ausente de otra fuerza opresora o impulsora. La conveccin libre tiene una velocidad de flujo menor que en la conveccin forzada, debido a la naturaleza en que esta se desarrolla as como la velocidad de flujo es pequea, las velocidades de transferencia de calor tambin lo son. La conveccin libre y su flujo se debe a las fuerzas horizontales dentro del fluido. Las fuerzas boyantes combinan los gradientes de velocidad y la fuerza gravitacional, o centrfuga en maquinaria de rotacin de fluidos. Podemos afirmar que los resultados experimentales no presentan gran desviacin de los calculados en correlaciones empricas, tanto para el agua dentro del tubo, como en la condensacin de vapores en la superficie externa vertical. En los resultados se encuentran desviaciones de los coeficientes convectivos para el primer ensayo, adems se anexan clculos adicionales de coeficientes peliculares por Colburn,Nusseltt y otras 2 ecuaciones encontradas por el Bird. Todos los datos son similares, para el segundo ensayo y tercer ensayo y no presentan desviaciones del experimental. Se pueden entrar a comparar los coeficientes para el agua y el vapor obtenindose grandes diferencias a comparacin con los resultados de conveccin forzada, esto es debido a que la naturaleza de la conveccin libre no incluye la mezcla de fluido, aumentando la resistencia al paso del calor. Para la condensacin del vapor se asume que loo realiza en pelcula, estos coeficientes tambin fueron calculados p RE<2100 , ya que ara para que una condensacin sea turbulenta el tubo debe tener un gran dimetro. Es claro en las graficas que la conveccin libre es una funcin mas fuerte del numero de Gr que del numero de Re. El numero de Gr es el

responsable de las fuerzas boyantes, que al mismo tiempo son las que rigen la conveccin libre.

CUESTIONARIO 1. CUANDO OCURRE LA CONVECCION LIBRE? Cuando Una superficie slida entra en contacto con un fluido de temperatura distinta a la de aquella superficie. Las diferencias de densidad proporcionan la fuerza del cuerpo que se requiere para desplazar el fluido; opuesta a la conveccin forzada la conveccin libre no utiliza una fuerza externa para el desplazamiento del fluido. Para que exista este fenmeno de conveccin libre, debe existir una diferencia de temperaturas critica que permita que la fuerza boyante supere la influencia de las fuerzas viscosas y as el fluido de temperatura mayor circulara en direccin en que disminuye la temperatura al mismo tiempo que s densidad u aumentara y el fluido fro circulara en direccin del aumento de temperatura al mismo tiempo que su densidad disminuye.

2. En que circunstancias es mas comn su desarrollo? En la conveccin libre las velocidades de flujo son generalmente mucho menores que las asociadas con la conveccin forzada. Por lo tanto las correspondientes velocidades de transferencia de calor son tambin mas pequeas que la forzada. En muchos sistemas que estn envueltos por muchas formas de transferencia de calor, la conveccin libre provee la mayor resistencia a la transferencia de calor y por lo tanto juega un importante rol en el diseo y desarrollo de sistemas en los que se desea minimizar la velocidad de transferencia de calor o minimizar los costos de operacin.

3. Pueden darse simultneamente la forzada y la libre? Si pueden presentarse. De la misma forma que hay una regin de transicin y no un punto de separacin entre el flujo turbulento y laminar debe haber

tambin una regin de transicin entre la conveccin libre de un fluido que permanece estacionario y el flujo laminar.

4. Cuando se presentan las dos? A bajas velocidades lineales ambas son indudablemente operacionales. La influencia de las corrientes de conveccin libre se disipan en la regin de transicin y de turbulencia debido a lo delicado de la naturaleza de las corrientes de conveccin libre. Los efectos de las fuerzas boyantes (particulares de la conveccin libre ) sobre los flujos forzados depende de la direccin de la fuerza boyante relativa al flujo forzado. Se estudian 3 casos. a) Fuerza boyante y flujo forzado en la misma direccin b) Direcciones opuestas c) Direcciones perpendiculares

El Nu de conveccin mezclada viene dado por: Nun = Nun f +- Nun n n=3, 7/2, 4 y el signo depende del flujo y de la direccin, expuesta en los tres casos. 5. Que determina el desarrollo de la conveccin libre? El desarrollo de la conveccin libre la determina la diferencia en temperaturas y densidad, el efecto neto se denomina fuerza boyante. En la mayora de los casos el gradiente de densidad esta unido al gradiente de temperaturas y la fuerza corporal o del cuerpo del fluido a la fuerza gravitacional.

Nota: Hay muchas aplicaciones, entre las cuales estn la transferencia de calor de calentadores basados en tableros elctricos o de radiadores de vapor hacia el aire de habitaciones y en la disipacin de calor de los anillos y serpentines de las unidades de refrigeracin hacia el aire circundante. Es tambin importante en las ciencias ambientales donde es el responsable de los movimientos ocenicos y atmosfricos Graficas
GRAFICAS:

Gr vs. Nu
2500000 2450000 Gr 2400000 2350000 2300000 1 2 Nu 3

Gr vs. Nu
2500000 2450000 Gr 2400000 2350000 2300000 1 2 Nu 3

BIBLIOGRAFA BIRD, Byron R y otros. Fenmenos de transporte, primera edicin, editorial Revert, 1978,Espaa. INCROPERA,Frank P y DEWITT David, Fundamentals of heat and mass transfer, cuarta edicin,Editorial Wiley,1996,Estados Unidos de Norteamrica. PERRY,Robert y otros. Manual del ingeniero qumico,tercera edicin en espaol,editorial McGrawHill,1992, Estados Unidos de Norteamrica. KERN,Donald Q. Procesos de transferencia de calor, primera e dicion, editorial Cecsa, 1997,Mexico.