Anexo A

TABLA DE CONTENIDO
Transferencia convectiva 2
Coeficientes convectivos 3
Números adimensionales 6
Flujos difusionales 6
Flujos convectivos: 9
Correlaciones para estimar coeficientes convectivos 12
Flujo externo 12
Placa Plana Horizontal 12
Flujo Transversal a Cilindros 21
Esferas 22
Lechos Empacados 23
Flujo interno 25
Pierce 25
Flujo Laminar en Tubos 26
Flujo turbulento en tubos circulares lisos 27
Ecuación de Colburn 27
Ecuación de Sieder y Tate, Whitaker, Petukov 28
Tuberías rugosas 28
Flujo Turbulento en conductos lisos no circulares 29
Factor de fricción 29
Flujo Turbulento de Metales Líquidos Dentro de Tubos Lisos 30
Flujo interno con cambio de temperatura 31
Convección natural 46
Placa vertical 46
Cilindros verticales 47
Cilindro horizontal 48
Placas Horizontales: 48
Convección natural desde esferas: 49
Convección natural y convección forzada combinadas 49
Influencia de la convección natural 50
Conductos horizontales 51
Flujo de fluidos sobre bloques de tubos 51
Flujo sobre bloques de tubos con tabiques deflectores 53
Condensación tipo película de Vapores Puros 54
Condensacion al interior de tuberias 56
Intercambio de calor radiante entre dos superficies negras 56
El factor de visión 69
Intercambio de calor entre dos superficies grises 71
Superficies rerradiantes 72
Radiación desde gases 73
Efectos de la forma geométrica 77
EJERCICIOS PROPUESTOS 81
Solución de la ecuación diferencial ordinaria lineal, de segundo orden, no homogénea: 83
2 TRANSFERENCIA MOLECULAR DE CALOR, MASA Y/O CANTIDAD DE MOVIMIENTO
ANEXO A
CONVECCIÓN Y RADIACIÓN
Transferencia convectiva
La transferencia de calor convectiva ocurre cuando un gas o un líquido fluyen sobre una superficie
sólida cuya temperatura es diferente a la del fluido. Se distinguen dos categorías de transferencia de
calor por convección, la convección forzada y la convección natural (o libre). En la convección
forzada, el movimiento del fluido es causado por un agente externo como una bomba o un ventilador.
En la convección natural, el movimiento del fluido es el resultado de las fuerzas de flotación
originadas en diferencias de densidad creadas por las diferencias de temperatura dentro del fluido.
A diferencia de la transferencia de calor por conducción, los problemas de transferencia de calor por
convección generalmente se resuelven por medio de correlaciones empíricas derivadas de datos
experimentales y análisis dimensional. La razón es que para resolver un problema de convección a
partir de los principios básicos, uno debe resolver la ecuación de cantidad de movimiento del fluido
junto con la ecuación del balance de energía. Aunque se han obtenido muchos resultados importantes
al resolver las ecuaciones fundamentales para los problemas de convección en los que el flujo es
laminar, no se ha conseguido lo mismo para el flujo turbulento.
Las correlaciones empíricas generalmente se expresan en términos de un coeficiente de transferencia

de calor, h, que se define por la relación:
Q  hAT (1.11)
En esta ecuación, Q es la velocidad de transferencia de calor entre la superficie sólida y el fluido, A

es el área sobre la que se produce la transferencia de calor, y T es una diferencia de temperatura
característica entre el sólido y el fluido. La ecuación (1.11) a menudo se denomina Ley de Newton del
enfriamiento, aunque Newton tuvo poco que ver con su desarrollo, y no es realmente una ley física.
h h 2 2
Es simplemente una definición de la cantidad . Las unidades de son W/m K o Btu/hr.pie .°F .
La ecuación (1.11) puede parecer similar a la Ley de Fourier de conducción de calor. Sin embargo, el
coeficiente, h , es un tipo de entidad completamente diferente de la conductividad térmica, k, que
aparece como la constante de proporcionalidad en la Ley de Fourier. En particular, h no es una
propiedad material. Depende no solo de la temperatura y la presión, sino también de factores tales
como la geometría, el régimen fluidodinámico (laminar o turbulento) y, en el flujo turbulento, la
intensidad de la turbulencia y la rugosidad de la superficie sólida. Por lo tanto, el coeficiente de
transferencia de calor no debe considerarse como una cantidad fundamental, sino simplemente como
un símbolo a través del cual se implementan los métodos empíricos.
Desde el punto de vista de la transferencia de calor, el flujo turbulento es altamente deseable. En

general, las tasas de transferencia de calor se pueden ordenar de acuerdo con el mecanismo de
transferencia de calor de la siguiente manera:
conducción < convección natural < convección forzada laminar < convección forzada turbulenta
La razón por la cual el flujo turbulento es tan efectivo en la transferencia de calor es que los
remolinos turbulentos pueden transportar rápidamente fluido de un punto a otro. Cuando esto ocurre,
la energía térmica del fluido se transporta junto con el fluido en sí. Este mecanismo de transporte de
2
3 COEFICIENTES CONVECTIVOS
remolino es mucho más rápido (típicamente, aproximadamente dos órdenes de magnitud) que el
mecanismo de transporte molecular de la conducción de calor.
Las correlaciones de transferencia de calor presentadas en este capítulo son válidas para la mayoría de
los fluidos newtonianos comunes. No son válidos para metales líquidos o fluidos no newtonianos. Se
requieren correlaciones especiales para este tipo de fluidos.
Coeficientes convectivos, definición: Estos coeficientes están asociados con el flujo en una interface.
Analizaremos el caso entre una interface sólido-fluido. En primer lugar debemos recordar que el flujo
total de cualquier entidad comprende la suma de la cantidad difusiva y la cantidad convectiva o de
arrastre. Así, refiriéndonos a la figura A.1 para el transporte de energía térmica en la dirección " z " se
debe cumplir, para cualquier elemento de superficie  xy  que:
dT
qz z 0
 k   CPT Vz z 0
dz (A.1)
En la superficie sólida, debido a la ley de no resbalamiento, la velocidad del fluido respecto a la

superficie sólida es cero, por lo que esta densidad de flujo solo tendrá el término difusivo. Para
resolverla es necesario conocer el perfil de temperaturas que usualmente está ligado al perfil de
velocidades y por tanto solo en condiciones muy limitadas se podrá resolver analíticamente. Por esta
razón se introduce el valor empírico del coeficiente convectivo de la ley de Newton del enfriamiento:
dT
qz z 0
 k  h TS  TF 
dz z 0 (A.2)
Figura A.1
Como el gradiente en el origen no se conoce, h tampoco. Se postula entonces la existencia de un

perfil lineal equivalente de temperaturas que se extiende desde la interface hasta la corriente principal.
Esta abstracción sirve para definir la teoría pelicular, se usa para definir el espesor z F de una película
de fluido efectiva que representaría la misma resistencia a la transferencia de calor que la presentada
por el sistema actual. Refiriéndonos a la figura A.1, se observa que la derivada al perfil de
temperaturas en la superficie se puede representar por
T TF  TS

z z 0 zF (A.3)
3
La corriente de fluido tiene velocidad característica V , temperatura TF , concentración c AF . La

superficie tiene temperatura TS y concentración c As .
Reemplazando en (A.2) se obtiene h  k / z F correlacionando el espesor de la película imaginaria con

el coeficiente convectivo. Por esta razón estos coeficientes también se denominan coeficientes
peliculares. Es de anotar que z F tampoco se conoce, y cualquiera de los dos deberá encontrarse
experimentalmente.
La cantidad h es una función de las propiedades del fluido, la geometría y la rugosidad de la

superficie, y el patrón de flujo del fluido. Existen varios métodos para evaluar h: los métodos
analíticos se pueden utilizar para algunos sistemas en flujo laminar; pero para los sistemas turbulentos
se usan métodos integrales, teoría de longitud de mezcla, y el análisis dimensional y claro está, la
experimentación.
Figura A.2
Transferencia de masa convectiva. Vamos a ilustrar este término haciendo uso de la Figura A.2. Esta
figura representa el flujo turbulento de un gas o líquido más allá de un límite líquido o sólido que se
muestra sombreado a la izquierda. Ese límite puede ser la pared de un conducto, o la interfaz que
separa dos fases. Se supone que la transferencia de masa ocurre desde una concentración c AS del
límite a una concentración más baja c AF en la mayor parte del fluido que fluye. Esto puede ocurrir si
el límite consiste en una substancia soluble o si un líquido volátil se evapora en una corriente de gas
que fluye, así como los procesos inversos de condensación y cristalización. Tenemos pues como
ejemplo cuatro tipos de transferencia de masa de una sola película: (a) evaporación; (b) condensación;
(c) disolución; (d) cristalización.
En los cuatro casos mencionados, las concentraciones y las presiones parciales en la fase fluida están
en equilibrio con la fase condensada vecina. Esta condición se denota con un asterisco. Por lo tanto,
p* es la presión de vapor de equilibrio del líquido, y C* es la solubilidad en equilibrio del sólido.
Evaporación Condensación Disolución Cristalización
4
La transferencia de masa tiene lugar inicialmente a través de una subcapa laminar, o capa límite, que
se encuentra inmediatamente adyacente a la interfaz. La transferencia a través de esta región, también
conocida como una "capa sin agitación" en aplicaciones biológicas, es relativamente lenta y
constituye la parte preponderante de la resistencia al transporte masivo. A esta capa le sigue una zona
de transición donde el flujo cambia gradualmente a las condiciones turbulentas que prevalecen en la
mayor parte del fluido. En el cuerpo principal del fluido vemos paquetes macroscópicos de fluido o
remolinos que se mueven rápidamente de una posición a otra, incluida la dirección hacia y desde el
límite. La transferencia de masa tanto en la zona de transición como en la región completamente
turbulenta es relativamente rápida y contribuye mucho menos a la resistencia global al transporte que
la subcapa laminar.
Se nota además que con un aumento en la velocidad del fluido hay un aumento concomitante en el
grado de turbulencia y los remolinos pueden penetrar más profundamente en las capas de transición y
límite. En consecuencia, el último disminuye en espesor y la velocidad de transporte experimenta un
aumento correspondiente en la magnitud. Por lo tanto, altos caudales significan un mayor grado de
turbulencia y, por lo tanto, una transferencia de masa más rápida.
Las concentraciones en el régimen turbulento generalmente fluctúan alrededor de un valor medio.

Estas fluctuaciones no se pueden cuantificar fácilmente. Esto se puede superar postulando la
existencia de un perfil de concentración lineal equivalente que se extiende desde el límite hacia la
masa de fluido principal. Este postulado está resumido en el concepto conocido como teoría de la
película, y la dimensión de la película en cuestión se denomina espesor efectivo de la película,
denotado como zF en la figura A.2.
Veamos cómo se puede cuantificar este concepto en una ley de flujo. Comenzamos con la ley de Fick
aplicándola al gradiente constante de la teoría de película obtenemos:
dc A  c  cAF   k c  c
N A   DAB  DAB AS c  AS AF 
dz zF
La relación entre la difusividad y el espesor de la película DAB / z F se fusiona en un único término

denominado coeficiente de transferencia de masa kc .
Es de interés tener una idea del orden de magnitud del coeficiente de transferencia de masa y del
espesor de película asociado. Esto parecería ser una tarea imposible, dada la amplia gama de
condiciones de flujo, configuraciones geométricas y propiedades físicas encontradas en la práctica.
Sin embargo, podemos llegar a algunas estimaciones razonables de los límites superiores e inferiores
a pesar de esta diversidad. Esto se debe a tres factores: Primero, es una práctica de ingeniería común
asociar el límite superior del flujo turbulento normal con velocidades del orden de 1 m/s en el caso de
los líquidos y de 10 m/s para los gases. Esto se aplica a los sistemas industriales (flujo de tuberías y
ductos), así como dentro de un contexto ambiental (viento, flujo del río) y se mantiene incluso en
casos extremos. Los vientos con fuerza de huracán, por ejemplo, pueden alcanzar hasta 30 m/s, pero
todavía están dentro del orden de magnitud citado.
En segundo lugar, las difusividades para una amplia gama de substancias son sorprendentemente
constantes. Se agrupan, en el caso de los gases, alrededor de un valor de 10-5 m2/s y para líquidos de
alrededor de 10-9 m2/s. En tercer lugar, si nos limitamos al flujo sobre un plano como una
configuración representativa, se encontrará que los coeficientes de transferencia de masa varían
5
inversamente con la potencia 2/3 del número de Schmidt, Sc, que también es bastante constante. Es
del orden 1 para gases, y unos tres órdenes de magnitud mayor para el transporte de solutos de
tamaño modesto en líquidos. Basándose en la correlación para coeficiente promedio en flujo
turbulento sobre placas planas:
k Re L  106
St  c  0.036 Re L1/5 Sc 2/3
V
Y tomando nota de la dependencia extremadamente débil del número de Reynolds, Re, obtenemos las
siguientes estimaciones de orden de magnitud para la transferencia de masa de flujo turbulento:
2
Para gases: kc : 10 m/s
5
Para líquidos: kc : 10 m/s
Para el espesor efectivo de la película, z FM  DAB / kc , los valores correspondientes son:

Para gases: z FM : 1 mm
Para líquidos: z FM : 0.1 mm¨
Los valores de kc dados aquí representan el orden de magnitud de los límites externos de lo que se
puede lograr, es decir, la tasa máxima de transferencia de masa obtenible o el tiempo mínimo
requerido para lograr la transferencia de una masa dada hacia o desde una superficie plana.
La transferencia de masa por difusión molecular reside en el otro extremo del espectro. Produce la
menor tasa de transferencia de masa posible y establece un límite superior en los requisitos de tiempo.
El derrame de disolvente considerado en la evaporación de un solvente, por ejemplo, requiere varios
días para completar la evaporación por difusión en aire estancado. Los mismos datos aplicados al
flujo de aire turbulento a la misma temperatura arrojan una estimación de varios minutos, inferior en
tres órdenes de magnitud. Por lo tanto, se puede ver el factor de 1000 como separando los dos
extremos de la transferencia de masa difusiva y turbulenta.
Problema de práctica, evaporación de un derrame de solvente: Se derrama una carga de solvente

sobre un área grande, dando como resultado una capa de líquido con una densidad de área de 1 kg/m2.
El solvente tiene una presión de vapor de 10 kPa y una masa molar promedio de 100. La difusividad
se estima en 10-5 m2/s. (a) Calcule el tiempo de evaporación en ausencia de corrientes de aire. Este es
el máximo tiempo esperado. Compare para un flujo de aire de 5 m/s paralelo a una longitud de 1 m.
Respuesta (a): 136 h
NÚMEROS ADIMENSIONALES
Para planificar exitosamente la fase experimental se parte inicialmente del análisis dimensional, lo
que reduce el número de variables pues bastará con variar las cantidades adimensionales y no todas
las variables involucradas. Los principales números adimensionales que surgen en este análisis se
pueden obtener de las definiciones de flujos convectivos y difusionales.
Flujos difusionales: Hasta ahora se han usado los mismos modelos básicos para desarrollar las
siguientes leyes de flujo (ecuaciones constitutivas o expresiones fenomenológicas) para el transporte
de energía, masa y cantidad de movimiento:
Energía Ley de Fourier.
6
T
qz   k
z
Materia c A Ley de Fick
J Az   DAB
z
Cantidad de movimiento V Ley de Newton de la
unidireccional  zx    x viscosidad.
z
En cada caso las ecuaciones toman la forma:
Densidad de Flujo = (Propiedad de Transporte)x(Gradiente de Potencial)

Donde k , DAB y se llaman las “propiedades de transporte moleculares”, y T , c A y Vx son los
“potenciales”. La temperatura y la concentración son escalares; la velocidad es un vector.
Aunque estas ecuaciones son similares, ellas no son completamente análogas debido a que las
propiedades de transporte tienen unidades diferentes. Notando que las dimensiones de la difusividad
másica son [longitud al cuadrado/tiempo], podemos definir difusividades para calor y cantidad de
movimiento como
Difusividad Térmica:   k /  CP , donde CP es capacidad calorífica a presión constante.

Difusividad de cantidad de movimiento:    /   llamada también viscosidad cinemática.
Suponiendo
 y C constantes, reescribimos las leyes de flujo como:
P
Tabla A.1
Energía Térmica  (  C PT ) Fourier A.4a
qz  
z
Masa de A c A Fick A.4b
J Az   D AB
z
Cantidad de Movimiento  ( v x ) Newton A.4c
 zx  
z
 A
j Az   D AB
En el caso de usar concentraciones másicas la ley de Fick se escribe z
Transporte de cantidad de movimiento x en la dirección y por difusión
7
Transferencia de calor por conducción Transferencia de masa por difusión
Notemos que   CPT  tiene unidades de energía por unidad de volumen o concentración de energía
por analogía con c A (moles de A por unidad de volumen) o  A (masa de A por unidad de volumen).
Además  Vx  tiene dimensiones de cantidad de movimiento por unidad de volumen y puede
interpretarse como concentración de cantidad de movimiento. Aquí las leyes de flujo están escritas en
la forma difusional:
Flujo =  Difusividad x Gradiente de Concentración.
Podemos generalizar esta similitud y escribirla como  mz      / z  , siendo  mz cualquiera de

las densidades de flujo difusivas,  cualquiera de las difusividades y  cualquiera de las
concentraciones. Aunque en este capítulo no vimos transferencia convectiva neta de calor ni de
cantidad de movimiento, observamos que cuando hay transferencia de masa, este hecho puede
generar una densidad de flujo convectiva neta que representamos por  AVz . Podemos extender la
analogía para cubrir esta situación diciendo entonces que las densidades de flujo pueden tener
componente difusiva  mz y componente convectiva  cz  Vz . Esto lo podemos extender para un
sistema multidireccional y es de mucha utilidad en los balances generalizados.
Como las difusividades poseen las mismas dimensiones, su relación dará cantidades adimensionales:
   CP  CP
Pr   
Número de Prandtl:   k k
 
Sc  
Número de Schmidt: DAB  DAB
 k Sc
Le   
Número de Lewis: DAB  CP DAB Pr
Estas cantidades aparecen en situaciones donde hay transporte simultáneo de calor e impulso; masa e
impulso; o calor y masa respectivamente.
Hemos visto que las leyes físicas que rigen el transporte de masa, energía e impulso, así como la de la
electricidad, se basan en la noción de que el flujo de estas entidades es inducido por un potencial de
conducción. Esta fuerza impulsora se puede expresar de dos maneras. En el caso más general, se toma
como gradiente o derivada de ese potencial en la dirección del flujo. En la Tabla A.1 aparece una lista
de algunas leyes de tasas basadas en dichos gradientes. En el segundo caso, más especializado, el
8
gradiente se considera constante. La fuerza motriz se convierte simplemente en la diferencia de

potencial sobre la distancia recorrida. Es el caso de los coeficientes convectivos.
Los flujos convectivos:

Cantidad de movimiento:  s  12 fV  V  0  (A.5a)
Calor h (A.5b)
qs    CPTs   CPT 
 CP
Masa N As  kc  c As  c A  (A.5c)
Estos flujos de interface están dados como [coeficiente (m/s) por concentración]. Por tanto, la
relación de los coeficientes es adimensional así (se exceptúa f que ya es adimensional):
h
St H 
Número de Stanton para transferencia de calor  CPV (A.6a)
kc
StM 
Número de Stanton para transferencia de masa V (A.6b)
La relación entre los flujos convectivos y difusionales también genera cantidades adimensionales. Los
flujos difusionales, ecuaciones (A.4) pueden entenderse como
Flujo molecular = (Difusividad)(concentración)/(Longitud característica) (ecuaciones A.4)
Flujo de interface = (coeficiente)(concentración) (ecuaciones A.5)
La relación (A.5)/(A.4) dará (coeficiente)(longitud característica)/(difusividad), obteniendo
Transferencia de calor, número de Nusselt
hLcar
Nu 
k (A.7a)
Transferencia de masa, número de Sherwood
kc Lcar
Sh 
DAB (A.7b)
Debe tenerse presente que la Lcar se refiere a una cantidad característica para el flujo, que en tubos es
el diámetro interno si el flujo es interno y el diámetro externo si el flujo es externo al flujo. Depende
de cómo el fluido "ve" al sólido. Así mismo, para esferas es el diámetro de la esfera, inclusive en el
caso de lechos empacados. Solamente para flujo paralelo a placas coincide con la longitud de la placa
en la dirección de flujo. También es importante tener en cuenta que las propiedades de transporte, k y
DAB son las del fluido, de tal manera que tanto Sh como Nu (y los St) pueden considerarse
coeficientes convectivos adimensionales, muy diferentes del número de Biot que es una relación de
resistencias en dos medios distintos.
9
Como la mayoría de los problemas no tienen solución teórica, gran parte del esfuerzo de ingeniería se
ha hecho en la experimentación, cuyos resultados se presentan en forma de correlaciones. Estas
correlaciones no son generales sino específicas a geometrías, flujos y propiedades de las sustancias.
Además tienen un margen de error que puede ser del orden de ±15% o inclusive más dependiendo de
las condiciones de la experimentación y la apropiada similitud con el proceso al que se quiera aplicar.
Estas correlaciones están dadas en términos de cantidades adimensionales. Por ejemplo, las
correlaciones usadas para determinar el factor de fricción, el coeficiente de convección para
transferencia de masa o de calor para convección forzada se expresan generalmente en la forma:
f  f  Re 
Nu  f1  Re, Pr 
Sh  f 2  Re, Sc 
Re representa el número de Reynolds, relación de fuerzas inerciales y viscosas en un flujo:
Vcar Lcar
Re 
 (A.7c)
La velocidad característica Vcar y la longitud característica Lcar se definen en cada geometría.
Para convección natural, donde el movimiento se genera por gradientes de densidad originados a su
vez por gradientes de temperatura o de concentración, el número de Reynolds es sustituido por el
número de Grashof:
 g  2 L3  T 
Gr 
2 para transferencia de calor (A.8a)
 es el coeficiente de expansión volumétrica para el fluido:
1  V 
  
V  T  P
V
= volumen específico 
V  1  d  V    dV  Vd   0 . Por lo tanto:
1    1     
    
  T  P  T  T 
Si este fluido puede considerarse gas ideal entonces   1 / T
Para transferencia de masa definimos un “coeficiente de expansión” originado en los gradientes de

concentración:
10
1     
 AB  
   A   A 
Consecuentemente definimos un número de Grashof para transferencia de masa:
gL3  2   g      S  L3c
GrAB   1 =
 2 S  S 2 (A.8b)
Las correlaciones adimensionales toman la forma
Nu  f 3  Gr , Pr 
Sh  f 4  GrAB , Sc 
Frecuentemente se usa el número de Rayleigh:
g  TS  T  L3car
Ra  Gr Pr 
 (A.9c)
Recordar que   k /  CP y    /  .
De los resultados experimentales se puede observar que en general, para la misma geometría y
régimen fluidodinámico, f1 y f 2 son similares de tal forma que se puede obtener por ejemplo una
correlación para Sh de una correlación experimental para Nu intercambiando en la misma Nu por Sh
y Pr por Sc, o a la inversa. Ocurre igual con los resultados para f 3 y f 4 . Esto da origen a las llamadas
analogías entre transferencia de masa y calor, ocasionalmente siendo posible extenderlas también a
transferencia de cantidad de movimiento.
Es importante destacar la relación existente entre diferentes grupos adimensionales:
Nu
St H 
Re Pr (A.6c)
Sh
StM 
Re Sc (A.6d)
Esto, aunado a la similitud encontrada en los resultados experimentales da origen a una correlación
importante y sencilla conocida como la analogía de Chilton y Colburn:
J H  St H  Pr   StM  Sc 
2/3 2/3
 JM (A:9)
Estos productos se conocen como coeficientes tipo J . Cuando la analogía se puede extender a
transferencia de cantidad de movimiento (no debe ocurrir separación de la capa límite),
J H  JM  f / 2 = fD / 8 (A.9a)
11
El factor de fricción de Darcy f D . En la literatura Europea puede encontrarse un factor de fricción

tipo J .
CORRELACIONES PARA ESTIMAR COEFICIENTES CONVECTIVOS
A continuación, se escriben correlaciones útiles para encontrar los coeficientes convectivos en

algunas de las geometrías y flujos más corrientemente encontrados. Las ecuaciones presentadas en
esta sección, así como las correlaciones similares que aparecen en la literatura, no son muy precisas.
En general, no se debe esperar que la precisión de los valores calculados sea mejor que  25% a 
30% cuando se usan estas ecuaciones. Esta limitación debe tenerse en cuenta al interpretar los
resultados de los cálculos de transferencia de calor.
FLUJO EXTERNO
Muchos problemas de interés de la ingeniería implican la transferencia de calor a los fluidos que
fluyen sobre objetos como tuberías, tanques, edificios u otras estructuras. En esta sección, se
presentan las correlaciones para varios de estos sistemas a fin de ilustrar el método de cálculo.
Placa Plana Horizontal
El flujo sobre una placa plana se ilustra en la Figura A.3. La velocidad del fluido no perturbado y la
temperatura aguas arriba de la placa son V y T , respectivamente. La temperatura de la superficie
de la placa es TS y L es la longitud de la placa en la dirección del flujo. El fluido puede fluir sobre
uno o ambos lados de la placa.
Figura A.3: Flujo en convección forzada sobre una placa plana
La situación se describe en la figura A.4; la longitud de la placa en la dirección de flujo es L y su

ancho, perpendicular al plano del dibujo, es w. El número de Re se define en función de la distancia
desde el borde de ataque donde x = 0 y cambia con la distancia desde allí:
V x
Re x 

Para el flujo de estas características, cerca de la superficie de la placa se desarrolla una capa límite
que según el valor de Rex puede presentar características de flujo laminar o turbulento. El valor crítico
para Reynolds se acepta comúnmente en Rec  V xc /  = 5x105.
12
Flujo laminar: Si la longitud total de la placa es menor que la longitud crítica, xc, toda la placa estará
en régimen laminar y se pueden aplicar las correlaciones para régimen laminar:
Coeficiente local de fricción
0.664 L 1.328
fx  f L  L1  f x dx 
Re x ; 0 Re L ;
Coeficiente local de transferencia de calor
hx x 1
Nu x   0.332 Re x2 Pr 3
1
k (A-10a)
Figura A.4: Capa límite fluidodinámica sobre una superficie plana horizontal
Coeficiente local de transferencia de Masa
k x x 1
Shx   0.332 Re x2 Sc
1
3
DAB (A-10b)
Flujo turbulento: En la zona donde el régimen es turbulento se aplican las correlaciones:
Factor de ficción local
fx = 0.058(Rex)0.2 ; fL = 0.072(ReL)0.2
Coeficiente de transferencia de calor local
hx x
Nu x   0.0292 Re0.8
1
x Pr
3
k (A-10c)
Coeficiente de transferencia de masa local

1
Shx  0.0292 Re 0x.8 Sc 3
(A-10d)
13
Los coeficientes de fricción f x y convectivos hx y kcx cambian con la posición, y si se quieren

utilizar para conocer la fuerza total o la cantidad total de calor o de masa transferidos desde la
superficie, deberán promediarse sobre la misma, teniendo en cuenta que un valor medio se define
como
x2
ym  x2  x1    y  x  dx
x1
Figura A.5
Si tenemos presente que la integral describe el área bajo la curva, ym se selecciona para que estas áreas
sean iguales, bajo la curva y(x) y bajo la recta en ym.
El coeficiente promedio de transferencia de calor se obtiene de las siguientes correlaciones:
L
hL  L1  hx dx
0
Régimen laminar: Si toda la placa está en flujo laminar, es decir si L  xc , entonces de (A-10a):
L 1/ 2
1 k  V x 
hL 
1
 0.332   Pr dx
3
L 0 x   
hL L 1 Re L  5 105 Pr  0.6
Nu L   0.664 Re L2 Pr 3
1
k (A.10e)
Flujo mixto: Si la longitud de la placa es tal que L  xc , una parte de la placa se hallará en régimen
laminar y la otra en régimen turbulento. Podemos encontrar un coeficiente promedio para este flujo
en régimen mixto así:
1  xc
     1/3
L
hm   0.332 k Re0.5 dx   0.0292 k Re0.8
x dx  Pr
L 
 0 x x xc x  (A.10f)
14
Nu m 
hm L
k
  1
 0.664 Re C0.5  0.0365 Re 0L.8  Re C0.8 Pr 3 
(A.10g)
0.6  Pr  60; 5x105  ReL  108.
Para escribir en EES: NuL=(0,664*Rec^0,5+0,0365*(Re^0,8-Rec^0,8))*Pr^(1/3)
Donde ReC es el Reynolds crítico; Re L es el Reynolds para toda la placa; si L  xC use solamente la
primera parte del paréntesis cuadrado (ecuación A.10e) tal como se desprende de las correlaciones ya
mostradas para capa límite en régimen laminar. Para números de Pr grandes el coeficiente 0.664 se
convierte en 0.678. En la literatura se encuentra una expresión simplificada útil solamente si el
Reynolds crítico toma el valor ReC = 500000:
Num   0.0365Re0.8
L  853  Pr
1/3
(A-10h)
En estas ecuaciones, todas las propiedades de los fluidos se evalúan a la temperatura de la película,
T f , definida por:
Tf 
Ts  T 
2
Los coeficientes de transferencia de calor calculados a partir de las ecuaciones (A.10e) hasta (A.10i)
son valores promedio para toda la placa. Por lo tanto, el flujo de calor entre la placa y el fluido está
dado por:
Q  hA Ts  T (A.12)
Aquí A es la superficie total para transferencia de calor por convección.
Para transferencia de masa cambiamos Nu por Sh y Pr por Sc, a saber:
kc L
 0.664 ReC0.5  0.0365  Re0.8
L  ReC   Sc
1
ShL  0.8
 3
DAB  (A.10i)
0.6  Sc  3000, 5x105  ReL  108.
Para flujo laminar de metales líquidos sobre placas planas:
Nu x  0.564 Pex0.5 Pr  0.05, Pex  100. (A.13)
Pe = Re.Pr es el número de Peclet.
Cuando el calentamiento no comienza en el borde de ataque sino a una distancia x0, el coeficiente
local adimensional Nux para flujo laminar sobre la placa está dado por:
15
Nu x x0  0
Nu x  1/ 3
  x0  3 / 4 
1   x  
    (A.14a)
y para flujo turbulento:
Nu x x0  0
Nu x  1/ 9
  x0  9 / 10 
1   x  
    (A.14b)
También es posible tener flujo uniforme de calor en la pared, en lugar de temperatura uniforme. Para
flujo laminar (Kays, W. M. y M. E. Crawford, Convective Heat and Mass Transfer, Mac Graw – Hill,
New York, 1980):
Nu x  0.453 Re 0x.5 Pr 1 / 3 Pr  0.6 (A.15a)
Flujo turbulento
Nu x  0.0308 Re 0x.8 Pr 1 / 3 0.6  Pr  60 (A.15b)
La presencia de una zona sin calentar puede nuevamente predecirse usando las ecuaciones (A.14).
Si se conoce el flujo de calor podemos determinar la variación de la temperatura con la posición.

Aquí parece importante determinar una temperatura promedia, más que un coeficiente promedio:
L
 L
qS  x q L
1
(TS  T ) m  (TS  T )dx   dx  S
L L  kNu x kNu L
 0
0 (A.15c)
NOTAS PARA TRANSFERENCIA DE MASA CONVECTIVA
Las expresiones obtenidas para los coeficientes de transferencia de masa aplicando la similitud con
los coeficientes para transferencia de calor tienen validez para Vs = 0, es decir, cuando la componente
“y” de la velocidad en la superficie de la placa vale cero. Sin embargo, si hay transferencia de masa,
algún componente de velocidad habrá allí, pues, analizando la ley de Fick se observa que:
nAs  j As   AsVs  j As  wAs  nAs  nBs 
Si nBs = 0, entonces Vs  nAs /   As / wAs   nAs /  s
16
Para tener en cuenta este componente “y” de la velocidad en la superficie, el coeficiente 0.332 de la
ecuación para Shx se modifica según sea el parámetro (VS / V ) Re x :
Tabla A-2
(VS / V ) Re x 0.60 0.50 0.25 0.00 -2.50
(Coeficiente) mS 0.01 0.06 0.17 0.332 1.64
El valor del coeficiente se puede aproximar en el intervalo por:
mS  0.332 – 0.528 (VS / V ) Re x
Para (VS / V ) Re x > 0.62 la capa límite se desprende y dejan de ser aplicables las ecuaciones
basadas en Vs = 0. La componente de la velocidad Vs aumenta el espesor de la capa límite
decreciendo nAs y  s .
Lienhard (2011, p 653) propone corregir los coeficientes obtenidos de las correlaciones para
transferencia de calor multiplicándolos por el "factor de soplado" de acuerdo al siguiente raciocinio:
nAs  j As   AsVs (i)
Suponiendo conocidos los perfiles de concentración en la capa límite el desarrollo difusional viene
dado por:
d  As  As
nAs   DAB   nAs  nBs 
dz 
En términos de fracciones másicas y haciendo nAs  nBs  n
dwA
nAs    DAB  wA ns
dz z 0 (ii)
con nAs constante, separamos variables e integramos
zF
 
wA
dwA
 dz    DAB 
0 wAs nAs  nwA
Llegando a:
 DAB RA  wA
n ln
zF RA  wAs (iii)
Aquí RA es la relación entre la densidad de flujo de A y la densidad de flujo total nA / n . Se anota la

similitud con el desarrollo que lleva a la ecuación (3.1) del Treybal.
17
Lienhard usa la ecuación (i) para definir un coeficiente de transferencia de masa para el caso corriente
en el cual se desconocen con precisión los perfiles de concentración en la corriente convectiva:
nAs  k    As   A    AsVs (iv)
En términos de fracciones másicas
nAs  wAs n  k    wAs  wA  (iva)
Dividiendo y multiplicando el lado izquierdo por n observamos que:
nAs  wAs n   RA  wAs  n ; despejando de (iva) se obtiene:
 w  wA 
n  k    As   k  B
 RA  wAs  (v)
Se le ha dado la abreviatura B al paréntesis cuadrado. Si el término de arrastre es despreciable, de (iv)

se obtiene
nAs  k    As   A  = k    wAs  wA  (vi)
y de (ii) obtenemos
 DAB
nAs   wAs  wA 
zF (vii)
Por lo tanto k   DAB / z F . Si al numerador de la fracción en (iii) le sumamos y restamos wAs se

obtiene:
 DAB  wAs  wA   ln 1  B  

n ln 1     k  B
zF  RA  wAs   B  (iiia)
Al dividir la ecuación (iiia) por la ecuación (v) se observa que los coeficientes convectivos obtenidos
a partir de analogías con transferencia de calor, donde el arrastre es despreciable deben multiplicarse
por un factor de corrección dado por el paréntesis recto de la ecuación (iiia):
ln 1  B 
BF 
B (A-11a)
Este símbolo se le da porque este término se conoce como "factor de soplado" (Blowing Factor en
inglés)
wA  wAS
B
wAS  RA (A-11b)
18
El término RA se obtiene de información adicional, no difusiva. Por ejemplo, en el caso de A

difundiendo en B estancado, RA = 1. Si B > 0 el espesor de la capa límite aumenta y el coeficiente de
transferencia disminuye. Para el caso contrario (succión) esta disminuye y el coeficiente de
transferencia se ve favorecido. Si estas fracciones másicas son pequeñas, B  0 y BF  1. Por
ejemplo, si B = 0.06, BF = 0.97 y la corrección es del orden de un 3 % lo que está dentro de la
tolerancia de los cálculos de ingeniería.
EJEMPLO A.1: Una corriente de aire a 600 mm Hg se usa para secar una serie de placas fotográficas
cada una de longitud Li = 0.25 m en la dirección de flujo de aire. El aire está seco y a temperatura
igual a la de las placas, T = TS = 50 C. La velocidad del aire es V = 9.1 m/s.
(a) ¿Cómo varía el coeficiente local de transferencia de masa con la distancia desde el borde inicial?
(b) ¿Cuál de las placas se secará primero? Calcule la velocidad de secado por metro de ancho para
esta placa en kg/s.m.
(c) ¿A qué velocidad deberá suministrarse calor a la placa que se seca más rápidamente para
mantenerla a TS = 50 C constante durante el proceso de secado?
Figura 1 ejemplo A.1
Datos: Para aire a 50 °C,  = 18.2x10-6 m2/s; a partir de tablas de vapor para agua a 50 °C, su
concentración de saturación es (inverso del volumen específico como gas)  As = 0.082 kg/m3; h fg =
2383 kJ/kg; a 25 °C, y 1 atm, DAB = 0.26x10-4 m2/s; a 50 °C y 600 mmHg es:
DAB = 2.6x10-5x(323/298)1.75(760/600) = 3.79x10-5 m2/s
La viscosidad cinemática es inversamente proporcional a la presión, por lo que se debe corregir:
 760  
  1.82  105    2.31  10
5
 600  m2/s Sc = /DAB = 0.61
Solución: Considerando Reynolds crítico Rec = 5x105, la longitud de transición es
Los coeficientes de transferencia locales varían con la distancia desde el borde de ataque así:
19
Figura 2 ejemplo A.1. Variación cualitativa de los coeficientes locales sobre una placa plana horizontal
Según esto los coeficientes mayores (y por tanto la mayor velocidad de evaporación) se deben
presentar sobre la primera o sobre la quinta filas de fotografías. Para la primera fila:
Como toda la zona está en régimen laminar adaptamos la ecuación (A.10e):
kc L kcm1
ShL   0.664 Re0.5
L Sc
1/3
DAB = 0.030 m/s
Para la quinta fila, toda en régimen turbulento, hallamos el coeficiente promedio a partir de (A.10d):
x5 0.8
1  DAB  V x 
kcm 5   0.0292    Sc1/3 dx
x  x4 x   
En esta ecuación x = 0.25 m, x4 = 1.0 m y x5 =1.25 m. Se obtiene:
kcm 5 x
 0.0365  Re50.8  Re0.8
4  Sc
1/3
DAB
Se obtienen los siguientes resultados:
Re5  625000 Re 4  500000 kcm 5  0.0333 m/s
Este valor es mayor que el coeficiente obtenido para la primera fila. Por lo tanto, la quinta fila es la
que se seca primero. La velocidad de evaporación es:
nA5  kcm 5  x  1  As  0  = 6.83x10-4 kg/s.m
El calor que se debe suministrar para mantener esta placa a 50 °C es el de evaporación:
Q5  nA5  h fg = 1626 W/m
Es de notar que la concentración másica en la interface se calcula normalmente a partir de la presión

TS
de vapor del líquido a la temperatura del mismo, PA , así:
20
TS
P MA
 AS  A
RT
Como el coeficiente se obtuvo de una correlación para transferencia de calor es importante verificar el
factor de soplado. Usando las ecuaciones (A.11) con wA  0 porque el aire está seco, y la densidad
del aire a 50 °C y 600 mmHg (80 kPa),  = 0.8624 kg/m3, B=0,1051 y BF=0,9509. La evaporación
disminuye a nAS = 6.5x10-4 kg/s.m, y el suministro necesario de calor a Q5 =1547 W/m, 5% menores.
Código para obtener kcm 5 :

kc5*dx/Dab=0,0365*(Re5^0,8-Re4^0,8)*Sc^0,33
dx=0,25: Dab=3,79e-5: Sc=0,48: Re4=5e5: Re5=V*L5/nu
nu=1,82e-5: V=9,1: L5=1,25
Flujo Transversal a Cilindros
El número de Nusselt promedio para convección forzada entre un cilindro de diámetro D y un fluido
que se desplaza perpendicularmente al eje del cilindro está correlacionado con los números de
Reynolds y Prandtl según Hilpert por
Figura A.6
(A.16)
hD VD
NuD  Re 
k 
Los valores de C y m dependen del número de Reynolds como se indica en la Tabla A-3.
Tabla A-3: Constantes para la ecuación (A.16) con barras cilíndricas
21
Rango de ReD C m
0.4 — 4 0.989 0.330
4 — 40 0.911 0.385
40 — 4000 0.683 0.466
4000 — 40000 0.193 0.618
40000 — 400000 0.027 0.805
Todas las propiedades deben calcularse a la temperatura de película.
Churchill y Bernstein proponen una ecuación que cubre todo el rango de ReD, Pr > 0.2:
4/5
0.62 Re 0D.5 Pr 1 / 3   Re D  
5/8
Nu D  0.3  1    

1  0.4 / Pr  
2 / 3 1/ 4
  282000   (A.17)
Para escribir en EES:
Nusselt=0,3+(0,62*Re^0,5*Pr^(1/3))/(1+(0,4/Pr)^(2/3))^0,25*(1+(Re/282000)^(5/8))^(4/5)
Esta correlación tiene validez para Re Pr  0.2 y las propiedades del fluido se evalúan a la
temperatura de película. La tasa de transferencia de calor viene dada por la ecuación (A.12) con A =
DL, donde L es la longitud del cilindro.
Para un gas que fluye perpendicularmente a un cilindro que tiene una sección transversal cuadrada
(como un conducto de aire), se aplica la siguiente ecuación [Serth p 63]:
Nu  0.102 Re0.675 Pr1/3 5000  Re  105
Los números de Nusselt y Reynolds se calculan como para un cilindro circular pero con el diámetro,
D, reemplazado por la longitud de un lado de la sección transversal cuadrada. Todas las propiedades
del fluido se evalúan a la temperatura de la película.
Esferas
Para el flujo sobre una esfera de diámetro D, se han propuesto numerosas correlaciones. Whitaker, S.
(AIChEJ., 18, 361, 1972) recomienda la siguiente expresión:
 
Nu D  2  0.4 Re 0D.5  0.06 Re 2D/ 3 Pr 0.4  /  S 
0.25
(A.18)
0.71  Pr  380 3.5  ReD  7.6x104 1.0  (  /  S)  3.2

Nuss=2+(0,4*Re^0,5+0,06*Re^(2/3))*Pr^0,4*(mu/mu_s)^0,25
Todas las propiedades menos  s se evalúan a T ; la viscosidad  s se evalúa a la temperatura
superficial, Ts, y los resultados pueden aplicarse a problemas de transferencia de masa simplemente
reemplazando NuD y Pr por ShD y Sc respectivamente.
22
Un caso especial de transferencia convectiva de calor y/o masa desde esferas está relacionado a la
transferencia desde gotas líquidas descendiendo libremente y la correlación de Ranz, W. y W.
Marshall (Chem. Eng. Prog., 48, 141, 1952) se usa frecuentemente:
Nu D  2  0.6 Re 0D.5 Pr 1 / 3 (A.19a
Para esta ecuación las propiedades se calculan a T . En la literatura no aparece rango de validez.
Solamente que se aplica para gotas que caen libremente en aire.
En el límite cuando ReD  0, estas dos ecuaciones tienden a NuD = 2, valor que corresponde a la
transferencia de calor por conducción (en ausencia de convección natural) desde una superficie
esférica.
Para transferencia de masa a baja velocidad, la ecuación de Ranz y Marshall se puede adaptar
intercambiando Nu por Sh y Pr por Sc:
Sh  2  0.6 Re0.5 Sc1/3 2  Re  200 0.6  Sc  2.7 (A.19b)
Frossling (1938) propuso la siguiente correlación:
Sh  2  0.552 Re0.5 Sc1/3 2  Re  800 0.6  Sc  2.7
Otras correlaciones se encuentran en la literatura. Por ejemplo, la tabla 6.3, p 276 del Skelland aporta
otras doce correlaciones.
Lechos Empacados
El flujo de gases (o líquidos) a través de lechos empacados de partículas sólidas tiene importancia en
muchos procesos industriales que incluyen transferencia y almacenamiento de energía térmica,
reacción catalítica heterogénea y secado.
Los lechos compactos de partículas sólidas, como el lecho de bolas que se muestra en la figura, se
usan a menudo como intercambiadores de calor o bien, para almacenar energía térmica. Las partículas
se calientan haciendo pasar un fluido caliente a través del lecho, y la energía almacenada se extrae
posteriormente haciendo pasar un fluido frio. El objetivo del lecho puede ser simplemente transferir
calor de una corriente fluida a otra: en ese caso se dice que el lecho es un intercambiador de calor
23
regenerativo. Otra opción posible es usar el lecho para almacenar energía térmica durante un tiempo
determinado. Por ejemplo, en un sistema de calefacción doméstica por energía solar, el lecho se
calienta durante el día, cuando la energía solar está disponible, y el calor se extrae por la noche para
calentar la casa. Los lechos compactos también se usan como intercambiadores de materia en
procesos de separación, donde se encuentran partículas de muchas formas. La fracción de vacío  de
un lecho compacto se define como
 = (Volumen del lecho - Volumen total de las partículas)/Volumen del lecho
Los valores típicos de , varían entre 0.3 y 0.5. El área de transferencia depende del tamaño y la
forma de las partículas, así como de la fracción de vacío. El área específica a [m-1] de un lecho
compacto es el área de transferencia, por unidad de volumen del lecho:
a = Área total de las partículas/Volumen del lecho
Para un lecho de partículas de volumen, cada una, Vp y área Ap,
6 1   
a
dp
Por ejemplo, para una partícula esférica de diámetro dp, Ap/Vp = 6/dp y a = 6(l-)/dp. En general,
podemos calcular el área específica si conocemos la geometría de las partículas y la fracción de vacío
del lecho. La velocidad superficial V se define como la velocidad del fluido en el lecho si no hubiese
partículas:
V=m'/Ac
donde m' [kg/s] es el flujo másico,  es la densidad del fluido y Ac es el área de la sección transversal
del lecho. La velocidad superficial también puede verse como la velocidad del fluido un instante antes
de entrar al lecho o un instante después de salir de él. Puesto que el área de la sección transversal
disponible para el flujo es solo Ac, la velocidad intersticial es V/.
En la literatura se encuentran muchas correlaciones útiles para estimar el coeficiente convectivo para
flujo en lechos empacados:
Para gases, Gupta y Thodos, 1963 - 1964 (Treybal 1981 p 74; Hines y Maddox p 206 1987)
 jH   j AB  2.06 Re D0.575 Pr = 0.7, Sc = 0.6 90  ReD  4000 (A.20a)
El número de Reynolds se define con base en la velocidad que existiría si el lecho estuviera vacío y la
dimensión característica de las partículas dP = 6Vp/Ap que para una partícula esférica es el diámetro.
Para geometrías no esféricas se puede expresar en términos de su área superficial como
d P  0.567 AP donde AP es el área de la superficie de la partícula. Para empaques diferentes a las
esferas, la correlación anterior usando los factores de corrección siguientes: cilindro de longitud igual
al diámetro multiplicar por 0.79; para cubos multiplicar por 0.71.
Wilson y Geankoplis (1966) recomiendan para flujo de líquidos:
24
 jH   j AB  1.09 Re 2/3
D 0.0016  ReD  55 165  Sc  70600 (A.20b)
 jH   j AB  0.250 Re D0.31 0.55  ReD  1500 165  Sc  10690 (A.20c)
En ambos casos 0.35    0.75
Específicamente para Transferencia de masa, Dwivedi y Upadhyay, 1977 (Tosun 2002 p 105),
proponen una correlación válida para gases y líquidos en lechos empacados y fluidizados:
0.765 0.365
 JM   0.386
Re0.82
D
Re D 0.01 < ReD <15000 (A.20d)
No se recomienda para gases con Re < 10.
FLUJO INTERNO
La transferencia de calor a los fluidos que fluyen dentro de tuberías y conductos es un tema de gran
importancia práctica. En consecuencia, se han propuesto muchas correlaciones para calcular los
coeficientes de transferencia de calor en este tipo de flujo. Sin embargo, las correlaciones presentadas
a continuación son adecuadas para manejar la gran mayoría de las aplicaciones de transferencia de
calor de procesos. Las correlaciones se dan para cada uno de los tres regímenes de flujo; turbulento,
de transición y laminar. La situación física descrita por las correlaciones se muestra en la Figura A.7.
El fluido entra a la tubería a una temperatura promedio Tbl y sale a una temperatura promedio Tb2. (Tb
significa temperatura de mezcla (es la temperatura del fluido promediada sobre el área de flujo). La
temperatura de la pared de la tubería es Twl en la entrada y Tw2 en la salida. El fluido se calienta o se
enfría a medida que fluye a través de la tubería, dependiendo de si la temperatura de la pared es
mayor o menor que la temperatura del fluido.
Figura A.7
Pierce, (Chemical Engineering 86, 27, p 113, 1979), propone una sola ecuación que correlaciona el
factor de Colburn con todos los números de Reynolds. Esta ecuación produce una curva continua y
suave conveniente para utilizar en los análisis y diseños ayudados por computador o calculadora:
25
1
  6 8
3 / 2 1 

12
0.14
 1  Re 1.6
 1.986  10     b 
JH   9.36    
 
     
 S
14
Re   7.831  10  Re     
   
  (A.21))
 h  2/3
J H  St H Pr 2/3    Pr
Es aplicable para 0.7  Pr  16700.   CPV  y todos los rangos de Reynolds.
Para copiar al código en EES:
J=(1/Re^9,36+((Re^1,6/7,831e-14+(1,986e6/Re)^8)^1,5)^(-1))^(1/12)
La extensión a transferencia de masa es obvia haciendo JH = JAB = StAB Sc2/3 = (kC/vm)Sc2/3

y haciendo igual a la unidad la corrección por viscosidad. Otras expresiones que tradicionalmente se
han usado se mencionan a continuación.
Flujo Laminar en Tubos
Para el flujo laminar en tuberías circulares (Re <2100), una correlación de datos experimentales fue
hecha por F. N. Sieder y G. E. Tate (Ind. Eng. Chem., 1429, 1936), y generalmente tiene una
precisión de hasta  25%:
1/ 3 0.14
 D  b 
Nu D  1.86 Pe   
 L  S  (A.22a)
La cantidad Pe  Re Pr se conoce como número de Peclet.
Todas las propiedades excepto  S se evalúan a la temperatura media global. El último paréntesis de
la ecuación (18) tiene en cuenta la distorsión del perfil de velocidades por el gradiente de temperatura
 
entre le pared y el fluido. Válida para 0.48  Pr  16700; TS constante; 0.0044  ( / S )  9.75.
Whitaker recomienda su uso cuando Nu  3.72, pues para valores menores se presentan condiciones
de perfiles completamente desarrollados donde Nu = 3.66, valor teórico para un tubo muy largo con
temperatura de pared constante. Además para números de Reynolds muy pequeños los efectos de la
convección natural deben tenerse en cuenta.
Para transferencia de masa la anterior correlación toma la forma
1/3 1/3
 D  D
Sh  1.86  Re Sc   Re Sc 
 L  L >2 (A.22b)
Hausen, H. (1943) presentó la siguiente fórmula del número medio de Nusselt en la región térmica de
entrada en el caso perfil parabólico de velocidad y pared de temperatura constante:
0.14
 C1 ( D / z ) Re Pr   b 
Nu   Nu   
n 

 1  C 2 ( D / z ) Re Pr    S  (A.23)
26
Según los valores de C1, C2 y n, la ecuación (19) tiene diferentes aplicaciones:
Tabla A-4
Condición en la Perfil de Pr Nu Nu C1 C2 n
superficie velocidad
Ts constante Parabólico Todos Medio 3.66 0.0668 0.04 2/3
Ts constante Plano 0.7 Medio 3.66 0.1040 0.016 0.8
qs constante Parabólico Todos Local 4.36 0.023 0.0012 1.0
qs constante Plano 0.7 Local 4.36 0.036 0.0011 1.0
La extensión a transferencia de masa es sencilla si tenemos presente que JD = JH.
Los cálculos de transferencia de calor en la región de transición están sujetos a un mayor grado de
incertidumbre que aquellos en los regímenes turbulentos laminar y completamente desarrollados.
Aunque el equipo industrial a veces está diseñado para operar en la región de transición,
generalmente se recomienda evitar trabajar en este régimen de flujo si es posible.
Flujo turbulento en tubos circulares lisos
Con respecto a la transferencia de calor en tuberías circulares, se logra un flujo turbulento

completamente desarrollado con un número de Reynolds de aproximadamente 104.
Para tubos lisos y longitudes superiores a sesenta diámetros, por experimentación se encuentra que el
coeficiente convectivo de transferencia de calor h    D, k , CP ,  ,  , V  . A partir del análisis
dimensional encontramos el siguiente agrupamiento de tales variables:
hD  vD C P  
   , 
k   k  (A.24)
Nusselt experimentó con tres gases (aire, gas carbónico y gas natural). Al representar en coordenadas
doble logarítmicas h contra la velocidad másica superficial  V  , encontró que para cada uno de los
gases, los resultados se alineaban sobre una recta por encima de un cierto valor crítico de  V  ,
resultando tres rectas paralelas. Estos resultados se explican si la función  de la ecuación (A.24)
puede expresarse de manera potencial, a saber: Nu  C Re Pr . A partir de la pendiente de las rectas
a b
en el gráfico doble logarítmico encontramos el valor del exponente del número de Reynolds como a
= 0.8. El hecho de que para valores de la velocidad másica superficial inferiores a la crítica fueran
superiores a los correspondientes por extrapolación de las rectas de pendiente 0.8, fue explicado por
Nusselt como debido a la contribución de la convección natural a los flujos gaseosos a baja velocidad.
De experimentaciones posteriores con otros gases y líquidos se obtuvo el valor de 0.023 para C, y se
encontró que b estaría entre 0.3 y 0.4, siendo más apropiado el valor 0.3 cuando el fluido se enfriaba y
0.4 cuando se calentaba.
Para este régimen de flujo (Re > 104), las siguientes correlaciones son ampliamente utilizadas:
Ecuación de Dittus Boelter (1930)
Nu  0.023Re0.8 Pr n 0.7  Pr  100 L / D  60 (A.25)
27
con n = 0.4 para calefacción y n = 0.3 para enfriamiento. Las propiedades del fluido deben calcularse
a la temperatura media aritmética de las másicas extremas del fluido.
Esta ecuación hace una corrección ambigua al efecto de la distorsión del perfil de velocidades por los
gradientes de temperatura radiales, pues mientras los gases aumentan su viscosidad con el aumento de
temperatura, los líquidos disminuyen su viscosidad con la temperatura.
Ecuación de Colburn (1933)
Nu  0.023Re0.8 Pr1/3 (A.26a)
J H  St H Pr 2/3  0.023Re 0.2 (A.26b)
Lo anterior recordando que St H  Nu / Re Pr , confirmando que los coeficientes tipo J son solamente
función de Reynolds.
Las propiedades del fluido excepto CP en el número de Stanton ( St H  h /  CPV ), se evalúan al valor
medio de la temperatura de película Tf. Re >10000; 0.7  Pr  160; L/D > 60.
Ecuación de Sieder y Tate (1936)
Nu  0.027 Re0.8 Pr1/3  m /  S 

0.14
(A.27a)
J H  m / s 
0.14
 0.027 Re 0.2 (A.27b)
Se recomienda esta relación para la transferencia de calor en fluidos cuya viscosidad cambie
marcadamente con la temperatura y es aplicable para 0.7  Pr  16700; L/D > 60 y Re > 10000.
Todas las propiedades excepto  s se determinan a la temperatura media del fluido.
Las siguientes correlaciones son de investigaciones más recientes:
Whitaker (1972) propone la correlación:
Nu  0.015Re0.83 Pr 0.42   /  s 
0.14
(A27c)
Válida para 2300 < Re < 105 0.48 < Pr < 592 0.44 < m /  s < 2.5
Todas las propiedades físicas excepto  s se evalúan a la temperatura promedia. La viscosidad  s se

evalúa a la temperatura de la superficie.
Tuberías rugosas
Una correlación que se atribuye a Petukhov (1970) es de la forma
28
( f / 2) Re D Pr
Nu D 

1.07  12.7 f / 2 Pr 2 / 3  1  (A.28)
0.5  Pr  2000; 10000  Re  5x106. Para obtener concordancia con valores en la zona de transición,
Gnielinski (1976) modificó la expresión así:
( f / 2)(Re D  1000) Pr   D  
2/3
NuD  1    
1  12.7 f / 2  Pr 2/3  1   L   (A.29)
0.5  Pr  2000; 3000  Re  5x106. Las ecuaciones (A.28) y (A.29) se aplican tanto a temperatura
constante como a flujo constante en la pared. Esta última predice los coeficientes en general dentro de
un  20%. Las propiedades se evalúan a la Temperatura global media. Aquí f es el factor de fricción
de Fanning que puede obtenerse del diagrama de Moody o de la ecuación (A.31).
Flujo Turbulento en conductos lisos no circulares
Estas ecuaciones para flujo turbulento, se pueden aplicar para conductos de sección transversal no
circular si el diámetro se reemplaza por el diámetro equivalente Deq :
Deq = 4 x radio hidráulico = 4 x área de flujo del fluido/perímetro
Observe que el área de flujo del fluido es diferente al área de flujo para el calor.
El caso que con más frecuencia se encuentra de flujo laminar en conductos no circulares es el flujo en
el espacio anular de un intercambiador de tubo en tubo. Para este caso se recomienda la siguiente
ecuación de Gnielinski:
0.8
D 
0.8
0.19 1  0.14 D2 / D1   Re Pr Deq / L 
Nu  3.66  1.2  2   0.467
 D1  1  0.117  Re Pr Deq / L  (A.30)
D1 es el diámetro exterior del tubo interno

D2 es el diámetro interior del tubo externo
Deq es el diámetro equivalente del espacio anular: D2 - D1
Tanto el número de Nusselt como el de Reynolds en esta ecuación están basados en el Deq .
Todas las correlacione para flujo interno dan el coeficiente promedio para toda la longitud del tubo.
Factor de fricción
Combinando ecuaciones para tuberías rugosas en regímenes de flujo laminar y turbulento Churchill
(1977) desarrolló una sola ecuación para el factor de fricción en flujo laminar, de transición o
turbulento, y para tuberías tanto lisas como rugosas fanning:
29
1 / 12
f  8   1  
12 3/ 2
     
2  Re   A  B   (A.31)
16 16
  1    37530 
A  2.457 ln   B 
 7 / Re   0.27( / D)    Re 
0.9

Las ecuaciones son más convenientes que las tablas o las correlaciones gráficas en diseño y operación
ayudadas por computador. Prueba y error son necesarios si en lugar de la velocidad de flujo se
especifica la caída de presión.
La ecuación (A.31) no sólo reproduce el gráfico del factor de fricción sino que evita interpolaciones y
da valores únicos en la región de transición. Estos valores están, naturalmente, sujetos a alguna
incertidumbre, debido a la inestabilidad física inherente en esta región.
Una correlación simplificando Colebrook para tuberías lisas es (White) Fanning y Darcy:
= 4 f   0.782ln Re 1.51

2
fD
(A.32)
Otra ecuación que se encuentra en la literatura es la ecuación de Notter y Sleicher (Notter, R.H. and
C.A. Sleicher, "A solution to the turbulent Graetz problem III. Fully developed and entry region heat
transfer rates," Chem. Eng. Sci., Vol. 27, (1972), pp. 2073-2093)
Nu = 5 + 0.016 Rea Prb; a = 0.88  0.24/(4 + Pr); b = 0.33 + 0.5e0.6Pr (A.33)
0.1  Pr  104; 104  Re  106 ; L/D > 25.
Región de entrada térmica
Nu = 0.036 Re0.8 Pr1/3 (D/L)0.055; 10  L/D  400 (A.34)
Las propiedades del fluido se toman a la temperatura media.
Casi se puede decir que existen tantas correlaciones como investigadores.
Flujo Turbulento de Metales Líquidos Dentro de Tubos Lisos
Para flujo de calor constante en la pared Skupinski et al. (1965) recomendaron la siguiente
correlación:
NuD = 4.82 + 0.0185 Pe0.827 (A.35)
3.6x103  Re  9.05x105 ; 100  Pe  104
Para temperatura constante en la pared Seban y Shimazaki (1951) recomiendan:
Nu = 5.0 + 0.025 Pe0.8 ; Pe > 100 (A.36)
30
Para transferencia de masa, flujo al interior de tubos, en Welty et al. se encuentra la correlación
(Gilliland, Sherwood, Linton 1950). Las restricciones según tabla 3.3 de Treybal.
Sh  0.023Re0.83 Sc1/3 2000  Re  70000 0.6  Sc  3000 (A.26c)
Para transferencia de masa Gilliland y Treybal (1934) encontraron experimentando en columnas de

pared húmeda (gases) la siguiente correlación
Sh  0.023Re0.83 Sc 0.44 (A.26d)
Válida para 2000< Re <35000 y 0.6 < Sc < 2.5
FLUJO INTERNO CON CAMBIO DE TEMPERATURA
A diferencia del flujo externo, la fuerza guía, es decir, la diferencia de temperaturas entre la pared y la
corriente de fluido cambia a lo largo del tubo. Por esto se debe establecer un balance de energía
térmica para una longitud diferencial del conducto donde la temperatura del fluido es T.
Las condiciones térmicas en la superficie generalmente se pueden aproximar para tener una temperatura
superficial constante (Ts = constante) o un flujo de calor superficial constante (qs = constante). Por
ejemplo, la condición de temperatura de superficie constante se obtiene cuando se produce un proceso de
cambio de fase, como ebullición o condensación, en la superficie exterior de un tubo. La condición de
flujo de calor superficial constante ocurre cuando el tubo se somete a radiación o calentamiento por
resistencia eléctrica uniformemente desde todas las direcciones.
Figura A.8: Esquema de flujos de calor en un elemento de volumen al interior de un conducto
Temperatura de la superficie es constante: Como la temperatura T del fluido dentro del tubo cambia
longitudinalmente, la fuerza guía, T  TS  también lo hace. Analizando una longitud
dz
de tubo, el
fluido tiene un cambio dT en su temperatura. Igualando el cambio de energía interna del fluido con el
calor intercambiado con el medio tenemos:
m ' CP dT  h TS  T  D  dz (i)
Separando variables en (i):
31
dT h  D 
 dz
TS  T  m ' CP (ii)
Integrando entre T  T1 para x  0 y T  T2 para z  L se obtiene:
T2  TS h  D  L
ln 
T1  TS m ' CP (iii)
T2  TS  h  D  L 
 exp   
T1  TS  m ' CP  (iiia)
m '  V  D 2 / 4  
Reconociendo que , se puede modificar el exponente así:
T2  TS  L h   L 
 exp  4   exp  4 St H 
T1  TS  D  CPV   D  (iiib)
Se reconoce que la última fracción es el número de Stanton. Esta relación también se puede usar para
determinar la temperatura del fluido en cualquier posición z reemplazando L por z. Note que la
diferencia de temperatura entre el fluido y la superficie decae exponencialmente en la dirección del
flujo, y la tasa de decaimiento depende de la magnitud del exponente h  D  L / m ' CP , como se
muestra en la Figura. Este parámetro adimensional se denomina número de unidades de transferencia,
denotado por NUT, y es una medida de la eficacia de los sistemas de transferencia de calor. Para
NUT = 5, la temperatura de salida del fluido se vuelve casi igual a la temperatura de la superficie, T2
= Ts:
NUT 0.01 0.05 0.10 0.50 1.00 5.00 10.00

% 1 1.125 9.5 39.38 63.25 99.38 100
Un NUT de aproximadamente 5 indica que se alcanzó el límite para la transferencia de calor, y la

transferencia de calor no aumentará sin importar cuánto extendamos la longitud del tubo.
32
Escribiendo (iii) como
 T T    DhL
m ' CP ln  1 S    DhL m ' CP 
 T2  TS   T T 
ln  1 S 
 T2  TS  (iv)
El calor transferido a lo largo del tubo se calcula por un balance global alrededor de todo el tubo:
Q  m ' CP (T2  T1 ) (v)
Reconociendo que TS  T1   TS  T2   T2  T1  y reemplazando m'CP de (iv) en (v) se obtiene:
Q  h  DL  T  ML (vi)
El promedio logarítmico de las fuerzas guía (diferencias de temperatura) a la entrada y salida del
intercambiador es:
33
TS  T1   TS  T2 
 T ML 
 T  T  
ln  S 1 
 TS  T2   (vii)
La expresión (vi) permite escribir el flujo de calor en forma semejante a la ley de Newton del
enfriamiento. En el caso presente, la temperatura de la superficie se puede suponer constante y
conocida, pero si el coeficiente externo no fuera tan grande como lo es en la condensación de vapor
de agua, y la resistencia de la pared del tubo no fuera despreciable, en lugar de h debemos usar un
coeficiente tipo U que incluya las diferentes resistencias en serie.
Flujo de calor qS constante en la superficie: El análisis en este caso es matemáticamente simple

porque el calor transferido al fluido es:
Q  qS   D  L
El balance en una longitud diferencial es:
m ' CP dT  qS  D  dz
La variación de la temperatura con la posición vendrá dada por
qS   D 
T  T1  z
m ' CP
La temperatura del fluido cambia linealmente con la distancia. Como el coeficiente convectivo al
interior de conductos de sección transversal constaante es constante, de la ley del enfriamiento se
concluye que la diferencia de temperaturas entre el fluido y la superficie dl conducto también debe ser
constante:
qS  h TS  T   qS
TS  T  
h
34
EJEMPLO A.2: (Tosun pp 91) Vapor que condensa en la superficie exterior de un tubo circular de
pared delgada de 65 mm de diámetro mantiene una temperatura superficial uniforme de 100 °C. Fluye
aceite por el interior del tubo con velocidad media de 1 m/s. Determinar la longitud del tubo necesaria
para aumentar la temperatura del aceite de 40 °C a 60 °C. Las propiedades físicas del aceite a 50 °C
son las siguientes:
μ = 12.4×10−3 kg/m-s
ν = 4.28×10−5 m2/s
Pr = 143
A 100 ◦C: μ = 9.3×10−3 kg/m-s.
Solución: En este caso se conocen las temperaturas de entrada y de salida pero se desconoce la
longitud del conducto. Es un típico problema de dimensionamiento. La expresión útil en esta
situación es la (iii) donde se puede dejar la longitud L de forma explícita. Las propiedades del fluido
deben ser determinadas a la temperatura media aritmética entre entrada y salida, a 50 °C. Es necesario
calcular el coeficiente convectivo h . Para seleccionar entre las ecuaciones para flujo interno en
tuberías lisas se debe conocer si el flujo es laminar o turbulento:
VD
Re   1519

Flujo laminar. Si se usa la ecuación (A.22a), encontramos que esta da el coeficiente convectivo en
función de la longitud que es la incógnita:
1/3 0.14
 D   
Nu D  1.86  Re Pr   b 
 L   S 
Puede hacerse una solución por prueba y error. Sin embargo, recordando que:
h Nu
St H  
 CPV Re Pr
Reemplazando el valor de Nusselt sin la corrección de viscosidad, se escribe la ecuación (iii) = (iv):
L 1.86  Re Pr D / L 
1/3
T T
ln 2 s  4
T1  Ts D Re Pr
Se obtiene L / D en forma explícita:
35
3/ 2
L  1  T  T 
 Re Pr  ln  s 1  
D  4  1.86  Ts  T2   =2602
La longitud del tubo debe ser entonces L = 169 m.
El resultado usando EES:
NOTA: Debe tenerse en cuenta que el software toma Nu como la letra griega  . Se introduce la
subrutina "function" para que el mismo programa calcule casos de flujo laminar o turbulento:
"Coeficiente convectivo para Flujo en tuberías"

$UnitSystem SI mass C Pa Rad J
function Nuh(Re;Pr;D;L;mu;mu_s)
if Re < 2100 Then
Nuh = 1,86*(Re*Pr*D/L)^(1/3)*(mu/mu_s)^0,14 "laminar"
Else
Nuh = 0,015*Re^0,83*Pr^0,42*(mu/mu_s)^0,14 "turbulento"
endif
end
"Definiciones"
Re=rho*V*D/mu
"Pr=mu*C_p/k"
"Datos"
mu=12,4e-3: mu/rho=4,28e-5: Pr=143: mu_s=9,3e-3
D=65e-3: T_o=40: T_s=100: T_L=60: V=1
"Ecuaciones"
Nuh=Nuh(Re;Pr;D;L;mu;mu_s) {se llama a la "function"}
ln((T_L-T_s)/(T_o-T_s))+4*L/D*Nuh/Re/Pr=0
Resultados: L=169,1 [m] Nuh=8,464 Re=1519 (laminar)
El uso de la ecuación (A.21) puede hacer la solución más directa sin necesidad de incluir la
"function".
EJEMPLO A.3: Entra aire a un tubo de pared delgada de 1.5 cm de diámetro interno a 20 °C, con
velocidad de 50 m/s. La temperatura de la superficie se puede considerar constante en 150 °C gracias
a que en el exterior condensa vapor de agua (steam) a esa temperatura. Determine la longitud de
tubería necesaria para aumentar la temperatura del aire hasta 90 °C.
Datos: El aire a 20 °C tiene  = 1.2047 kg/m3. Para aire a 55 °C la viscosidad es  = 19.8x106

kg/m.s, a 150 °C,  s = 23.86 x106 kg/m.s; Pr = 0.707
Solución: se usa el mismo programa teniendo en cuenta que el aire es un fluido compresible pero que
si el caudal másico m '  VA es constante, con A constante, V es constante.
36
Resultado: L=1,005 [m] Nuh=93,05 Re=45633 (turbulento)
NOTA: En el caso en el cual las resistencias externa y de la pared del tubo no son despreciables, pero
la temperatura del fluido externo se puede considerar constante T , el balance es similar, pero en la
ecuación (ii) se reemplaza h, el coeficiente convectivo interno por U, un coeficiente global que tiene
en cuenta las tres resistencias en serie:
U  Per 
T
 dT
2
L
 
T1 T  T  V   CP  Ax
0
dx
T1  T UA
ln 
T2  T m ' CP
Similarmente la ecuación (vii) se escribe
Q  UA  T  ML (viia)
con
T  T1   T  T2 
 T ML 
 T  T  
ln   1 
 T  T2   (viiia)
Recordando la definición del coeficiente global:
1

1

ln R2 R1

1  Resistencias
UA 2 R1 Lh1 2 kT L 2 R2 Lh 2 =
La ecuación (viia) se escribe
 T ML
Q
 Resistencias (viib)
Ocasionalmente puede convenir seleccionar el coeficiente global referido al área interna del tubo. Si
así se hace, para un conducto circular, Per   D1 ; Ax   D1 / 4 ; Per / Ax  4 / D1 . Integrando se
2
obtiene:
T2  T U  Di  L L Ui
ln  i  4
T1  T m ' CP Di  CPV
1 1 R1 ln R2 R1 R1
  
U i h1 kT R2 h 2
37
Ejemplo A-4: (20.38 y 20.39 p 321 Welty et al. modificados) La válvula de un tubo de agua caliente,
dispuesto horizontalmente, se abre justamente lo suficiente para permitir un flujo de 0.02 m/s. El agua
entra a 80 °C y el tubo, de 1/2 plg nominal, cédula 40 (Schedule), que conduce el agua, está en aire
estancado a 20 °C. Las características de esta tubería son: diámetro interior 0,622 plg; diámetro
exterior 0.840 plg (dimensiones estándar de tuberías, Apéndice M, Welty et al.). La conductividad
térmica del material a la temperatura de operación es 50 W/m.K.
a) ¿Cuál es la pérdida total de calor en 5 m de la tubería de agua bajo tales condiciones? ¿Cuál es la
temperatura de salida del agua?
b) Cuando la válvula en la línea de agua mencionada, se abre totalmente, la velocidad del agua es de 3
m/s. ¿Cuál es la pérdida de calor en 5 m de tubería en este caso, si las temperaturas inicial del agua y
la del aire son las mismas que se especificaron anteriormente? ¿Cuál la nueva temperatura de salida?
Solución: Como la temperatura T del líquido dentro del tubo cambia longitudinalmente, la fuerza
guía, T  T  también lo hace. Analizando una longitud
dx dT
de tubo, el fluido tiene un cambio en
su temperatura. Igualando el cambio de energía interna del fluido con el calor intercambiado con el
medio tenemos: m ' CP dT  qs Pdx . Teniendo en cuenta que qs incluye varias resistencias en serie:
m ' CP dT  U o  Do T  T  dx (i)
1 1 ln  Do / Di  1
  
U o  Do L  hi  Di L  2 kT L ho  Do L  (ii)
Separando variables en (i):
dT U  Do 
 o dx
T  T  m ' CP
Integrando entre T  T1 para x  0 y T  T2 para x  L se obtiene:
38
T2  T U  Do  L
ln  o
T1  T m ' CP (iii)
T2  T  U  Do  L 
 exp   o 
T1  T  m ' CP  (iiia)
Solución aproximada: A partir de la ecuación (ii) se evidencia que el mayor esfuerzo está en el
cálculo de los coeficientes convectivos interno y externo.
a) Coeficiente convectivo interno: Para las propiedades físicas es necesario usar la temperatura media
aritmética del fluido. Como se desconoce la temperatura de salida habrá que suponer algún valor. La
temperatura de salida T2 estará entre T1 = 80 °C y T = 20 °C. A 320 K en la Tabla A.6 de Incropera
 
et al, 4a ed., = 989 kg/m3, CP = 4180 J/kg.K, = 577x10-6 kg/m.s, k = 0.640 W/m.K, Pr = 3.77.
Se calcula el criterio Reynolds para conocer el régimen de flujo:
VDi 989  0.02  0.622  25.4  103 

Re 
 = 577  106 = 542 < 2100
El régimen es laminar. Se usa la ecuación (A.22a)
1/3
D   
0.14
 0.016 
1/3
hi Di
 1.86  542  3.77 
Nu D  1.86  Re Pr   b   ki
 L   S   5 
= = 3.13 =
En este cálculo preliminar omitimos la corrección por viscosidades. Se obtiene hi = 125 W/m2.K.
Coeficiente convectivo externo. La transferencia es con aire estancado, por lo cual es transferencia de
calor por convección natural. Se usa una correlación simplificada de la tabla A.7 asumiendo Rayleigh
laminar:
ho  1.32 Tso  T  / Do  = 1.32  50  20  / 0.021 = 8.1 W/m2.K

1/ 4 1/ 4
Se asumió una temperatura promedia para la superficie externa del tubo la misma temperatura
promedio con la que se leyeron las propiedades del fluido interno. Retomando la ecuación (ii)
tenemos:
1 1 ln  0.021 / 0.016  1  UA
  
UA 125  0.016  5  2 50  5 8.1 0.021  5  = 0.406 K/W = 2.46 W/K
Retomamos la ecuación (iiia)
 UA      T
T2  T 2.46 2
 exp    exp   
T1  T  VAi CP   989  0.02   0.016 / 4  4180 
2
= = 0.86 = 71.8 °C
39
Qlaminar  974  0.02    0.0162 / 4  4195  80  71.8   131.4 W
Se observa que el producto UA es constante e igual al inverso de la suma de resistencias. Observe

que Ai se refiere al área para flujo de fluidos y A es el área para transferencia de calor y cambia si se
refiere al interior o el exterior del tubo. Con este resultado, se deduce que se deben hacer nuevos
cálculos con las propiedades leídas a 350 K y la temperatura promedia de la superficie externa del
tubo estimada por un balance térmico.
b) Se hacen los mismos pasos anteriores.

Coeficiente interno. Partimos de propiedades para el agua a 350 K según tabla A.6:  = 974 kg/m3,
CP = 4195 J/kg.K,  = 365x10-6 kg/m.s, k = 0.668 W/m.K, Pr = 2.29. Se observa que la propiedad
más sensible al cambio de temperatura es la viscosidad afectando principalmente Re y Pr. Se calcula
el criterio Reynolds para conocer el régimen de flujo:
VDi 974  3  0.016 

Re 
 = 365  106 = 1.28x105 > 2100
Régimen turbulento. Usamos Sieder y Tate, (A.27a), sin corrección de viscosidad:
Nu  0.027 1.28  105  2.291/3

0.8
hi Di  hi
= 433.6 = ki = 18100 W/m2.K.
Coeficiente externo. Se usa la misma correlación que en el caso anterior:
ho 1.32  75  20  / 0.0211/ 4
=   = 9.4 W/m2.K
Tomando el coeficiente convectivo global referido al área interna del conducto:
1 1 R ln R2 R1 R 1 1 0.008ln 0.0210.016 0.016  U

  1  1    i
U i hi kT U
R2 ho = i 18100 50 0.021  9.4 =0.081 = 12.32 W/m2.K
T2  T  L Ui  T2  20  5 12.32   T2
 exp  4   exp  4 
T1  T  Di  CPV  = 80  20  0.016  974  4195  3  = 79.9 °C
(T2-20)/60=EXP(-20/0.016*12,32/974/3/4195)
Qturbulento  974  3    0.0162 / 4  4195  80  79.9   141.6 W
Las propiedades del agua se calculan a la temperatura promedio aritmética, Tma  T1  T2  / 2 , (no la
temperatura calórica). El coeficiente convectivo externo, por convección natural, se estima a la
temperatura exterior promedia del tubo, Tso .Estos valores deben estimarse lo que implica prueba y
error. La secuencia sugerida es la siguiente: se supone la temperatura de salida del agua, T2 ; se
calcula la temperatura promedia del fluido, Tma  T1  T2  / 2 , se calcula la temperatura promedio
para la superficie externa del tubo Tso haciendo el siguiente estimativo:
40
Tma  Tso  1.32 Tso  T 

1.25

 1 R1 ln R2
R1  Do0.25
  
 h1 kT 
Debido a las características físicas planteadas las primeras suposiciones pueden ser T2
aproximadamente igual a T1 , TS aproximadamente igual a Tm y se itera hasta convergencia.
El coeficiente convectivo interno se calcula para tubería lisa. La expresión de Pierce (A-21):
1
  8 3/ 2  
1 12
 1.986  106    
0.14
 1  Re1.6  b 
J H   9.36          

Re  

7.831  1014  Re      S 
    (A.21))
 h  2/3
J H  St H Pr 2/3    Pr
Es aplicable para 0.7  Pr  16700.   CPV  y todos los rangos de Re.
J=(1/Re^9,36+((Re^1,6/7,831e-14+(1,986e6/Re)^8)^1,5)^(-1))^(1/12)*(mu_b/mu_s)^0,14
El código propuesto en EES es:
$UnitSystem SI Mass C J Pa Deg

"Datos"
V=0,02: T1=80: Tf=20: L=5: Q=m*Cp*(T1-T2)
ID=0,622*0,0254: OD=0,84*0,0254
k_T=50:
"Ecuaciones"
(T2-Tf)/(T1-Tf)=exp(-4*L/ID*U/rho/Cp/V)
1/U=1/hi+ID/k_T/2*ln(OD/ID)+ID/OD/ho
ho=1,32/OD^0,25*(abs(Ts-Tf))^0,25
Tm=(T1+T2)/2
(Tm-Ts)/Res=1,32/OD^0,25*(abs(Ts-Tf))^1,2
Res=1/hi+ID/2/k_T*ln(OD/ID)
J=(1/Re^9,36+((Re^1,6/7,831e-14+(1,986e6/Re)^8)^1,5)^(-1))^(1/12)*phi
phi=(mu/mu_s)^0,14
J=St*Pr^2/3
St=hi/rho/Cp/V
Re=rho*V*ID/mu
mu=Viscosity(Water;T=T;P=P)
mu_s=Viscosity(Water;T=Ts;P=P)
rho=Density(Water;T=T;P=P)
Cp=Cp(Water;T=T;P=P)
Pr=Prandtl(Water;T=T;P=P)
T=75: Ts=70: P=101325
m=V*A*rho: A=pi*ID^2/4
Resultados para flujo laminar interno

A=0,000196 [m^2] Cp=4190 [J/kg-K]
hi=222,3 [W/m^2-K] "coeficiente convectivo interno"
ho=9,323 [W/m^2-K] "coeficiente convectivo externo"
41
ID=0,0158 [m]
J=0,005359 [-]
k_T=50 [W/m-K] L=5 [m]
m=0,00381 [kg/s] "caudal másico, calculado a la entrada"
mu=0,0003782 [kg/m-s] "viscosidad a Tm "
mu_s=0,0003876 [kg/m-s] "viscosidad a TS "
OD=0,02134 [m] P=101325 [Pa]
phi=0,9966 [-] "factor de corrección por distorsión del perfil de velocidad"
Pr=2,427 [-] "número de Prandtl"
Q=161,9 [W] "calor disipado por convección"
Re=811,9 [-] "Re < 2100, flujo laminar"
Res=0,004546 [K/W] "resistencia interna y pared tubo"
rho=971,8 [kg/m^3] "densidad a la temperatura de entrada"
St=0,00273 [-] "número de Stanton"
T=74,93 [C] "para cálculo de propiedades T = Tm "
T1=80 [C] "temperatura de entrada"
T2=69,86 [C] "temperatura de salida"
Tf=20 [C] "temperatura del aire"
Tm=74,93 [C] "temperatura media aritmética"
Ts=73,09 [C] "temperatura promedia de la superficie exterior"
T_s=73,09 [C] "temperatura promedia de la superficie exterior"
U=11,91 [W/m^2-K] "coeficiente global"
V=0,02 [m/s] "velocidad de entrada del agua"
Con V=3 m/s

hi=15580 [W/m^2-K] "coeficiente convectivo interno"
ho=9,609 [W/m^2-K] "coeficiente convectivo externo"
J=0,002183 [-]
m=0,5715 [kg/s] "caudal másico, calculado a la entrada"
mu=0,0003547 [kg/m-s] "viscosidad a Tm "
mu_s=0,0003549 [kg/m-s] "viscosidad a TS "
phi=0,9999 [-] "factor de corrección por distorsión del perfil de velocidad"
Pr=2,267 [-] "número de Prandtl"
Q=192,8 [W] "calor disipado por convección"
Re=129871 [-] "Re muy alto para flujo interno"
Res=0,0001117 [K/W] "resistencia interna y pared tubo"
rho=971,8 [kg/m^3] "densidad a la temperatura de entrada"
St=0,001274 [-] "número de Stanton"
T=79,96 [C] T1=80 [C] T2=79,92 [C] Tf=20 [C] Tm=79,96 [C] Ts=79,91 [C] T_s=79,91 [C]
U=12,96 [W/m^2-K]
V=3 [m/s] "máxima velocidad recomendada para flujo en tuberías"
Es de notar que aunque, con flujo turbulento el cambio de temperatura es mucho menor que con flujo
laminar, el calor transferido es mayor. La razón principal es el aumento del caudal, pues aunque el
coeficiente convectivo interno aumenta notablemente, el efecto sobre el coeficiente global no es
mayor pues la resistencia dominante está en el exterior.
42
EJEMPLO A.5: Considere el flujo de un gas de densidad  y caudal másico m’ (es decir, con caudal
volumétrico Q’= m’/) por el interior de un tubo cuya superficie interna está cubierta con un líquido o
moldeada en un sólido sublimable que garantiza una concentración uniforme de interfase As. Derive
la siguiente expresión para la variación de concentración de la corriente gaseosa A a lo largo del
conducto, donde A1 es la concentración del soluto a la entrada y la concentración a la salida:
 A   As  P x  Per
 exp   er k  
 A1   As  m'  siendo el perímetro.
Demuestre que la transferencia total de vapor para un tubo de longitud L se puede escribir como
m A'  k  PL  AS   A ML   A1   As     A2   As 

  As   A ML 
     As  
ln  A1 
donde    A 2   As  
Nota: Observe que por un balance global también puede decir que mA  Q '   A 2   A1 
'
SOLUCIÓN
Figura Ejemplo A.5
La fuerza guía a la entrada, As A1) = (A)1, difiere de la fuerza guía a la salida (AsA2) = (A)2,
y varía constantemente a lo largo del tubo siendo (As A) = () en cualquier punto x. Por tanto
debemos hacer un balance sobre una longitud diferencial, notando que el incremento en contenido de
soluto por parte de la corriente fluida se debe a la densidad de flujo difusiva nA en la interfase:
Q ' d  A  k    As   A  Per.dx
Para baja velocidad de transferencia, el caudal volumétrico puede considerarse constante, y con
perímetro, coeficiente y concentración superficial constantes, luego podemos separar variables e
integrar así:
Q' 
 A2
d A L
 ln
 AS   A2   k  PL
k P 
  dx  AS   A1  Q'
  AS   A  0
 A1  , es decir
43
 AS   A2   exp  k PL    A 2
 AS   A1   Q'   A 1
 (i)
Esta es la ecuación pedida, evaluada para x = L.
Si deseamos encontrar una fuerza guía promedio para toda la longitud, partimos de un balance global
así:
mA'  Q'  A 2   A1   k  PL A M
Notando que  AS   A1    AS   A 2    A 2   A1    A 1   A 2 y con ayuda de (i):
Q'  A 1   A 2   A 1   A 2 

 A M  
 A 1
  A ML
k ln
 A 2
Por definición el promedio logarítmico de dos cantidades cualesquiera a y b es:
ab ML  a ab  b ba

ln ln
b a , siempre positivo.
Por lo tanto, la fuerza guía promedio en este caso, es el promedio logarítmico de las fuerzas guías en
los dos extremos del conducto. El caso es análogo a la transferencia de calor desde un tubo con
temperatura superficial constante hacia un fluido con velocidad completamente desarrollada. Estas
dos situaciones corresponden a perfiles de velocidad completamente desarrollados y perfil de
concentraciones o de temperaturas en desarrollo. Se debe tener presente que en cualquier sección
transversal del tubo, existen gradientes radiales, por lo cual se debe entender que  A o T
corresponden a valores promediados con respecto a r, los cuales cambian solamente con x.
EJEMPLO A.6: (Problema 28_16 del WWW 5a ed) Un sistema de aireación de membrana de tubo
concéntrico sin burbujeo se usa para transferir oxígeno gaseoso a agua líquida. Agua libre de oxígeno
disuelto se añade al lado del tubo, en la entrada. Gas oxígeno puro, mantenido a una presión constante
de 1,0 atm, fluye a través del espacio anular. El diámetro interior del tubo es de 1,0 cm. La pared del
tubo está hecha de silicona, un polímero que es altamente permeable al gas O2, pero no al vapor de
agua de manera que la concentración de O2 disuelto en el agua en el interior de la pared se establece
por la presión parcial de O2 en el espacio anular. La longitud del tubo es de 500 cm, y la velocidad
media de líquido a través del tubo es de 50 cm/s. El proceso se mantiene a 258 °C. La constante H de
la ley de Henry en el agua a estas condiciones es 0.78 atm/(molg/m3) y el coeficiente de difusión para
el oxígeno disuelto en agua es de 2.1x10-5 cm2/s, y la viscosidad cinemática del agua es 9.12x10-3
cm2/s. Estime la concentración de oxígeno a la salida del tubo y cuál debería ser la longitud del
mismo si se quisiera que el agua saliera 60% saturada.
Solución: El análisis es similar al del ejemplo A.5. La concentración de oxígeno en la superficie

interior del tubo interior hace que la película de agua en contacto con esa superficie se mantenga
44
saturada a la presión y temperatura existentes de acuerdo con la ley de Henry: c As  PO2 / H . La

concentración de mezcla en la fase líquida pasa de c A1 a la entrada hasta c A2 a la salida:
c A2  c As  k P L
 exp   c er 
c A1  c As  VS  (i)
(c_As-c_A1)/(c_As-c_A2)=exp(k_c*pi*D*L/Q) "Ecuación (i)"

c_A1=0: P_o=1: D=1e-2: L=500e-2: V=50e-2: P_o=H*c_As: H=0,78
D_AB=2,1e-9: nu=9,12e-7
"c_A2=0,6*C_As"
Sh=0,023*Re^0,83*Sc^(1/3) "Flujo turbulento en tubos ecuación (A.26c)"
Sh=k_c*D/D_AB "Sherwood, coeficiente adimensional transferencia de masa"
Sc=nu/D_AB "Número de Schmidt"
Re=V*D/nu "Número de Reynolds"
Q=V*A "Caudal volumétrico"
A=pi*D^2/4 "Área para flujo de fluidos"
Resultados
c_A1=0 molg/m3 Concentración de entrada
c_A2=0,1712 molg/m3 Concentración después de L = 5 m (13.35 % de saturación)
c_As=1,282 molg/m3 Concentración de saturación
D=0,01 m D_AB=2,100E-09 m /s H=0,78 atm/(molg/m3) V=0,5 m/s; A=0,00007854 m2
2
L=5 m; nu=9,120E-07 m2/s; P_o=1 atm; Q = 0,00003927 m3/s
Re=5482 Reynolds
Sc=434,3 Schmidt
Sh=221 [-] Sherwood
k_c=0,0000464 m/s Coeficiente convectivo
c_As=1,282 molg/m3 Concentración de saturación
c_A2=0,7692 molg/m3 Concentración 60% de saturación
L=24,68 m Longitud necesaria para 60% de saturación
EJEMPLO A.7: El aire pulmonar se inhala a través de la tráquea a razón de 10 L/min. La tráquea
promedio tiene diámetro interno de 20 mm y longitud 125 mm. Asumiendo superficie saturada con
agua y temperatura del sistema de 37 °C, calcule el coeficiente de transferencia de masa y el agua que
se evapora diariamente por el paso del aire en la tráquea. Como el aire entra y sale, tome el doble de
45
la longitud de la tráquea como longitud efectiva, suponiendo que no hay humidificación adicional en
el pulmón.
Figura ejemplo A.7
Datos: Para aire a 310 K  = 1.128 kg/m3;  = 1.893x10-5 Pa.s

DAB = 0.26x10-4 m2/s a 298 K; Para el agua L = 993 kg/m3; presión de vapor 6221 Pa.
Solución: Como el aire cambia de composición constantemente a lo largo del conducto mientras la
concentración en la superficie permanece constane se hace un balance en una longitud diferencial en
la cual el aire sufre un cambio en el contenido de soluto debido a la transferencia desde la superficie
por convección:
V ' dc A  kc  c As  c A   Ddx (i)
En esta ecuación V ' es caudal volumétrico [m3/s]; D el diámetro interno del conducto [m]; c As es la
concentración de saturación, en la superficie interna del conducto, constante para toda la longitud.
Separando variables e integrando con c A1 a la entrada y c A 2 a la salida se obtiene
 cAS  cA2   exp   kc DL 

 V ' 
 cAS  cA1   (ii)
mA  V '  c A 2  c A1   kc DL  c A  ML (iii)
Con c A1 = 0 suponiendo el aire seco a la entrada y usando la ley de gases ideales:
PA2  k  DL 
 1  exp   c
P T
A  V '  (iv)
Esta expresión nos da la humedad relativa del aire al terminar su recorrido. Para encontrar el
coeficiente convectivo de transferencia de masa comenzamos por conocer el régimen fluidodinámico
calculando el número de Reynolds:
VD 4  V '
Re  
  D
V ' = 10/(1000x60) = 1.667x10-4 m3/s. El metro cúbico es una unidad de volumen grande por lo que es
conveniente escribir al menos tres cifras decimales para evitar la pérdida de litros.
46
Con las propiedades del aire que es el fluido que avanza por el tubo Re = 632 < 2100, significa flujo
laminar en tuberías. La longitud de entrada es L = 0.05xRexD = 0.632 m por lo que debemos
considerar perfil en desarrollo. Empleamos la ecuación (17)
1/ 3 0.14
 D  b 
Nu D  1.86 Pe   
 L  S 
Adaptada a transferencia de masa con el último paréntesis igual a 1 por ser sistema isotérmico:
1/3
 D
ShD  1.86  Pe 
 L (17a)
Para transferencia de masa el número de Peclet es Pe = Re.Sc; la difusividad a 310 K:
1.75
 310 
4
DAB  0.26 10  
 298 
Por lo tanto
 1.893 105
Sc    0.60
 DAB 1.128  2.786 105
Reemplazando en (17a) con L = 0.250 m, Sh = 5.8 > 3.72; 5.8 = kcD/DAB

kc = 8.09  10-3 m/s; calculando para valores promedio en el conducto, B = 0.03; BF = 0.985 valor que
hace innecesaria la corrección.
Reemplazando en (v) con L = 0.250 m, Sh = 5.8 > 3.72; 5.8 = kcD/DAB

kc = 8.09  10-3 m/s; calculando para valores promedio en el conducto, B = 0.03; BF = 0.985 valor que
hace innecesaria la corrección.
Reemplazando en (iv) la presión parcial de vapor de agua en el aire de salida son 3318 Pa, es decir cA2
= 1.2875x10-3 kgmol/m3. Observando la ecuación (iii) no es necesario usar el promedio logarítmico de
la fuerza guía pues es más sencillo el balance global: mA = 1.667x10-4x1.2875x10-3 = 2.146x10-7
kgmol/s. En un día serán:
'
mA = 1.854x10-3 kgmol/s o mA = 333.8 g/s, aproximadamente una tercera parte de un litro.
EJEMPLO A.8: Se transfiere acetato de etilo desde una corriente líquida hacia una corriente de aire
en una columna de pared húmeda de 5 cm de diámetro y 10 m de longitud. El sistema es isotérmico a
27 °C y 1 atm, a estas condiciones la presión de vapor del acetato de etilo es 0.080 atm. El aire entra a
la columna a velocidad promedia de 0.2 m/s.
a) Encuentre la concentración de vapor de acetato de etilo en la corriente de salida y la cantidad
evaporada en la columna en kgmol/h si el aire entrante está libre de estos vapores.
b) Si en lugar de esta columna de pared húmeda, se usa una columna empacada, ¿Qué altura de
empaque será necesaria para obtener el mismo resultado? El diámetro de las partículas es 0.571 cm y
la fracción de vacío en la columna empacada es 0.75.
47
La viscosidad cinemática de la fase gaseosa se puede tomar como  = 1.569x10-5 m2/s, la difusividad
del vapor de acetato de etilo en aire, a 1 atm y 273 K es 0.0709x10-4 m2/s.
Solución: Considerando bajas concentraciones el caudal volumétrico se mantiene aproximadamente

constante. Un balance para una longitud diferencial de la columna de pared húmeda es:
Q ' dc A  kc  c As  c A  Pdz P D

perímetro interno de la columna
Separando variables e integrando entre entrada y salida
 cAS  cA2   exp   kc PL   exp  4 L kc 

   DV 
 cAS  cA1   Q'   
Para calcular el coeficiente convectivo o pelicular, determinamos Reynolds
VD
Re 
 = 637
Es decir, laminar. Se usa la ecuación 17 p 407 RBG:
1/3
 D
Sh  1.86  Re Sc 
 L
DAB = 8.36x10-6 m2/s; Sc = 1.877; Sh = 3.376; kC = 5.64x10-4 m/s; cA2 = 2.91x10-3 kgmol/m3
m' = 4.1x10-3 kgmol/hr
Para una columna empacada el balance es similar pero se debe tener en cuenta el área específica del
empaque:
 D2  D2
V dc A  kc  c As  c A  a dz
4 4
V  c  cA1   L
ln AS
kc a  c AS  c A 2 
El coeficiente debe calcularse con alguna correlación adecuada para flujo en empaques:
jH  j AB  2.06 Re D0.575 Pr = 0.7 Sc = 0.6 90  ReD  4000 (16)
Específicamente para transferencia de masa, Dwivedi y Upadhyay (Tosun 2002 p 105) proponen una
correlación válida para gases y líquidos en lechos empacados y fluidizados:
0.765 0.365
 JM   0.386
Re0.82
D
Re D
Validez para 0.01 < ReD <15000
48
6 1   
a
dP =262.7 m-1
kc 2/3
JM  Sc
Re = 0.2x0.00571/1.569x10 = 72.79; JM = 0.0925;
-5
V ; kc = 0.0162 m/s
L = 0.106 m
El desarrollo en EES es:
$UnitSystem SI Mass J Pa K rad

"Datos"
D=5e-2 [m]
T=300 [K]
L=10 [m]
Pv=0,08*convert(atm;Pa)
V=0,2 [m/s]
Dp=0,571e-2 [m]: epsilon=0,75
nu=1,569e-5 [m^2/s]: D_273=7,09e-6 [m^2/s]
D_300=D_273*(300/273)^1,75
(Cs-C2)/(Cs-C1)=exp(-4*L/D*kc/V)
Cs=Pv/R/T: C1=0
R=R#: Sh= kc*D/D_300
Sh=1,86*(Re*Sc*D/L)^(1/3)
Re=V*D/nu: Sc=nu/D_300
alpha=kf/rho/Cp
kf=conductivity(Air_ha;T=T;P=P)
rho=density(Air_ha;T=T;P=P)
Cp=cp(Air_ha;T=T;P=P)
P=101325
m_A=Q*C2: Q=V*A: A=pi*D^2/4
Le=V/kce/ae*ln((Cs-C1)/(Cs-C2))
epsilon*Jm=0,765/Rep^0,82+0,365/Rep^0,386
ae=6*(1-epsilon)/Dp
Rep=V*Dp/nu
jm=kce/V*Sc^(2/3)
m_Ah=m_A*3600 [kgmol/hr]
Pr=nu/alpha
Resultados
A=0,001963 [m^2] área de flujo
ae=262,7 [1/m] área específica
alpha=0,00002227 [m^2/s] difusividad térmica
C1=0 [kgmol/m^3] concentraión de entrada
C2=0,00291 [kgmol/m^3] concentración de salida
Cp=1007 [J/kg-K] calor específico
Cs=0,00325 [kgmol/m^3] concentración saturación interface
D=0,05 [m] diámetro conducto
Dp=0,00571 [m] diámetro partículas
D_273=0,00000709 [m^2/s] difusividad a 273 K
D_300=0,000008362 [m^2/s] difusividad a 300K
epsilon=0,75 [-] fracción vacío empaque
49
Jm=0,1233 [-] coeficiente adimensional lecho empacado

kc=0,0005646 [m/s] coeficiente conducto circular
kce=0,01621 [m/s] coeficiente lecho empacado
kf=0,02638 [W/m-K] conductividad térmica aire
L=10 [m] longitud columna pared húmeda
Le=0,1061 [m] longitud lecho empacado
m_A=0,000001143 [kgmol/s] evaporación
m_Ah=0,004114 [kgmol/hr] evaporación horaria
nu=0,00001569 [m^2/s] viscosidad cinemática
P=101325 [Pa] presión total
Pr=0,7045 número de Prandtl
Pv=8106 [Pa] presión de vapor del acetato
Q=0,0003927 [m^3/s] caudal olumétrico
R=8314 [J/kmol-K] constante universal de los gases
Re=637,3 [-] Reynolds en tubo
Rep=72,79 [-] Reynolds lecho empacado
rho=1,177 [kg/m^3] densidad aire
Sc=1,876 [-] número de Schmidt
Sh=3,376 [-] coeficiente adimensional tubo
T=300 [K] V=0,2 [m/s]
Transferencia en lecho empacado
8 - (Ejemplo 6.5 p 206 Hines y Maddox modificado) Aire a 75 °C fluye por un lecho de 0.05 m de
largo de partículas cilíndricas de naftaleno. Tanto la longitud como el diámetro de las partículas vale
0.01 m. Para empaques diferentes a esferas se pueden usar la correlación (A2-20a) con el factor de
corrección. La velocidad intersticial de 0.60 m/s y la porosidad (fracción de vacío)  = 0.5. Determine
la concentración promedio del naftaleno en el aire de salida y el porcentaje de saturación del mismo.
El aire de entrada está libre de naftaleno. La presión de vapor del naftaleno a 75 °C es 5.61 mm de
Hg. Considere el sistema isotérmico.
Figura Ejemplo A.5
La composición de la corriente gaseosa cambia constantemente porque es flujo interno. Si la

transferencia es solamente con la pared del tubo, el balance es
Q ' d  A  k    As   A  Pdx
Para baja velocidad de transferencia, el caudal volumétrico puede considerarse constante, y con
perímetro, coeficiente y concentración superficial constantes, luego podemos separar variables e
integrar así:
50
Q' 
 A2
d A L
 ln
 AS   A2   k  PL
   dx
k  P   AS   A  0
 AS   A1  Q'
 A1  , es decir
 AS   A2   exp  k PL    A 2
 AS   A1   Q'   A 1
 (i)
así:
mA'  Q'  A 2   A1   k  PL A M
así:
mA'  Q'  A 2   A1   k  PL A M
Notando que  AS   A1    AS   A 2    A 2   A1    A 1   A 2 y con ayuda de (i):
Q'  A 1   A 2   A 1   A 2 

 A M  
 A 1
  A ML
k ln
 A 2
Por definición el promedio logarítmico de dos cantidades cualesquiera a y b es:
ab ML  a ab  b ba

ln ln
b a , siempre positivo.
Por lo tanto, la fuerza guía promedio en este caso, es el promedio logarítmico de las fuerzas guías en
los dos extremos del conducto. El caso es análogo a la transferencia de calor desde un tubo con
temperatura superficial constante hacia un fluido con velocidad completamente desarrollada. Estas
dos situaciones corresponden a perfiles de velocidad completamente desarrollados y perfil de
concentraciones o de temperaturas en desarrollo. Se debe tener presente que en cualquier sección
transversal del tubo, existen gradientes radiales, por lo cual se debe entender que  A o T
corresponden a valores promediados con respecto a r, los cuales cambian solamente con x.
LECHO EMPACADO
FRACCIÓN DE VACÍO
VT volumen del tanque
Vtodasp Volumen de todas las partículas
(VT - Vtodasp) Volumen vacío en el tanque
(VT - Vtodasp)/VT =  fracción de vacío en el tanque
51
(1 - )VT = [1 - (VT - Vtodasp)/VT]VT = Vtodasp
  VT  Vtodasp  
1   VT  1     VT  VT  VT  Vtodasp
  VT 
Entonces, el número de partículas que hay en el tanque, si VP   d P / 6 es el volumen de una

3
partícula de diámetro d P , será Num = Vtodasp/VP. La superficie disponible para transferencia de masa,
la superficie de todas las partículas en el tanque es entonces Num   d P . En cualquier otro tanque que
2
tenga la misma disposición de partículas podemos encontrar esta superficie si usamos el área
específica del empaque a  Num   d P / VT
2
Para una columna empacada el balance es similar al de flujo en tubos, pero se debe tener en cuenta el
área específica del empaque que como se acaba de estudiar es
a = Área total de las partículas/Volumen del lecho = Vtodasp/VP*Ap/VT = (1 - )VT/VP*Ap/VT
Así, para un lecho de partículas de volumen, cada una, Vp y área Ap, a =(Ap/Vp)(1-)
La fracción de vacío  de un lecho compacto se define como

 = (Volumen del lecho - Volumen total de las partículas)/Volumen del lecho
Los valores típicos de , varían entre 0.3 y 0.5. El área de transferencia depende del tamaño y la
forma de las partículas, así como de la fracción de vacío. El área específica a [m-1] de un lecho
compacto es el área de transferencia, por unidad de volumen del lecho:
Por ejemplo, para una partícula esférica de diámetro dp, Ap/Vp = 6/dp y a = 6(l-)/dp. En general,
podemos calcular el área específica si conocemos la geometría de las partículas y la fracción de vacío
del lecho.
BALANCE DIFERENCIAL LECHO EMPACADO
Cambio de composición diferencial de la corriente = transferencia interfacial
 D2  D2
V dc A  kc  c As  c A  a dz
4 4
V  c  cA1   L
ln AS
kc a  c AS  c A 2 
 cAS  cA2   exp   kc aL 

 V 
 cAS  cA1    (i)
V=V' = 0.5*0.6=0.3 m/s velocidad superficial
 jH   j AB  2.06 Re D0.575 Pr = 0.7, Sc = 0.6 90  ReD  4000 (A.20a)
52
Para calcular la longitud característica se determina el área de la partícula: AP  2 D / 4   DD ; por

2
lo tanto d P  0.567 AP = 1.5D = 0.0123 m. El área específica es entonces:
6 1    Vd P
a Re 
dP = 244 m2/m3 
A estas condiciones  = 2.072E-5 m2/s y Re = 178 Sc = /DAB. De la tabla 14-2 p 781 Cengel, la
difusividad de vapor de naftaleno en aire a 300 K y 1 ata es 0.65E-5 m2/s. A Tf = 348 K es:
1.75
 348 
DAB  0.65  105  
 300  = 8.43E-6 m2/s Sc = 2.46
kc 2/3 2.06
Sc 
V 1780.575
 0.5 teniendo en cuenta el factor 0.79 se obtiene kc = 0.0274 m/s
cA1 = 0; PAS = 5.61/760x101325 = 748 Pa cAS = PAS/RT = 2.59E-4 kgmol/m3
Reemplazando en (i) se obtiene cA2 = 1.74E-4 kgmol/m3. El porcentaje de saturación es

cA2/cASx100 = 67.2 %
CONVECCIÓN NATURAL
Placa vertical: Las correlaciones para la transferencia de calor por convección libre son de naturaleza
similar a las correlaciones para convección forzada en flujos externos. Sin embargo, en los problemas
de convección natural no existe una velocidad de flujo, V , sobre la cual basar un número de
Reynolds. El grupo adimensional que toma el lugar del número de Reynolds en la caracterización de
la convección libre se llama el número de Grashof, Gr, y se define por:
g 2 L3 TS  T
GrL 
2 ; Para el código en EES: Gr=g*beta*L^3*(Ts-Tinf)/nu^2

1  V  V
 
V  1  
d V   dV  Vd   0
V  T  P ; = volumen específico   . Por lo tanto:
1    1     s 
    
  T  P  s T  Ts 
Si este fluido puede considerarse gas ideal entonces, reemplazando

  PM / RT ,   1/ T .
53
Figura A.9
La analogía con transferencia de masa puede mantenerse para geometrías similares si definimos un
“coeficiente de expansión” originado en los gradientes de concentración:
1     s 
 AB  
 s   A   As 
y consecuentemente definimos un número de Grashof para transferencia de masa:
g  AB  2 L3c   AS   A  gL3c     S gx 3  
GrABL =  GrABx  1
2  S 2  2 S
La transferencia de calor por convección libre desde una placa vertical calentada se ilustra en la
Figura A.9. El fluido adyacente a la superficie se calienta y se eleva en virtud de su menor densidad
con respecto al volumen del fluido. Esta capa de fluido ascendente arrastra el fluido frio desde la
región casi inactiva más allá de la superficie calentada, de modo que el caudal másico del fluido
ascendente aumenta con la distancia a lo largo de la placa. Cerca del fondo de la placa, el flujo es
laminar, pero en algún punto, la transición a flujo turbulento puede ocurrir si la placa es lo
suficientemente larga.
Los coeficientes de transferencia de calor para la convección libre de una placa vertical calentada o
enfriada de altura L se pueden obtener a partir de la correlación [Churchill y Chu en Incropera, 9.26].
 
2
0.1  Gr Pr  1012
0.387  Gr Pr 
1/6
 
Nu P  0.825  
9/16 8/ 27
 1   0.492 / Pr   
    (A.36)
Para EES: Nu_p=(0,825+0,387*(Gr*Pr)^(1/6)/(1+(0,492/Pr)^(9/16))^(8/27))^2
54
Precisa en  30%. La longitud característica utilizada en los números de Nusselt y Grashof es la

altura de la placa, L. Las propiedades del fluido, incluido  , se evalúan a la temperatura de la
película.
Cilindros verticales: La ecuación (A.36) también se puede usar para la convección libre de cilindros
verticales si la relación de diámetro a altura es mayor que (35/Gr0.25). Para valores menores de D/L, los
efectos de la curvatura de la superficie se vuelven significativos. Sin embargo, se puede aplicar un
factor de corrección de la siguiente manera [White]
  L  
0.9
NuC  Nu P 1  1.43  0.25  

  DGr   (A.37)
Donde Nu P es el número de Nusselt calculado a partir de la ecuación (A.36).
Para temperatura uniforme en la superficie McAdams recomienda, para placa o cilindro vertical (si el
radio es mucho mayor que el espesor de la capa límite) las siguientes relaciones:
Nu m  C GrL Pr  hm L
n
Nu m 
; k
Donde: L es la altura de la placa o cilindro vertical, Gr es el número de Grashof para transferencia de

calor y  es el coeficiente de expansión volumétrica para el fluido:
Las propiedades del fluido se calculan a la temperatura de película. Para valores del producto GrLPr
entre 104 y 1013 C y n se leen en la Tabla A - 5:
Tabla A-5
Tipo de flujo Rango de GrLPr C n
Laminar 104 a 109 0.59 ¼
Turbulento 109 a 1013 0.10 1/3
En el caso de metales líquidos (Pr < 0.03) para placa plana isotérmica vertical
Num = 0.68(GrLPr2)1/4
Cilindro horizontal: Para la convección libre de un cilindro horizontal, se puede usar una correlación
similar a la ecuación (A.36) [Churchill y Chu].
 
2
105  Gr Pr  1012
0.387  Gr Pr 
1/6
 
NuCH  0.60  8/ 27 
 1   0.559 / Pr 9/16  
    (A.38)
Nuss=(0,6+0,387*Ra^(1/6)/(1+(0,559/Pr_L)^(9/16))^(8/27))^2
Ra=g*beta_L*(Tsm-Tinf)*Do^3/nu_L/alpha_L
En este caso, sin embargo, la longitud característica utilizada en los números de Nusselt y Grashof es
el diámetro del cilindro. Las propiedades del fluido se evalúan a la temperatura de la película.
55
Por McAdams
Nu m  C GrD Pr  hm D g 2 D 3 TS  T
n
Nu m  GrD 
k 2
D es el diámetro externo del cilindro. C y n de la Tabla A - 5:
Tabla A-6
Tipo de flujo Rango de GrLPr C n
Laminar 104 a 109 0.53 ¼
Turbulento 109 a 1012 0.13 1/3
Para cilindros horizontales de temperatura superficial constante en metales líquidos

Num = 0.53(GrDPr2)1/4
Placas Horizontales:
Superficie superior caliente o superficie inferior fría
Nu  0.54  Gr Pr  104  Gr Pr  107

1/ 4
(A.39a)
Nu  0.15  Gr Pr  10  Gr Pr  10
1/3 7 11
(A.39b)
Cara superior fría o cara inferior caliente
Nu  0.27  Gr Pr  105  Gr Pr  1011

1/ 4
(A.39c)
La longitud característica utilizada en los números de Nusselt y de Grashof es el área dividida por el
perímetro. Propiedades a temperatura de película.
Convección natural desde esferas:
1/12
 
0.598  Gr Pr 
1/ 4
 7.44  108 Gr Pr 
Nu  2  4/9 1  
9/16 16/9
1   0.469 / Pr 9/16   1   0.469 / Pr   
      (A.40)
Para el código en EES:
Nu=2+0,598*Ra^0,25/(1+(0,469/Pr)^(9/16))^(4/9)
La longitud característica es el diámetro de la esfera y la correlación es válida para Gr Pr  10 y

13
0.6  Pr  100 . Para valores mayores del número de Prandtl, una mejor coincidencia con los datos
experimentales es obtenida eliminando el término entre corchetes elevado a la potencia 1/12. Las
propiedades del fluido se evalúan a la temperatura de la película.
56
Para convección libre hacia aire a presión atmosférica y temperaturas moderadas (37 C a 81.5 C)
los valores se toman de la Tabla A – 7 [Necati]
Tabla A-7
Geometría L Tipo de flujo Rango hm hm
de GrL Pr [ Btu / hr. pie . F ]
2
[W / m 2 .K ]
Placas Altura Laminar 104 a 109 0.29  T / L  1.42  T / L 
1/ 4 1/ 4
verticales Turbulento 109 a 1013

0.19  T  1.31 T 
1/3 1/3
Cilindros Diámetro Laminar 104 a 109 0.27  T / L  1.32  T / L 

1/ 4 1/ 4
horizontales externo Turbulento 109 a 1012

0.18  T  1.24  T 
1/3 1/3
Para esferas pueden usarse las mismas correlaciones que para los cilindros horizontales, pero la
longitud característica es el radio. Para la última columna se usan unidades del SI y para la penúltima
del sistema inglés. Se pueden adaptar para presiones diferentes a la atmosférica así: unidades inglesas
y flujo laminar, multiplicar por (P/14.7)0.5 para flujo turbulento por (P/14.7)2/3; P en psia. SI y flujo
laminar multiplicar por (P/1.0132)0.5, para flujo turbulento multiplicar por (P/1.0132)2/3. P en bar.
Convección natural y convección forzada combinadas
Cabe señalar que los efectos de convección libre a veces pueden ser importantes en los problemas de
convección forzada, especialmente en el caso del flujo laminar a bajos números de Reynolds. La
importancia relativa de los efectos de convección libre y forzada se puede determinar mediante la
comparación de los números de Reynolds y Grashof de la siguiente manera:
La convección libre es despreciable si Gr/Re2 << 1 y la convección forzada es despreciable cuando

Gr/Re2 >>1. La convección combinada se puede presentar cuando Gr/Re2  1. Se acostumbra
correlacionar los coeficientes para convección combinada, tanto para flujo interno como para externo
por una expresión de la forma Nu  Nu F  Nu N . Para la geometría específica que interese, los
n n n
valores de NuF y NuN se determinan a partir de las correspondientes expresiones para convección
forzada y natural. El signo (+) se utiliza para flujo asistido y para flujo transversal, mientras que el
signo () se usa para flujo opuesto. La mejor correlación de los datos se obtiene para n = 3, aunque
valores como 3.5 y 4 se adaptan mejor para flujo transversal involucrando placas y cilindros o esferas,
respectivamente.
Además, en situaciones donde predomina la convección libre, la transferencia de calor por radiación a
menudo es significativa, incluso a temperaturas relativamente bajas. Esto es una consecuencia del
hecho de que las tasas de transferencia de calor en la convección libre son típicamente mucho
menores que en la convección forzada.
Influencia de la convección natural. La convección natural puede jugar papel importante en los
casos de transmisión de calor durante el flujo interno y laminar de los fluidos, si los perfiles de
velocidad y temperatura se alteran suficientemente por este tipo de convección.
Conducto vertical. Temperatura de pared constante.
57
Martinelli y Boelter (Univ. Calif. – Berkeley – Publis. Engr., 5, 23 (1942)), investigaron la

calefacción de fluidos en flujo interno vertical ascendente y el enfriamiento de los mismos en flujo
interno vertical descendente. En ambos casos las velocidades de los fluidos en las proximidades de las
paredes se incrementan por la convección natural originada en las variaciones de densidad surgidas
por los gradientes de densidad por diferencia de temperaturas. Proponen la siguiente ecuación [12 N°
4]:
Nu = 1.75F1[Gz+0.0722F2(GrSPrSD/L)0.84]1/3 (A41)
Aquí, F1 = (TS T)ML/(TS  T)MA; F2 prevé la disminución de la fuerza convectiva cuando la

temperatura del fluido se aproxima a la de la pared, y se puede aproximar por la correlación:
F2 = 1  0.4145(Nu/Gz) para (Nu/Gz)  1.93
F2 = 5.7143  2.857(Nu/Gz) si (Nu/Gz) > 1.93
(Nu/Gz) = (T2  T1)/(TS  T)MA.
Donde Gz es el número de Graetz: Gz = m’Cp/kz = (/2)PeR/z, siendo Pe = RePr el número de

Péclet y R el radio del conducto.
El número de Grashof, característico en convección natural: Gr = g2D3T/µ2, se basa en el diámetro

de los tubos y la diferencia entre la temperatura de sus paredes y la temperatura másica o media del
fluido.  es el coeficiente de expansión térmica del sistema.
La ecuación (31) se aplica para las siguientes condiciones:
1) Gz se evalúa con las propiedades del fluido correspondientes a la temperatura másica media del
mismo (T1 + T2)/2.
2) GrS y PrS se evalúan con las propiedades del fluido correspondientes a la temperatura TS de la
pared. Para GrS: T = TS T1.
3) Calefacción en flujo ascendente o enfriamiento en flujo descendente.
4) Para calefacción y flujo descendente o enfriamiento y flujo ascendente se cambia +0.0722 por
0.0722.
5) El efecto de la convección natural en tuberías verticales u horizontales con régimen laminar, se

vuelve importante cuando se cumple que:
Re  0.1334(Gr.Pr.d/L)0.875 y Gr.Pr.(d/L)  1x105
Re < 177.8(Gr.Pr.d/L)0.25 y Gr.Pr.d/L > 1x105
Ambas condiciones deben además cumplir que 102 > Pr(D/L) >1.
Conductos horizontales: En el flujo interno de fluidos por tubos horizontales con temperaturas de
pared TS constante, el aumento de temperatura del fluido en las vecindades de la pared determina una
circulación periférica ascendente complementada con una descendente por la parte central. Estos
58
efectos superpuestos al flujo forzado ocasionan el avance en espiral del fluido, con la consiguiente
mezcla adicional y aumento del coeficiente de transmisión de calor individual. Para este flujo Eubank
y Proctor (Thesis, Dep. of Chem. Eng. MIT 1951) proponen:
Nu = 1.75[Gz + 0.04(Gr.Pr.D/L)0.75]1/3[µ/µS]0.14 (A42)
Gz, Gr y Pr se evalúan con las propiedades del fluido correspondientes a la temperatura másica media
del mismo (T1 + T2)/2. Para Gr, T = (TS T)MA.
Para el flujo laminar de gases por el interior de tubos cilíndricos, cualquiera que sea su posición,
puede utilizarse con suficiente aproximación la simplificación propuesta por Cholette (Chem. Eng.
Prog., 44, 81 (1948)) y Kroll, (“Heat Transfer and Presure Drop for Air Flowing in Small Tubes”,
Thesis Dep. Of Chem. Eng. M.I.T. (1951)): Nu = 1.5 Gz0.4. Las propiedades se evalúan a la
temperatura másica media.
FLUJO DE FLUIDOS SOBRE BLOQUES DE TUBOS
Flujo perpendicular a bloques de tubos sin tabiques deflectores (p. ej. en intercambiadores de
coraza y tubos)
Las dos disposiciones usuales de estos bloques son o alineadas o al tresbolillo. En 1933 Colburn [12]
con los datos experimentales disponibles hasta entonces, dedujo la ecuación:
0.6 1
 hDo   Do Gmáx   C P   3
   0 . 33   
 k T f   T f  k T f (A-43)
Condiciones:
1. 10  Re  40000
2. Velocidad másica superficial máxima Vmax  Gmax , basada en la mínima sección transversal de
flujo.
3. Propiedades del fluido evaluadas a Tf = 0.5(T + T0).
4. Disposición de los tubos alternada al tresbolillo.
5. Bloques con al menos diez filas de tubos.
En el caso de bloques de tubos alineados en paso rectangular, Colburn, teniendo en cuenta otros datos
de la bibliografía y los suyos propios, indica que debe sustituirse la constante de proporcionalidad
(factor de forma) 0,33 por 0,26.
En 1937, Grimison, a partir de datos experimentales, propuso una ecuación concreta para el caso de
que el gas fuera aire, que generalizaba para cualquier otro gas de la siguiente forma:
n 1
 hDo   D G  C   3
   c o máx   P 
 k T f   T f  k T f (A-44)
Condiciones:
59
1. 2000  Reo  40000.

2. Velocidad másica superficial máxima Gmáx, basada en la mínima sección transversal de flujo.
3. Propiedades del fluido evaluadas a Tf = 0.5(T + T0).
4. Bloques con al menos 10 filas de tubos.
5. Valores de c y n tomados de la Tabla A - 7
Tabla A-7
Disposición b/Do 1.25 1.50 2.00 3.00
c n c n c n c n
a/Do
Alternada 0.6 0.241 0.636
0.9 0.504 0.271 0.453 0.581
1.0 0.561 0.558
1.125 0.540 0.565 0.585 0.560
1.250 0.585 0.556 0.570 0.554 0.586 0.556 0.589 0.562
1.50 0.509 0.568 0.519 0.562 0.511 0.568 0.591 0.568
2.0 0.456 0.572 0.470 0.568 0.544 0.556 0.507 0.570
3.0 0.350 0.592 0.402 0.580 0.497 0.562 0.476 0.574
Alineada 1.25 0.393 0.592 0.310 0.608 0.113 0.704 0.071 0.552
1.50 0.414 0.586 0.282 0.620 0.114 0.702 0.076 0.755
2.0 0.472 0.570 0.337 0.620 0.259 0.632 0.224 0.648
3.0 0.327 0.601 0.403 0.584 0.423 0.581 0.323 0.608
Schneyder midió coeficientes de transmisión de calor individuales medios de tubos correspondientes

a distintas filas de bloques de ellos alternados, sobre los que perpendicularmente fluía aire. En líneas
generales, sus resultados coincidieron con los alcanzados previamente por otros investigadores al
determinar los coeficientes locales sobre la periferia de los indicados tubos de filas distintas del
bloque. Los números de Nusselt medios aumentan hasta la tercera fila, decreciendo luego ligeramente
y permaneciendo prácticamente constantes a partir de la quinta fila. Pudo expresar sus resultados
mediante la misma ecuación, pero con los valores de c y n que se indican en la Tabla A - 8, para 8000
 Re0  20000, a/Do = 1.8 = b/Do; disposición alternada al tresbolillo.
Los valores de h calculados con la ecuación de Colburn resultan algo menores que los calculados con
la ecuación de Grimison que se considera de mayor precisión. Los valores de h que se encuentran con
esta última ecuación para las filas 5 a 10 de un bloque de tubos son aproximadamente 80 por 100
superiores a los coeficientes medios de los tubos de las primeras filas. En el caso de bloques de menos
de diez filas de tubos, los valores de h calculados con las ecuaciones anteriores deberán multiplicarse
por los factores que se indican en la tabla A-9.
Tabla A-8
Fila número c n
1 0.226 0.588
2 0.215 0.608
3 0.183 0.638
4 0.191 0.640
5 0.169 0.654
6 0.201 0.640
7 0.234 0.625
8 0.240 0.624
9 0.244 0.622
10 0.252 0.620
60
Tabla A-9
Número de filas 1 2 3 4 5 6 7 8 9 10
Tresbolillo 0.68 0.75 0.83 0.89 0.92 0.95 0.97 0.98 0.99 1.0
Alineadas 0.64 0.80 0.87 0.90 0.92 0.94 0.96 0.98 0.99 1.0
Flujo sobre bloques de tubos con tabiques deflectores
Se da este caso en el flujo de fluidos, por la parte de la carcasa, externamente a los tubos de los
cambiadores de calor multitubulares.
El estudio más completo sobre pérdidas de carga y transmisión de calor en este tipo de flujo se debe a
Tinker. Donohue, basándose en sus datos y en los de otros investigadores, propuso las siguientes
correlaciones para los cambiadores multitubulares:
- Parte de la carcasa de un cambiador multitubular sin tabique deflector alguno:
0 .6 0.33 0.14
hDO  D G  C     
 0.00825De0.6  o c   P   
k     k   o  (A-45)
Condiciones:
1. 2.5  De (DoGc/)  500
2. De= 4(área sección transversal/perímetro mojado) = 4(ab  Do/4)/Do expresado en metros,
teniendo los parámetros a, b el significado del paso vertical y horizontal entre tubos para el bloque de
tubos en el interior de la carcasa, y Do el diámetro exterior de los tubos.
3. Gc velocidad másica del fluido en la carcasa, libre de tabiques deflectores, referida a la superficie
libre de flujo.
4. Propiedades físicas del fluido evaluadas a su temperatura másica, excepto 0 (a la temperatura de
las paredes de los tubos).
- Parte de la carcasa de un cambiador multitubular con tabiques deflectores de discos-coronas:
0.6 0.33 0.14

hDO  D G  C     
 0.005842 De0.6  o e   P   
k     k   o  (A-46)
Condiciones:
1. 0,5  De (DoGe/)  700.
2. De como antes.
3. Ge = (GcorGcarc)0.5·siendo Gcor la velocidad másica referida a la abertura de las coronas y Gcarc la
velocidad másica referida a la sección transversal de la carcasa.
4. Propiedades físicas del fluido a su temperatura másica, excepto 0.
- Parte de la carcasa de un cambiador multitubular con tabiques deflectores de segmentos circulares:
61
0.6 0.33 0.14

hDO  D G  C     
 0.25 o e   P   
k     k   o  (A-47)
Condiciones
1. 10  (DoGe/)  5000
2. Ge = (GsegGcarc)0.5 siendo Gseg la velocidad másica referida a la abertura del segmento; Gcarc la
velocidad másica referida a la sección transversal de la carcasa.
3. Propiedades físicas del fluido a su temperatura másica, excepto 0
William y Katz (Trans. ASME, 74, 1307 (1952)) llegaron a una ecuación similar a la anterior con
valores distintos de la constante de proporcionalidad o factor de forma, que precisan para cada caso
concreto, incluso cuando las paredes de los tubos del cambiador están extendidas con aletas.
Condensación tipo película de Vapores Puros
Placas inclinadas y tubos verticales: A partir de la teoría de Nusselt modificada por McAdams
(Necati p 387), el coeficiente promedio de transferencia de calor para la condensación sobre placas
inclinadas un ángulo  con la horizontal se obtiene de:
1/ 4
  l g sin  ( l   v )h fg kl3 
hm  1.13  
 l (Tv  Tw ) L  (A-48)
También, en función del número de Reynolds para la película, despreciando v respecto a l
1/3
  2f 
1/3
 4   4 
hm  3 2   1.76   
 k f  f g sin    f   f 
    Ref =  1800 (A-49)
Observe que las unidades del coeficiente convectivo dependen de las unidades (consistentes) en que
se reemplacen los diferentes parámetros de los paréntesis. Si la relación entre el diámetro del tubo y el
espesor de la película de condensado son grandes (> 10) estas ecuaciones también pueden usarse para
condensación en tubos verticales.
La cantidad  es el caudal másico de condensado por unidad de ancho (ver películas descendentes, p
72)
62
Tubos horizontales: La ecuación para el coeficiente promedio de transferencia de calor para

condensación en el exterior de tubos horizontales (7.13 en Manrique; 10.14 en Karlekar y Desmond;
14.18 en Necati Ösizik; 21.30 en Welty Wicks y Wilson; 10.41 en Incropera DeWitts, etc.) es:
1/ 4
  l g ( l   v )h fg kl3 
hm  0.725  2/3 
 l (Tv  Tw ) DN  (A-50)
Aquí D es el diámetro exterior del tubo, Tv, Tw son las temperaturas del vapor y de la superficie de la
pared respectivamente; hfg es el calor latente de condensación. Los subíndices “l” indican que la
propiedad es para la fase líquida. N es el número de tubos dispuestos verticalmente en un arreglo
multitubular.
Si se desprecia la densidad del vapor frente a la del líquido, es decir (l v) l podemos escribir la
ecuación para un solo tubo como:
1/3 1/3
 2   4 
hm  3 l 2   1.514  
 kl  i g   l  (A-51)
 es caudal másico de condensado por unidad de longitud del tubo, dado por:
hm DL Tv  Tw 

h fg L
CONDENSACION AL INTERIOR DE TUBERÍAS: Para determinar el coeficiente convectivo de

transferencia de calor por condensación en película en el interior de tuberías horizontales, Kern
(página 269 - primera edición) recomienda utilizar la misma correlación usada para determinar el
coeficiente pelicular para condensación en el exterior de tubos horizontales tomando un caudal
másico de condensado por unidad de longitud,  doble del real. Esto equivale a escribir la ecuación
(12 - 40) del texto en cuestión de la siguiente manera:
1/3 1/3
 l2   8 
hm  3 2   1.514  
 kl  l g   l  (A-52)
63
Aquí  es el caudal másico actual por unidad de longitud. Demuestre que esto es equivalente a
escribir la ecuación para condensación en el exterior de tubos horizontales en la siguiente forma:
1/ 4
  l g ( l   v )h fg kl3 
hm  0.725  
 2 l (Tv  Tw ) D  (A-53)
RADIACIÓN: Todo cuerpo a una temperatura absoluta finita emite radiación electromagnética. Esta
radiación, cuando está en el rango de longitud de onda comprendido entre los 0.2 m y los 100 m se
denomina térmica. Cualitativamente puede explicarse su origen a variaciones en los estados
electrónico, vibracional y rotacional de átomos o moléculas. Conforma solo una pequeña parte de
todo el espectro de radiación e incluye parte de la radiación ultravioleta, la radiación visible (0.35 a
0.78 m) y parte del infrarrojo.
Figura 1.7: Espectro Electromagnético Parcial
Aunque nos vamos a centrar en la radiación desde superficies sólidas, ésta puede ocurrir también en
líquidos y gases.
64
Figura 1.7
El cuerpo negro: Se denomina así una superficie que en todas las longitudes de onda emite y absorbe
la máxima cantidad posible de radiación. A partir de la segunda ley de la termodinámica Max Planck
(1900) dedujo que la potencia con la que un cuerpo negro emite energía en una longitud de onda y a
una temperatura dadas, en el vacío, viene dada por
C
Eb  1
 e
5 C 2 / T

1 (1.12)
Eb = Potencia emisiva espectral o monocromática de un cuerpo negro a una temperatura y longitud de
onda dadas, W/m2.m,
 = Longitud de onda de la energía emitida, m,
T = Temperatura absoluta de la superficie, K,
C1 = 3.7405x1016 W.m2,
C2 = 0.0143879 m.K.
En la figura 1.7 se muestra una línea que une los máximos a diferentes temperaturas. Diferenciando la
ley de Planck e igualando a cero se obtiene su ecuación conocida como ley del desplazamiento de
Wien:
maxT  2897.8  m.K  5215.6  m. R
La potencia a la cual una superficie emite energía radiante en todas las longitudes de onda, por unidad
de área se denomina potencia emisiva superficial E. Existe un límite superior para esta potencia
emisiva, que está prescrito por la ley de Stefan Boltzmann:
65

Eb   Eb d  TS4
0 (1.13)
Donde TS es la temperatura absoluta de la superficie y  es la constante de Stefan Boltzmann

(5.6697x108 W/m2.K4; 0.1714x108 Btu/hr.pie2.R4). Este tipo de superficies se denominan radiadores
ideales o cuerpos negros.
En algunas ocasiones es necesario calcular la energía emitida en una banda específica de longitudes
de onda lo que se consigue integrando la ecuación para Eb entre los respectivos límites. La fracción
de la energía total emitida en estas franjas de longitudes de onda se obtiene así:
2

C
2  1
d
 e  1
 C2 /  T ( T ) 2
 Eb d   5
 
1 1  C1d (T )
 
 (T )5  eC2 / T  1
 
T 4

 Eb d  
( T )1
0 (1.14)
Esta expresión sólo depende de  T  .
El flujo emitido por una superficie real es menor que el emitido por una superficie ideal a la misma
temperatura y está dado por
E  TS4 (1.13a)
Aquí  es una propiedad radiativa de la superficie denominada emisividad que tiene valores en el
rango 0    1, esta propiedad nos da una medida de la eficiencia con que la superficie emite
energía radiante. Depende fuertemente del material y de su terminado. Si   1 la superficie se llama
gris.
La radiación también puede incidir sobre una superficie desde sus

alrededores. La radiación puede originarse desde una fuente
especial, tal como el sol, o desde otra superficie hacia la cual la
superficie de interés está expuesta. Indistintamente de cual sea la
fuente nosotros designamos la velocidad a la cual esta radiación
incide por unidad de área de la superficie como la irradiación G.
Esta irradiación puede ser absorbida, transmitida y/o reflejada en

proporciones  ,  ,  , respectivamente. Por la ley de Kirchhoff
 =  si las temperaturas de absorbedor y emisor son del mismo
orden de magnitud. Si la fracción  (transmisividad) transmitida
es cero, el cuerpo se considera opaco; si  = 1, entonces será
transparente a la irradiación. Si la absortividad  = 1 será negro.
Una porción o toda la irradiación puede ser absorbida por la

superficie aumentando la energía térmica del material. La velocidad a la cual la energía radiante es
absorbida por unidad de área superficial puede evaluarse conociendo la absortividad  . Esto es, Gabs=
 G con 0    1. Si  < 1, y la superficie es opaca, la otra porción de la irradiación es reflejada.
Si la superficie es semitransparente porciones de la irradiación pueden ser transmitidos. Sin embargo,
66
mientras que la radiación absorbida o emitida aumenta o reduce, respectivamente, la energía térmica
de la materia, la radiación reflejada y transmitida no tiene efecto en ésta energía. Nótese que el valor
de  depende de la naturaleza de la irradiación así como de la superficie misma. Por ejemplo, la
absortividad de una superficie a la radiación solar puede diferir de su absortividad a la radiación
emitida por las paredes de un horno.
Error de Termocupla
Una termocupla o cualquier otro dispositivo usado para medir la temperatura de un gas en un
conducto puede dar una lectura significativamente diferente de la temperatura del fluido si las paredes
tienen otra temperatura. Esta diferencia ocurre debido a que el elemento detector estabiliza su
temperatura intercambiando calor con el fluido por convección y con las paredes por radiación (y
eventualmente por conducción).
Figura 1.9 error de termocupla
Si, por ejemplo, la temperatura de la pared Ts es menor que la temperatura del gas Tg, la termocupla
indicará algún valor intermedio Tt. El balance de calor será
hAt Tg  Tt    At Tt 4  Ts4 
Este balance sigue siendo válido para el caso contrario, cuando el gas está más frío que la superficie,
pues bastará multiplicar por 1 ambos lados de la ecuación. La temperatura mostrada por el detector
de temperatura Tt siempre flotará entre los otros dos valores Tg y Ts .
Cuando lo que se requiere es conocer el valor de la temperatura del gas esta se puede despejar de
manera explícita de la ecuación, pero si se trata de la temperatura del termómetro será necesaria
prueba y error o algún método numérico.
EJEMPLO 1.2: Fluye aire a través de un conducto cuyas superficie interior tiene temperatura de
500 °F. Una termocupla alojada en un pozo de acero oxidado, de emisividad  = 0.9, colocado
perpendicular a la corriente gaseosa, indica una temperatura de 300 °F. La velocidad másica
superficial del aire es de G = 3600 lb/hr-pie2, y el diámetro externo del pozo del termómetro es 1/2
plg. (a) Estime la verdadera temperatura del aire en el conducto. Suponga que la termocupla y el pozo
están a la misma temperatura. Para aire en flujo perpendicular a cilindros entre temperaturas de 100
°F a 300 °F, el coeficiente convectivo de transferencia de calor puede ser estimado con la siguiente
ecuación dimensional (Bennet y Myers p 445):
67
0.6
G
h  0.026
D 0.4 [Btu/hr-pie2-°F] ; G [lb/hr-pie2]; D [pie] (i)
El coeficiente h se obtiene en Btu/hr-pie2-°F; G es la velocidad másica del aire en lb/hr-pie2 y D es el

diámetro exterior del cilindro en pie. Compare el resultado obtenido con el que se consigue usando las
ecuaciones correspondientes del anexo A.
(b) Si se pule la superficie del portatermómetro de tal manera que su emisividad térmica disminuya
hasta  = 0.1, ¿cuánto marcará la termocupla?
Solución: La temperatura indicada por el termómetro, Tt , está influenciada por la radiación recibida
desde la superficie y por la conducción a lo largo de las paredes del portatermómetro. Su temperatura
fluctuará entre la de la superficie Ts y la del gas Tg :
Tg  Tt  Ts
(a) En estado estable, desestimando el efecto conductivo se establece el balance:
 Ts4  Tt 4   h Tt  Tg 
(ii)
Las temperaturas deben ser absolutas por tener radiación. Como alternativa a la ecuación que propone
el problema se puede estimar el coeficiente convectivo h utilizando la ecuación propuesta por
Churchill y Bernstein (A.17):
4/5
0.62 Re 0D.5 Pr 1 / 3   Re D  
5/8
Nu D  0.3  1    

1  0.4 / Pr 
2 / 3 1/ 4

  282000  
Usando la ecuación propuesta en el enunciado se encuentra
 3600 
0.6
h  0.026
 0.5 / 12 
0.4
= 12,61 [Btu/hr-pie2-°F]
La ecuación (ii) nos da de manera directa el valor de Tg :
Tg  760  1.712 x109 x0.9  9604  7604  / 12.61  697 °R
(b) En este caso, la incógnita es Tt que no puede obtenerse explícito de la ecuación (ii). Se sugiere la
siguiente técnica de solución: a partir de (ii) se escribe Tt en función de Tt .
4
Tt  Tg   Ts4  Tt 4  / h
(iii)
68
Las iteraciones se hacen suponiendo un valor para Tt que se reemplaza en el lado derecho de la
ecuación (iii) calculando un valor resolviendo la misma (iii), y así hasta convergencia. Esta se puede
acelerar usando cualquier técnica numérica adecuada.
Usamos el software EES:
"error de termocupla"
$UnitSystem Eng Mass R Btu Deg atm "Se fijan las unidades"
"Datos"
sigma=sigma#: epsilon=0,9: T_s=500+460 [R]: T_t=300+460 [R]
"Correlaciones"
sigma*epsilon*(T_s^4-T_t^4)=h*(T_t-T_g) "Ecuación (i)"
G=3600 [lb/hr-ft^2]: D=0,5/12 [ft]
h=0,026*G^0,6/D^0,4 "Ecuación (ii)"
"Coeficiente usando la ecuación de Churchill y Bernstein (A.17)"

hb*D/k=0,3+0,62*Re^0,5*Pr^(1/3)/(1+(0,4/Pr)^(2/3))^(1/4)*(1+(Re/282000)^(5/8))^0,8
Re=G*D/mu
mu=Viscosity(Air;T=T_g)
k=Conductivity(Air;T=T_g)
Pr=Prandtl(Air;T=T_g)
Resultados para la temperatura del gas en el primer caso ( = 0.9)

D=0,04167 [ft]; h=12,61 [Btu/hr-pie2-°F]
sigma=1,712E-09 [Btu/hr-ft^2-R^4]
T_g=697 [R] = 237 °F; T_s=960 [R] = 500 °F; T_t=760 [R] = 300°F; T_t –T_g = 63 [°F]
Usando la correlación (A.16) se obtiene un resultado que concuerda muy bien con la ecuación
propuesta:
hb=11,95 [Btu/hr-pie2-°F]; k=0,01846 [Btu/hr-ft-R]; mu=0,05418 [lb_m/ft-hr]; Pr=0,7093; Re=2769
Resultados para temperatura del termómetro en el segundo caso ( = 0.1) y T_g=697 [R] = 237 °F:
T_s=960 [R] = 500 °F; T_t=705 [R] = 245 °F; T_t –T_g = 8 [°F]
El ejemplo 1.2 muestra que el error en la lectura puede disminuirse cambiando la superficie del sensor
para que tenga baja emisividad y absortividad. Un método alterno es usar escudos de radiación, por
ejemplo cubrir la termocupla con un tubo más pequeño, orientado paralelamente al flujo del gas, que
permita a este pasar a través de él y sobre la termocupla pero protege de la radiación a la termocupla.
El escudo intercambia calor con el gas y con la superficie y alcanza alguna temperatura intermedia.
La termocupla intercambia con el gas y con el escudo y da un valor más aproximado al verdadero
para el gas. Otro factor de error se origina en bajas velocidades del gas lo que influye para bajo
coeficiente convectivo. De otra parte, si el gas contiene moléculas de tres o más átomos, estas pueden
absorber y emitir cantidades importantes de radiación, pero esta circunstancia favorece la lectura del
termómetro. Puede además introducirse errores por conducción a lo largo del vástago del termómetro
o de los alambres de la termocupla o a lo largo de la pared del portatermómetro que actuaría como
una aleta hueca.
EJEMPLO 1.3: El Sol puede suponerse que actúa como un cuerpo negro a 5800 K. Calcular:
a) El poder emisor total.
b) La longitud de onda a la que se consigue el poder emisor máximo
69
c) El poder emisor monocromático máximo.

d) Porcentaje de energía total emitida que corresponde a longitudes de onda del espectro visible.
Solución
a) El poder emisivo total de una superficie negra está dado por la ley de Stefan - Boltzmann:
Eb = T4 = (5,67x108)(5800)4 = 6.42x107 J/s.m2
b) Según la ley del desplazamiento de Wien:
maxT = 2,8978x103 m.K  max = 2,8978x103/5800 = 5.0x107m = 0.50 m.
c) De acuerdo con la ecuación (1.12), el poder emisor monocromático máximo es:
C1 3.7405 x1016 J
Eb max    8.41x1013
5
e C2 /  max T
 1  5.0 x10   e
7 5
 s.m 2
3
max
0.0143879/ 2.8978 x10
1
d) El porcentaje de energía total emitida dentro del espectro visible (0.35 - 0.75 m) teniendo en
cuenta la ecuación (1.14) podrá expresarse como:
 T   4.3510 
 3.7405  1016  d  T 
3
 2 C1d  T  
 
 T 1   T   e  1  5.6697  108   T   e0.0143879/ T  1
5 C2 /  T
 5
=  2.03103 
= 0.4692
Es decir que el 46.92% de la energía radiada por el sol se emite en el rango visible del espectro. La
integral se evaluó numéricamente usando el método de Romberg (ver anexo al final del capítulo).
Analíticamente puede integrarse la ecuación (1.14) realizando el cambio de variable C2/x = z, con x =
T. Así, dx =  (C2/z2)dz. Reemplazando, agrupando constantes, aplicando desarrollo en series de
Taylor y adecuando los valores límites se obtiene una expresión de la forma
z2 z2
  z
3
z d ( z) z 3e  z d ( z )
 z e 1  e  e  e 3 z   dz
2
C   C   C  3 z z 2 z


z1
 e z  1 

z1
1  e  z1
z 
Cada una de las integrales resultantes se integra por partes. También pueden usarse las tablas
tabuladas por Dunkle (1954) y reproducidas en varios textos de transferencia de calor como
“funciones de radiación del cuerpo negro”.
Usando el EES obtenemos
F=Integral(C_1/(sigma*x^5*(exp(C_2/x)-1));x;2,03e-3;4,35e-3) "Ecuación (1.14)"

C_1=3,7405e-16 [W-m^2]
C_2=0,0143879 [m-K]
sigma=5,67e-8 [W/m^2-K^4]
Solución F=0,46923 Fracción de energía emitida en el rango visible.
70
EJEMPLO 1.4: Un invernadero está construido con un vidrio que posee una tramitancia de 0.92 para
todas las longitudes de onda comprendidas entre 0,35x106 m y 3x106 m y es opaco al resto de las
radiaciones. El Sol se comporta como un cuerpo negro a 5800 K, mientras que el interior del
invernadero, que también puede considerarse negro, se encuentra a 300 K. Determinar:
a) Tramitancia total del vidrio del invernadero a la radiación solar.
b) Tramitancia total del vidrio a la radiación interior del mismo.
c) Flujo neto de calor por radiación que recibe el invernadero.
Suponga que el invernadero recibe del sol un flujo de energía radiante de 1350 J/s.m2.
Solución
a) De forma paralela a la ecuación (1.14) que nos determina la fracción de la energía total emitida en
un rango de longitudes de onda a una temperatura dada por un cuerpo negro, podemos obtener la
fracción transmitida, absorbida o emitida por un cuerpo gris conociendo ,  o . Así, la tramitancia
total de esta superficie vendrá dada por
1T   0.35  106   5800   2.03  103
2T   3  106   5800   17.4  103
( 1T ) (17.4 x103 )

 
C1d (T ) C1d (T )
      0.92 
 0.8347
 ( T )  e  1  (T )5  eC2 / T  1
5 C2 /  T 
 
 
( 1T ) (2.03 x103 )
En esta ocasión  es independiente de  y puede sacarse de la integral.
b) Cuando la fuente de radiación es el interior del invernadero, análogamente al apartado anterior se

llega a:
2T = (3x106)(300) = 0.9x103 m.K: 1T = (0.35x106)(300) = 0.105x103 m.K
( 2T ) (0.9 x103 )

 
C1d (T ) C1d (T )
      
 8.00  105
 ( T )  e  1  ( T )  e  1
5 C2 /  T  5 C2 /  T
 
 
( 1T ) (0.105 x103 )
c) El flujo neto de calor por radiación que recibe el invernadero será entonces
qneto = (0.835)(1350) x)(5.67x108)(300)4 = 1127.2 J/s.m2
Puede apreciarse como el vidrio del invernadero transmite el 83,5 % de la energía del Sol al interior
del mismo, mientras que apenas transmite el 0.008 % de la energía radiante del invernadero hacia el
exterior.
EJEMPLO 1.5: Una bombilla incandescente es un aparato económico pero muy ineficiente que
convierte la energía eléctrica en luz. Transforma alrededor de 10% de la energía eléctrica que
consume en luz, mientras que convierte el 90% restante en calor. El bulbo de vidrio de la lámpara
71
se calienta con mucha rapidez, como resultado de absorber todo ese calor y disiparlo hacia los
alrededores por convección y radiación. Considere una bombilla de 60 W, de 8 cm de diámetro,
en un cuarto a 25 °C. La emisividad del vidrio es de 0.9. Suponiendo que 10% de la energía pasa
a través del bulbo de vidrio como luz, con absorción despreciable, y que el resto de ella es
absorbida y disipada por el propio bulbo por convección natural y radiación, determine la
temperatura equilibrio de dicho bulbo de vidrio. Suponga que las superficies interiores del cuarto
se encuentran a la temperatura ambiente.
Solución: Para convección natural alrededor de esferas, Churchill, S. W. (“Free Convection

Around Inmersed Bodies”) en Incropera et al. 6ª ed. p583, ch 9, recomienda para fluidos con Pr
 0.7 y Ra  1011, la correlación:
0.589 Ra1/ 4
Nu  2  4/9
1   0.469 / Pr 9/16 
 
con Raleygh dado por
g  TS  T D 3
Ra  Pr
2
Todas las propiedades se calculan a la temperatura de película T f  TS  T  / 2 . Por tanto, para
gases que puedan considerarse ideales
1

Tf
Como la temperatura de la superficie es desconocida se requiere asumir algún valor para calcular
las propiedades y hacer una primera iteración. Asumiendo un valor de 80 ºC, la temperatura de
película sería 326 K. En la tabla A4 del Incropera et al, hay valores a 300 K y 350 K. Parece fácil
interpolar linealmente para 325 K así:
k = 0.02815 W/m.K; Pr = 0.704;  = 18.41x10-6 m2/s;  = 1/325 = 3.08x10-3 K-1.
En estado estable, toda la energía suministrada a la bombilla se entrega al medio, parte como
energía luminosa, parte como calor. El balance de energía térmica se reduce a:
0.90  Potencia = hS TS  T    S TS4  Talrr

4

Todos los valores menos TS se conocen. La forma geométrica se puede aproximar a la de una
esfera con diámetro 8 cm. La solución de la ecuación requiere prueba y error o un Newton
Raphson.
Usando el software EES aprovechando la base de datos de este para calcular las propiedades
requeridas. Para facilitar el cálculo, en Options, Información de variables se da un guess para la
TS y la de película de acuerdo con lo esperado según la experiencia y el sentido común. Se
obtiene, por el software temperatura de 169.4 °C, bastante diferente del valor supuesto. Las
72
modificaciones en las propiedades del aire las modifica en las iteraciones internas el mismo
software.
Como en la práctica, la temperatura de la superficie de este objeto es del orden de 80 °C, se

concluye que el vidrio del que está hecho el bulbo de la bombilla es transparente a la mayor parte
de la radiación térmica emitida por el filamento de tungsteno. Como solo el porcentaje absorbido
afecta la temperatura del objeto, ¿cuál será el valor de esa absortividad?
En la solución solamente se cambia el valor de la absortividad igual a la emisividad. Se observa
que la transmisividad de la bombilla es del orden del 86 % a la energía térmica emitida por el
filamento.
$UnitSystem SI Mass J K Deg Pa

"Problema 9-47 Cengel: bombilla"
"Datos"
Q=alpha*60 [W]
D=0,08 [m]
T_infinity=298 [K]
epsilon=alpha
"Propiedades"
Fluid$='air'
k=Conductivity(Fluid$; T=T_film)
Pr=mu*Cp/k
Cp=1008
rho=Density(Fluid$; T=T_film; P=101300)
mu=Viscosity(Fluid$; T=T_film)
nu=mu/rho
beta=1/T_film
T_film=1/2*(T_s+T_infinity)
g=9,807 [m/s^2]
sigma=5,67E-8 [W/m^2-K^4] "constante de Stefan-Boltzmann"
T_sc=80
"Correlaciones"
Ra=(g*beta*abs(T_s-T_infinity)*D^3)/nu^2*Pr
Nusselt=(2+(0,589*Ra^(1/4))/(1+(0,469/Pr)^(9/16))^(4/9))
h=k/D*Nusselt
A=pi*D^2
Q=h*A*(T_s-T_infinity)+epsilon*A*sigma*((T_s)^4-(T_surr)^4)
T_surr=T_infinity: T_sc=T_s-273
Resultados: A=0,02011 m2; beta=0,0027 K-1; Cp = 1008 J/kg-K; h= 7.84 W/m2K; k=0,03074 W/m.K
mu=0,00002169 kg/m.s; nu=0,00002275 m2/s; rho= 0.9533 kg/m3
Nusselt=20.41; Pr=0,71; Ra=2.692E+06 Q=54 W
T_film=370.2 K; T_infinity=298 K; T_s = 442.4 [K] (169.4 ºC); T_surr=298 K.
Como comprobación de los valores dados por la tabla A6 recalculamos:
1/ 4
 T 
h  1.32  
 L  = 10.23 W/m2.K
73
Este resultado es 30 % mayor que el valor predicho por la correlación de Churchill, pero obtenido con
significativamente menos esfuerzo, lo que lo hace útil para cálculos preliminares.
Para el segundo análisis, la incógnita es la absortividad. Se obtienen los siguientes valores:

alpha=0,1394; beta=0,003072 K-1; Cp=1008 J/kg-K; epsilon=0,1394; rho=1,084 kg/m3.
h=6,467 W/m2.K; k=0,02753 W/m-K; mu=0,00001974 kg/m-s; nu=0,00001821 m2/s
Nusselt=18,8; Pr=0,7229; Q=8,366 W (calor absorbido por el vidrio); Ra=1,850E+06.
T_film=325,5 K; T_infinity=298 [K]; T_s=353 K; T_sc=80 °C
EJEMPLO 1.6: Calentamiento de una placa por energía solar (Cengel Ejemplo 1-13 p 34
modificado). Una placa metálica delgada está aislada en la parte posterior y expuesta a la radiación
solar en la superficie del frente. La superficie expuesta de la placa tiene una absortividad de 0.6, para
la radiación solar. Si la radiación solar incide sobre la placa a una rapidez de 700 W/m2 y la
temperatura del aire circundante es de 25 °C, determine la temperatura de la superficie de la placa
cuando la pérdida de calor por convección y radiación es igual a la energía absorbida por la propia
placa. Suponga que el coeficiente de transferencia de calor por convección es de 30 W/m2.°C.
Figura 1 ejemplo 1.6
Solución: Del calor que incide, solo afecta la temperatura de la superficie el que se absorbe. La
superficie pierde hacia los alrededores por convección y radiación. En estado estable lo absorbido es
igual a lo perdido. Para evaluar el efecto radiante separamos la convección y la radiación:
 qA  hA T  T f    A T 4  T f4 
La incógnita T , la temperatura de la superficie, no puede expresarse explícitamente a partir de esta

ecuación. Debe usarse prueba y error o algún método numérico. Lo más práctico es dejar T en
función de T 4 y hacer prueba y error usando valores lógicos ( T debe tener algún valor mayor al del
ambiente):
T  T f   q / h   T 4  T f4  / h
74
Como se usa la ecuación para transferencia de calor por radiación las temperaturas deben ser
absolutas. En esta ocasión se transfiere calor hacia alrededores que tienen la misma temperatura del
ambiente. El software EES nos da la solución ddirectamente:
Usando EES, el balance hacia y desde la superficie se escribe simplemente:
$UnitSystem SI Mass Pa J K
alpha*A*q=h*A*(T-Tf)+sigma*epsilon*A*(T^4-Tf^4) "balance de energía"
"datos"
A=1 [m^2]: alpha=0,6 [-]: epsilon=0,6 [-]: h=30 [W/m^2-K]: q=700 [W/m^2]: Tf=298 [K]
sigma=sigma#: Tc=converttemp(K;C;T)
hm=hr+h "coeficiente combinado"
hr=sigma*epsilon*(T^2+Tf^2)*(T+Tf) "coeficiente por radiaación"
Si en la pestaña Windows hacemos click en Formatted Equations aparece
Resultados
hm=33,83; hr=3,832 [W/m^2-K]; sigma=5,670E-08 [W/m^2-K^4]; T=310,4 [K]; Tc=37,26 [C]
Se coloca "h=30" entre comillas para poder explorar como cambia la temperatura de la superficie al
cambiar el coeficiente convectivo creando una tabla paramétrica. En Tables, New Parametric Table
para variar entre 10 y 50 en intervalos de 5 K nos da 8 intervalos o 9 valores por lo que cambiamos
los "runs" de 10 que aparece por defecto a 9 y con Add pasamos a la derecha las variables que
queremos en la tabla que son h, hm, hr y Tc. Damos 0k y en la pestañita de h abrimos una ventana de
diálogo que nos permite introducir el primer valor, 10 y el último 50 con lo que la tabla se llena
automáticamente. Luego vamos a Calculate, Solve Table y obtenemos:
h hm hr Tc
[W/m^2-K] [W/m^2-K] [W/m^2-K] [C]
10 14,17 4,174 54,48
15 19,02 4,021 46,93
20 23,93 3,932 42,4
25 28,87 3,873 39,4
30 33,83 3,832 37,26
35 38,8 3,802 35,67
40 43,78 3,779 34,44
45 48,76 3,76 33,46
50 53,75 3,745 32,66
Para copiarla seleccionamos el primer cuadro (10) y sosteniendo la tecla para mayúsculas hacemos
click en la última casilla (305,8). Haciendo click con el botón secundario del ratón, seleccionamos
copiar con headers (encabezados) y lo pegamos en el documento word. La tabla aparece con 4
columnas y 11 filas. Para ocupar menos espacio se copia en Excel, se selecciona un espacio con 4
filas y 11 columnas y se da el comando: "=transponer(xx:yy)" donde los términos en el paréntesis
corresponden con las coordenadas del primer y último recuadro de la tabla copiada (donde están h y
75
32,66 en este caso), y se acciona ctrl+shift+enter obteniéndose la tabla tal como se muestra,
horizontal.
h [W/m^2-K] 10 15 20 25 30 35 40 45 50
hm [W/m^2-K] 14,17 19,02 23,93 28,87 33,83 38,8 43,78 48,76 53,75
hr [W/m^2-K] 4,174 4,021 3,932 3,873 3,832 3,802 3,779 3,76 3,745
Tc [C] 54,48 46,93 42,4 39,4 37,26 35,67 34,44 33,46 32,66
Para hacer el gráfico vamos a Plots XY se especifica qué variable deseamos en abscisas y cuál(es) en
ordenadas y listo. Podemos explorar otros detalles como rango, grids etc.
50
h [W/m2.K]
40
30
20
70
65
60
55
50
Tc
45
40
35
Figuras 2 y 3 ejemplo 1.6
Otra opción está en crear un informe en PDF para lo cual es necesario guardar el trabajo con nombre.
Luego en File Create Latex/PDF Document.
EJEMPLO 1.7: (Cengel p 90) Se deja una plancha sobre una tabla de planchar con su base expuesta al
aire. La placa base tiene espesor L = 0.5 cm y área para transferencia de calor 0.02 m2. La
conductividad térmica de esta placa es constante de 15 W/m.K y su emisividad es 0.6. La resistencia
provee una potencia de 1250 W de los cuales el 80 % efectivamente pasan a través de la placa. El aire
circundante y los alrededores se encuentran a 20 ºC, y el coeficiente convectivo se estima en 35
W/m2.K. Estime la temperatura de las dos superficies de la placa de la plancha.
76
Figura ejemplo 1.7
Solución: La temperatura de la superficie externa se obtiene sabiendo que toda la energía convertida
en energía térmica se disipa por convección y radiación desde la superficie externa a temperatura T2:
Q   hconv  hrad  S T2  Talr 
La temperatura de la otra superficie de la base se halla usando la ecuación (1.10) por conducción:
kS T1  T2 
Qtotal 
L
Resultados
Q_conv=457,8 W; Q_rad=542,2 W; Pot = 1000 [W]; Talr =293 K; T2 = 947 K = 674 ºC
T1 = 690.7 °C.
Código en EES
$UnitSystem SI Mass Pa Deg J K
"Problema 1-101 Cengel p 55"
"Datos"
Pot=0,8*1250 [W]
h=35 [W/m^2-K]
epsilon=0,6
A=0,02 [m^2]
sigma=5,67e-8 [W/m^2-K^4]
T_alr=20+273 [K]
"La energía se disipa por convección y radiación"
Pot=Q_conv+Q_rad
Q_conv=h*A*(T_iron-T_alr)
Q_rad=epsilon*A*sigma*(T_iron^4-T_alr^4)
T_s=T_iron-273
INTERCAMBIO DE CALOR RADIANTE ENTRE DOS SUPERFICIES NEGRAS
EL FACTOR DE VISIÓN
El flujo de calor radiante emitido por la superficie “1” e interceptado por la superficie “2” es
Q12  A1 F12T14
77
El flujo de calor radiante emitido por la superficie “2” e interceptado por la superficie “1” es
Q21  A2 F21T24
La cantidad Fij es la fracción de la energía radiante que sale de la superficie i y es interceptada por la
superficie j y se conoce como factor de visión.
El intercambio neto de calor por radiación entre las dos superficies será
  
Q  A1 F12 T14  T24  A2 F21 T14  T24  (A-54)
A1F12 = A2F21 pues depende solo de la geometría del sistema y es válido tanto para superficies negras
como para superficies grises.
El factor de visión presenta un “álgebra” que permite en muchas ocasiones determinarlo fácilmente.
Algunas de sus reglas son
i) A1F12 = A2F21 Reciprocidad

ii) En recinto cerrado ∑ Fij = 1
iii) ∑ A1 F1j = A1
iv) Para superficie plana o convexa Fii = 0 (no puede verse a sí misma).
Figura 1.7 Factor de visión
Una ilustración sencilla es determinar los factores de visión para un sistema compuesto por una
superficie hemisférica y su base. Llamando esta última la superficie “1” y el hemisferio la superficie
“2”, tendremos:
F11 = 0 porque es plana  F12 = 1 por (ii). Ahora como

A1F12 = A2F21  F21 = (1)(A1/A2) = (R2)/(2R2) = ½ y F22 = ½ nuevamente por (ii).
Aquí determinamos 22 = 4 factores de visión. En general en un recinto con N zonas se deben

determinar N2 factores de visión. Por las reglas de reciprocidad y simetría se reducen los cálculos a
N(N-1)/2.
El flujo de calor neto entre dos superficies que emiten y absorben de forma perfecta se puede calcular,
según la ecuación (A-54) como:
78
Q12  A1 F12  Eb1  Eb 2  

 Eb1  Eb 2 
1
A1 F12
Recordar que A1 F12  A2 F21 (ley de la reciprocidad) y que Eb   T .

4
El intercambio de calor entre superficies negras donde     1 es sencillo. Para superficies reales
puede ser complejo. El desarrollo que aquí se presenta se limitará a superficies grises difusas y con
temperatura uniforme.
Superficies grises    
La transferencia de calor radiante neta desde una superficie gris se determina por la diferencia entre la
radiación que deja la superficie (emitida mas reflejada) y la radiación G que incide sobre ella. Si
llamamos J (radiosidad) el calor radiante por unidad de área que deja la superficie este será
J   G   Eb  1    G   Eb (A-57)
Intercambio de calor entre dos superficies grises
El flujo neto de energía que sale de la superficie es la diferencia de la radiosidad y la irradiación:
Q   J  G A
Despejando de (A-37) se obtiene:
J   Eb
G
1 
A  Eb  J   Eb  J 
Q 
1  1  
 
 A  (A-57a)
79
Entonces el intercambio neto de calor entre dos superficies grises puede darse en función de sus
radiosidades y de la visión relativa que tengan entre sí:
J1  J 2
Q  A1 F12  J1  J 2   A2 F21  J1  J 2  
1
A1 F12 (A-58)
Comparando las Ecuaciones (A-57a) y (A-58) se desprende que la red de radiación puede construirse
conectando una “resistencia superficial” a cada una de las superficies y una “resistencia espacial”
entre ambos potenciales de radiosidad como se ilustra en el esquema de la analogía eléctrica. En
consecuencia, la Ecuación (A-58) también puede escribirse en términos de las potencias emisivas de
cuerpo negro como:
Eb1  Eb 2  Ts4  Talrr

4

Q12  
1  1 1 1   2 1  1 1 1 2
   
A11 A1 F12 A2 2 A11 A1 F12 A2 2 (A-59)
En el caso límite de una superficie A1 localizada dentro de una envolvente muy grande tal que F12 = 1

y A1 / A2 0, la expresión se simplifica a:
Qrad  AS (TS4  Talr4 ) (A-60)
Otras situaciones en las que F12 = 1 y de simetrías simples que pueden analizarse fácilmente con la
ecuación (A-59) son dos placas planas paralelas infinitas, dos cilidros concéntricos y dos
esferasconcéntricas.
Superficies rerradiantes
En la práctica se puede encontrar que las superficies que intercambian calor radiante estén unidas por
superficies aisladas o adiabáticas que reemiten toda la energía que absorben en forma difusa, creando
un flujo de calor adicional que puede expresarse modificando el factor de visión. Para el caso de una
superficie refractaria R conectando las superficies radiantes “1” y “2”, el factor de visión se
modificará así
 FR 2 
F12  F12  F1R  
 FR1  FR 2  (A-60)
F12 es la suma de la fracción radiada directamente más la fracción rerradiada y se conoce como factor
de intercambio. Teniendo en cuenta que la superficie adiabática puede verse a sí misma y ver las otras
dos superficies, el álgebra del factor de visión en este caso será
80
A1F1R = ARFR1 ; ARFR2 = A2F2R. Multiplicando la ecuación para F12 por A1 y reorganizando
1 1 1
A1 F12  A1 F12   A1 F12   A1 F12 
FR1 FR 2 1 1 1 1
  
A1 F1R FR 2 A1 F1R FR 2 AR FR 2 A1 F1R A2 F2 R A1 F1R
Así el intercambio neto de calor radiante entre dos superficies en presencia de superficies adiabáticas
será
  
Q  A1 F12 T14  T24  A2 F21 T14  T24  (A-61)
Se demuestra que entre dos superficies grises entre las cuales hay flujo neto de calor, conectadas por
cualquier número de superficies rerradiantes, el flujo neto de calor radiante puede representarse por la
ecuación

Q  A1  12 T14  T24  (A-62)
 1
1 1   A  1 
   1   1    1
Donde 12 =
F12   1   A2    2  (A-62a)
Esta ecuación se aplica muy fácilmente a un número de sistemas simples pero útiles, pues cuando no
existen superficies reradiantes F12 = F12 y las ecuaciones (A-62) se reducen a la ecuación (A-59).
Por ejemplo, cuando la energía radiante se intercambia entre dos placas planas paralelas,
F12  1
1 1
 1
= F12 = 1 y A1 = A2, por lo que 12 =
1  2 (A-60b)
Radiación desde gases

Mientras que líquidos y sólidos emiten radiación sobre un espectro continuo, los gases no lo hacen.
De hecho, la mayoría de los gases monoatómicos y diatómicos, tales como los que están presentes en
la atmósfera, escasamente emiten radiación. Otros, como el vapor de agua y el dióxido de carbono
emiten radiación solo en bandas espectrales específicas. Solamente los gases con momento dipolar,
triatómicos y de más átomos emiten radiación en cantidades apreciables. Los gases que emiten
energía también la absorben, pero en las mismas bandas de radiación en las que emiten. Los gases que
no emiten radiación tampoco la absorben.
El ancho de las bandas de radiación para emitir (absorber) un gas en particular depende tanto de la
temperatura como de la presión. Una representación aproximada para vapor de agua nos da bandas
entre 2.2 a 3.3, 4.8 a 8.5, y 12 a 25 m. Para el dióxido de carbono las bandas son entre 2.4 a 3.0, 4.0
a 4.8 y 12.5 a 16.5 m.
81
El intercambio neto de calor radiante entre una masa de gas a temperatura TG y una superficie negra
que lo rodea A, a temperatura TS se puede determinar como

Q  A  GTG4   GTS4  (A-61)
La emisividad de la masa de gas G se obtiene de un diagrama de emisividad a la temperatura del gas

y la absortividad G se obtiene del mismo diagrama pero a la temperatura TS de la pared, pues esta
sería la temperatura de la superficie absorbente si se encontrara en equilibrio térmico con el emisor.
Estos diagramas son función de la temperatura, de la presión parcial del gas y de una trayectoria
equivalente que depende de la forma geométrica del sistema.
Figura 1.10a Emisividad del CO2 a 1 atm.
Las figuras 1.10 a 1.12 son tomados del trabajo de Hottel (Hottel, R. C.: Radiant Heat Transmision,
en Costa Novella Ed. Alhambra, 1986).
82
Figura 1.10b Corrección a la emisividad del CO2 para presiones diferentes a 1 atm.
Figura 1.11a Emisividad de vapor de H2O a1 atm
La fig. 1.10a da la emisividad co2 del dióxido de carbono en una masa de gas a una presión total de 1
atm, en función de la temperatura del gas Tg para varios productos PCO2L, en donde Pco2 es la presión
parcial del dióxido de carbono en atmósferas y L el radio, en m, de la masa hemisférica de gas. La fig.
1.10b da un factor de corrección aproximado Cco2 para ajustar los valores de la emisividad de la fig.
1.10a a presiones totales diferentes de una atmósfera; es decir, al multiplicar la emisividad del dióxido
de carbono co2 leída de la fig. 1.10a por el factor Cco2 del diagrama de corrección, fig. 1.10b, se ajusta
la emisividad a otras presiones totales.
La fig. 1.11a muestra la emisividad H2O del vapor de agua en una masa de gas a una presión total de 1
83
atm, en función de la temperatura del gas Tg para diferentes valores del producto PH2OL. La fig. 1.11b
da un factor de corrección aproximado CH2O para convertir la emisividad del vapor de agua a valores
de PH2O distintos a los cercanos a cero y a valores de presión total PT diferentes de 1 atm; es decir, la
emisividad del vapor de agua H2O de la fig. 1.11a se debe multiplicar por el factor de corrección CH2O
de la fig. 1.11b para ajustar la emisividad a valores PH2O y PT diferentes de 0 y 1 atm respectivamente.
Las absortividades del CO2 o del vapor de agua en las condiciones ya descritas se estiman de la
misma manera como se hizo con las emisividades, pero a la temperatura de la superficie y cambiando
los PGL por PGL(TS/TG). Los valores de emisividades así obtenidos se multiplican por (TG/TS)0.65 para el
CO2 o por (TG/TS)0.45 para el vapor de agua para obtener las correspondientes absortividades.
Figura 1.11b Corrección para la emisividad del vapor de agua
Figura 1.12 Corrección de la emisividad cuando se encuentran presentes vapor de agua y CO2.
Cuando en una masa de gas hay vapor de agua y dióxido de carbono simultáneamente, la emisividad
de la mezcla  m es algo menor, en una cantidad  , de la determinada al sumar separadamente las
emisividades  c y  w del dióxido de carbono y del vapor de agua calculados en las figs. 1.10 y 1.11
respectivamente. La razón es que se presenta una absorción mutua de la radiación. La emisividad  m
y la absortividad  m de la mezcla de dióxido de carbono y del vapor de agua se puede determinar
84
aproximadamente por la relación
 m    c   w 1  Cse  ;  m   c   w 1  Cse  (A-62)
en donde  c y  w son las emisividades del dióxido de carbono y del vapor de agua dadas en las figs.
1.10 y 1.11 y la corrección Cse de la emisividad por absorción mutua se puede obtener de la fig. 1.12.
En esta figura Pco2 y PH2O representan respectivamente las presiones parciales del dióxido de carbono y
del vapor de agua.
Efectos de la forma geométrica

Los diagramas de emisividad anteriores son para masas de gases que tienen una geometría
hemisférica y la longitud de trayectoria L. se refiere al radio del hemisferio. Para generalizar estos
diagramas a masas de gases que tengan otra forma se puede usar una longitud de trayectoria
equivalente L en el producto PL que aparece en los diagramas de emisión de las figuras 1.10 a 12. La
tabla 1-1 da las longitudes de trayectoria equivalente en el caso de varias formas geométricas simples.
TABLA 1-1
Geometría L
Esfera 0.60D
Cilindro infinito 0.90D
Placas paralelas grandes 2.0Distancia entre placas
Banco de tubos en triángulos equiláteros con paso D 2.8D
Banco de tubos en triángulos equiláteros con paso 2D 2D
Banco de tubos en cuadrados equiláteros con paso D D
EJEMPLO A.9: Determinar la emisividad de una masa gaseosa a Tg = 1000 °R y presión total PT =
1.5 atm que contiene 10% de vapor de agua, si la forma geométrica de la masa gaseosa tiene una
longitud de trayectoria equivalente L = 2 pies = 0.61 m.
Solución La presión parcial Pw del vapor de agua en la masa de gas es Pw = 1.5x0.1 = 0.15 atm. Luego
PWL = 0.15x0.61 = 0.091 m.atm. De la fig. 1.11a para Tg = 1000 °R = 555.6 K y PwL = 0.091 m.atm
se encuentra que la emisividad del vapor de agua a presión de 1 atm es g = 0.165. Este resultado debe
corregirse con respecto a la presión total y la presión parcial. De la figura 1.11b para (Pw + PT)/2 =
(0.15 + 1.5)/2 = 0.825 y PWL = 0.091 m.atm, se encuentra que el factor de corrección es Cw = 1.4.
Entonces la emisividad corregida de la masa de gas es  corregida = 0.165x1.4 = 0.23.
EJEMPLO A.10: Una chimenea de 3 pies de diámetro contiene un gas con 5% de CO2 a 2000 F y
una atm de presión total. El coeficiente convectivo entre el gas y la superficie refractaria es 1.5
Btu/hr.pie2.F. Si la superficie está a 1900 F y tiene emisividad 0.8, calcule el calor transferido por
radiación y convección entre el gas y la superficie.
La ecuación (i) tiene validez para el intercambio de calor entre un gas participante no luminoso y una
superficie negra que lo rodea. Si la superficie envolvente es gris de emisividad s, la emisividad
efectiva de la superficie envolvente estará comprendida entre la propia s y 1.0 en una forma más o
menos compleja. Sin embargo para este caso podemos hacer una aproximación que implica un error
85
aceptable en aras de la simplicidad suponiendo una emisividad efectiva igual al promedio aritmético
entre las dos mencionadas (Costa Novella N° 4 p 495). Por tanto, la ecuación (i) se modifica así:
s 1
Q A
2
  GTG4   GTS4 
(ii)
Para determinar el intercambio de calor radiante entre el gas con trazas de gas carbónico usamos el
gráfico 1.4a con PCL = 0.041 y Tg = 1366.7 K obteniendo C = 0.07. Como la presión total es 1 atm no
es necesario hacer correcciones.
Para la absortividad del gas a la temperatura de la superficie usamos PCL(TS/TG) = 0.0395 leyendo del
gráfico 1.4a  = 0.08. Se corrige multiplicando por (TG/TS)0.65 = 1.027 obteniendo G = 0.082. Por
tanto, reemplazando en la ecuación (ii) por unidad de área se obtiene
qradiación = 98.15 W (334.9 Btu/hr)

qconvección = 473.6 W (1616 Btu/hr)
qTotal = 571.75 W (1951 Btu/hr)
h=8,517 [W/m^2-C]
La diferencia de temperatura entre los gases y la superficie de la chimenea es pequeña por lo que el
aporte del intercambio radiante es del 17 % del total. La inexactitud en la lectura de los gráficos hace
que las correcciones por relaciones de temperatura para calcular la absortividad sean insignificantes.
EJEMPLO A.11: Una de las paredes laterales de una cava de refrigeración mide 5 [m] de largo por
2 [m] de alto y está construida de ladrillo aislante (k = 0.4 W/m.K) de 10 cm de espesor. El ambiente
interior es aire estancado a 10 °C bajo cero, mientras que sobre el exterior se estima un viento
horizontal paralelo de 8 m/s y 30 °C de temperatura, la presión atmosférica es de 101 kPa. Calcule la
temperatura de las superficies de la pared y el flujo de calor a través de la misma.
SOLUCIÓN: Se trata de de una pared plana con efectos convectivos en sus superficies. En una de
ellas se presenta convección natural y en la otra convección forzada. Como las propiedades físicas
deben evaluarse a temperatura de película es preciso asumir valores a las temperaturas superficiales
que se comprobarán iterativamente.
a) Ecuación para el flujo de calor:
Q
Ti  To 
1 / hi S  L / kS  1 / ho S
Luego de los cálculos se obtiene un valor de Q = -706.3 W. El signo indica flujo hacia el interior de la
cava.
b) Determinación del coeficiente convectivo interno hi . Se trata de convección natural sobre una
placa plana vertical. Es adecuada la ecuación (A.36)
86
 
2
0.1  Gr Pr  1012
0.387  Gr Pr 
1/6
 
Nu P  0.825  
9/16 8/ 27
 1   0.492 / Pr   
    (A.36)
Se usa el número de Rayleigh:
g  TS  T  L3car
Ra  Gr Pr 
 (A.9c)
   /  ,   k /  CP , para gases ideales   1/ T
La longitud característica en este caso es la altura de la pared. La temperatura intermedia entre el

interior y el exterior de la cava es 10 °C. El coeficiente por convección forzada es mayor que el
coeficiente por convección natural. Se asume Tso = 20 °C y Tsi = 10 °C. Así, la temperatura promedia
para evaluar las propiedades al interior de la cava serán 0 °C y al exterior 25 °C.
Al realizar los cálculos se obtiene

Tsi = 8.47 °C, Tmi = -0,8 °C, Ra = 2.15x1010, Nup = 324.6, hi = 3.826 W/m2.K
c) Determinación del coeficiente convectivo externo ho . Es convección forzada, flujo externo, fluido
pasando sobre una superficie plana. Se asume Reynolds crítico igual a 500000. El número de
Reynolds se calcula con base en la longitud en la dirección del flujo, 5 m.
V L
Re L 
 = 2.518x106 > 5x105
Como este valor es superior al crítico el flujo sobre la pared es mixto, parte laminar y parte
turbulento. Se usa la ecuación
Nu m 
hm L
k
   1
 0.664 Re C0.5  0.0365 Re 0L.8  Re C0.8 Pr 3
= 3540 (A.10g)
87
0.6  Pr  60; ReL  108.
Los resultados son
Tso = 26.12 °C, Tmo = 27.95 °C, Num = 3540, ho =18.22 W/m2/K.
El código en el software EES es el siguiente:
$UnitSystem SI J Pa K Rad
"Pared de cava"
"Flujo de calor"
Q=(Ti-To)/Suma
Suma=Resi+Resk+Reso "Resistencias en serie"
Resi=1/hi/S: Reso=1/ho/S: Resk=esp/k/S
Ti=263: To=303: S=5*2: k=0,4: esp=0,1 "Datos conocidos"
"Coeficiente convectivo interno"
Nup=(0,825+0,387*Ra^(1/6)/(1+(0,492/Pr_i)^(9/16))^(8/27))^2
Ra=g*beta*abs(Tsi-Ti)*Lci^3/nu_i/alpha "Rayleigh"
g=g#: beta=1/(Tmi): Tmi=(Ti+Tsi)/2: Lci=2 "Tmi se actualiza en cada iteración"
nu_i= mu_i/rho_i
mu_i=Viscosity(Air;T=Tmi)
rho_i=density(air;T=Tmi;P=P1)
P1=101325: Pr_i=0,71
alpha=k_i/rho_i/Cp_i
k_i=conductivity(air;T=Tmi)
Cp_i=1006: Nup=hi*Lci/k_i
"Coeficiente convectivo externo"
Re=rho_o*V*Lco/mu_o
Num=(0,664*Re_c^0,5+0,0365*(Re^0,8-Re_c^0,8))*Pr_o^(1/3)
Num=ho*Lco/k_o
mu_o=viscosity(air;T=Tmo)
rho_o=density(air;T=Tmo;P=P1)
Pr_o=0,71
Tmo=(To+Tso)/2: Re_c=5e5 "Tmo se actualiza en cada iteración"
V=8: Lco=5
k_o=conductivity(air;T=Tmo)
"Temperaturas superficiales"
Q=(Ti-Tsi)/Resi
Q=(Tso-To)/Reso
Tsic=Tsi-273
Tsoc=Tso-273
Resultados
a) Coeficientes y propiedades finales
alpha=0,00001806 [m^2/s] beta=0,003673 [1/K] Cp_i=1006 [J/kg-K]
Ra=2,215E+10 [-] Nup=324,6 [-] hi=3,825 [W/m^2-K]
Re=2,518E+06 [-] Num=3540 [-] ho=18,22 [W/m^2-K]
k_i=0,02357 [W/m-K] k_o=0,02573 [W/m-K]
mu_i=0,00001725 [kg/m-s] mu_o=0,00001862 [kg/m-s]
rho_i=1,297 [kg/m^3] rho_o=1,173 [kg/m^3]
b) Resistencias y flujo de calor
Resi=0,02614 [K/W] Resk=0,025 [K/W] Reso=0,00549 [K/W]
Suma=0,05663 [K/W] Q=-706,3 [W]
88
c) Temperaturas
Ti=263 [K] Tsi=281,5 [K] Tsic=8,47 [C] Tmi=272,2 [K] = -0,8 °C
To=303 [K] Tso=299.1 [K] Tsoc=26.12 [C] Tmo=301.1 [K] = 28 °C
Este ejercicio está parcialmente resuelto en los documentos siguientes usando Excel. Obtienen como
resultados aproximados temperaturas de 0 °C y 10 °C y flujo de calor 665 W.
Parte 1 https://www.youtube.com/watch?v=59rGhJuul8I
Parte 2 https://www.youtube.com/watch?v=KUjtvYtJ3iU
Parte 3 https://www.youtube.com/watch?v=v7lwDYC9k58
EJERCICIOS PROPUESTOS
1. Es una práctica corriente recuperar el calor residual de un horno a gas o aceite, usando los gases de
salida para precalentar el aire de combustión. Un dispositivo comúnmente usado para éste propósito
es un arreglo de tubos concéntricos en el cual los gases de salida pasan a través del tubo interior,
mientras que el aire frío fluye en paralelo a través del espacio anular. Consideremos la situación para
la cual existe una transferencia uniforme de calor por unidad de longitud, Q/L = 1.25x105 W/m, desde
los gases calientes hacia la superficie interior del anillo, mientras fluye aire a través del pasaje anular
a razón de m’ = 2.1 kg/s. El tubo interior de pared delgada (resistencia térmica despreciable) tiene
diámetro Di = 2 m, mientras que el tubo exterior, el cual se halla perfectamente aislado de los
alrededores, tiene diámetro Do = 2.05 m. Las propiedades del aire pueden tomarse como: CP = 1030
J/kg. K,  = 270x10-7 N.s/m2, k = 0.041 W/m.K, y Pr = 0.68.
(a) Si el aire entra a T1 = 300 K y L = 7 m, cuál es la temperatura de salida T2 ?
(b) Si el flujo de aire está completamente desarrollado en la región anular (coeficiente convectivo
constante), cuál será la temperatura del tubo interior a la entrada (TS1) y a la salida (TS2), y cuál la
temperatura del tubo exterior a la entrada (TSO1) ?
Nota: Recuerde que el coeficiente convectivo para conductos anulares se calcula con las mismas
ecuaciones que para los conductos circulares pero usando el diámetro equivalente De adecuado a
transferencia de calor tal como lo recomienda Kern.
2. Considere una delgada capa de agua reposando sobre el suelo en una noche clara con temperatura
ambiente 20 °C. La temperatura efectiva del firmamento es -50 oC, la emisividad de la superficie del
agua es 1.0 y el coeficiente convectivo de transferencia de calor entre aire y agua es 30 W/m2.K. (a)
Estime el coeficiente de transferencia de masa para la evaporación del agua. (b) Calcule la
temperatura del aire por debajo de la cual el agua se congelará. Asuma que no hay transferencia de
calor por conducción hacia el piso, y que el aire está seco. La presión de vapor del agua a 0 oC es 611
Pa y su calor latente de evaporación es 2.502x106 J/kg (Incropera-DeWitt Tabla A.6).
3. Un procedimiento común para aumentar el contenido de humedad del aire es burbujearlo a través
de una columna de agua. Suponga que las burbujas de aire son esferas de radio R = 1 mm y se
encuentran en equilibrio térmico con el agua a 25 C. ¿Cuanto tiempo deberán permanecer las
burbujas en el agua para alcanzar una concentración de vapor en el centro que sea 99% de la máxima
posible (saturación)?. El aire se encuentra seco cuando entra al agua. Suponiendo que las burbujas
alcanzan su velocidad terminal de forma instantánea, qué distancia habrán recorrido en este tiempo?
Algunos datos que pueden (o pueden no) serle útiles son: DAB = 0.26x10-4 m2/s, fase gaseosa, 298 K;
 aire = 1.1614 kg/m3;  agua = 1000 kg/m3;  agua = 855x10-6 Pa.s. La presión de vapor del agua puede
considerarse, para el intervalo de interés, dada por:
PA = a.exp [b/(T + 460)] en mm Hg si T en °F, con a = 1.1832x109 y b = -9524.86
89
4. Los gases de salida de un horno de procesado de alambre se descargan a una chimenea alta, y las
temperaturas del gas y de la superficie interna a la salida de la chimenea se deben estimar. El
conocimiento de la temperatura de la salida del gas es útil para predecir la dispersión de efluentes en
la corriente térmica, mientras que el conocimiento de la superficie interna de la chimenea a la salida
indica si puede ocurrir condensación de productos gaseosos. La chimenea cilíndrica tiene 0.5 m de
diámetro interno y 6 m de altura. Su pared esta construida por ladrillo refractario de 10 cm de espesor
cuya conductividad térmica es 1.21 W/m.K. El caudal de gases es 0.5 kg/s, y su temperatura de
entrada es 600 °C. La temperatura ambiente y la velocidad del viento son 4 °C y 5 m/s,
respectivamente. Tomando las propiedades termofísicas del gas como las del aire, estime la
temperatura de salida del gas y de la superficie interna de la chimenea para las condiciones dadas.
5. Estime la altura necesaria para que aire puro a 2 atm y 25°C, con flujo másico 710 kg/h, resulte
70% saturado en los equipos siguientes: a) Una columna de pared húmeda constituida por una
película de agua que desciende por la pared interior de un tubo vertical de 15 cm de diámetro interno.
b) Un lecho empacado con esferas de 1 cm de diámetro saturadas de agua. En este caso la porosidad
del lecho puede tomarse como  = 0.44
6. Considere un cilindro vertical de 3 plg de diámetro y 5 pie de longitud, componente de un aparato

de humidificación, por cuya superficie exterior desciende una delgada película de agua. En ángulo
recto con el cilindro circula aire seco a 110 F y 600 mm Hg con velocidad de 22 pie/s. El agua no
debe gotear por abajo del cilindro.
(a) Determine la temperatura del agua.
(b) Determine el caudal al cual debe suministrarse el agua en la parte superior del tubo. Para agilizar,
tome como constantes las siguientes propiedades del aire: viscosidad dinámica  = 0.01876 cP =
0.0454 lb/pie.h; DAB = 1.021 pie2/h. ¿A qué presión y temperatura deberán calcularse estas
propiedades?
7. Una placa porosa colocada paralelamente a una corriente de aire puro que fluye con velocidad de
30 pie/s, está inmersa por su parte inferior en un recipiente con tetracloruro de carbono, cuya
temperatura se mantiene igual a la de la interfase del plato y el aire. Esta interfase permanece
humedecida con CCl4 gracias a la acción capilar a través de la placa. Si la presión del aire es 460 mm
Hg y la temperatura de la corriente de aire es 262 °F, estime la temperatura de la superficie, la
velocidad de evaporación en un punto ubicado a 6 pulgadas del borde de ataque y la velocidad de
evaporación total si la placa mide 5x5 pie y la transición a flujo turbulento se estima para ReC =
3.5x105.
Use las correlaciones adecuadas para estimar la viscosidad del vapor de CCl4 y para estimar la
viscosidad de la mezcla aire vapor, la conductividad térmica, el calor específico, así como la
difusividad del sistema a las condiciones de la interfase y la corriente principal de aire para poder
obtener valores promedio adecuados para la capa límite. No asuma baja velocidad de transferencia.
8. A través de un lecho empacado con esferas porosas de 5 mm de diámetro humedecidas con agua
fluye aire seco a una atm. y velocidad superficial de 1.0 pie/s. Suponga que la evaporación ocurre en
la superficie de las esferas húmedas a temperatura uniforme y constante que es esencialmente la de
bulbo húmedo. Suponga flujo pistón y estime la altura de empaque requerido para que el aire llegue a
estar 90% saturado con agua (Nos referimos a la temperatura del liquido). Datos: Puede usarse  y µ
del aire seco a 20 °C, con DAB = cm2/s. La fracción vacía en el lecho de esferas es  = 0.41.
9. Por un conducto vertical de tres cm de diámetro interno asciende glicerina con caudal másico de
100 kg/hr, cuya temperatura inicial es de 308 K. Calcular la temperatura de la glicerina en función de
la distancia recorrida en el tramo inicial de la conducción si la temperatura de la superficie se
mantiene a temperatura constante de 373 K. Las propiedades físicas de la glicerina son: CP = 2.5
90
kJ/kg.K; k = 0.286 W/m.K;  = 0.54x10-3 K-1;  = 1217 kg/m3. La viscosidad varía según la
siguiente tabla:
 (kg/m.s) 10.715 2.083 0.893 0.270 0.150

T (K) 272 288 300 308 311
SOLUCIÓN DE LA ECUACIÓN DIFERENCIAL ORDINARIA LINEAL, DE SEGUNDO

ORDEN, NO HOMOGÉNEA:
d2y
 Ay  B
dx 2 A, B son constantes (i)
Tiene como solución la suma de la solución homogénea y la solución particular.
La solución de la homogénea, dependiendo de si el signo es + o – da una solución en términos de

funciones trigonométricas o de funciones hiperbólicas:
yh = C1sin(mx) + C2cos(mx) para el signo positivo.
yh = C1sinh(mx) + C2cosh(mx) para el signo negativo.
m2 = A
La solución particular será, en este caso, una constante, que se determina haciéndola satisfacer la
ecuación (i), es decir si yp = K, entonces K = B/A.
Las constantes de integración C1 e C2 se determinan aplicando las condiciones límite del problema a
la solución general y = yh + yp.
Un caso más general es el siguiente: Ecuación diferencial ordinaria lineal de segundo orden, con
coeficientes constantes, homogénea:
d2y dy
a 2
 b  cy  0
dx dx (ii)
Tiene como ecuación característica
aD 2  bD  c  0
Esta ecuación tiene raíces D1 y D2 . Las soluciones a (ii) tienen la siguiente forma:
D1 , D2 Solución
1. Distintas y reales y  C1 exp  D1 x   C2 exp  D2 x  = C1 sinh  D1 x   C2 cosh  D2 x 
2. Iguales y reales y  C1 exp  Dx   C2 x exp  Dx 
3. Imaginarias D1,2 = bi y  C1 cos  bx   C2 sin  bx 
4. Complejas conjugadas y  exp  ax   C1 cos  bx   C2 sin  bx  
D1,2  a  bi
91
Correlaciones para predecir las propiedades de transporte se pueden estudiar en el capítulo 3.
TRANSFERENCIA DE CALOR O DE MASA SUPERPUESTA A UN CAMPO DE FLUJO
Determinación de coeficientes convectivos en flujo laminar y geometrías simples
INTRODUCCIÓN
Para analizar un problema del mundo real de manera científica se debe modelar matemáticamente a
través de conceptos básicos. Estos conceptos son la conservación de las especies químicas, la
conservación de la materia, la conservación de la cantidad de movimiento y la conservación de la
energía. Para cada una de estas entidades que se conservan se puede escribir un balance de flujos que
describa la transformación de las mismas. Este balance se expresa en forma matemática, tanto a nivel
macroscópico como microscópico. La desigualdad entrópica también es un concepto básico, pero
esencialmente indica si un proceso es factible o no y no es base de una ecuación de balance.
Para simular un fenómeno físico tal como flujo de fluidos, transferencia de calor, transferencia de
masa, los principios de conservación se expresan en términos de ecuaciones diferenciales parciales.
Estas ecuaciones se conocen como ecuaciones de conservación.
Estas ecuaciones se desarrollan en términos de “concentraciones” es decir, la entidad conservada por

unidad de volumen. Por ejemplo, la ecuación de cantidad de movimiento que expresa la conservación
de la cantidad de movimiento lineal en términos de la cantidad de movimiento por unidad de volumen
es decir el producto de densidad por velocidad V . La ecuación para la conservación de la energía
expresa la conservación de la energía por unidad de volumen, producto de densidad por capacidad
calorífica por temperatura  CPT . La conservación de las especies químicas será indicada por
concentraciones másicas, por ejemplo, densidad por fracción másica  wA   A , o concentraciones

molares c A  cx A . Esta concentración de cualquiera de las tres entidades la representaremos por .
Usaremos un volumen de control de aristas xyz como se muestra en la figura.
El principio de conservación establece que la acumulación de  con el tiempo en el volumen de

control es igual a la entrada neta de  al volumen de control más la generación neta de  dentro
del volumen de control.
La entrada neta de  al volumen de control se establece por dos mecanismos primarios: difusión,
debida a movimientos moleculares y convección debida al movimiento del fluido.
92
La densidad de flujo difusivo x se escribe

 dif x   
x (1)
 2
Aquí es una propiedad del sistema, una difusividad, con dimensiones de longitud /tiempo .
La densidad de flujo convectivo x se escribe
 convx  Vx (2)
La densidad de flujo total en x es:
 x   convx   dif x (3)
El flujo neto de entrada en la dirección x es   x   x x  y z  . Teniendo presente que expresiones

y z
similares se escriben para las direcciones y , ordenando términos y dividiendo por xy z t ,
el balance es:
 t t   t   x   x x    y   y y    z   z z 

   
t x y z .
Tomando el límite cuando x, y, z , t  0 se obtiene
  x  y  z
   
t x y z (4)
Escrito de manera abreviada (notación tensorial vectorial) este balance generalizado toma la forma:
     / t   .
93
La aplicación del operador  (operador nabla) a un vector se llama la divergencia de ese vector, o
más simplemente el producto punto y el resultado es un escalar. Recordamos que el gradiente (un
vector por un escalar) da un vector. El operador  tiene la ventaja adicional de que puede ser
expresado en coordenadas curvilíneas. Esto es de gran ayuda cuando las ecuaciones de balance se
expresen en sistemas coordenados alternos. Esto es posible ya que una ecuación vectorial escrita en
notación vectorial se aplica a todos los sistemas coordenados. O sea que una ecuación escrita en
notación vectorial puede transformarse a cualquier sistema coordenado.
Reemplazando la densidad de flujo como la suma de términos moleculares y convectivos, para el

problema tridimensional, la ecuación vectorial correspondiente es la divergencia de la densidad de
flujo, así:
       dif     conv            V  .
Incluyendo esta expresión en el balance generalizado, obtenemos:

   V         
t .
Balance generalizado para fluido incompresible y propiedades de transporte constantes
En este caso el balance anterior, se simplifica así:
constante e isotrópica, nos aparece al lado derecho el término       . El


Para difusividad
producto punto      operando sobre un escalar ocurre tan frecuentemente que se le ha dado un
símbolo especial  , y un nombre especial, el de operador Laplaciano.
2
De otra parte, el término que incluye la velocidad se analiza de la siguiente manera:
   V       V   V    .
Un balance global de materia, donde no habrá término de generación (la masa ni se crea ni se
destruye a no ser en reacciones nucleares), ni gradientes, por lo que  es la densidad del sistema, 
, masa total por unidad de volumen,
 
   V         V   V      0
t t ,
Esta es denominada la ecuación de continuidad. Claramente si la densidad  es constante, se

concluye que:
 V   0 .
Teniendo los análisis anteriores en cuenta, el balance toma la forma:
94

V         2   
t .
En sistemas de conducción en sólidos, y sistemas difusivos, donde el primer término (o de arrastre) se

pueda despreciar o sea efectivamente nulo, el modelo matemático se simplifica aún más:

    2   
t
La Derivada Substancial: La concentración  (de masa, energía calorífica o cantidad de

movimiento) puede, en un caso general, ser función de la posición y/o del tiempo, lo que en
coordenadas rectangulares nos lleva a:
d  dx  dy  dz 
      
dt x dt y dt z dt t
El significado físico de esta derivada puede entenderse si imaginamos un instrumento que lea  y
pueda moverse en el fluido. Las cantidades dx/dt, dy/dt, dz/dt son las componentes de velocidad del
elemento detector y d/dt será una derivada total. Si el detector está fijo en el espacio entonces d/dt
= /t. Si el elemento detector sigue el flujo, sus componentes de velocidad serán las mismas del
fluido y tenemos:
D    
 Vx  Vy  Vz 
Dt x y z t
Esta nueva derivada se denomina la “derivada substancial o material”, o en ocasiones la derivada de

Lagrange. Representa la velocidad de cambio de  para una masa constante de fluido que se mueve;
el volumen de esta masa de fluido puede cambiar con el tiempo mientras se mueve. Se observa que
D/Dt = (V) + /t.
Transferencia de Masa en una Película Liquida Descendente.
La situación a analizarse consiste en una película líquida de propiedades constantes que desciende
con flujo estacionario, por una superficie sólida plana y lisa. La superficie libre de la película que
desciende es adyacente a una fase gaseosa. Puede ocurrir transferencia de masa, tanto entre la
superficie sólida y la película como entre la película y la fase gaseosa. Cada posibilidad se examinará
a su turno, para condiciones de baja concentración de soluto y bajas velocidades de transferencia de
masa, en forma tal que la velocidad de arrastre perpendicular a la placa, debida a la difusión, pueda
despreciarse.
Transferencia de masa entre una fase gaseosa y una película liquida descendente.
Esta es una situación que se presenta por ejemplo en absorción gaseosa en una columna de pared
húmeda, para un sistema en el cual la resistencia controlante se encuentra en la fase líquida. Para éste
caso se puede asumir que la concentración de soluto en la superficie es constante en As, la solubilidad
del gas en el líquido (establecida por ejemplo por la ley de Henry). Se supone además que las
concentraciones de A son bajas de tal forma que las velocidades perpendiculares a la pared son
95
efectivamente cero, y que la difusión en la dirección x es despreciable comparada al transporte por

movimiento convectivo.
En tres dimensiones, el balance para la especie A es (ver ecuación (4)):
n Ax n Ay n Az  A
    A
x y z t (2.28)
nAx  j Ax   AVx nAy  j Ay   AVy (ver ecuación (3))
Figura 2.22
n AZ  0 (no hay gradientes ni flujo en la dirección z)
 A
0
t (estado estable)
 A  0 (no hay reacción química homogénea)
 AVy  0 ; Componente de velocidad y despreciable.
J Ax   AVx ; El transporte convectivo de A en la dirección x es mucho mayor que el transporte

difusional en tal dirección (Pex >> 1).
Reemplazando en (2.28) obtenemos:

 A 2 A
Vx  DA B 0
x y 2 (2.60)
96
Esta ecuación diferencial parcial no puede resolverse sin resolver la ecuación de cantidad de
movimiento para encontrar el perfil de velocidades Vx que es función de y. Es conveniente ubicar el
origen del sistema coordenado en la superficie de la película descendente como se indica en la figura
2.22.
Para determinar el perfil de velocidades se realiza se realiza un balance de cantidad de movimiento x

en la dirección y, en un elemento de volumen de espesor dy:
Salida menos entrada de cantidad de movimiento difusivo en y y convectivo en x es igual a la suma de

fuerzas de presión y gravitacionales que actúan sobre el elemento de volumen:
( yx S y   yx S y )  ( VxVx xL
 VxVx x 0
)(Wdy )  ( P0  PL )(Wdy )   g y S y dy
y  dy y
Como la velocidad en la dirección x no cambia sobre la misma línea de corriente (fluido

incompresible), el segundo término del lado izquierdo, correspondiente al cambio de cantidad de
movimiento convectivo en la dirección x es cero. Dividimos por S y dy :
d P  PL
 yx   0   gx
dy L (2.23)
Al observar el paralelismo entre las ecuaciones (1.20) o (4) y (2.23), nos lleva forzosamente a
concluir que el término de generación de cantidad de movimiento es la suma de fuerzas, es decir:
P0  PL
M    gx
L (2.23a)
que tiene dimensiones de cantidad de movimiento por unidad de tiempo y unidad de volumen (o
fuerza por unidad de volumen). Las condiciones límite para (2.23) son: y = 0, yx = 0; y = , Vx = 0
por ley de no resbalamiento. La diferencia de presiones, P0 - PL se debe a la columna hidrostática de la
capa de aire exterior que por la diferencia de densidades entre gas y líquido se puede despreciar.
Reemplazando la ley de Newton de la viscosidad obtenemos:
d 2Vx
   gx  0
dy (i)
De acuerdo a la disposición de los ejes, g x , la componente de la gravedad en la dirección de flujo

g x  g sin   , siendo  el ángulo que la placa hace con la horizontal.
La misma expresión (i) se obtiene partiendo del balance generalizado de cantidad de movimiento para
un fluido incompresible escrita en coordenadas cartesianas:
DV  V 
   (V  )V    2V  P   g
Dt  t  (2.45)
Esta expresión se conoce como la ecuación de Navier – Stokes y al expandirla resultan 27 términos.
Suele interpretarse como una forma de la segunda ley de Newton del movimiento: masa por
97
aceleración = Suma de fuerzas (viscosas, de presión y gravitatorias) que actúan sobre el sistema. Todo
por unidad de volumen. La componente x se escribe:
 V V V V    2V  2V  2V  P
 Vx x  Vy x  Vz x  x     2x  2x  2x     gx
 x y z t   x y z  x
Para un sistema incompresible en estado estable el primer término se hace cero. Además, en nuestro
caso, Vx es función solamente de y.
Resolvemos (i):
CL 1: y = 0 ; dVx/dy = 0 CL 2: y =  ; Vx = 0
Integrando (i) una vez:
dVx g y
  x  C1
dy 
CL 1 : C1 = 0
Integrando nuevamente
 g x 2
Vx    C2
2
 g x 2
 C2
CL 2 : 2
Por lo tanto:
  y  
2
  gx  2   g x 2
Vx      y     1   
2

 2   2       (2.61)
98
Se observa que la velocidad máxima viene dada por:
 g x 2
Vmax 
2
El caudal líquido en la película será:


W  g x 3
Q '  W  Vx dy 
0
3 donde W es el ancho de la placa en dirección z.
O sea, que el espesor de la película es:
1/3
 3 Q ' 
  
W  g x 
La velocidad promedio en la película, caudal volumétrico sobre área de flujo:
Q '  g x 2 2
Vm    Vmax
W 3 3
Definimos un número de Re para la película como:
 Vm Deq 4 AF 4 W
Re f  Deq  4 RH    4
 ; PH W
es decir:
4  Q ' 4
Re f  
W 
Donde  es el caudal másico por unidad de ancho de la película. El número de Reynolds es un

indicador de las condiciones de flujo. Para Ref  30 la película es laminar y libre de ondas. Hasta Ref
= 1800 se forman rizos u ondas sobre la película, y a Ref  1800 la transición de flujo laminar a
turbulento es completa.
Reemplazando (2.61) en (2.60), la ecuación diferencial que describe este proceso es:
99
  y  2   2  A
Vmax 1     A  DAB
     x y 2 (ii)
Con las siguientes condiciones límite:
x0 0 yb  A   AO concentración uniforme a la entrada
para todo x, y = 0,  A   AS la solubilidad, en la interfase existe equilibrio
 A
0
y = , todo x, y la placa es impermeable al soluto A.
Johnston y Pigford (Trans. A.I.Ch.E. 38, 25-51, 1942 "Distillation in a wetted wall column")
encuentran la solución en la serie siguiente. Aquí c AmL es la concentración de mezcla en la posición x
= L. Observe que esta ecuación no depende de y sino exclusivamente de x:
c As  c AmL
 0.7857 exp  5.1213Fo   0.1001exp  39.318 Fo   0.03599 exp  105.64 Fo   ⋯
c As  c Ao
El número de Fourier que es un tiempo de exposición adimensional está dado por:
DAB L
Fo 
 2Vmax
La fuerza guía para la absorción cambia desde  c As  c Ao  en la parte superior de la película líquida,
hasta  c As  c AmL  en la parte inferior de la misma y tiene un valor cambiante  c As  c Amx  en cualquier
punto de la película. Se hace un balance diferencial en una sección de altura dx a lo largo de la cual la
composición de la película cambia un dc Amx . El balance es:
Vm W   dc Amx  kc  c As  c Amx Wdx 
Separando variables e integrando se obtiene:
Vm  c As  c Ao 
kc  ln
L  cAs  cAmL  (iii)
Para tiempos donde Fo > 0.1, despreciando todos los términos de la serie menos el primero, con 
como longitud característica se obtiene el coeficiente adimensional Sherwood:
Vm 2  exp  5.1213  Vm

2
Sh  ln    0.241  5.213Fo 
DAB L  0.7857  = DAB L (iv)
Si Fo > 1.0, despreciando el primer término del paréntesis y reemplazando Fo y Vm por su valor se
obtiene Sh = 3.414. Estas correlaciones se recomiendan para Ref < 100.
100
Estos coeficientes se usan para calcular la absorción en toda la película usando una diferencia de
concentraciones promedio adecuada  c As  c Amx  M :
Vm W   c AmL  c Ao   kc  c As  c Amx  M WL  (v)
Reemplazando (iii) en (v) se obtiene la diferencia de concentraciones promedia a usar:
 cAs  cAo    cAs  cAmL 

 cAs  cAmx M 
 c  c Ao 
ln  As 
 c As  c AmL 
Este es el promedio logarítmico de las fuerzas guías en los dos extremos de la película líquida.
Tiempos cortos de exposición.
Cuando el tiempo de exposición texp = L/Vmax es corto, la sustancia A penetrará poco en la película, y
podemos asumir perfil de velocidad uniforme e igual a Vmax. Además, como no alcanza a enterarse de
la existencia de la pared, la película será, desde el punto de vista del soluto A, de espesor infinito. Por
tanto:
Figura 2.22a
 A 2  A
Vmáx  DAB
x y 2 (2.61)
 g  2  3
Vx  Vmáx   x   Vmed  Constante
 2  2
Condiciones límite:
x  0 A = Ao Condición de entrada
y=0 A = As Superficie libre.
y A = Ao Poca penetración.
101
Haciendo un cambio de variable, la ecuación (2.61) se convierte en una ecuación diferencial

ordinaria:
y
n
Asumimos A = f(n) con 2 DAB x / Vmax
 A f dn  A f dn 2 A 2 f  dn 
2
   2  
x n dx y n dy x 2 n  dx 
entonces ; ;
dn n dn n
 
dy y ; dx 2x
Reemplazando en (2.61)
2
n  n 
D AB   f " n   vmax   f ' n   0
 y  2x 
y 2 vmax
 2n 2
Reconociendo que 2 xDAB obtenemos
f " n   2nf ' n   0
Con las condiciones límite:
CL1: para y = 0 ; n= 0 ; A = As
CL2: para y   o para x = 0 ; n   , A = Ao
Substituyendo f' por  se llega a una ecuación de primer orden de variables separables que se puede
resolver para dar:
d

 2ndn  ln  n 2  ln C1   
df
dn
 C1 exp  n 2  
n
 
f n    A  C1  exp  n 2 dn  C2
0 (2.62)
Aquí se elige arbitrariamente el límite inferior de la integral indefinida, que no puede resolverse en
forma no acotada. Si se cambiase el límite inferior por otro valor se alteraría simplemente el valor de
las constantes de integración, no determinadas aún. Aplicando las condiciones límite se obtiene:
CL 1 : C2 = AS
102
 A0   AS
C1   
 exp n dn
2
CL 2: 0
Haciendo u = n2 ; dn = (1/2) u-½ du
 

 exp  n dn 
2
 1  12 1
2 0
u exp u du
0
por definición de la función Gama:

    u  1 exp u du
0
entonces:


 exp n dn     
2 1 1
2 2
0
2
Por tanto:
n
 A   AS
exp n 2 dn
2
 A0   AS
 
 0 (2.63)
La relación:
x
 exp  u  du
2
x
2
  exp  u  du  erf  x 
0 2


 exp  u  du
2 0
Se denomina función error, y se abrevia "erf" (ver Anexo D).
103
La función complementaria de error es: [1  erf(n)] = erfc(n)
Podemos expresar la solución como:
 A   A0  y  x
 erfc   texp 
 As   A0 Vmax
 2 DAB texp 
(2.64)
Método de Transformada de Laplace
La ecuación (2.61) se puede resolver también usando el método de la transformada de Laplace:
 A  2A t  x / Vmax
 DAB
 t  y2 (2.61)
Definiendo Y   A   Ao , con  Ao constante obtenemos:
 Y  2Y
 DAB
 t  y2 , (1.73)
con las siguientes condiciones límite:

x= 0, t  0, Y  0   A0   A0 , 0  y   comienza la película
0 < x < L, t  0 a lo largo de la película
Y  YS  (  As   Ao ), y  0 superficie de la película
Y  0    Ao   Ao  , y poca penetración
Recordemos las transformadas estándar implicadas:
  f  y, t  
Lt    p f  y, p   f  y,0 
 t 
  n f  y, t    n f  y, p 
Lt  
  y  yn
n
 ,
Usándolas la transformada de (1.73) se escribe como:
 dY   d 2Y 
Lt    Lt  DAB 
dy 2  .
 dt  
Reemplazando y organizando:
104
d 2Y pY
2
 0
dy DAB
Esta es la ecuación de la aleta. Las condiciones límite con una condición en el infinito sugieren la
solución:
 
1
Y  C1 exp   qy   C2 exp  qy  ; q  p
2
DAB
Las condiciones límite también se transforman:
YS
Y  , p   0; Y  0, p  
p .
Aplicando la primera condición límite C2  0 . Aplicando la otra condición límite: C1  YS / p . La

solución será:
Y exp   qy 

YS p
Aplicando la transformada inversa, numeral 30 de la tabla de trasformadas de Arpaci obtenemos:
 A   Ao  y   y 
 erfc    1  erf  
 As   Ao  2 DAB t   2 DAB t  . (1.72a)
El método de cambio de variable es adecuado pero el de transformada de Laplace es más directo.
La función de error es una función monótona creciente cuyo intervalo va de cero a uno y alcanza el
valor de 0.99 cuando el argumento vale 2. Podemos utilizar este hecho para definir una profundidad
de penetración del soluto, ypenetración, como la distancia y para la cual (A  Ao) ha disminuido hasta un
valor de (0.005)(As.-Ao), es decir, erfc(n) = 0.005  n = 2 y la profundidad penetrada por el soluto
será:
1
D x
2
ypenetración  4  AB 
 Vmáx  (2.65)
Esta distancia es una medida de la cantidad de masa de A que ha penetrado en la película de fluido.
Observe que es proporcional a la raíz cuadrada del tiempo de exposición.
Para valores típicos de tiempo de exposición de 2 segundos y de difusividades en líquidos del orden
de 105 cm2/s, el soluto, en la parte inferior de la placa habrá entrado hasta capas donde la velocidad
del fluido es todavía el 99% de la velocidad máxima.
105
La densidad de flujo de materia a través de la interfase en un punto x:
 A
n Ay   DAB
y 0 y y 0 ; partiendo de la ecuación (2.64) obtenemos:
 A n  2   As   A0    exp   n  
2
 A
   
y y 0
n y     2 DAB texp  y  0
1
D V 
2
nAy   AB máx    As   A0 
y 0
 x 
Idéntico resultado se obtiene a partir de la ecuación (2.63). Esta es una densidad de flujo local que
varía inversamente con x0.5, es decir cuando aumenta el valor penetrado que es una medida del camino
de difusión. Un valor promedio de la densidad de flujo lo encontramos integrando nAy a lo largo de x y
dividiendo por la longitud L de la placa:
    A0   DABVmáx 
1 1
L L
1  DABVmáx 
2 2
nAS   nAy dx  As    0 x dx  2   L 
 12
  As   A0 
L0 L  
O en términos de la velocidad promedio:

1
 6D V 
2
nAS   AB med    As   A0 
 L  (2.66)
La masa total de A, absorbida por la película será el producto de nAS por el área de la superficie.
Usando la definición del coeficiente convectivo, nAS = k(AS  A0), nos aparece definido un
coeficiente de transferencia de materia que usa como fuerza guía un valor constante e igual a la
máxima diferencia de concentraciones entre la superficie y el fluido y que además es proporcional a
DAB½. Esto es característico de la "Teoría de la Penetración" propuesta por Higbie en 1935.
1
 6D V 
2
k    AB med 
 L  (2.67)
Este coeficiente tiene subíndice  porque actúa sobre diferencias de concentración expresadas en
unidades de concentración másicas volumétricas y, como puede comprobarse sus dimensiones son las
de una velocidad (L/t). Es idéntico a kc.
En términos adimensionales corresponde a:
1/ 2
kc  6 2Vm  
Sh    0.691  Re1/f 2 Sc1/ 2
DAB  DAB L L
106
Esta ecuación se recomienda para 1200 < Ref < 100
Al escribir los resultados en términos de coeficientes adimensionales facilita extenderlos a

transferencia de calor (o al contrario) usando la analogía.
La ecuación (2.67) escrita como:
 4D 
1
2 L
texp 
k    AB  Vmax
  texp 
Se aplica para la absorción gaseosa hacia chorros líquidos y transferencia de masa desde burbujas
ascendentes si la penetración del soluto es pequeña.
La transferencia para películas líquidas descendiendo por el interior o exterior de tubos (columnas de
pared húmeda), de radio R, pueden tratarse con las mismas expresiones que la película líquida
descendente plana, cuando son hidrodinámicamente equivalentes, lo que ocurre para números de
Goucher mayores que tres (Tallmadge y Gutfinger en Skelland p 135):
0.5
 g 
Go  R   3
 2  (2.68)
Donde  es la tensión superficial (dina/cm = g/s2).
Los resultados anteriores para películas líquidas descendentes dan resultados bajos comparados con la
experimentación. Esto se origina en la presencia de ondulaciones en la superficie aún para Ref del
orden de 25. Estas ondulaciones aumentan la superficie ideal de la película e inducen algún grado de
turbulencia. Stirba y Hurt (A. I. Ch. E. J., 1, 178-184 1955) proponen el uso de una difusividad
aparente o efectiva basados en experimentos en películas verticales. Sus resultados pueden
correlacionarse como:
DABef  7.253  1012 Re 2f  8.590  108 Re f  8.919  105 [pie 2 / hr] 100 < Ref < 5000
D_ef=7,253e-12*Re_f^2+8,59e-8*Re_f+8,919e-5
Transferencia simultanea de calor y cantidad de movimiento
Transmisión de calor desde una pared a una película descendente: tiempos de contacto cortos.
Una película de un líquido frío que desciende por una pared sólida vertical o inclinada, tal como se
indica en la figura 2.23, ejerce un considerable efecto de refrigeración sobre la superficie del sólido.
Estimar la velocidad de transmisión de calor desde la pared al fluido para tiempos de contacto tan
cortos que la temperatura del fluido sólo varía apreciablemente en la región inmediata a la pared.
A partir de las ecuaciones (1) a (4) en estado estable, sin generación, y asumiendo que en la dirección

y prima la densidad de flujo difusiva dif y y que en la dirección z lo hace la densidad de flujo
convectiva  convz con concentración de energía    CPT , el balance de energía se reduce a:
107
 dif y  convz
 0
y z
     C PT   
   Vz  CPT   0
y y  z
Asumiendo propiedades constantes con   k /  CP y sabiendo que Vz es función de y:
T  2T
 CPVz k 2
z y
Esta ecuación requiere el balance de cantidad de movimiento para conocer el perfil de velocidades en
la película. Suponiendo flujo laminar y que el perfil corresponde al perfil isotérmico se obtiene:
 2 y  y 2   2 g
Vz  Vmax      Vmax 
2
      siendo
Figura 2.23
Como se asume baja penetración, cerca a la superficie de la pared el perfil de velocidades puede
asumirse lineal con pendiente igual a la del perfil parabólico en el origen:
 dV  y
 g
Vz   z y y
 dy  
 y 0  =
Combinando los resultados anteriores se obtiene la ecuación diferencial aproximada:
108
T  2T
y  2
z y (2.73)
Con las siguientes condiciones límite :
CL1: z=0 y0 T = T0

CL2: z>0 y=0 T = TS
CL3: Todo z y T = T0
Esta ecuación puede transformarse a una ecuación diferencial ordinaria introduciendo el cambio de
variable:
1/ 3
  
n  y 
 9 z  y suponiendo T = f(n).
Así:
T f n n 1n
 ; 
z n z z 3z
T f n  2T  2 f  n 
2
n n
 ;  
y n y y 2 n 2  y  y y
  ;
Reemplazando en (2.72) obtenemos:
2
 yn n
 f ' n      f ' ' n 
3z  y
Observando que
  
y3    3n 3
 3z 
y reorganizando:
f ' '3n 2 f '  0 (2.74)
donde las comillas indican derivadas con respecto a n.
Las tres condiciones límites originales se reducen a dos (la condición límite 1 y 3 se hacen iguales):
CL1: n   T  T0 z  0 ó y  
CL2: n  0 T  TS  y  0
109
d
f '   f ''
Haciendo dn la ecuación (2.74) queda:
d
 3n 2  0
dn
d ln   n 3  ln C1
 3n 2 dn
dn , entonces:
  C1 exp n 3   
dT
;  C1 exp  n 3
dn
T n
 dT  C1  exp  n dn
3
 
TS 0 (2.75)
Para evaluar la constante de integración, C1, usamos las condiciones límite así:
T0 
 dT  C1  exp  n dn
3
 
TS 0 (2.75a)


I   exp  n3 dn 
Llamemos 0
u  n3 ; n  u1 / 3 1
dn  u  2 / 3 du
Sea entonces y así 3
Si n  0, u  0 ; si n   , u  ; luego :

1
I   u  2 / 3 exp u  du
30
Por definición la función gamma (ver anexo E):

    u 1 exp u  du
0
En consecuencia:

1 1 1
I   u 1 / 31 exp u du   
30 3 3
110
Utilizando integración por partes se demuestra la siguiente propiedad de la función gamma:
 p  1  p p 
por esto:
1 1 1  4
I       1   
3 3 3  3

4
 
    exp  x 3 dx  0.8930
3 0
Por esto, reemplazando en (2.75)
T0  TS 
C1 
0.8930
de donde, (2.74) se transforma en:
T  TS n

1

T0  TS 0.8930 0

exp  n 3 dn 
(2.76)
La integral de esta expresión puede evaluarse mediante el uso de series u otro método numérico y
tiende rápidamente a su valor asintótico variando entre 0 y 0.89295, para n variando entre 0 y 2.
Podemos usar el perfil de temperaturas para determinar el coeficiente de transferencia de calor en

función de z:
dT
q S  hTS  T0   k
dy y 0
para z constante:
T T n

y n y
1/ 3
n   

y  9z 
T T0  TS

n 0.8930
exp  n 3  
Pues la derivada de una integral con respecto al límite superior es el integrando.
111
1
T  T  TS   
 
3
 0 exp  n 3   9z 
y y 0  0.8930   y 0
1/ 3
k   
hz 
0.8930  9z 
Este es el coeficiente convectivo de transferencia de calor a una distancia z del sitio donde comienza
el calentamiento. Como ß es la pendiente del perfil de velocidades en la superficie, y suponiendo que
el perfil de velocidades está completamente desarrollado, o se obtiene:
k TS  T0     3 hm Ts  To 
1
q y , med   
y 0 2  43   9  L  =
Este análisis es la base para el trabajo de Lévêque en tubos circulares cortos. Una de las soluciones
más sencillas para el coeficiente de transferencia de calor con flujo laminar en tubos circulares es la
de Lévêque. Este análisis se aplica directamente a la transferencia de calor laminar sobre una placa
plana, pero estos resultados pueden aplicarse fácilmente a tubos circulares. Esta solución se desarrolla
matemáticamente asumiendo que el perfil de temperaturas está confinado a una estrecha zona cercana
a la pared del tubo, como en los casos donde el caudal másico es alto y la longitud del tubo no es
grande, lo que ocurre para m’CP/kL > 20. Como criterio tenemos que:
El perfil de velocidad se halla completamente desarrollado a una distancia
LeM
 0.05 Re,
D
El perfil de temperatura se desarrolla completamente a:
LeT
 0.05 Re Pr
D ;
y el perfil de concentraciones a:
Lec
 0.05 Re Sc
D .
Consideramos, por ejemplo el caso de agua fluyendo a través de un tubo a Re = 100 y 68 °F. A esta
temperatura Pr = 7. El número de Schmidt para varios materiales en agua es substancialmente
dependiente del soluto, pero un valor típico es del orden de 1200.
Entonces las longitudes de entrada de concentraciones, térmica e hidrodinámica son respectivamente:

Lec
 6000
D
112
LeT
 35
D
Le
5
D
Esto nos indica que el perfil de concentraciones completamente desarrollado puede no obtenerse en
fluidos con un número de Schmidt grande y cuando la zona de calentamiento empieza luego de que se
ha desarrollado el perfil de velocidades podemos encontrar esta situación válida en transferencia de
calor con flujos grandes y poca longitud de conducto.
Para este caso, y debido a la poca penetración de la honda térmica podemos asimilar el sistema al de
una placa plana. Consideramos un fluido fluyendo sobre una superficie plana bajo las siguientes
condiciones:
- Las propiedades del fluido son constantes.
- La temperatura en la superficie es uniforme TS.
- La temperatura del fluido fuera de la capa límite es T0 (igual a la temperatura uniforme a la entrada
del conducto).
-La transferencia en la dirección y (radial) ocurre por difusión exclusivamente.
-Se asume distribución lineal de velocidades o sea :
Vz = ß(R. r) = ßy, donde ß es el gradiente de velocidad en la pared.
El perfil de velocidades para el tubo es (2.72):
  r 2  Vz 4Vm
Vz  2Vm 1     
  R   r rR R
rR y Vz V
 z
con r y
O sea:
1/3 1
k  8Vm   Vm  3
hz 
0.8930  9 Dz  = 1.076 k   zD 
o en modo adimensional
hz D
 1.076 Re 3 Pr 3 D z  3
1 1 1
Nu z 
k
113
Esta ecuación da los mismos valores que los obtenidos por Graetz para flujo parabólico desarrollado
en el rango RePr(D/z) > 100. Como la ecuación de Lévêque no es válida en la región más allá de
donde el perfil de temperaturas alcanza el centro de la tubería, la región donde coinciden ambas
soluciones puede tomarse como una medida de la longitud de la entrada térmica.
La suposición del gradiente de velocidad ß constante no es válida para la entrada inmediata al tubo
donde el gradiente es infinito; pero sin embargo, los gradientes de velocidad junto a la pared se
desarrollan rápidamente y pronto se aproximan al utilizado en el desarrollo anterior.
El coeficiente convectivo promedio se halla integrando sobre z y dividiendo por la longitud L del
conducto:
L
 hz dz 1 1
k  8v m  L  13  v 
3 3
hm  0
    z dz  1.614 k  m 
L 0.8930 L  9 D  0  DL 
El valor promedio para el coeficiente adimensional es:
hm D
 1.614 Re 3 Pr 3 D / L  3
1 1 1
Nu m 
k (2.77)
La ecuación (2.77) es la base de la correlación de Sieder y Tate (1936) quienes, para transferencia de
calor con flujo laminar en tubos propusieron:
1/ 3 0.14
D  b 
Nu  1.86 Re 1/ 3
Pr 1/ 3
   
L  s  (2.78)
donde µ es la viscosidad del fluido a la temperatura promedio global y µs es la viscosidad a la

temperatura de pared. El termino (µb/µs)0.14 es una corrección empírica para la distorsión del perfil de
velocidad resultante del efecto de la temperatura en la viscosidad.
El coeficiente obtenido con esta ecuación, usa como fuerza guía el promedio logarítmico de las
diferencias de temperatura en los extremos. Es conservativa según McAdams, por suponer perfil
desarrollado en tubos cortos.
ANEXO
"coeficiente condensado en tubos horizontales"
$UnitSystem SI Mass J kPa C rad
hs=0,725*(rho_cond^2*g*hfg*ks^3/(2*mus*DELTAT*ID))^0,25 "coeficientecondensación"
rho_cond=Density(R134a;T=T;x=x2)
ks=Conductivity(R134a;T=T;x=x2)
h1=Enthalpy(R134a;T=T;x=x1)
g=g#
DELTAT=5
h2=Enthalpy(R134a;T=T;x=x2)
T=39,4: x1=1: x2=0
hfg=(h1-h2)
mus=Viscosity(R134a;T=T;x=x2)
114
ID=0,008411
"Resistencias"
1/UA1=1/hs/Ai/Fs+1/ho/A_T/eta_T/Fs
Ai=0,1057: ho=44,68: A_T=1,874: eta_T=0,9559 "de problema anterior"
{Fs=0,7 "supuesto"}
"NUT"
m_hat_R=0,0028: m_hat_c=0,07224: Cp_c=1007: Tci=20
Qs=m_hat_R*hfg
Cc=m_hat_c*Cp_c: Cmins=Cc*Fs "es Cmin"
Qmaxs=Cmins*(T-Tci)
epsilon=Qs/Qmaxs
NUT=-ln(1-epsilon)
NUT=UA/Cmins
UA=UA1
Variables in Main
Ai=0,1057 [m^2]
A_T=1,874 [m^2]
Cc=72,75 [W/K]
Cmins=42,66 [W/K]
Cp_c=1007 [J/kg-K]
DELTAT=5 [C]
epsilon=0,5538 [-]
eta_T=0,9559 [-]
Fs=0,5864 [-]
g=9,807 [m/s^2]
h1=271052 [J/kg]
h2=107385 [J/kg]
hfg=163668 [J/kg]
ho=44,68 [W/m^2-K]
hs=2084 [W/m^2-K]
ID=0,008411 [m]
ks=0,07603 [W/m-K]
mus=0,0001624 [kg/m-s]
m_hat_c=0,07224 [kg/s]
m_hat_R=0,0028 [kg/s]
NUT=0,807 [-]
Qmaxs=827,5 [W]
Qs=458,3 [W]
rho_cond=1149 [kg/m^3]
T=39,4 [C]
Tci=20 [C]
UA=34,42 [W/K]
UA1=34,42 [W/K]
x1=1 [-]
x2=0 [-]
RBG Marzo 2021
115

Anexo A

Cargado por

Copyright:

Formatos disponibles

Anexo A

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Anexo A

Cargado por

Copyright:

Formatos disponibles

TABLA DE CONTENIDO

Las correlaciones empíricas generalmente se expresan en términos de un coeficiente de transferencia

En esta ecuación, Q es la velocidad de transferencia de calor entre la superficie sólida y el fluido, A

Desde el punto de vista de la transferencia de calor, el flujo turbulento es altamente deseable. En

En la superficie sólida, debido a la ley de no resbalamiento, la velocidad del fluido respecto a la

Como el gradiente en el origen no se conoce, h tampoco. Se postula entonces la existencia de un

La corriente de fluido tiene velocidad característica V , temperatura TF , concentración c AF . La

Reemplazando en (A.2) se obtiene h  k / z F correlacionando el espesor de la película imaginaria con

La cantidad h es una función de las propiedades del fluido, la geometría y la rugosidad de la

Evaporación Condensación Disolución Cristalización

Las concentraciones en el régimen turbulento generalmente fluctúan alrededor de un valor medio.

La relación entre la difusividad y el espesor de la película DAB / z F se fusiona en un único término

Para el espesor efectivo de la película, z FM  DAB / kc , los valores correspondientes son:

Problema de práctica, evaporación de un derrame de solvente: Se derrama una carga de solvente

Energía Ley de Fourier.

En cada caso las ecuaciones toman la forma:

Densidad de Flujo = (Propiedad de Transporte)x(Gradiente de Potencial)

Difusividad Térmica:   k /  CP , donde CP es capacidad calorífica a presión constante.

Transporte de cantidad de movimiento x en la dirección y por difusión

Transferencia de calor por conducción Transferencia de masa por difusión

Flujo =  Difusividad x Gradiente de Concentración.

Podemos generalizar esta similitud y escribirla como  mz      / z  , siendo  mz cualquiera de

gradiente se considera constante. La fuerza motriz se convierte simplemente en la diferencia de

Los flujos convectivos:

Flujo molecular = (Difusividad)(concentración)/(Longitud característica) (ecuaciones A.4)

Flujo de interface = (coeficiente)(concentración) (ecuaciones A.5)

La relación (A.5)/(A.4) dará (coeficiente)(longitud característica)/(difusividad), obteniendo

Transferencia de calor, número de Nusselt

Transferencia de masa, número de Sherwood

Re representa el número de Reynolds, relación de fuerzas inerciales y viscosas en un flujo:

La velocidad característica Vcar y la longitud característica Lcar se definen en cada geometría.

 es el coeficiente de expansión volumétrica para el fluido:

Si este fluido puede considerarse gas ideal entonces   1 / T

Para transferencia de masa definimos un “coeficiente de expansión” originado en los gradientes de

Consecuentemente definimos un número de Grashof para transferencia de masa:

Las correlaciones adimensionales toman la forma

Frecuentemente se usa el número de Rayleigh:

Es importante destacar la relación existente entre diferentes grupos adimensionales:

El factor de fricción de Darcy f D . En la literatura Europea puede encontrarse un factor de fricción

CORRELACIONES PARA ESTIMAR COEFICIENTES CONVECTIVOS

A continuación, se escriben correlaciones útiles para encontrar los coeficientes convectivos en

Placa Plana Horizontal

Figura A.3: Flujo en convección forzada sobre una placa plana

La situación se describe en la figura A.4; la longitud de la placa en la dirección de flujo es L y su

Coeficiente local de fricción

Coeficiente local de transferencia de calor

Coeficiente local de transferencia de Masa

Flujo turbulento: En la zona donde el régimen es turbulento se aplican las correlaciones:

Factor de ficción local

Coeficiente de transferencia de calor local

Coeficiente de transferencia de masa local

Los coeficientes de fricción f x y convectivos hx y kcx cambian con la posición, y si se quieren

El coeficiente promedio de transferencia de calor se obtiene de las siguientes correlaciones:

0.6  Pr  60; 5x105  ReL  108.

Para escribir en EES: NuL=(0,664*Rec^0,5+0,0365*(Re^0,8-Rec^0,8))*Pr^(1/3)

Aquí A es la superficie total para transferencia de calor por convección.

Para transferencia de masa cambiamos Nu por Sh y Pr por Sc, a saber:

0.6  Sc  3000, 5x105  ReL  108.

Para flujo laminar de metales líquidos sobre placas planas:

Nu x  0.564 Pex0.5 Pr  0.05, Pex  100. (A.13)

Para escribir en EES: NuL=(0,664Rec^0,5+0,0365(Re^0,8-Rec^0,8))*Pr^(1/3)