Trnaferencia de Calor 5.
Trnaferencia de Calor 5.
Trnaferencia de Calor 5.
Ingeniería Mecánica
6.-Sobre una esfera de catalizador de 1 cm de diámetro tiene lugar la reacción en fase gaseosa A ?
2B. Cuando se coloca la esfera en un recinto de grandes dimensiones, lleno del gas A, al alcanzarse
régimen estacionario se observa que la velocidad de desaparición de componente A es de 1.5 10-6
mol/s, siendo su composición sobre la superficie de la esfera del 90% en moles. Sabiendo que la
concentración total es de 40 mol/m3, calcular la difusividad de la mezcla A+B, supuesta constante.
Partir de las ecuaciones de continuidad y deducir la ecuación correspondiente indicando el
significado de todas las variables.
12.- / / / / / / / / / / / / / / / / / / / / / / /
13.- Se evapora agua de una presa de 10Km2 de superficie, la difusión (DAB=1ft2/h) a través de
una capa de .15cm de aire estancado a 20°C y 1 atm. Estime la caída del nivel de agua en la presa
si el aire fuera de la capa tiene una humedad correspondiente al .25% de la humedad de
saturación (.0189lb de H2O/lb de aire) Cuantos metros cubicos de agua se evaporan diariamente?
15.- La pared de un horno está constituida por una capa de ladrillo refractario y otra, de doble
espesor, de ladrillo corriente. Si se sabe que la conductividad calorífica del ladrillo aislante es
menor que la del ladrillo corriente, ¿cuál de los siguientes perfiles de temperatura será el que
previsiblemente se establezca en la pared al alcanzarse régimen estacionario?
5.1.1 INTRODUCCION
INTRODUCCIÓN
La convección es una de las tres formas de transferencia se caracteriza porque se
produce por intermedio de un fluido (aire, agua) que transporta el calor
Capa Limite Hidrodinámica El movimiento de un fluido con relación a una superficie solida o
liquida determina la formación de una zona cercana a la superficie en la que se dan gradientes
fuertes de velocidad. Comprobado experimentalmente. Se indica que en esta zona el esfuerzo
cortante es elevado. Se conoce esta zona como la capa limite hidrodinámica:
Capa Limite Hidrodinámica Capa Limite Hidrodinámica Si un fluido se aproxima con una velocidad
v0 a una superficie solida en dirección paralela a la misma, la capa de fluido en contacto con la
superficie se adhiere a ella. En cualquier posición x a lo largo de la superficie de la placa habrá
capas del fluido colocados a cierta distancia de la misma, que “ no se han enterado” de la
existencia de la misma y conservaran su velocidad original v0 Capa Limite Hidrodinámica Entre
ambas capas fluidas, la velocidad nula adherida a superficie solida y la velocidad v0 se darán todos
los valores intermedios de la misma. Se define como capa limite la zona entre Vx=0 a Vx=0.99V0 .
Al aumentar el valor de x, aumenta el espesor de la capa hasta un punto en que el numero de
Reynolds local, Rex= xV0ρ/μ aprox 1.5-2 X 105, en este punto empieza aparecer turbulencia. Capa
Limite Hidrodinámica El criterio para saber el tipo de capa límite presente es la magnitud del
número de Reynolds, Rex, conocido como número local de Reynolds, basado en la distancia x del
borde de ataque. El número local de Reynolds se define como: Para un flujo que pasa por una
placa plana, los datos experimentales indican lo siguiente: Capa Limite Hidrodinámica (a) Rex < 2 X
105 la capa límite es laminar (b) 2 x 105
CONVECCIÓN FORZADA
1. Salvo indicaciones expresas en algunas correlaciones, en general se evalúan las propiedades del
fluido a la temperatura media de película, definida como:
2
2. El número de Reynolds (Re) es el parámetro que permite identificar el régimen de flujo en
convección forzada.
FLUJO EXTERNO
Flujo laminar
Churchill y Ozoe
/ ///
0.3387
1 0.0468/
Pr cualquiera // Pex > 100
Flujo turbulento
Analogía de Prandtl
.
Analogía de Chilton-Colburn
0.0292 . /
Zhukauskas
.
Jakob
Correlación aplicable a gases sobre geometrías no circulares.
/
Propiedades del fluido evaluadas a Tm El
valor de las constantes C y n se toman de la siguiente tabla:
-----------------------------------------------------Esfera -----------------------------------------------------
Whitaker
0,71 < Pr < 380 // 3,5 < Re < 7,6·104 // 1 < (μ/μs) < 3,2
Propiedades evaluadas a Tm excepto μs a Ts.
umax es la velocidad de paso máxima, que modela la velocidad que tendrá el fluido a su paso entre
los tubos; mayor que la del fluido libre a la entrada debido a la disminución de la sección de paso.
Asociado a esta velocidad se considera un número de Reynolds según:
Grimison
/
1.13 ,
Flujo laminar
.
• Flujo de temperatura superficial constante en zona de flujo completamente desarrollado
3.66 .
Hausen
Se emplea poco porque supone una determinada longitud de entrada térmica.
0.0668 /
3.66 1 0.04 / /
Sieder-Tate
/ . 1.86
2
Propiedades del fluido evaluadas a Tm excepto μs a Ts
Mijeev-Mijeeva
. . 1.4
(L/D) < 0.067·Pr5/6·ReD // 0.06 < (Pr/Prs) < 10
.
4
Flujo turbulento
Las siguientes correlaciones se adecúan a las condiciones de flujo de calor o de temperatura
superficial constante
Colburn
0.023/ .
Dittus-Boelter
.
0.023
0.3 si enfriamiento (Ts < Tm)
0.4 si calentamiento (Ts > Tm)
0,7 < Pr < 160 // Re > 10000 // (L/D) > 10
Propiedades del fluido evaluadas a Tm
Solo para diferencias de T moderadas:
| | 6 º para líquidos
| | 60 º para gases
Esta restricción permite utilizar la expresión como herramienta aproximada, evaluando las
propiedades a la temperatura del fluido sin perturbar, lo cual resulta de interés en aquéllos casos
en los que se desconozca la temperatura de la superficie interior del conducto.
Sieder-Tate
.
/ .
0.027
Petukhov
/8
. /
1
1.07 12.7 /8
0.5 < Pr < 2000 // 104 < ReD < 5·106
Propiedades del fluido evaluadas a Tm
El factor de fricción se obtiene del diagrama de Moody para tuberías rugosas y de la siguiente
expresión para tuberías lisas o de rugosidad suave:
0.79 1.64
3000 < ReD < 5·106
a. Para Pr > 0.7 y flujo turbulento pueden utilizarse las correlaciones de conductos circulares
teniendo en cuenta que:
b. La longitud característica es el diámetro hidráulico definido como:
4
Los miembros derechos de las ecuaciones (8-47) y (8-47a) se encuentran dentro de un rango de
3%. Se puede atribuir esta diferencia a la técnica de aproximación para la integral de la cual se
obtuvieron estas ecuaciones. Si ignoramos esta pequeña diferencia, podemos escribir
A la ecuación (8-48) se le conoce como la analogía de Colburn y demuestra que la pérdida de calor
debida a la fricción en el fluido y el coeficiente de transferencia de calor se encuentran
relacionados. Para muchas situaciones de transferencia de calor por convección, incluyendo
diferentes geometrías y circundando flujo turbulento, Colburn (referencia 12) demostró que esta
analogía es válida para un amplio rango de número de Prandtl (0.5 a 50), siempre y cuando no
exista separación de capa frontera. Sin embargo, se debe insistir en que existe un buen número de
situaciones en que falla la analogía de Colburn; por ejemplo, transferencia de calor en metales
líquidos (Pr ~ 0.01) y aceites pesados (Pr > 10³), y flujo de calor unitario uniforme para el flujo
laminar sobre una placa plana. Para Pr = 1, la analogía de Colburn se reduce a la analogia de
Reynolds que discutimos en la sección 8-4.
introducción
los coeficientes de transferencia de calor por convección generalmente son determinados a través
de experimentos que envuelven el uso de complejos instrumentos y mediciones difíciles. un
método alternativo para obtener el coeficiente de transferencia de calor por convección, o en
general estudiar procesos donde existe transferencia de calor, es conducir los experimentos a
través de la transferencia de masa. los métodos de medición basados en la transferencia de masa
en general reportan menores índices de incertidumbre y son relativamente más simples para el
montaje y para el control de los parámetros del experimento.
para expresar las ecuaciones de una forma compacta y general se utiliza la notación tensorial en
coordenadas cartesianas con parámetros adimensionales. los términos con unidades de longitud
son adimensional izados a través de una longitud de referencia (l0); las velocidades por medio de
una velocidad de referencia (n0); los términos de tiempo según un tiempo de referencia (l/n0) y
las presiones según la presión de referencia dada por (rv20). esto es:
Con esta notación las ecuaciones de la continuidad y de Navier-Stokes toman la siguiente forma:
(2)
Se debe notar que la ecuación (3) en realidad representa tres ecuaciones, una para cada dirección
cartesiana. El número de Reynolds Re0 que aparece en (3) está definido sobre la base de los
parámetros de referencia, i.e.:
Como un flujo con propiedades constantes es independiente del hecho que exista o no exista
transferencia de calor o de masa, se puede superponer a dicho campo de velocidades un proceso
de transferencia de calor o de masa sin producir ninguna alteración sobre él.
(4)
El número de Prandtl (Pr) expresa la relación entre la viscosidad cinemática y la difusividad térmica
del fluido, y de igual forma el número de Prandtl (Prt) expresa la relación entre la difusividad
turbulenta de momentum ε y la difusividad turbulenta
de calor εh . El número de Prandtl turbulento es dado por una relación del
tipo determinada por consideraciones de semejanza y por evidencias
experimentales.
En la ecuación (4) se puede observar que el campo de temperatura obtenido como solución tendrá
la forma:
(5)
Entonces, el número de Nusselt puede ser usado para describir la relación funcional (5) expresada
anteriormente; de ello se obtiene:
(7)
(8)
definida como:
(9)
;
Un análisis de las ecuaciones (4) y (8) manifiesta que las dos son similares; la ecuación (4) describe
el campo de temperatura de la misma forma que la ecuación (8) describe el campo de
concentración, como puede ser verificado de las soluciones adimensionales encontradas
para Nu y Sh. Mas para que exista analogía plena entre las dos ecuaciones y en consecuencia en
las soluciones, aún debe ser satisfecha la condición: Prt = Sct.
Un = cRe mPr n
De acuerdo con la analogía calor-masa, los resultados de transferencia de masa pueden ser
correlacionados en la misma forma:
Sh = cRe mPr n
La razón de estas dos ecuaciones indica que los resultados de la transferencia de masa pueden ser
convertidos en sus análogos de transferencia de calor:
(10)
(11)
Donde:
Stm y Stc son los números de Stanton para transferencia de calor y de masa respectivamente.
Sustituyendo estos en la ecuación (10) se obtiene:
jm = jc (12)
con y
Las expresiones de la analogía fueron presentadas en términos de números adimensionales
locales, pero pueden ser integradas y resultar expresiones semejantes
en términos de valores globales. La expresión de la analogía es válida para flujo laminar
o turbulento, siempre que Prt = Sct.
Lewis mostró que los coeficientes de transferencia de calor y de masa pueden ser relacionados de
forma más precisa siguiendo un tratamiento más estricto. Usando una expresión derivada del
perfil universal de velocidad de la capa límite turbulenta y asumiendo iguales la difusividad
turbulenta de calor y de cantidad de movimiento (Prt = 1), puede ser obtenido el coeficiente de
transferencia de calor de la capa límite turbulenta en términos del coeficiente de fricción y del
número de Prandtl.
Siguiendo esta línea de análisis fue demostrado que para flujos turbulentos simples, tales como el
flujo sobre placa plana o flujo turbulento totalmente desarrollado en tubo, existen discrepancias
con la analogía de Colburn. Según Colburn, la relación Φ = (St/Stm) debe tener un valor constante
dado que Pr y Sc se mantienen prácticamente constantes, sin embargo Lewis encontró que en
general el factor Φ disminuye con el incremento del número de Reynolds. Para ilustrar lo
expresado hasta aquí se presentan a continuación dos ejemplos:
Placa plana
Usando la correlación analítica de Von Karman para St se tiene:
(13)
(14)
(15)
(16)
A partir de estas relaciones se puede encontrar Φ como función del número de Reynolds, cuyo
comportamiento se muestra en la figura 1. Nótese de la figura (construida para flujo de aire Pr =
0,71 e intercambiando masa con una superficie de naftaleno Sc = 2,28) que Φ de Colburn es
constante e igual a 2,03 para todo Re, sin embargo el Φ hallado según Von Karman va
disminuyendo con Re, de forma tal que para Re = 107 la diferencia alcanza hasta un 25% de
diferencia con relación al valor de Colburn.
(17)
(18)
A partir de estas relaciones se hace el mismo análisis que para la placa plana y se obtienen los
resultados mostrados en la figura 2. Obsérvese en la figura que Φ se comporta nuevamente como
en el caso de la placa plana, solo que la disminución con el número de Reynolds es aún mayor,
alcanzando una diferencia de hasta 30 % para
Re = 106.
5.1.5 FORMULAS EMPÍRICA: DUCTOS, CILINDROS, ESFERAS, BANCOS DE TUBO Y EN METALES
LIQUIDO
5.2 CONVECCIÓN NATURAL
CONVECCIÓN NATURAL
1. La convección natural es inducida solamente por fuerzas de flotación, y por tanto, en ella no
existe una velocidad de convección forzada bien definida.
2. El número de Grashof (Gr) juega en la convección libre o natural el mismo papel del número
de Reynolds (Re) en convección forzada.
5.2.2 convección natural en : placa plana vertical, tubo vertical, placa inclinadas , placa
horizontales, cilindro horizontal y esferas.
• Mc Adams
Churchill y Chu
/
0.387
0.825 1 0.492/ //
Para régimen laminar se puede obtener mayor precisión utilizando la siguiente correlación:
0.67 /
0.68 1 0.492/ //
Ra ≤ 109
• Condición de Gebhart para cilindros verticales
Las correlaciones anteriores para placas planas verticales son aplicables a cilindros verticales de
altura L si se satisface la siguiente condición:
35
a. Para geometrías planas horizontales la dimensión característica que se debe emplear es:
í
b. En placas inclinadas (y en particular en el caso extremo de placas horizontales) las fuerzas de
flotación favorecen o impiden la transferencia de calor por convección natural en función de la
disposición de las superficies (inferior o superior) y de su temperatura (caliente o fría).
Por ello, es fundamental tener en cuenta estos aspectos a la hora de aplicar las correlaciones.
• Mc Adams
Ra C m
Cara superior caliente o inferior fría
104- 107 0.54 1/4
• Churchill y Chu
Para ciertos casos con placas inclinadas, es posible emplear la correlación de Churchill y Chu para
placas planas verticales modificando el valor de la constante de la gravedad por el de su componente
paralela a la superficie (g·cosΦ), siendo Φ el ángulo que forma la superficie con la vertical.
/
0.387
0.825 1 0.492/ //
• Morgan
RaD C n
-----------------------------------------------Esferas ----------------------------------------------
• Yuge
Nu 2 0.43 Ra/
1 < Ra < 105
• Churchill
0,589 · Ra/
2
Nu 1 0.469/Pr / /
Pr ≥ 0.7 // Ra ≤ 1011
---------------------- 5.2.3Convección natural en espacios cerrados ----------------------
a. La longitud característica es el espesor de la capa de fluido, L, entre las dos superficies
involucradas, que se encuentran a temperaturas T1 y T2.
b. La altura de la cavidad se denota como H
c. Los números adimensionales se definen según:
T* = (T1+T2)/2
Jakob
- Para cavidades verticales:
Tipo de fluido Pr Ra C n
Gas / líquido - ≤ 1700 1 0
Gas 0.5 - 2 1700 - 7000 0.059 0.4
Gas 0.5 - 2 0.21 1/4
7000 – 3.2·10
Gas 0.5 - 2 > 3.2·10 0.061 1/3
Líquido 1 - 5000 1700 - 6000 0.012 0.6
Líquido 1 - 5000 6000 – 3.7·104 0.375 0.2
Líquido 1 - 20 3.7·104 - 108 0.13 0.3
Líquido 1 - 20 > 108 0.057 1/3
/ /
Tipo de fluido Pr Ra C n
Gas / líquido 1 - 5000 6000 - 106 0.11 0.29
Gas / líquido 1 - 5000 106 - 108 0.40 0.2
0.228.
Gas o líquido // 0.7 < Pr < 4000 // 100 < Ra < 109
5.2.4 CONVECCIÓN MIXTA
El proceso de convección mixta que ocurre en una cavidad semiabierta, que contiene
un fluido que circula por esa cavidad, puede ser analizada como el resultado de dos
formas básicas de convección: la convección natural causada por fuerzas de empuje y,
la convección forzada debido a la circulación del fluido frío que entra en la cavidad con
una determinada distribución de velocidades.
La convección forzada ocurre cuando un fluido está sometido a una cierta distribución
de velocidades, forzando de esta forma la transferencia de calor entre el fluido y el
equipamiento o dispositivo en cuestión. La intensidad de la convección forzada es
caracterizada por el número de Reynolds (Re).
la convección mixta causada por las fuerzas de empuje y por la circulación del fluido
que ingresa en la cavidad rectangular La cavidad tiene una relación geométrica A=2,
siendo H la altura de la cavidad y L su ancho. La distancia adimensional de la región de
entrada y salida, es D=0,2. El fluido entra frío en con un perfil de velocidades
parabólico en la parte inferior de la cavidad. El modelo numérico considera que el
fluido sale en adiabáticamente (en la parte superior de la cavidad) con un perfil de
velocidades también parabólico. Los resultados numéricos fueron obtenidos utilizando
el método de elementos finitos, a partir de un programa computacional desarrollado
por los autores, en lenguaje FORTRAN. Para el presente estudio fue utilizada una malla
uniforme con 2970 elementos triangulares lineal. El programa computacional calcula el
número de Nusselt medio en la superficie caliente ( ) en función de diversos
parámetros geométricos y térmicos. Son presentados resultados de las distribuciones
de la función corriente adimensional ( ) y temperatura adimensional ().
Fig. 1: Esquema de la cavidad estudiada
Fig. 2: Malla de elementos finitos utilizada
DESARROLLO TEÓRICO
Adimensionalización
(1)
Ecuaciones Fundamentales
a) régimen no permanente;
b) flujo bidimensional y laminar;
c) flujo incompresible;
d) función de disipación viscosa no considerada;
e) las propiedades físicas del fluido ( , , cp, ), son constantes, excepto la densidad
en los términos de empuje;
f) sin generación interna de calor;
g) los efectos de compresibilidad no son considerados.
(2)
(3)
(4)
, y (5)
y los números adimensionales de Prandtl (Pr), Grashof (Gr), Reynolds (Re), Rayleigh
(Ra) y Richardson (Ri) están dados, respectivamente, por:
, ,
y (6)
El cálculo del número de Nusselt medio en la superficie caliente ( ) está dado por: