PROBLEMAS RESUELTOS DE

TRANSMISIÓN DE CALOR

Mª Isabel Lamas Galdo

PROBLEMAS RESUELTOS DE TRANSMISIÓN DE CALOR
Autoría de: Mª Isabel Lamas Galdo

Idioma de la obra: Español

Nº de páginas: 105 pp.
Índice: ii

Depósito legal: C 287-2022

ISBN: 978-84-9749-835-7

DOI: https://doi.org/10.17979/spudc.9788497498357

Edición: Universidade da Coruña, Servizo de Publicacións <http://www.udc.gal/publicacions>

© de los textos, sus autores
© de las imágenes, sus autores
© 2021 de la edición, Universidade da Coruña

Esta obra se publica bajo una licencia Creative Commons

Atribución-NoComercial-CompartirIgual 4.0 Internacional
(CC BY-NC-SA 4.0)

i
ÍNDICE

Prólogo ............................................................................................................................................................ iii

Nomenclatura ................................................................................................................................................. iv

Tema 1. Introducción a la transmisión de calor ......................................................................................... 1

Tema 2. Conducción de calor unidimensional en estado estacionario .................................................. 8

Tema 3. Conducción de calor unidimensional en régimen transitorio...................................................27

Tema 4. Convección forzada. Flujo exterior ..............................................................................................42

Tema 5. Convección forzada. Flujo interior ...............................................................................................55

Tema 6. Convección libre .............................................................................................................................62

Tema 7. Condensación y ebullición ............................................................................................................75

Tema 8. Radiación térmica ..........................................................................................................................81

Tema 9. Intercambiadores de calor ............................................................................................................86 

ii
PRÓLOGO

El presente texto tiene como objetivo servir de apoyo a las asignaturas Transmisión de Calor y
Frío Industrial/Refrigeración de los grados en Ingeniería Naval y Oceánica, Ingeniería
Mecánica e Ingeniería en Tecnologías Industriales que se imparten en la Escuela Politécnica
Superior de la Universidade da Coruña. Su contenido también puede ser usado en asignaturas
relacionadas con la transmisión de calor en otros grados y universidades.

El libro consta básicamente de problemas resueltos procedentes del material utilizado para
docencia, así como de exámenes realizados a lo largo de varios años. Esta colección de
problemas aborda una visión global de la transmisión de calor, abarcando desde los fundamentos
hasta los modos de transferencia de calor: conducción, convección y radiación, así como
aplicaciones prácticas.

iii
NOMENCLATURA

Símbolo Descripción Unidades

c Calor específico J/kg∙K

C Capacidad térmica W/K
A Área m2
D Diámetro m
e Espesor m
E Energía J
g Aceleración de la gravedad o energía libre de Gibbs m/s2 o J/kg
h Coeficiente de transferencia de calor por convección o entalpía W/m2∙K o kJ/kg
h Coeficiente de transferencia de calor por convección promedio W/m2∙K
hlv Calor latente de evaporación J/kg∙K
I Intensidad A
k Conductividad térmica W/m∙K
m Masa kg

m Caudal másico kg/s
p Presión Pa
p Perímetro m
Q Calor J
Q Tasa de transferencia de calor W
q ' Tasa de transferencia de calor por unidad de longitud W/m
q '' Tasa de transferencia de calor por unidad de área W/m2
R Radio o resistencia térmica m o W/K
t Tiempo o espesor de aleta som
T Temperatura K o ºC
T Temperatura promedio K o ºC
u Velocidad m/s
W Trabajo J

W Potencia W

iv
Parámetros adimensionales
Bi Número de Biot
Fo Número de Fourier
Gr Número de Grashof
Gz Número de Graetz
Ja Número de Jacob
Nu Número de Nusselt
Nu Número de Nusselt promedio
Pr Número de Prandtl
Ra Número de Rayleigh
Re Número de Reynolds
Pe Número de Péclèt

Símbolos griegos

α Difusividad térmica o absortividad m2/s o -

ε  Emisividad -
ƞ  Eficiencia -
θ  Temperatura adimensional -
λ  Autovalor rad
μ Viscosidad dinámica Pa∙s
ρ Densidad kg/m3
σ Tensión superficial o constante de Stephan-Boltzmann N/m∙K o W/m2∙K4
ν Viscosidad cinemática m2/s

Subíndices
alr Alrededores
cond Conducción
conv Convección
D Diámetro
f Película (film)
gen Generación
i Interno
L Lontigud
m Promedio
max Máximo
min Mínimo
rad Radiación

v
 TEMA 1

INTRODUCCIÓN A LA
TRANSMISIÓN DE CALOR

1
Problema 1.1
Por el interior de un conducto caliente fluyen 0,3 kg/s de aire, el cual entra a 25 ºC y sale a 32 ºC.
Para simplificar los cálculos asumir que el calor específico del aire es de cp = 1 kJ/kg∙K. Calcular la
tasa de transferencia de calor que recibe el aire desde la entrada hasta la salida del conducto en
estado estacionario.

Solución:

Aplicando balance de energía al sistema abierto que constituye el conducto:

dE 
dt
QW 
  m
e
 e (he  v e2 / 2  gz e )   m
s
 s (hs  v s2 / 2  gz s )

En estado estacionario, con potencia nula y despreciando variaciones de energía cinética y

potencial:

dE
dt
 W
Q e  e (he  v e2 / 2  gze )  s m s (hs  v 2s / 2  gz s )
  m

  m(h
0Q  e  hs )

Asumiendo el aire como gas ideal con cp constante, he – hs = cp(Te - Ts), por tanto:

  mc
0Q  p (Te  Ts )

  mc
Q  p (Ts  Te )  0,3  1 (32  25)  2,1 kW

Problema 1.2
En una etapa de un proceso de recocido, una hoja de acero inoxidable 304 (ρ = 7900 kg/m3, c =
578 J/kgꞏK) se lleva de 300 K a 1250 K conforme pasa a una velocidad de 10 mm/s a través de un
horno calentado eléctricamente. El espesor y ancho de la hoja son 8 mm y 2 m (dimensión
perpendicular al papel), respectivamente. Calcular la potencia eléctrica que se requiere suministrar
al horno asumiendo que en estado estacionario un 15% de esa potencia se pierde como
disipación de calor hacia el exterior.

2
Solución:

Aplicando balance de energía a este sistema abierto:

dE 
dt
QW 
  m
e
 e (he  v e2 / 2  gz e )   m
s
 s (hs  v s2 / 2  gz s )

En estado estacionario, con potencia nula y despreciando variaciones de energía cinética y

potencial:

dE
dt
 W
Q 
  m
e
 e (he  v e2 / 2  gz e )   m
s
 s (hs  v 2s / 2  gz s )

  m(h
0Q  e  hs )

Si las pérdidas representan el 15% de la potencia total resulta:


0Q  
elec  0,15Qelec  m(he  hs )

Como el calor específico es constante, se cumple que he – hs = c(Te - Ts). Por tanto:


0Q  
elec  0,15Q elec  mc(Te  Ts )


0 = 0,85Q elec + ρVewc(Te - Ts )

siendo ρ la densidad, V la velocidad, e el espesor y w el ancho. Sustituyendo valores resulta:

 7900  0,01 0,008  2  578  (1250  300)

Q elec   816544 W
0,85

Problema 1.3
Por un cable desnudo de cobre de propiedades k = 395 W/mꞏK y α = 10-4 m2/s, diámetro D = 1,4
mm y resistencia eléctrica por unidad de longitud R´ = 0,01 Ω/m se transmite una intensidad de I =
150 A debido a una avería eléctrica. El ambiente y alrededores se encuentran a 20 °C y el
coeficiente de transferencia de calor combinado convección-radiación linealizado es 25 W/m2ꞏK.
Asumir que la temperatura del cable es uniforme para cada instante de tiempo y que inicialmente
se encuentra a 20 ºC. Sabiendo que la temperatura de fusión del cobre es de 1080 ºC, calcular el
tiempo que tarda en fundirse el cable.
NOTA: Por efecto Joule q 'Joule  I2R ' (W/m).

Solución:

Aplicando balance de energía al sistema cerrado que constituye el cable:

3
 dE  
QW
dt

Con potencia nula:

dE  
Q W
dt

La energía total es la suma de la contribución de la energía interna, cinética, potencial y otras que
no se tratarán en el presente texto. Despreciando variaciones de energía cinética y potencial, la
energía total resulta:

E  me  m(u  v 2 / 2  gz )  mu

Por tanto, el balance de energía queda:

d(mu) 
Q
dt

Teniendo en cuenta el calor por radiación, convección y la generación debida a efecto Joule,
resulta:

d(mu)
 hconv rad A (T  T )  I2R 'L
dt

Nótese que en la expresión anterior se ha tenido en cuenta el coeficiente de transferencia de

calor combinado convección-radiación linealizado, hconv-rad, y que el calor debido a convección +
radiación se ha puesto negativo porque es calor disipado por el sistema. Haciendo operaciones:

d(ρcT)
 hconv rad πDL(T  T )  I2R'L
dt

siendo ρ la densidad y el  volumen.

dT
ρc  hconv rad πDL(T  T )  I2R 'L
dt

πD2 L dT
ρ c  hconv rad πD L (T  T )  I2R' L
4 dt

4
 πD2
ρ c
4 dT
1
hconv rad πD(T  T )  I R dt
2 '

πD2
ρ c t
4
hconv rad πD(T  T )  I2R '
dT  0 dt  t  40,46 s

Al resolver la ecuación anterior tener en cuenta que α = k/(ρc).

Problema 1.4
Un secador de pelo se puede idealizar como un conducto circular por el que un ventilador dirige
aire ambiente que se calienta al pasar sobre una resistencia eléctrica. En estado estacionario, un
secador opera con una potencia eléctrica de 500 W y calienta aire desde una temperatura
ambiente de 20 ºC hasta 45 ºC.
a) Calcular el flujo másico de aire. Asumir que para el aire cp = 1000 J/kg∙K y despreciar las
pérdidas de calor hacia el exterior.
b) Considerar un secador con un conducto de 150 mm de largo y 70 mm de diámetro externo con
una emisividad superficial de 0,8. Si el coeficiente de transferencia de calor con el aire ambiente
exterior es de 4 W/m2∙K y tanto la temperatura del aire exterior como de los alrededores es de 20
ºC, calcular las pérdidas de calor hacia el exterior del conducto y confirmar que son efectivamente
despreciables. Para simplificar los cálculos, asumir que la superficie exterior del conducto en
contacto con el exterior se encuentra a una temperatura uniforme de 30 ºC.

Solución apartado a:

Aplicando balance de energía a este sistema abierto:

dE 
dt
QW 
  m
e
 e (he  v e2 / 2  gz e )   m
s
 s (hs  v s2 / 2  gz s )

En estado estacionario y despreciando variaciones de energía cinética y potencial:

dE
dt
 W
Q 
  m
e
 e (he  v e2 / 2  gz e )   m
s
 s (hs  v 2s / 2  gz s )

 W 
 Q  
siendo Q resistencia  Qpérdidas  Wventilador . Según el enunciado, en el apartado a se

 W
desprecian las pérdidas de calor hacia el exterior, con lo cual resulta Q   500 W, por tanto:

 e  hs )
0  500  m(h

Como cp es constante:

 p (Te  Ts )
0  500  mc
5
   1000  (20  45)
0  500  m    0,02 kg/s
m

Solución apartado b:

Se procederá a calcular las pérdidas de calor y verificar que son considerablemente menores que
500 W:

  
pérdidas  Qconv  Qrad  hA (Ts  T )  εAσ (Ts  Talr )  hπDL (Ts  T )  επDLσ (Ts  Talr ) 
4 4 4 4
Q

 4  π  0, 07  0, 15  (Ts  T )  0, 8  π  0, 07  0,15  5, 67  10 8 (303 4  293 4 ) 

 1, 31 1, 58  2, 89 W  500 W

Nótese que en la expresión anterior, se han indicado las temperaturas en grados kelvin al calcular
Ts4  Talr4 .

Como se puede observar, al calcular las pérdidas de calor resulta 2,89 W, considerablemente
menor que 500 W. Por tanto, efectivamente las pérdidas de calor son despreciables.

Problema 1.5
Por un conducto de aire acondicionado de 15 m de longitud, sección transversal exterior 200ˣ350
mm y emisividad superficial 0,5 fluyen 0,3 kg/s de aire. Exteriormente al conducto hay aire
ambiente cuyo coeficiente de transferencia de calor por convección es de 4 W/m2∙K y tanto la
temperatura del aire ambiente exterior como de los alrededores es de 24 ºC. Para simplificar los
cálculos asumir que el calor específico a presión constante del aire es de 1 kJ/kg∙K y que en
estado estacionario la superficie exterior del conducto en contacto con el exterior se encuentra a
una temperatura promedio de 22 ºC que, para simplificar los cálculos, se puede asumir uniforme.
Sabiendo que el aire entra en el conducto a 12,5 ºC, calcular la temperatura del aire a la salida del
conducto.

6
Solución:

Llamando 1 a la entrada, 2 a la salida y aplicando balance de energía a este sistema abierto:

dE 
dt
QW e  e (he  v e2 / 2  gz e )  s m s (hs  v s2 / 2  gzs )
  m

dE 
dt
QW 
  m
1
 1 (h1  v12 / 2  gz1 )   m
2
 2 (h2  v 22 / 2  gz 2 )

En estado estacionario, potencia nula y despreciando variaciones de energía cinética y potencial:

dE
dt
 W
Q 
  m
1
 1 (h1  v12 / 2  gz1 )   m
2
 2 (h2  v 22 / 2  gz 2 )


0Q  
conv  Qrad  m(h1  h2 )

Asumiendo el aire como gas ideal con cp constante, h1 – h2 = cp(T1 – T2), por tanto:


0Q  
conv  Qrad  mc p ( T1  T2 )

Sustituyendo las expresiones correspondientes para convección y radiación resulta:

 p (T1  T2 )
0  hA(T  Ts )  σεA(T4  Ts4 )  mc  T2  13,26 ºC

siendo el área de transferencia de calor: A = 2∙(0,2 + 0,35) ∙15 = 16,5 m2.

Nótese que en la expresión anterior se ha puesto como negativo el calor disipado por el sistema.

7
 TEMA 2

CONDUCCIÓN DE CALOR
UNIDIMENSIONAL EN
ESTADO ESTACIONARIO

8
Problema 2.1
Una placa de acero de espesor 3 mm y conductividad térmica 61 W/m∙K se sitúa sobre una placa
de aluminio de 8 mm de espesor y conductividad térmica 238 W/m∙K. La resistencia de contacto
entre ambas placas es de 10-5 m2∙K/W. La superficie superior del acero se encuentra a 35 ºC,
mientras que la superficie inferior del aluminio a 30 ºC. Calcular el salto térmico a través de la
resistencia de contacto.

Solución:

Se resolverá aplicando analogía eléctrica. El circuito eléctrico análogo es el siguiente:

La tasa de transferencia de calor por unidad de área viene dada por:

Ts  Ti
q '' 
R ''
conds  R ''tc  Rcondi
''

A su vez, la tasa de transferencia de calor por unidad de área también se puede expresar como:

T
q ''  '' tc
R tc

Igualando se obtiene que q ''  53883 W/m2 y el salto térmico a través de la resistencia de contacto
Ttc  0,54 ºC.

Problema 2.2
Una pared plana se compone de dos materiales. El material A tiene una generación de calor
uniforme de 106 W/m3, conductividad térmica 50 W/mꞏK y espesor 50 mm. El material B no tiene
generación y su conductividad y espesor son 20 W/mꞏK y 10 mm, respectivamente. La superficie
interior del material A está bien aislada, mientras que la superficie exterior del material B se expone
a un ambiente convectivo a T∞ = 25 ºC y h = 500 W/m2ꞏK. En estado estacionario y despreciando la
transferencia de calor por radiación, determinar las temperaturas T (x = 0), T (x = 50 mm) y T (x =
60 mm).
9
Solución:

Solamente hay conducción de calor a lo largo del eje x.

Cálculo de T(x = 60 mm = 0,06 m):

  h[T ( x  0, 06)  T ]
Balance de energía alrededor del material B: q
''

 ''  q '''gen e A  50000 W / m 2

Balance de energía alrededor del material A: q

Por tanto:

h[T ( x  0, 06 )  T ]  q '''gen e A  T ( x  0, 06 )  125 º C

Cálculo de T(x = 50 mm = 0,05 m):

T ( x  0, 05 )  T ( x  0, 06)
q ''  k B  T ( x  0, 05 )  150 º C
eB

Otra forma de calcular T(x = 0,05) sería aplicando analogía eléctrica a partir del material B, nótese
que en el material A no se puede aplicar analogía eléctrica debido a que hay generación de calor:

T ( x  0, 05 )  T
q ''   T ( x  0, 05 )  150 º C
eB 1

kB h

Cálculo de T(x = 0):

Como no se puede aplicar analogía eléctrica al material A, se calcula T(x = 0) aplicando la

ecuación general de conducción del calor. El subíndice A de TA se indica para referirse a la
temperatura del sólido A:

10
 q '''gen T
 TA 
2

c t

En estado estacionario, coordenadas cartesianas y conducción de calor solamente a lo largo del

eje x queda:

'''
d2 TA q gen
  0
dt 2 c

Teniendo en cuenta que α = k/(ρc):

'''
d2 TA q gen
 0
dt 2 kA

Resolviendo la ecuación diferencial queda:

q '''gen x 2
TA    C1x  C 2
kA 2

Aplicando la condición de contorno en x = 0 se obtiene que C1 = 0:

dTA
x  0: 0  C1  0
dx x 0

Otra forma de llegar a la misma conclusión de que C1 = 0 es mediante la siguiente condición de

contorno en x = 0,05 m:

dTA dTB
x  0, 05 : k A  k B  C1  0
dx x  0 ,05 dx x  0 ,05

Ya que al principio de la solución se ha calculado que T(x = 0,05) 150ºC y se desprecia que haya
resistencia térmica de contacto entre los materiales A y B, para calcular C2 se puede aplicar la
siguiente condición de contorno en x = 0,05 m:

q '''gen e 2A
x  0, 05 : TA  150 º C  C 2  TA ( x  0, 05 )   175
kA 2
Sustituyendo valores, la temperatura de A resulta:

q '''gen x 2
TA    C1x  C 2  10000x 2  175
kA 2

Por tanto, la temperatura de A en x = 0 resulta:

TA ( x  0 )  175 º C

11
Problema 2.3
Del tubo de la figura se conoce la conductividad térmica, k = 0,16 W/mꞏK y los radios interior y
exterior, ri = 10 mm y re = 12 mm. El tubo experimenta transferencia de calor por convección con un
fluido que circula por el interior, a T∞i = 200 ºC y hi = 1000 W/m2ꞏK y con un fluido que circula por el
exterior a he = 50 W/m2ꞏK y T∞e = 15 ºC. Calcular el flujo de calor por unidad de longitud de tubo
mediante los siguientes procedimientos.
a) Aplicando la ecuación general de conducción de calor.
b) Aplicando analogía eléctrica.

Solución apartado a:

Ecuación general de conducción de calor:

q '''gen T
 T 
2

c t

En coordenadas cilíndricas:

 1   T  1   T   2 T  q '''gen T
 r  2   2  
 r r  r  r     z  c t

En estado estacionario, sin generación de energía y con conducción de calor solamente a lo largo
de la coordenada radial queda:

 1   T  1   T   2 T  q '''gen T
 r  2    
     z c t
2
 r r  r  r  

1 d  dT 
 r 0
r dr  dr 

d  dT 
r 0
dr  dr 

Resolviendo resulta:

T  C1 ln(r)  C 2

Las constantes de integración C1 y C2 se determinan aplicando condiciones de contorno:

r  ri :  kA i  T r r  hA i [Ti  T(r  ri )] 
i

r  re :  kA e  T r r  hA e [Te  T(r  re )] 
e 

12
 C1 
k  hi [Ti  C1 ln(ri )  C 2 ] 
ri 

C
k 1  h e [C1 ln(re )  C 2  Te ] 
re 

C1 
0,16  1000[200  C1 ln(0,01)  C 2 ] 
0,01 

C1
0,16  50[C1 ln(0,012)  C 2  15] 
0,012 

Resolviendo este sistema de dos ecuaciones y dos incógnitas resulta C1 = -397,9 y C2 = -1639.
Por tanto, la distribución de temperatura viene dada por la siguiente expresión:

T  C1 ln(r)  C 2  397,9 ln(r)  1639

El flujo de calor se obtiene mediante la ley de Fourier:

  k  T  k2 rL  T
Q


Q
q '   k2 r T
L

Sustituyendo valores se obtiene el flujo de calor en r = ri o en r = re que, lógicamente, tienen que

resultar idénticos al no existir generación de energía:

dT
q r' ri  k2 ri  T r r  k2 ri  400 W / m 2
i
dr r  ri

dT
q r' re  k2 re  T r r  k2 re  400 W / m 2
e
dr r  re

Solución apartado b:

El circuito eléctrico análogo es el siguiente:

El flujo de calor viene dado por:

Ti  Te
q '   400 W / m 2
1 ln(re / ri ) 1
 
2 rh
i i 2 k 2 reh e

13
Problema 2.4
Por un cable de cobre (kc = 395 W/mꞏK) de diámetro D = 1,4 mm y resistencia eléctrica por unidad
de longitud R´ = 0,01 Ω/m se transmite una intensidad de I = 10 A. Dicho cable se cubre con un
aislante de espesor ea = 1 mm y conductividad térmica ka = 0,2 W/mꞏK. Se mide la temperatura de
la superficie externa del aislante y resulta 28 ºC.
a) Calcular la temperatura de la superficie del cable.
b) Calcular la temperatura del centro del cable.
  I2R ' (W/m).
NOTA: Por efecto Joule se disipa q'

Solución apartado a:

Aplicando analogía eléctrica:

Tsc  Tsa
q ' 
R 'cond

 '  I2R '  1 W / m . Por tanto:

siendo q

Tsc  28
1  Tsc  28,71 º C
ln(r2 / r1 )
2 k a

Solución apartado b:

En el cobre hay generación de energía, por tanto no se puede aplicar analogía eléctrica a este
material. Por tanto, el procedimiento para resolver este apartado es mediante la ecuación general
de conducción de calor:

q '''gen Tc
 c  Tc 
2

c c c t
Nótese que se ha puesto el subíndice c refiriéndose al cobre.

14
  gen  q ' q '
El calor por unidad de volumen es q
'''
  649612 W / m 3
A D / 42

En coordenadas cilíndricas la ecuación general de conducción de calor resulta:

 1   Tc  1   Tc   2 Tc  q gen

'''
T
c  r  2   2 
  c
 r r  r  r     z   c c c t

En estado estacionario y con conducción de calor solamente a lo largo de la coordenada radial

queda:

 1   T  1   Tc   Tc  q gen
2 '''
Tc
c  r
c
     2 
 
 r r  r  r     z  c c c t
2

'''
1 d  dTc  q gen
c r  0
r dr  dr  c c c

Teniendo en cuenta que αc = kc/(ρccc):

'''
1 d  dT  q gen
r  0
r dr  dr  k c

Resolviendo se obtiene la distribución de temperaturas en función de las constantes de

integración C1 y C2:

q '''genr 2
Tc    C1 ln(r)  C 2
4k c

Las constantes de integración C1 y C2 se determinan aplicando condiciones de contorno:

dT
r 0:  0  C1  0
dr
D
r  r1  1 : T  Tsc  28,71 º C  C2  28,71 º C
2

Por tanto, la distribución de temperatura (en grados centígrados) viene dada por la siguiente
expresión:

q '''genr 2
Tc    C1 ln(r)  C 2  411,14r 2  28,71
4k c

Particularmente, para r = 0 se corresponde con la temperatura del centro del cable:

Tc (r  0)  28,71 º C

15
Como se puede observar, esta temperatura resulta igual que la temperatura de la superficie, lo
cual no es de extrañar debido a la elevada conductividad térmica del cobre y al reducido diámetro
del cable.

Problema 2.5
Una aleta cilíndrica de aluminio (k = 236 W/mꞏK) de longitud 100 mm y diámetro 5 mm cuya base
se encuentra a 50 ºC se expone al aire ambiente a T∞ = 20 ºC y h = 30 W/m2ꞏK. Asumir
transferencia de calor unidimensional. Representar gráficamente la temperatura a lo largo de la
longitud.

Solución:

La distribución de temperatura a lo largo del eje x se obtiene a partir de la siguiente expresión:

 T  T cosh[m(L  x)]  (h / mk)senh[m(L  x)]
 
b Tb  T cosh(mL)  (h / mk)senh(mL)

siendo:

hP hD 4h
m    10,08 m1
kA c kD / 4
2
kD

Por tanto, la representación gráfica de temperatura a lo largo de la longitud resulta:

50

40

30
T (ºC)

20

10

0
0 0,02 0,04 0,06 0,08 0,1
x (m)

Problema 2.6
Calcular la pérdida de calor a través de las aletas siguientes:
a) Una aleta cilíndrica de aluminio (k = 236 W/mꞏK) de longitud 100 mm y diámetro 5 mm cuya
base se encuentra a 50 ºC se expone al aire ambiente a T∞ = 20ºC y h = 30 W/m2ꞏK.
b) Una aleta cuadrada de aluminio (k = 236 W/mꞏK) de longitud 100 mm y lado 5 mm cuya base
se encuentra a 50 ºC se expone al aire ambiente a T∞ = 20ºC y h = 30 W/m2ꞏK.
c) Una aleta anular de aluminio (k = 236 W/mꞏK) de radio externo 100 mm, radio interno 50 mm y
espesor 5 mm cuya base se encuentra a 50ºC se expone al aire ambiente a T∞ = 20 ºC y h = 30
W/m2ꞏK.

16
Solución apartado a:

Una forma de resolverlo es mediante la siguiente expresión:

   sen h(mL)  (h / mk)cosh(mL)

Q a
cosh(mL)  (h / mk)senh(mL)

siendo:

hP hD 4h
m    10,08 m1
kA c kD2 / 4 kD

D2
  b hPkA c  (Tb  T ) hPkA c  (Tb  T ) hDk  1,401 W
4

Por tanto:

h
sen h(mL)  cosh(mL)

Qa   mk 
h
cosh(mL)  senh(mL)
mk

30
senh(10,08  0,1)  cosh(10,08  0,1)
10,08  236
 1,401  1,08 W
30
cosh(10,08  0,1)  senh(10,08  0,1)
10,08  236

Al calcular el la pérdida de calor con la expresión anterior, tener en cuenta que la calculadora debe
de estar en radianes.

Otra forma de calcular la pérdida de calor es utilizando la siguiente expresión, correspondiente a

aleta con extremo adiabático:

   tan h(mL )
Q a c

Siendo Lc la longitud corregida. Hay que utilizar la longitud corregida en lugar de la longitud real
para poder aplicar la expresión anterior, correspondiente a aleta con extremo adiabático a una
aleta con extremo que intercambia calor por convección. Para esta aleta cilíndrica, la longitud
corregida viene dada por:

R D
Lc  L   L   0,10125 m
2 4

Por tanto:

   tan h(mL )  1, 401  tanh(10,08  0,10125)  1,08 W

Q a c

Otra forma de calcular la pérdida de calor es utilizando el concepto de eficiencia de aleta, según el
cual la pérdida de calor a través de la aleta es igual a la eficiencia de aleta multiplicada por el calor
perdido por la aleta ideal:

 Q
Q 
a a aideal
17
La aleta ideal es aquella que se encuentra íntegramente a la temperatura de la base, por tanto:


Q aideal  hA a (Tb  T )

La eficiencia de aleta para esta aleta cilíndrica viene dada por:

tanh(mL c )
a   0,7544
mL c

Nótese que la anterior definición de eficiencia se basa en el concepto de longitud corregida. Por
tanto, sustituyendo estas expresiones, la pérdida de calor a través de la aleta resulta:

  Q
Q a a aideal  a hA a (Tb  T )  a h DL c (Tb  T )  1,08 W

Solución apartado b:

Una forma de resolverlo es mediante la siguiente expresión:

h
senh(mL) 
cosh(mL)
 
Q mk 
a
h
cosh(mL)  senh(mL)
mk

siendo:

hP h4t 4h
m    10,08 m1
kA c kt 2 kt

  b hPkA c  (Tb  T ) hPkA c  (Tb  T ) h4tkt 2  1,785 W

Por tanto:

   sen h(mL)  (h / mk)cosh(mL)  1,375 W

Q a
cosh(mL)  (h / mk)senh(mL)

Al calcular el la pérdida de calor con la expresión anterior, tener en cuenta que la calculadora debe
de estar en radianes.

Otra forma de calcular la pérdida de calor es utilizando la siguiente expresión, correspondiente a

aleta con extremo adiabático:

   tan h(mL )
Q a c

Siendo Lc la longitud corregida. Hay que utilizar la longitud corregida en lugar de la longitud real
para poder aplicar la expresión anterior, correspondiente a aleta con extremo adiabático a una
aleta con extremo que intercambia calor por convección. Para esta aleta cuadrada, la longitud
corregida viene dada por:

18
 t
Lc  L   0,10125 m
4
Por tanto:

   tan h(mL )  1,375 W

Q a c

Otra forma de calcular la pérdida de calor es utilizando el concepto de eficiencia de aleta, según el
cual la pérdida de calor a través de la aleta es igual a la eficiencia de aleta multiplicada por el calor
perdido por la aleta ideal:

  Q
Q 
a a aideal

La aleta ideal es aquella que se encuentra íntegramente a la temperatura de la base, por tanto:


Q aideal  hA a (Tb  T )

La eficiencia de aleta para esta aleta cuadrada viene dada por:

tanh(mL c )
a   0,7544
mL c

Nótese que la anterior definición de eficiencia se basa en el concepto de longitud corregida. Por
tanto, sustituyendo estas expresiones, la pérdida de calor a través de la aleta resulta:

 Q 
Q a a aideal  a hA a (Tb  T )  a h4tL c (Tb  T )  1,375 W

Solución apartado c:

Se procederá a calcular la pérdida de calor utilizando el concepto de eficiencia de aleta, según el

cual la pérdida de calor a través de la aleta es igual a la eficiencia de aleta multiplicada por el calor
perdido por la aleta ideal:

 Q
Q 
a a aideal

La aleta ideal es aquella que se encuentra íntegramente a la temperatura de la base, por tanto:


Q aideal  hA a (Tb  T )

Para determinar la eficiencia de esta aleta anular se utilizará la siguiente gráfica:

19
 siendo:

r1  0,05 m
r2  0,1 m
L  0,05 m
t 0,005
r2c  r2   0,1   0,1025 m
2 2
t 0,005
L c  L   0,05   0,0525 m
2 2
A p  L c t  0,0002625 m 2
h
L3c/ 2  0,2647
kA p
r2c / r1  2,05

Entrando en la gráfica con un valor de L3c/ 2 h

 0,2647 y r2c/r1 = 2,05 se obtiene que la
kA p
eficiencia de aleta resulta aproximadamente 92% (0,92 en base unitaria). Por tanto:

 Q 
aideal   a hA a (Tb  T )  a h2 (r2c  r1 )(Tb  T )  41,63 W
2 2
Q a a

Téngase en cuenta que el área de aleta es la que se indica en color azul en la figura (por arriba y
por abajo, que no se percibe en la figura), teniendo en cuenta que al utilizar el modelo de
eficiencia de aleta no se tiene en consideración el extremo indicado en color negro:

20
Problema 2.7
Una tarjeta de circuito de espesor despreciable y dimensiones 10 ˣ 10 mm disipa 1 W en estado
estacionario. Para simplificar los cálculos, este calor disipado se puede asumir uniformemente
distribuido. Para refrigerarla se dispone de 4 ˣ 4 aletas rectas de cobre, de sección circular,
diámetro 1,5 mm, longitud 5 mm y conductividad térmica 200 W/m∙K. El circuito se instala sobre
una base de 15 mm de espesor y conductividad térmica 1 W/m∙K. La resistencia térmica de
contacto circuito-base es de 10-4 m2∙K/W. Por arriba y abajo circula una corriente de aire que
ayuda a enfriar el circuito. Este aire se encuentra a 20 ºC y presenta un coeficiente de
transferencia de calor por convección de 40 W/m2∙K, tanto por la parte superior como por la
inferior.
a) Determinar la temperatura de la tarjeta de circuito.
b) Repetir el apartado anterior en caso de que no hubiera aletas.

Solución apartado a:

Este apartado se resolverá utilizando analogía eléctrica. El circuito eléctrico análogo es el que se
indica en la figura:

21
 
Del calor total que genera la tarjeta de circuitos, Q  ,y
, parte se disipa hacia la parte inferior, Q
total i

parte hacia la superior, Q s :

  
Q total  Q s  Q i

 Tc  T Tc  T
Q total  
1 R tc  R cond  R convi
1 1

R a R sl

siendo:


Q total  1 kW

R ''tc 10 4
R tc   1K/W
A 0,01  0,01

eb 0,015
R cond    150 K / W
k b A 1  0,01  0,01

1 1
R convi    250 K / W
hA 40  0,01  0,01

1 1
R sl    348,5 K / W
hA sl 40  (0,01  16    0,0015 2 / 4)
2

1
Ra 
hA a a

La eficiencia de aleta se calcula mediante la siguiente expresión:

tgh(mL c )
a 
mL c

Utilizando el modelo de longitud corregida, Lc = L + R/2 = 0,005375 m.

hP hD 4h
m    23,09 m1
kA c k aletas D / 4
2
k aletasD

A a  Naletas DLc  0,0004 m2

Por tanto, la resistencia que representa la superficie aleteada resulta:

1
Ra   62,82 K / W
hA a a

Y la temperatura de la tarjeta de circuito se obtiene a partir de la siguiente expresión:

22
  Tc  T Tc  T
Q total  
1 R tc  R cond  R convi
1 1

R a R sl

Tc  20 Tc  20
1   Tc  67,1 º C
1 1  150  250
1 1

62,82 348,5

Solución apartado b:

En ausencia de aletas el circuito eléctrico análogo es el indicado en la siguiente figura:


En este caso, también del calor total que genera la tarjeta de circuitos, Q , parte se disipa hacia
total
 , y parte hacia la superior, Q
la parte inferior, Q  :
i s

  
Q total  Q s  Q i

23
  Tc  T Tc  T
Q total  
R convs R tc  R cond  R convi

siendo:

1 1
R convs    250 K / W
hA 40  0,01  0,01

Por tanto:

Tc  20 Tc  20
1   Tc  174 º C
250 1  150  250

Como se puede observar, en ausencia de aletas la temperatura de la tarjeta de circuitos es

considerablemente mayor, como era de esperar ya que las aletas se colocan para disipar mejor el
calor generado por la tarjeta de chips.

Problema 2.8
Por el interior de una tubería aleteada circula agua a 60 ºC y coeficiente de transferencia de calor
por convección 2000 W/m2∙K. Exteriormente circula aire a temperatura 20 ºC y coeficiente de
transferencia de calor por convección 10 W/m2∙K. Los diámetros interno y externo de la tubería
son 11 y 12,7 mm respectivamente. Las aletas son anulares, de diámetro externo 32 mm y
espesor 1 mm. El espaciado entre aletas es tal que se colocan 235 aletas por cada metro de
longitud de tubo. Tanto el tubo como las aletas son de un cobre de conductividad térmica k = 400
W/m∙K. Calcular la pérdida de calor por metro de longitud de tubo. Despreciar la transferencia de
calor por radiación.

Solución:

Se realizará este problema aplicando analogía eléctrica. El circuito eléctrico análogo es el

siguiente:

24
siendo:

1
R convi   0,01448 K / W
D 0Lhi

ln(r1 / r0 ) ln(D1 / D 0 )
R cond    0,0005721 K / W
2 kL 2 kL

1 1
R sl    3,278 K / W
h e A sl 10  (1  235  0,001)    0,011

1
Ra 
a A a h e

Para calcular la eficiencia de aleta, ƞa, se utilizará la siguiente gráfica:

25
 siendo:

D1
r1   0,00635 m
2
D
r2  2  0,016 m
2
L  r2  r1  0,00965 m
t 0,001
r2c  r2   0,016   0,0165 m
2 2
t 0,001
L c  L   0,00965   0,01015 m
2 2
A p  L c t  10,15  10 6 m 2
h
L3c/ 2  0,05
kA p
r2c 0,0165
  2,6
r1 0,00635

h
Entrando en la gráfica con un valor de L c
3/2
 0,05 y r2c/r1 = 2,6 se obtiene que la eficiencia
kA p
de aleta resulta aproximadamente 98% (0,98 en base unitaria).

El área aleta es 235 veces la que se indica en color azul en la figura (por arriba y por abajo, que
no se percibe en la figura), teniendo en cuenta que al utilizar el modelo de eficiencia de aleta no se
tiene en consideración el extremo indicado en color negro:

A a  235  2    (r2c
2
 r12 )  0,3423 m 2

Sustituyendo valores:

1
Ra   0,2981 K / W
a A a h e

Por tanto, la pérdida de calor resulta:

  Ti  Te
Q  139 W
1
R convi  R cond 
1 1

R a R sl

26
 TEMA 3

CONDUCCIÓN DE CALOR
UNIDIMENSIONAL EN
RÉGIMEN TRANSITORIO

27
Problema 3.1
Una placa de aluminio (c = 903 J/kgꞏK, k = 236 W/mꞏK, ρ = 2702 kg/m3) de espesor 10 cm que
inicialmente se encuentra a 20 ºC se apoya sobre una superficie adiabática. El ambiente se
encuentra a 600 ºC y el coeficiente de transferencia de calor por convección es 20 W/m2ꞏK.
Calcular el tiempo que tarda la placa en alcanzar 350 ºC.

Solución:

Se procederá a comprobar si se puede aplicar el análisis simplificado. Para ello es necesario

calcular el número de Biot:

h( / A ) he
Bi    0, 00847
k k

Como el número de Biot resultó menor que 0,1 se puede aplicar el análisis simplificado, que
consiste en asumir la temperatura del sólido permanece uniforme para cada instante de tiempo.

Aplicando balance de energía a este sistema cerrado:

dE  
QW
dt

Con potencia nula:

dE  
Q W
dt

La energía total es la suma de la contribución de la energía interna, cinética, potencial y otras que
no se tratarán en el presente texto. Despreciando variaciones de energía cinética y potencial, la
energía total resulta:

E  me  m(u  v 2 / 2  gz )  mu

Por tanto, el balance de energía queda:

d(mu) 
Q
dt

Como se trata de un sólido con calor específico constante, se cumple que Δu = cΔT. Por tanto:

d(mcT ) 
Q
dt

Si la placa se está calentando, y asumiendo como positivo el calor recibido, resulta:

d(mcT )
 hA (T  T )
dt

Como la masa y el calor específico son constantes:

28
 dT
mc  hA(T  T)
dt

Indicando la masa como densidad por volumen:

dT
c  hA(T  T)
dt

Resolviendo la ecuación diferencial, resulta:

T
1 hA
t
 T  T  hA
T T  T dT  c 0 dt  ln   
 T  Ti  c
t
i

Gráficamente, la temperatura a lo largo del tiempo se comporta de la siguiente manera:

Sustituyendo valores se obtiene la solución:

 T  T  hA  T T  h  600  350  20
ln   t  ln   t  ln   t  t  10267 s
 T  Ti  c  T  Ti  ec  600  20  2702  0,1 903

Problema 3.2
Una resistencia eléctrica que tiene forma cilíndrica de diámetro 3,5 mm y longitud 10 mm genera
un calor de 0,1 W. Esta resistencia se monta junto a otros componentes de una placa electrónica
de modo que tanto ambiente como alrededores se encuentran a 30 ºC y el coeficiente de
transmisión de calor por convección es de h = 12 W/m2∙K. La emisividad superficial es de ε = 0,9.
Para simplificar los cálculos, asumir que la temperatura de la resistencia permanece uniforme para
cada instante de tiempo. Calcular la temperatura que alcanza en estado estacionario.

Solución:

En estado estacionario, el calor generado debe ser igual al disipado. En este caso se disipa por
convección y por radiación, por tanto:

  
Q gen  Q conv  Q rad


gen  hA ( T  T )  A ( T  Talr )
4 4
Q


gen  h DL ( T  T )  DL ( T  Talr )
4 4
Q

0, 1  12    0, 0035  0, 01 (T  303)  5, 67  10 8  0, 9    0, 0035  0, 01 (T 4  303 4 )

29
Nótese que la temperatura T∞ se ha indicado en kelvin. Resulta una ecuación de grado 4 con las
siguientes 4 soluciones, de las cuales la válida es la físicamente coherente, es decir, 350,5 K =
77,5 ºC.

Problema 3.3
Un cilindro de aluminio (k = 230 W/mꞏK, ρ = 2700 kg/m3, c = 1030 J/kgꞏK) de diámetro 0,05 m e
inicialmente a 25 ºC se expone a un ambiente convectivo a T∞ = 600 ºC y h = 100 W/m2ꞏK.
a) Calcular el tiempo que tarda el cilindro en alcanzar el 75% del almacenamiento máximo posible
de energía.
b) Calcular la temperatura del cilindro en ese momento.

Solución apartado a:

En primer lugar se procederá a comprobar si se puede aplicar el análisis simplificado. Para ello es
necesario calcular el número de Biot:

h( / A ) hR / 2
Bi    0, 0054
k k

Como el número de Biot resultó menor que 0,1 se puede aplicar el análisis simplificado, que
consiste en asumir la temperatura del sólido permanece uniforme para cada instante de tiempo.

La energía que almacena el cilindro en un tiempo t viene dada porla siguiente expresión:

t t
 
hA
t 
Q  Q 
 dt  hA (T  T )dt c(T  T )  1  e
 i  
c
 (J)
0 0  

La energía máxima se alcanza cuando t  ∞ en la expresión anterior, es decir:

Qmax  c(Ti  T )

El tiempo que se tarda en alcanzar el 75% del almacenamiento máximo posible de energía se
obtiene a partir de la siguiente expresión:

Q  0,75Qmax

 
hA
t
c(Ti  T )  1  e c   0,75c(Ti  T )
 
 

 
hA
t
c(Ti  T )  1  e c
  0,75 c(Ti  T )

 

hA
 t
1 e c
 0,75

2h
 t
1 e Rc
 0,75

30
 2h
 t
e Rc
 0,25

2h
 t  ln(0,25)
Rc

t  963,8 s

Solución apartado b:

Aplicando balance de energía:

dE  
QW
dt

Con potencia nula:

dE  
Q W
dt

La energía total es la suma de la contribución de la energía interna, cinética, potencial y otras que
no se tratarán en el presente texto. Despreciando variaciones de energía cinética y potencial, la
energía total resulta:

E  me  m(u  v 2 / 2  gz )  mu

Por tanto, el balance de energía queda:

d(mu) 
Q
dt

Como se trata de un sólido con calor específico constante, se cumple que Δu = cΔT. Por tanto:

d(mcT ) 
Q
dt

Si el cilindro se está calentando, y asumiendo como positivo el calor recibido, resulta:

d(mcT )
 hA (T  T )
dt

Como la masa y el calor específico son constantes:

dT
mc  hA(T  T)
dt

Indicando la masa como densidad por volumen:

dT
c  hA(T  T)
dt

31
Resolviendo la ecuación diferencial, resulta:

T
1 hA
t
 T  T  hA
T T  T dT  c 0 dt  ln   
 T  Ti  c
t
i

Gráficamente, la temperatura a lo largo del tiempo se comporta de la siguiente manera:

Sustituyendo valores se obtiene la solución:

 T  T  hA  T T  h  600  T  100
ln   t  ln   t  ln    963,8  T  107,1 ºC
 T  Ti  c  T  Ti  Rc  600  25  2700  0,05  1030

Problema 3.4
Una placa de propiedades k = 15 W/m∙K, c = 477 J/kg∙K, y ρ = 7900 kg/m3 tiene un espesor de
200 mm. Inicialmente se encuentra a 400 ºC y se apoya sobre una superficie adiabática. El
ambiente y alrededores se encuentran a 20 ºC y el coeficiente de transferencia de calor
combinado convección-radiación linealizado es de 250 W/m2∙K.
a) Calcular el tiempo que tarda la superficie inferior de la placa en alcanzar 300 ºC.
b) Calcular el tiempo que tarda la superficie inferior de la placa en alcanzar 375 ºC.

Solución apartado a:

Se procederá a comprobar si se puede aplicar el análisis simplificado. Para ello es necesario

calcular el número de Biot:

h( / A ) he 
Bi    3, 3
k k

Como el número de Biot resultó mayor que 0,1 no se puede aplicar el método simplificado, es
decir, la temperatura no se puede asumir uniforme. Por tanto, la solución viene dada por el
siguiente sumatorio:


θ*   Cn e λn Focos  λn x * 
2

n 1

32
 siendo:

x
x*  , en la superficie inferior x = 0 y por tanto x* = 0
e

T  T
* 
Ti  T

t
Fo 
R2

4 sin  λn 
Cn 
2λn  sin  2λn 

λn tg  λn   Bi

Se puede utilizar un único término del sumatorio si el número de Fourier resulta mayor que 0,2.
Como el número de Fourier depende del tiempo y es precisamente el tiempo lo que se pide
calcular en el enunciado, lo que se va a hacer es resolver el problema asumiendo que Fo > 0,2.
Una vez resuelto se comprobará si esta hipótesis resulta cierta.

Con un único término del sumatorio resulta:

θ*  C1e λ1 Fo cos  λ1x * 

2

siendo C1 = 1,22 y λ1 = 1,217 rad. Teniendo en cuenta además que x* = 0 resulta:

T  T 2 300  20 2
 C1e λ1 Focos(0)   1,22e1,217 Fo  Fo  0,334
Ti  T 400  20

Por tanto, como el número de Fourier ha resultado mayor de 0,2 el procedimiento de utilizar
solamente el primer término de la serie ha sido el correcto. A partir del número de Fourier se
obtiene el tiempo:

k
t
αt ρc
Fo  2  2  t  3359 s
e e

Solución apartado b:

Considerando un único término del sumatorio resulta:

T  T 2 375  20 2
 C1e λ1 Focos(0)   1,22e1,217 Fo  Fo  0,18
Ti  T 400  20

Como en este caso el número de Fourier resulta menor que 0,2 el error obtenido al utilizar
solamente el primer término de la serie sería relativamente considerable. Una alternativa a utilizar
el sumatorio es la siguiente gráfica, gráfica de Heisler:

33
Detalle para números de Fourier menores que 0,2:

T0  T
En este caso, los parámetros representativos son 1/Bi = h/(ek) = 0,3 y  0,93 . T0 es la
Ti  T
temperatura del plano medio de una placa de espesor 2e, que en este caso concreto sería la
temperatura de la superficie inferior (superficie adiabática) para la placa de espesor e. A partir de la
gráfica se obtiene el número de Fourier, aproximadamente 0,16, y a partir del mismo el tiempo:

αt
Fo   0,16  t  1608 s
e2

Problema 3.5
Un helado de forma cilíndrica se extrae de un congelador a una temperatura de -20 ºC. Calcular el
tiempo que tarda en empezar a derretirse (cuando algún punto empieza a alcanzar 0 ºC) en un
ambiente a 25 ºC con un coeficiente de convección de 50 W/m2∙K. Para simplificar los cálculos,
34
asumir transferencia de calor unidimensional, despreciar la transferencia de calor por radicación y
asumir el helado como un cilindro largo de hielo de radio 0,015 m. Propiedades del hielo: c = 1930
J/kg∙K, ρ = 913 kg/m3, k = 2,22 W/m∙K.

Solución:

Se procederá a comprobar si se puede aplicar el análisis simplificado. Para ello es necesario

calcular el número de Biot:

h( / A ) hR / 2
Bi    0, 17
k k

Como el número de Biot resultó mayor que 0,1 no se puede aplicar el método simplificado, es
decir, la temperatura no se puede asumir uniforme. Por tanto, la solución viene dada por el
siguiente sumatorio:


θ*   Cn e λn Fo J0  λnr * 
2

n 1

siendo:

r R
r*  , como el helado se empieza a derretir por los bordes resulta r *   1
R R

T  T
* 
Ti  T

t
Fo 
R2

2 J1  λn 
Cn 
λn J0  λn   J12  λn 
2

J1  λn 
λn  Bi
J0  λn 

Se puede utilizar un único término del sumatorio si el número de Fourier resulta mayor que 0,2.
Como el número de Fourier depende del tiempo y es precisamente el tiempo lo que se pide
calcular en el enunciado, lo que se va a hacer es resolver el problema asumiendo que Fo > 0,2.
Una vez resuelto se comprobará si esta hipótesis resulta cierta.
Con un único término del sumatorio resulta:

θ*  C1e λ1 Fo J0  λ1r * 
2

35
siendo C1 = 1,08 y λ1 = 0,79 rad. Teniendo en cuenta además que r* = 1 resulta:

T  T 0  25
 C1e λ1 Fo J0  λ1r *   1,08e 0,79 Fo J0  0,79 
2 2
  Fo  0,8
Ti  T 20  25

Por tanto, como el número de Fourier ha resultado mayor de 0,2 el procedimiento de utilizar
solamente el primer término de la serie ha sido el correcto. A partir del número de Fourier se
obtiene el tiempo:

k
t
αt ρc
Fo  2  2  t  144,3 s
R R

Problema 3.6
Un largo cilindro de diámetro 20 cm, de acero inoxidable (k = 14,9 W/mꞏK, α = 3,95ꞏ10-6 m2/s) se
retira de un horno a una temperatura uniforme de 600 ºC. El cilindro se enfría posteriormente en
un ambiente a 200 ºC con un coeficiente de transferencia de calor por convección de h = 80
W/m2∙K. Determinar la temperatura en el centro del cilindro después de 7 min de enfriamiento.

Solución:

Se procederá a comprobar si se puede aplicar el análisis simplificado. Para ello es necesario

calcular el número de Biot:

h( / A ) hR / 2
Bi    0, 26
k k

Como el número de Biot resultó mayor que 0,1 no se puede aplicar el método simplificado, es
decir, la temperatura no se puede asumir uniforme. Por tanto, la solución viene dada por el
siguiente sumatorio:


θ*   Cn e λn Fo J0  λnr * 
2

n 1

siendo:

r 0
r*  , lo cual en el centro del cilindro resulta r *   0
R R

T  T
* 
Ti  T

t 3,95  10 6  420 
Fo  2
, para este caso resulta Fo  2
 0,16
R 0,1

Téngase en cuenta de que en el cálculo del número de Fourier el tiempo se ha pasado a segundos,
es decir, 420 segundos en lugar de 7 minutos.

36
Se podría utilizar un único término del sumatorio si el número de Fourier resultara mayor que 0,2.
Como en este caso resulta menor que 0,2 el error obtenido al utilizar solamente el primer término
de la serie sería relativamente considerable. Una alternativa a utilizar el sumatorio es la siguiente
gráfica, gráfica de Heisler:

Detalle para números de Fourier menores que 0,2:


En este caso, los parámetros de la gráfica son Fo  0,16 y 1/Bi = h/(Rk) = 1,862, por tanto el
T0  T
cociente resulta aproximadamente 0,95. A partir de este valor se obtiene la temperatura
Ti  T
en el centro del cilindro:

T0  T T0  200
 0,95   0,95  T0  580 º C
Ti  T 600  200

37
Problema 3.7
Una placa de plástico de espesor 2 mm se encuentra sobre una superficie adiabática. Inicialmente
la placa presenta una temperatura uniforme de Ti = 200 ºC. Por encima de la placa circula una
corriente de aire a 25 ºC y h = 200 W/m2ꞏK. La conductividad y difusividad térmica del plástico son
k = 0,25 W/mꞏK y α = 1,2ꞏ10-7 m2/s, respectivamente. Despreciar la transferencia de calor por
radiación.
a) Representar gráficamente la temperatura adimensional de la superficie inferior, (To-T∞)/(Ti-T∞),
en función del módulo de Fourier para el intervalo 0 ≤ Fo ≤ 2. Comparar el resultado utilizando 1, 5
y 10 términos de sumatorio de cosenos de Fourier.
b) Representar gráficamente la temperatura de la superficie inferior en función del tiempo para el
intervalo de 0 a 50 segundos.

Solución apartado a:

Se procederá a comprobar si se puede aplicar el análisis simplificado. Para ello es necesario

calcular el número de Biot:

h( / A ) he
Bi    1, 6
k k

Como el número de Biot resultó mayor que 0,1 no se puede aplicar el método simplificado, es
decir, la temperatura no se puede asumir uniforme. Por tanto, la solución viene dada por el
siguiente sumatorio:


θ*   Cn e λn Fo cos  λn x * 
2

n 1

siendo:

x 0
x*  , lo cual en el la superficie inferior resulta x *   0
e e

T  T T T
*  , que en x = 0 resulta: 0*  0 
Ti  T Ti  T

t
Fo 
e2

4sen  λn 
Cn 
2λn  sen  2λn 

λn tan  λn   Bi

38
Los términos λn se denominan autovalores, gráficamente se representan en la siguiente figura:

Para un valor del número de Biot de 1,6 los 10 primeros autovalores resultan:

λ1 = 1,008 rad

λ2 = 3,564 rad

λ3 = 6,524 rad

λ4 = 9,590 rad

λ5 = 12,69 rad

λ6 = 15,81 rad

λ7 = 18,93 rad

λ8 = 22,06 rad

λ9 = 25,20 rad

λ10 = 28,33 rad

A partir de estos autovalores se obtienen los valores de Cn:

4sen  λ1 
C1   1,159
2λ1  sen  2λ1 

4sen  λ2 
C2   0,2081
2λ2  sen  2λ2 

4sen  λ3 
C3   0,07053
2λ3  sen  2λ3 

39
 4sen  λ 4 
C4   0,03375
2λ 4  sen  2λ 4 

4sen  λ5 
C5   0,01952
2λ5  sen  2λ5 

4sen  λ 6 
C6   0,01574
2λ6  sen  2λ 6 

4sen  λ7 
C7   0,008855
2λ7  sen  2λ7 

4sen  λ8 
C8   0,006535
2λ8  sen  2λ8 

4sen  λ9 
C9   0,005018
2λ9  sen  2λ9 

4sen  λ10 
C10   0,003973
2λ10  sen  2λ10 

Una vez obtenidos los primeros 10 λn y Cn ya se está en condiciones de representar gráficamente


el sumatorio θ0*   Cn e λn Fo para 1, 5 y 10 términos:
2

n 1

Como se puede observar, para números de Fourier mayores que 0,2 las soluciones utilizando 1, 5
y 10 términos son prácticamente similares. Sin embargo, para números de Fourier menores de 0,2
la solución para un término proporciona errores considerables ya que para un Fourier nulo la
temperatura adimensional debe resultar 1, como así se obtiene utilizando los primeros 5 y 10
términos de la serie.

40
Solución apartado b:

Lo que pide el enunciado es la misma gráfica que en el apartado anterior pero de forma
dimensional, es decir, en lugar de representar la temperatura adimensional, 0* , y el tiempo
adimensional, Fo, se representa directamente T0 frente a t:

Para realizar la gráfica anterior se han utilizado los 10 primeros términos del sumatorio y se ha
obtenido la temperatura y el tiempo a partir de los parámetros adimensionales:

T 0 T
0*   T0
Ti  T

t
Fo   t
e2

41
 TEMA 4

CONVECCIÓN FORZADA.
FLUJO EXTERIOR

42
Problema 4.1
Una placa electrónica de 10 ˣ 10 cm y reducido espesor se encuentra aislada por un lado,
mientras que por el otro se somete a una corriente de aire a 20 ºC y 10 m/s. Los alrededores
también se encuentran a 20 ºC. La superficie tiene una emisividad de ε = 0,60. En estado
estacionario se mide la temperatura de la placa y resulta 60 ºC, que para simplificar los cálculos
se puede asumir uniforme. Calcular el calor que está disipando la placa.

Solución:

El calor disipado será el debido a convección más radiación:

 Q
Q  Q 
rad conv

siendo:

  A(T 4  T 4 )  0,6  5,67  10 8  0,1  0,1  [(60  273)4  (20  273)4 ]  1,6 W
Q rad s 


Q conv  hL A(Ts  T )

Cálculo de hL :

  17  10 -6 m2 /s
T  T 60  20 
Se evalúan las propiedades necesarias: Tf  s   40 ºC Pr  0,72
2 2 k  0,027 W/mK


uL 10  0,1
ReL    58823  5  105  Flujo laminar
 17  10 6

NuL  0,664 Re1/L 2 Pr 1/ 3  144,34

hLL
NuL   hL  38,97 W/m2K
k

Por tanto:

Q conv  hL A(Ts  T )  15,58 W

 Q
 Q 
Q rad conv  1,67  15,58  17,25 W

43
Problema 4.2
Una placa de circuitos de 30 ˣ 30 cm disipa 50 W que, para simplificar los cálculos, se puede
asumir que disipa uniformemente. Dicha placa se enfría con aire que se encuentra a 20 ºC y a una
velocidad de 5 m/s. Despreciar pérdidas de calor por la parte inferior y despreciar también
transferencia de calor por radiación. Asumir el flujo turbulento ya que los componentes
electrónicos se suelen diseñar para actuar como turbuladores. Calcular:
a) Temperatura promedio de la placa.
b) Temperatura de la placa en el borde final.

Solución apartado a;

Para flujo turbulento sometido a flujo de calor uniforme, la expresión que proporciona la
temperatura promedio de la placa es la siguiente:

 /k
q''L
T s  T 
0,03696ReL4 / 5 Pr1/ 3

T s  T
Las propiedades necesarias se evalúan a la temperatura de película, Tf  , para la cual es
2
necesario disponer de la temperatura promedio de la superficie. Se realizará por tanto un
procedimiento iterativo.

- 1ª iteración:

En la primera iteración se evaluarán las propiedades a una temperatura de película de, por
ejemplo, 20ºC:

  15,14  10 -6 m 2 /s
T s  T 
Tf   20 º C Pr  0,7294
2 k  0,02513 W/mK


Una vez determinadas las propiedades necesarias ya se está en condiciones de determinar la

temperatura promedio de la superficie en esta primera iteración:

 /k
q''L  /k
q''L
T s  T  
0,03696ReL4 / 5 Pr 1/ 3 u L
4/5

0,03696    Pr1/ 3
  

44
 50
0,3 / 0,02513
0,3  0,3
T s  20  4/5
 5  0,3 
0,03696   6 
 0,72941/ 3
 15,14  10 

T s  20  20,08  T s  40,08 º C

Una vez obtenido un valor de temperatura promedio de superficie se procederá con la segunda
iteración.

- 2ª iteración:

  16,07  10 -6 m 2 /s
T s  T 40,08  20 
Tf    30,04 º C Pr  0,7269
2 2 k  0,02587 W/mK


Una vez determinadas las propiedades necesarias ya se está en condiciones de determinar la

temperatura promedio de la superficie en esta segunda iteración:

50
0,3 / 0,02587
0,3  0,3
T s  20  4/5
 5  0,3 
0,03696   6 
 0,72691/ 3
 16,07  10 

T s  20  20,5  T s  40,5 º C

Como se puede observar, la temperatura promedio de superficie obtenida en esta segunda

iteración es muy parecida a la obtenida en la primera iteración, por lo cual no se seguirá iterando y
se considerará como válida, aunque es importante aclarar que más iteraciones proporcionarían
una solución más precisa.

Solución apartado b:

En todo punto de la placa se verifica la ley del enfriamiento de Newton:

q ''  hx (Tsx  T )

Q 50 
Siendo el flujo de calor por unidad de área uniforme, q ''    555,5 W . Lógicamente,
A 0,3  0,3
en el borde final, x = L, también se cumple la ley del enfriamiento de Newton:

q ''  hL (TsL  T )

Por tanto, para determinar la temperatura de la placa en el borde final, TsL, es necesario calcular el
coeficiente de transferencia de calor por convección en el borde final, hL. Se calculará a partir del
número de Nusselt, el cual a su vez se calcula mediante la siguiente expresión:

4/5
 5  0,3 
NuL  0,0308ReL4 / 5 Pr1/ 3  0,0308   6 
 0,72691/ 3  262,08
 16,07  10 
45
El coeficiente de transferencia de calor por convección en el borde final resulta:

hLL kNuL 0,02587  262,08

NuL   hL    22,6 W / m2K
k L 0,3

Y, finalmente, la temperatura de la superficie en el borde final:


q ''  hL (TsL  T )  555,5  22,6  (TsL  20)  TsL  44,6 ºC

Problema 4.3
Una placa electrónica de 50 ˣ 50 cm disipa un flujo de calor uniforme de 200 W/m2. Por encima de
la placa circula una corriente de aire a 27 ºC, presión atmosférica y velocidad 2 m/s, mientras que
por debajo está aislada. La transferencia de calor por radiación es despreciable. El fabricante no
desea que en ningún punto de la placa la temperatura de operación exceda de 60 ºC con el fin de
alargar la vida útil. Indicar si se cumple el deseo del fabricante.

Solución:

T s  T
Es necesario evaluar propiedades a Tf  . La expresiones que proporcionan la
2
 /k
q''L
temperatura promedio de la superficie son T s  T  y
0,6795Re1/L 2 Pr1/ 3
 /k
q''L
T s  T  para flujo laminar y turbulento, respectivamente. Como para evaluar
0,03696ReL4 / 5 Pr1/ 3
las propiedades necesarias es necesaria la temperatura promedio se resolverá el problema
mediante un procedimiento iterativo.

1ª iteración:

En esta primera iteración se evaluarán propiedades por ejemplo a una temperatura de película de
T s  T
Tf  = 77ºC = 350 K, con lo cual las propiedades necesarian resultan:
2

  20,64  10 -6 m2 /s

T s T 
Tf   350 K Pr  0,7163
2 k  0,0293 W/mK


El número de Reynolds resulta:

u L 2  0,5
ReL    48442  5  105  Flujo laminar
 20,64  10 6

46
Como el número de Reynolds es menor de 5ꞏ105, se trata de flujo laminar y por tanto la expresión
que proporciona la temperatura de superficie promedio es la siguiente:

 /k
q''L 200  0,5 / 0,0293
T s  T    25,51  T s  52,51 º C
1/ 2
0,6795ReL Pr 1/ 3
0,6795  484421/ 2  0,71631/ 3

Este valor de temperatura de superficie promedio se utilizará como punto de partida en la segunda
iteración.

2ª iteración:

  16,98  10 -6 m2 /s
T s  T 52,51  27 
Tf    39,8 ºC  312,8 K Pr  0,7245
2 2 k  0,02659 W/mK


El número de Reynolds resulta:

uL 2  0,5
ReL   6
 58901  5  105  Flujo laminar
 16,98  10

Como el número de Reynolds sigue resultando menor de 5ꞏ105, se trata de flujo laminar y por
tanto la expresión que proporciona la temperatura de superficie promedio es también la siguiente:

 /k
q''L 200  0,5 / 0,0293
T s  T    25,39  T s  52,39 º C
1/ 2
0,6795ReL Pr 1/ 3
0,6795  589011/ 2  0,72451/ 3

Este valor de temperatura de superficie promedio es, lógicamente, más preciso que el obtenido en
la primera iteración. Si se continuaran haciendo iteraciones cada vez se obtendrían valores más
precisos. Sin embargo, dada la poca variación en el resultado obtenido entre la primera y la
segunda iteración, se adoptará este resultado como válido.

El siguiente paso es comprobar si algún punto de la placa excede 60ºC. La temperatura máxima
tiene lugar en el borde final, y se calcula a partir de la siguiente expresión:
  hL (TsL  T )
q''

Para despejar la temperatura del borde final a partir de la expresión anterior es necesario calcular
el coeficiente de transferencia de calor por convección en el borde final. Se calcula a partir del
número de Nusselt:

hLL NuLk
NuL   0,453Re x1/ 2Pr 1/3  98,74  hL   5,251 W / m2K
k L

Una vez que se tiene el valor de hL ya se está en condiciones de calcular la temperatura de

superficie en el borde final.

  hL (TsL  T )
q''  TsL  65,08 ºC

Como se puede observar, se superan los 60 ºC a los que se refiere el fabricante.

47
Problema 4.4
Un tubo largo de diámetro 10 cm y cuya temperatura superficial se mantiene a 110 ºC está
afectado por una corriente de viento a 10 ºC soplando a 8 m/s. Calcular la tasa de transferencia de
calor por unidad de longitud. Despreciar la tasa de transferencia de calor por radiación.

Solución:

La tasa de transferencia de calor será la debida a convección más radiación, despreciando

radiación el calor se disipa íntegramente por convección, y viene dado por la siguiente expresión:


Q conv  hA(Ts  T )

Cálculo de h :

  18,94  10 -6 m2 /s
T  T 
Se evalúan las propiedades necesarias: Tf  s  60 ºC Pr  0,7199
2 k  0,02807 W/mK


uD
PeD  ReD Pr  Pr  42194  0,7199  30375,5


Como el número de Péclèt resultó mayor de 0,2 se empleará la expresión de Churchill:

4/5
0,62ReD1/ 2 Pr1/ 3   ReD  
5/8
hD
NuD   0,3  1/ 4 
1     101,9
k 1   0,4 / Pr 2 / 3    282000  
 
hD NuDk
NuD   h  28,62 W / m2K
k D

Una vez que se tiene el coeficiente de transferencia de calor se está en condiciones de calcular la
tasa de transferencia de calor por unidad de longitud, la cual viene dada por:


Q conv  hA(Ts  T )  hDL(Ts  T )  q'conv  hD(Ts  T )  898,6 W / m

Problema 4.5
Por el interior de una tubería larga de acero (kacero = 80 W/m∙K) de diámetros interior y exterior 5 y
5,5 cm respectivamente fluyen 0,2 kg/s de agua líquida a una temperatura de 15 ºC, la cual para
simplificar los cálculos se puede considerar constante a lo largo de toda la longitud. La tubería se
cubre con un material aislante (kaislante = 0,05 W/m∙K) de 3 cm de espesor. El aire exterior se
encuentra a 25 ºC y circulando a 10 m/s debido al viento. Despreciar la transferencia de calor por
radiación. Calcular la temperatura de la superficie exterior del material aislante.

48
Solución:

Se procederá a resolver este problema mediante analogía eléctrica. Despreciando la transferencia

de calor por radiación, el circuito eléctrico análogo se muestra en la siguiente figura:

Los radios r1, r2 y r3 son los indicados en la siguiente figura, siendo r1 = 0,025 m, r2 = 0,0275 m y r3
= r2 + e = 0,0575 m:

La temperatura de la superficie del aislante, Ts, se puede obtener a partir de la siguiente expresión:

Ts  Ti T  Ts
 e
1 ln(r2 / r1 ) ln(r3 / r2 ) 1
 
2r1hi 2k acero 2k aislante 2r3he

Para despejar Ts de la expresión anterior es necesario calcular hi y he, lo cual se hace a

continuación.

Cálculo de hi:

 agua  0,002243 Pa  s

Propiedades de agua a 15ºC Pragua  8,29

k agua  0,58 W/mK

Para comprobar si el agua presenta flujo laminar o turbulento se calcula el número de Reynolds:

umaguaD1  agua
4m
ReDagua    13378  2300  flujo turbulento
 agua D1 agua

Para calcular el número de Nusselt se utilizará por ejemplo la correlación de Dittus-Boelter:

49
 hiD1 k aguaNuDagua
NuDagua   hi   1237 W / m 2K
k agua D1

Cálculo de he:

Para evaluar las propiedades necesarias del aire es necesaria la temperatura de la superficie del
aislante pero es precisamente esa temperatura la que se pide calcular. Se seguirá por tanto un
procedimiento iterativo.

- 1ª iteración:

En la primera iteración se asumirá por ejemplo una temperatura de película, Tf = (T∞e + Ts)/2, de
20ºC (lógicamente, esta temperatura debe ser mayor de 15ºC y menor de 25ºC). A esta
temperatura las propiedades necesarias son las siguientes:

 aire  15,14  10 -6 m2 /s
T e  Ts 
Propiedades de aire Tf   20 ºC Praire  0,7293
2 k  0,02513 W/mK
 aire

El número de Reynolds para el aire resulta:

uaireD3
ReDaire   75957
 aire

El número de Péclèt para el aire resulta:

uaireD3
P eDaire  R eDaire Praire  Praire  55396
 aire

A la vista del número de Péclèt, una correlación apropiada para determinar el número de Nusselt
para el aire es la de Churchill, dada por:

4/5
0,62Re1/Daire
2 1/ 3
Praire   ReDaire 5 / 8 
NuDaire  0,3  1/ 4 1      181,5
1   0,4 / Praire 2 / 3    282000  
 

A partir del número de Nusselt para el aire se obtiene el coeficiente de transferencia de calor por
convección:

heD3 k aire NuDaire

NuDaire   he   39,7 W / m2K
k aire D3

La temperatura de la superficie del aislante resulta:

Ts  Ti T  Ts
 e  Ts  24,7 º C
1 ln(r2 / r1 ) ln(r3 / r2 ) 1
 
2r1hi 2k acero 2k aislante 2r3he

50
- 2ª iteración:

En esta segunda iteración se utilizará el valor de Ts obtenido en la primera iteración, es decir, Ts =

24,71 ºC. Las propiedades del aire teniendo en cuenta este nuevo valor resultan:

 aire  15,59  10 -6 m2 /s
T e  Ts 25  24,7 
Propiedades de aire Tf    24,85 ºC Praire  0,728
2 2 k  0,0255 W/mK
 aire

El número de Reynolds para el aire resulta:

uaireD3
ReDaire   73779
 aire

El número de Péclèt para el aire resulta:

uaireD3
P eDaire  R eDaire Praire  Praire  53721
 aire

El número de Nusselt para el aire resulta:

4/5
0,62Re1/Daire
2 1/ 3
Praire   ReDaire 5 / 8 
NuDaire  0,3  1/ 4 1      178
1   0, 4 / Praire 2 / 3    282000  
 

A partir del número de Nusselt para el aire se obtiene el coeficiente de transferencia de calor por
convección:

heD3 k aire NuDaire

NuDaire   he   39,44 W / m2K
k aire D3

La temperatura de la superficie del aislante resulta:

Ts  Ti T  Ts
 e  Ts  24,7 º C
1 ln(r2 / r1 ) ln(r3 / r2 ) 1
 
2r1hi 2k acero 2k aislante 2r3he

Como se puede observar, en esta segunda iteración se ha obtenido una temperatura de superficie
de aislante igual que en la primera iteración, por lo que se considerará el valor obtenido como
válido y ya no es necesario seguir iterando.

Problema 4.6
Una línea de alta tensión de 25 mm de diámetro tiene una resistencia eléctrica de 10-4 Ω/m y su
temperatura superficial en estado estacionario es prácticamente uniforme e igual a 30 ºC. En flujo
cruzado sobre la línea circula aire ambiental a 10 ºC y 5 m/s. Despreciar la transferencia de calor
por radiación.
a) Calcular la corriente que está transmitiendo.
b) Si la línea se aproxima como una varilla sólida de cobre de conductividad térmica k = 396,1
W/mꞏK, ¿cuál es la temperatura del eje?
NOTA: Efecto Joule: q 'Joule  I2R ' (W/m).

51
Solución apartado a:

El calor generado se disipa por convección y radiación de la forma:

q 'Joule  q conv
'
 q rad
'

I2R '  hD(Ts  T )

Cálculo de h :

  15,14  10 -6 m2 /s
Ts  T 
Se evalúan las propiedades necesarias: Tf   20 ºC Pr  0,7293
2 k  0,02513 W/mK


uD
PeD  ReD Pr  Pr  8247  0,793  6021,3


Como el número de Péclèt resultó mayor de 0,2 se empleará la expresión de Churchill:

4/5
0,62ReD1/ 2 Pr 1/ 3   ReD  
5/8
hD
NuD   0,3  1/ 4 
1     48,76
k 1   0,4 / Pr 2 / 3    282000  
 
hD NuDk
NuD   h  49,03 W / m2K
k D

Una vez que se tiene el coeficiente de transferencia de calor la intensidad se calcula despejándola
de la siguiente ecuación:

I2R '  hD(Ts  T )  I  877,6 A

Solución apartado b:

Resolviendo la ecuación general de conducción de calor:

q '''gen T
 T 2

c t

En coordenadas cilíndricas:

 1   T  1   T   2 T  q gen
'''
T
 r  2   2  
 r r  r  r     z  c t

52
En estado estacionario y con conducción de calor solamente a lo largo de la coordenada radial
queda:

 1   T  1   T   2 T  q '''gen T
 r  2   2 
 
 r r  r  r     z  c t

'''
1 d  dT  q gen
 r  0
r dr  dr  c

'''
d  dT  q gen
r  r
dr  dr  k

Resolviendo resulta:

q '''gen
T r 2  C1 ln(r)  C 2
4k cobre

Las constantes de integración C1 y C2 se determinan aplicando condiciones de contorno:

dT 
r 0: 0  C1  0
dr 
C 2  30,02 º C
r  R : T  30 º C 

Sustituyendo valores, la distribución de temperatura viene dada por la siguiente expresión:

T  99,02r  30,02

En el eje, r = 0, la temperatura resulta 30,02 ºC.

Problema 4.7
Una esfera de 10 mm de diámetro se encuentra a 75 ºC. La esfera está inmersa en una corriente
de aire a presión atmosférica, temperatura 23 ºC y velocidad 10 m/s. Despreciar la tasa de
transferencia de calor por radiación. Calcular:
a) Coeficiente de transferencia de calor por convección.
b) Tasa de transferencia de calor.
DATO: Área de una esfera: 4πR2.

Solución apartado a:

53
Se utilizará la siguiente correlación:

1/ 4
hD  
NuD   
 2  0,4ReD1/ 2  0,06ReD 2 / 3 Pr0,4   
k  s 
si 0,71  Pr  3,2; 3,5  ReD  7,6  104

Previamente es necesario verificar que se cumpla el rango de aplicación. Pero antes de comprobar
esto se evalúan las propiedades que van a ser necesarias:

Ts  75 ºC   s  20,73  10 -6 Pa  s

   18,39  10 -6 Pa  s

  15, 42  10 m /s
-6 2
T  23 ºC  
k   0,02535 W/mK
Pr  0,7286
 

uD 10  0,01
ReD    6485,1
 15, 42  10 6

Se cumple el rango de aplicación de la correlación, por tanto se puede aplicar y el número de

Nusselt resulta:

1/ 4
hD  
NuD   
 2  0,4ReD1/ 2  0,06ReD 2 / 3 Pr0,4     47,38
k  s 

hD NuDk 
NuD   h  120,11 W / m2K
k D

Solución apartado b:

Despreciando la tasa de transferencia de calor por radiación, el calor total será el disipado por
convección:


conv  hA(Ts  T )  h4 R (Ts  T )  1,96 W
2
Q

54
 TEMA 5

CONVECCIÓN FORZADA.
FLUJO INTERIOR

55
Problema 5.1
Por un tubo circular de 60 mm de diámetro y longitud 7 m fluyen 0,01 kg/s de agua líquida, la cual
sale del tubo a 80 ºC. Debido a la acción solar, la superficie del tubo está sometida a un flujo de
calor de 2000 W/m2, el cual, para simplificar los cálculos, asumir uniforme. Calcular la
temperatura de la superficie a la salida del tubo.

Solución:

  355,1  10 -6 Ns/m2

Se evalúan las propiedades necesarias: Tm  80 ºC Pr  2,27
k  0,6561 W/mK


A continuación se calcula el número de Reynolds para comprobar si el flujo es laminar o turbulento:

umD 4m 
ReD    597,9  2300  flujo laminar
 D

A continuación se calcula el inverso del número de Graetz para comprobar si el flujo es

desarrollado o no desarrollado:

x /D
Gz 1   0,086  0,05  flujo desarrollado
ReD Pr

Como el flujo es laminar desarrollado y el flujo de calor se puede asumir uniforme, el número de
Nusselt resulta 4,36. A partir del mismo se obtiene el coeficiente de transferencia de calor a la
salida:

h sD NuDk
NuD   hs   47,68 W/m2K
k D

Por tanto, la temperatura de la superficie a la salida del tubo resulta:

q ''s  hs (Tss  Tms )  2000  47,68(Tss  80)  Tss  122 º C

Problema 5.2
Un aceite de propiedades k = 0,1443 W/m∙K; ν = 2,03∙10-3 m2/s y Pr = 23182 fluye a una velocidad
de 2 m/s por el interior de una tubería de 30 cm de diámetro. La tubería tiene una longitud de 200
mm y su superficie se mantiene a una temperatura prácticamente uniforme de 0 ºC. Calcular el
coeficiente de transferencia de calor para el aceite que circula por el interior de la tubería.

Solución:

Primeramente se procederá a calcular el número de Reynolds para ver si se trata de flujo laminar
o turbulento:

umD 2  0,3
ReD    294,8  2300  flujo laminar
 2,03  10 3

Dado que se trata de flujo laminar, es necesario si se trata de flujo desarrollado o no desarrollado.
Se calculará la inversa del número de Graetz a la salida del tubo, correspondiente a x = 0,2:

56
 x /D 0,2 / 0,3
Gz 1    9,75  10 8  0,05  flujo no desarrollado
ReD Pr 294,8  23182

Resulta por tanto flujo no desarrollado a la salida del tubo y por tanto para cualquier punto del tubo,
que tendría un valor de x todavía menor. Para flujo laminar no desarrollado con un número de
Prandtl tan elevado se utilizará la siguiente expresión para determinar el número de Nusselt
promedio desde la entrada hasta la salida del tubo:

hD 0,0668 D / L  ReD Pr
NuD   3,66   366,6
1  0,04 D / L  ReD Pr 
2/3
k

A partir del número de Nusselt promedio se obtiene el coeficiente de transferencia de calor

promedio desde la entrada hasta la salida del tubo:

hD kNuD
NuD   h  176,3 W / m2K
k D

Problema 5.3
Un tubo circular de 50 mm de diámetro está sometido a un flujo de calor constante 10000 W/m2.
Por el interior del tubo fluyen 0,5 kg/s de un líquido de propiedades k = 0,5 W/mꞏK, μ = 500ꞏ10-6
Nꞏs/m2, Pr = 3 y c = 4100 J/kgꞏK, el cual entra a 15 ºC y sale a 60 ºC. La pared externa del tubo
está bien aislada, de forma que todo este calor pasa al fluido. Calcular:
a) Longitud de tubo.
b) Temperatura de la superficie del tubo a la salida.

Solución apartado a:

Aplicando balance de energía al sistema abierto que constituye el conducto:

dE 
dt
QW e  e (he  v e2 / 2  gze )  s m s (hs  v s2 / 2  gzs )
  m

En estado estacionario, con potencia nula y despreciando variaciones de energía cinética y

potencial:

dE
dt
 W
Q 
  m
e
 e (he  v 2e / 2  gz e )   m
s
 s (hs  v s2 / 2  gz s )

  m(h
0Q  e  hs )

57
Como el calor específico es constante, la diferencia de entalpías resulta he – hs = c(Tme - Tms). Por
tanto:

  mc
0Q  p (Tme  Tms )

0  q ''s πDL  mc
 p (Tme  Tms )  L  58,7 m

Solución apartado b:

En todo punto se cumple que q  h(Ts  Tfluido ) . Particularmente, a la salida también se cumple:
''

q ''  hs (Tss  Tms )

Para calcular hs es necesario primeramente calcular el número de Reynolds para comprobar si el

flujo es laminar o turbulento:

umD 4m 
ReD    25477  2300  flujo turbulento
 D

Aplicando la correlación de Dittus-Boelter:

NuD  0,023ReD4 / 5 Pr n  119,53

Despejando h a partir del número de Nusselt resulta:

hD NuDk
NuD   h  1195,3 W / m 2K
k D

Por tanto, la temperatura de la superficie del tubo a la salida resulta:

q ''  hs (Tss  Tms )  Tss  68,36 ºC

Problema 5.4
0,05 kg/s de aire caliente fluyen por un conducto circular de diámetro 0,15 m y 5 m de largo. El
conducto se encuentra a una temperatura menor que la del aire, por lo que éste se enfría y
abandona el tubo a 77 ºC. Calcular el coeficiente de transferencia de calor por convección a la
salida del tubo.

Solución:

  20,82  10 -6 Ns/m2

Se evalúan las propiedades necesarias: Tm  77 ºC Pr  0,7163
k  0,0293 W/mK


A continuación se calcula el número de Reynolds para comprobar si el flujo es laminar o turbulento:

umD 4m 
ReD    20395  2300  flujo turbulento
 D

58
Al tratarse de flujo turbulento no es cecesario comprobar si es desarrollado o no desarrollado
porque en flujo turbulento la región de entrada es despreciable. Son de aplicación las correlaciones
de Dittus-Boelter, Sieder-Tate, Petukhov, Gnielinski, etc pero con los datos que aporta el enunciado
únicamente se puede utilizar la correlación de Dittus-Boelter:

NuD  0,023ReD4 / 5 Pr n  58,33

Nótese que el coeficiente n se ha tomado como 0,3 por tratarse de un enfriamiento. A partir del
número de Nusselt, el coeficiente de transferencia de calor a la salida del tubo resulta:

hD NuDk
NuD   h  11,38 W/m2K
k D

Problema 5.5
Un flujo de 0,15 kg/s de etileno glicol (cp = 2562 J/kg∙K, µ = 0,00522 N∙s/m2, k = 0,26 W/m∙K, Pr =
51,3) fluye por el interior de un tubo de pared delgada y diámetro 3 mm. El tubo se enrolla y se
sumerge en un baño de agua que mantiene la superficie del tubo a una temperatura
prácticamente constante de 10 ºC. El etileno glicol entra al tubo a 85 ºC y sale a 35 ºC. Despreciar
la transferencia de calor por radiación. Calcular la longitud del tubo.

Solución:

En este caso no es necesario consultar las propiedades en tablas puesto que se proporcionan en
el enunciado. A continuación se calcula el número de Reynolds para comprobar si el flujo es
laminar o turbulento:

umD 4m 
ReD    12202  2300  flujo turbulento
 D

Al tratarse de flujo turbulento no es cecesario comprobar si es desarrollado o no desarrollado

porque en flujo turbulento la región de entrada es despreciable. Son de aplicación las correlaciones
de Dittus-Boelter, Sieder-Tate, Petukhov, Gnielinski, etc pero con los datos que aporta el enunciado
únicamente se puede utilizar la correlación de Dittus-Boelter:

NuD  0,023ReD4 / 5 Pr n  139,28

Nótese que el coeficiente n se ha tomado como 0,3 por tratarse de un enfriamiento. A partir del
número de Nusselt, el coeficiente de transferencia de calor resulta:

hD NuDk
NuD   h  12071, 4 W/m 2K
k D

59
La longitud de tubo se calcula a partir de la siguiente expresión:

hD 12071,40,003
Ts  Tms 

mc
L
10  35  L
e p  e 0,152562
 L  3,71 m
Ts  Tme 10  85

Otra forma de resolver el problema es de la siguiente manera:

  mc
Q  p (Tms  Tme )  19215 W

  hAT  hDLT
Q  L  3,71 m
ml ml

(Ts  Tms )  (Ts  Tme )

siendo: Tml   45,51
 Ts  Tms 
ln  
 Ts  Tme 

Problema 5.6
Por el interior de un conducto rectangular de 1,5 m de largo y sección 50 ˣ 100 mm circulan 0,15
kg/s de aire. La superficie del conducto se mantiene a una temperatura prácticamente uniforme de
150 ºC y el aire abandona el mismo a 40 ºC y presión atmosférica. Calcular el coeficiente de
transferencia de calor a la salida del conducto.

Solución:

El primer paso consiste en evaluar las propiedades necesarias. Como el enunciado pide el
coeficiente de transferencia de calor a la salida de tubo, las propiedades del aire se evaluarán a la
salida del tubo, cuya temperatura se indica que es de 40ºC:

  19,17  10 -6 Ns/m2

Tm  77 ºC Pr  0,7245
k  0,0266 W/mK


Como no se trata de un conducto circular es necesario calcular el diámetro hidráulico:

4A t 4  0,05  0,1
Dh    0,0033 m
P 2  (0,05  0,1)

Para comprobar si el flujo es laminar o turbulento se calcula el número de Reynolds relativo al

diámetro hidráulico:

 h
umDh mD 0,15  0,0033
ReD     5164  2300  flujo turbulento
 A t  0,05  0,1  19,17  10 6

Al tratarse de flujo turbulento se pueden aplicar las correlaciones correspondientes a conductos

circulares. En flujo turbulento la región de entrada es despreciable y son de aplicación las
correlaciones de Dittus-Boelter, Sieder-Tate, Petukhov, Gnielinski, etc pero con los datos que
aporta el enunciado únicamente se puede utilizar la correlación de Dittus-Boelter:

NuD  0,023ReD4 / 5 Pr n  19, 49

60
Nótese que el coeficiente n se ha tomado como 0,4 por tratarse de un calentamiento. A partir del
número de Nusselt, el coeficiente de transferencia de calor a la salida del tubo resulta:

hDh NuDk
NuD   h  157,16 W/m2K
k Dh

61
 TEMA 6

CONVECCIÓN LIBRE

62
Problema 6.1
Una placa delgada de 0,5 ˣ 0,5 m se mantiene a una temperatura prácticamente uniforme de 80 ºC.
La parte inferior de la placa está perfectamente aislada, mientras que la parte superior se expone a
aire en reposo a 25 ºC. La emisividad superficial es unitaria y los alrededores también se
encuentran a 25 ºC. Calcular la transferencia de calor al ambiente.

Solución:

La transferencia de calor se debe a convección y a radiación:

 Q
Q  Q 
rad conv

El calor disipado por radiación viene dado por:

  A(T 4  T 4 )  5,67  108  1 0,5  0,5  (3534  2984 )  112,7 W

Qrad s alr

El calor disipado por convección viene dado por:


Q conv  hA(Ts  T )

Es necesario calcular el coeficiente de transferencia de calor por convección, h. Se trata de la

región superior de una superficie caliente. Es necesario calcular la longitud característica mediante
la siguiente expresión:

A s L2 L
LC     0,125 m
P 4L 4

Las propiedades necesarias se evalúan a la temperatura de película:

  18,2  10 -6 m2 / s

Pr  0,7216
Ts  T 25  80 
Tf    52,5 ºC k  0,02753 W / mK
2 2  
 25,22  10 -6 m2 / s
  1/ Tf  1/ 325,5  0,0031 K 1

El número de Rayleigh resulta:

g Ts  T L3C 9,81  0,0031  (80  25)  0,1253

RaLC   6 6
 7,1  10 6
 25,22  10  18,2  10

A la vista del numero de Rayleigh obtenido, se puede aplicar la siguiente correlación para el
número de Nusselt:

NuLC  0,54Ra1/LC4  27,9

A partir del número de Nusselt se obtiene el coeficiente de transferencia de calor por convección:
hL kNuLC
NuLC  C  h   6,14 W / m2K
k LC

63
Y el calor disipado por convección resulta:


Q conv  hA(Ts  T )  6,14  0,5  0,5  (80  25)  84,5 W

Finalmente, el calor total resulta:

  112,7  84,5  197,2 W

Q

Problema 6.2
Un flujo de aire caliente que circula a través de un largo conducto rectangular de sección 0,2 ˣ 0,6
m mantiene la superficie externa del conducto a una temperatura prácticamente uniforme de 40ºC.
El conducto se encuentra rodeado de aire en reposo a 20 ºC y alrededores también a 20 ºC. La
emisividad superficial es de ε = 0,6. Calcular la pérdida de calor por unidad de longitud.

Solución:

El flujo de calor disipado por el conducto es mediante radiación y convección:

q '  q rad
'
 q conv
'

El calor disipado por radiación, por unidad de longitud, viene dado por:

A
'
q rad   (Ts4  Talr
4
)  5,67  10 8  0,6  (2  0,2  2  0,6)  (313 4  293 4 )  121,2 W / m
L

El calor disipado por convección, en la parte externa del conducto es debido a convección libre. Se
debe a la aportación de las dos superficies laterales más la superior más la inferior, tal y como se
indica en la siguiente expresión:

q 'conv  2q lat

'
 q 'sup  q inf
'
 2  0,2  hlat (Ts  T )  0,6  hsup (Ts  T )  0,6  hinf (Ts  T )

Es necesario calcular los coeficientes de transmisión de calor por convección hlat , hsup y hinf . El
procedimiento se indica a continuación.

64
   16,06  10 -6 m2 / s

Pr  0,7269
Ts  T 
Evaluación de propiedades a Tf   30 ºC k  0,02587 W / mK
2  
 22,1  10 -6 m2 / s
  1/ Tf  1/ 303  0,0033 K 1

Superficies laterales:

Se trata de superficies verticales. Se comprobará si se trata de flujo laminar o turbulento mediante

el número de Rayleigh:

g Ts  T L3lat 9,81  0,0033  (40  20)  0,23

RaLlat    1, 46  107  109  flujo laminar
 22,1  10 6  16,06  10 6

Nótese que en la expresión anterior se ha puesto como longitud Llat la dimensión vertical 0,2 m. Al
tratarse de flujo laminar, la siguiente expresión es apropiada para calcular el número de Nusselt
para las superficies laterales:

0,67Ra1/Llat4
NuLlat  0,68  4/9
 32,5
1   0, 492 / Pr 9 /16 
 

A partir del Número de Nusselt se obtiene el coeficiente de transferencia de calor para las
superficies laterales:

h lat L kNuLlat 0,02587  32,5

NuLlat   h lat    4,2 W / m2K
k L 0,2

Superficie superior:

Se trata de analizar la región superior de una superficie caliente. Es necesario determinar la

longitud característica, la cual viene dada por:

As L  0,6 L  0,6
Lc     0,3 m
P 2  (L  0,6) 2 L

En la expresión anterior se ha considerado que la longitud del tubo es lo suficientemente elevada

como para considerar que L + 0,6 ≈ L.
El número de Rayleigh para la superficie superior resulta:

g Ts  T L3c 9,81  0,0033  (40  20)  0,33

RaLc sup    4,92  107
 22,1  10 6  16,06  10 6

Para el número de Rayleigh obtenido, se utilizará la siguiente correlación para calcular el número
de Nusselt para la superficie superior:

NuLC sup  0,15Ra1/LC3  54,95 W / m2K

65
A partir del Número de Nusselt se obtiene el coeficiente de transferencia de calor para la superficie
superior:
h L kNuLC sup
NuLC sup  sup C  h sup   4,74 W / m2K
k LC

Superficie inferior:

Se trata de analizar la región inferior de una superficie caliente. La longitud característica resulta la
misma que la que se ha calculado para la superficie superior, es decir, 0,3 m.

El número de Rayleigh también resulta el mismo que para la superficie superior, es decir, 4,92∙107.
Para este valor, se utilizará la siguiente correlación para calcular el número de Nusselt para la
superficie inferior:

NuLC inf  0,27Ra1/LC4  22,62

A partir del Número de Nusselt se obtiene el coeficiente de transferencia de calor para la superficie
inferior:

h inf L C kNuLC inf

NuLC inf   h inf   1,95 W / m2K
k LC

Por tanto, una vez calculados los coeficientes de transferencia de calor por convección para las
superficies laterales, superior e inferior, ya se está en condiciones de calcular el calor disipado por
convección, por unidad de longitud, el cual resulta:

'
q conv  2q lat
'
 q 'sup  q inf
'
 2  0,2  hlat (Ts  T )  0,8  hsup (Ts  T )  0,8  hinf (Ts  T ) 

 2  0,2  4,2  (40  20)  0,6  4,74  (40  20)  0,6  1,95  (40  20)  113,9 W / m

Y, finalmente, el calor total disipado por unidad de longitud resulta:

q '  q rad
'
 q conv
'
 151,6  125,7  265,5 W / m

Problema 6.3
La placa de la figura, de 0,8 ˣ 0,5 m, se mantiene a una temperatura prácticamente uniforme de 60
ºC. Tanto por la parte superior como inferior, la placa se expone a aire en reposo a 20ºC. Los
alrededores también se encuentran a 20 ºC. La emisividad superficial es de ε = 0,9. Calcular la
tasa de transferencia de calor por metro de longitud.

66
Solución:

El flujo de calor disipado por el conducto es mediante radiación y convección:

q '  q rad
'
 q conv
'

El calor disipado por radiación, por unidad de longitud, viene dado por:

A
'
q rad   (Ts4  Talr
4
)  5,67  10 8  0,6  (2  0,5  2  0,8)  (333 4  293 4 )  653,6 W / m
L

El calor disipado por convección, en la parte externa del conducto es debido a convección libre. Se
debe a la aportación de las dos superficies laterales más la superior más la inferior, tal y como se
indica en la siguiente expresión:

q 'conv  2q lat

'
 q 'sup  q inf
'
 2  0,5  hlat (Ts  T )  0,8  hsup (Ts  T )  0,8  hinf (Ts  T )

Es necesario calcular los coeficientes de transmisión de calor por convección hlat , hsup y hinf . El
procedimiento se indica a continuación.

  17  10 -6 m2 / s

Pr  0,7245
Ts  T 
Evaluación de propiedades a Tf   313 K k  0,02661 W / mK
2  
 23, 47  10 -6 m 2 / s
  1/ Tf  1/ 313  0,0032 K 1

Superficies laterales:

Se trata de superficies verticales. Se comprobará si se trata de flujo laminar o turbulento mediante

el número de Rayleigh:

g Ts  T L3lat 9,81  0,0032  (60  20)  0,53

RaLlat    3,93  108  109  flujo laminar
 17  10 6  23, 47  10 6

Nótese que en la expresión anterior se ha puesto como longitud Llat la dimensión vertical 0,5 m. Al
tratarse de flujo laminar, la siguiente expresión es apropiada para calcular el número de Nusselt
para las superficies laterales:

0,67Ra1/Llat4
NuLlat  0,68  4/9
 73,27
1   0, 492 / Pr 9 /16 
 

A partir del Número de Nusselt se obtiene el coeficiente de transferencia de calor para las
superficies laterales:

h lat L kNuLlat 0,02661  73,27

NuLlat   h lat    3,9 W / m2K
k L 0,5

67
Superficie superior:

Se trata de analizar la región superior de una superficie caliente. Es necesario determinar la

longitud característica, la cual viene dada por:

As L  0,8 L  0,8
Lc     0, 4 m
P 2  (L  0,8) 2 L

En la expresión anterior se ha considerado que la longitud del tubo es lo suficientemente elevada

como para considerar que L + 0,8 ≈ L.
El número de Rayleigh para la superficie superior resulta:

g Ts  T L3c 9,81  0,0032  (60  20)  0,83

RaLc sup    2,01  108
 17  10 6  23, 47  10 6

Para el número de Rayleigh obtenido, se utilizará la siguiente correlación para calcular el número
de Nusselt para la superficie superior:

NuLC sup  0,15Ra1/LC3  87,92 W / m2K

A partir del Número de Nusselt se obtiene el coeficiente de transferencia de calor para la superficie
superior:

h sup L C kNuLC sup

NuLC sup   h sup   5,85 W / m2K
k LC

Superficie inferior:

Se trata de analizar la región inferior de una superficie caliente. La longitud característica resulta la
misma que la que se ha calculado para la superficie superior, es decir, 0,4 m.

El número de Rayleigh también resulta el mismo que para la superficie superior, es decir, 2,01∙108.
Para este valor, se utilizará la siguiente correlación para calcular el número de Nusselt para la
superficie inferior:

NuLC inf  0,27Ra1/LC4  32,16

A partir del Número de Nusselt se obtiene el coeficiente de transferencia de calor para la superficie
inferior:

h inf L C kNuLC inf

NuLC inf   h inf   2,14 W / m2K
k LC

Por tanto, una vez calculados los coeficientes de transferencia de calor por convección para las
superficies laterales, superior e inferior, ya se está en condiciones de calcular el calor disipado por
convección, por unidad de longitud, el cual resulta:

q 'conv  2q lat

'
 q 'sup  q inf
'
 2  0,5  hlat (Ts  T )  0,8  hsup (Ts  T )  0,8  hinf (Ts  T ) 

 2  0,5  3,9  (60  20)  0,8  5,85  (60  20)  0,8  2,14  (60  20)  569,2 W / m

68
Y, finalmente, el calor total disipado por unidad de longitud resulta:

q '  q rad
'
 q conv
'
 653,6  569,2  1222,8 W / m

Problema 6.4
En una bombilla incandescente de 100 W se ha medido la temperatura superficial resultando 180
ºC. Asumir que esta temperatura superficial es prácticamente uniforme y que la bombilla es
esférica, con un diámetro de 10 cm. La emisividad de la superficie del cristal es 0,85 y el aire
ambiente se encuentra en reposo y a una temperatura de 20 ºC. Determinar el porcentaje de
energía que se pierde en forma de calor.
DATO: Área de la superficie de una esfera: 4πR2.

Solución:

La energía que se pierde en forma de calor es el disipado por radiación y por convección:

 Q
 Q 
conv  A(Ts  Talr )  hD A(Ts  T )   4 R (Ts  Talr )  hD 4 R (Ts  T )
4 4 2 4 4 2
Q rad

Es necesario calcular el coeficiente de transferencia de calor por convección, hD :

  23,04  10 -6 m2 / s

Pr  0,7199
Ts  T 
Se evalúan las propiedades necesarias: Tf   100 ºC k  0,031 W / mK
2  
 32,36  10 -6 m2 / s
  1/ Tf  1/ 373  0,002681 K 1

g Ts  T D3
Número de Rayleigh: RaD   5,643  10 6


El número de Nusselt se obtiene a partir de la siguiente correlación:

0,589Ra1/D 4
NuD  2  4/9
 22,76
1   0, 469 / Pr 9 /16 
 

A partir del número de Nusselt se obtiene el coeficiente de transferencia de calor por convección:

hD D kNuD
NuD   hD   7,041 W / m2K
k D

69
Por tanto, la energía que se pierde en forma de calor resulta:

 Q
 Q 
conv   4 R (Ts  Talr )  hD 4 R (Ts  T )  52,57  35,37  87,94 W
2 4 4 2
Q rad

Si la energía que se pierde por calor resulta 87,94 W, esto supone un 87,94% teniendo en cuenta
que la bombilla es de 100 W.

Es interesante notar de las importantes pérdidas de calor que presentan este tipo de bombillas.
Una bombilla incandescente produce luz mediante el calentamiento por efecto Joule de un
filamento de tungsteno hasta ponerlo al rojo blanco. Debido a su baja eficiencia este tipo de
bombillas están actualmente prohibidas ya que, por ejemplo en el caso de la bombilla de este
problema, un 87,94% de la electricidad se transforma en calor mientras que solamente el restante,
es decir, un 12,06%, se transforma en luz.

Problema 6.5
La condensación del vapor de agua del aire sobre una superficie se produce si la temperatura de
dicha superficie es inferior a la temperatura de rocío del aire. En la superficie interior de los
parabrisas de los automóviles la condensación puede evitarse impulsando aire caliente
paralelamente a dicha superficie. Calcular la temperatura mínima a la que debe de impulsarse aire
caliente a 6 m/s para desempañar el parabrisas de un automóvil que circula a 72 km/h cuando la
temperatura exterior es de 0 ºC y la de rocío se estima en 10 ºC para la humedad existente bajo
las condiciones que se desean analizar. Tomar como longitud del parabrisas 60 cm, espesor 6
mm y conductividad térmica del vidrio 0,95 W/m∙K. Asumir que tanto el aire exterior como el
interior fluyen en dirección paralela al parabrisas. Despreciar la transferencia de calor por
radiación y asumir que en el exterior del cristal no se produce condensación.

Solución:

Haciendo un balance de energía resulta:

q ''convi  q radi
''
 q cond
''
 q conve
''
 q rade
''

Para simplificar los cálculos, el enunciado sugiere despreciar la transferencia de calor por
radiación, por tanto:

q ''convi  q radi
''
 q ''cond  q conve
''
 q rade
''

70
 Tsi  Tse
hi (Ti  Tsi )  k  hi (Ti  Tsi )
e

Para simplificar los cálculos, las temperaturas interior y exterior del cristal, Tsi y Tse, se asumirán
uniformes. Lo que indica el problema es que Tsi ≥ 10 ºC. Se para el caso límite de Tsi = 10 ºC.

Son necesarias las temperaturas T∞i y Tse indicadas en la figura para evaluar algunas de las
propiedades necesarias para resolver el problema pero dichas temperaturas no se pueden
determinar sin antes resolver el problema. Se realizará por tanto un procedimiento iterativo.

1ª iteración:

En esta primera iteración se asumirá por ejemplo que T∞i = 30 ºC y Tse = 8ºC.

- Cálculo de hi (convección forzada):

 i  15,14  10 -6 m2 / s
T  Ti 30  10 
Tfi  si   20 ºC Pri  0,7293
2 2 k  0,02513 W / mK
 i

uiL
Rei   237723  5  105  flujo laminar
i

Como se indicó anteriormente, como aproximación que se está considerando que

las temperaturas en cada cara del parabrisas, Tse y Tsi, son uniformes. Por tanto,
para calcular el número de Nusselt se va a utilizar una correlación para temperatura
de superficie constante, concretamente la siguiente:

Nui  0,064 Re1/i 2 Pri1/ 3  294,1

hiL
Nui   hi  12,2 W/m2K
ki

- Cálculo de he (convección forzada):

 e  13,72  10-6 m2 / s
Tse  Te 8  0 
Tfe    4 ºC Pre  0,7333
2 2 k  0,02393 W / mK
 e

u eL
Ree   874863  5  105  flujo turbulento
e

Nue  (0,037Ree4 / 5  871,3)Pre1/ 3  1106

h eL
Nue   he  44,1 W/m2K
ke

71
Por tanto:

Tsi  Tse Ti  38,27 º C

hi (Ti  Tsi )  k  hi (Ti  Tsi ) 
e Tse  7,821 º C

Con estos valores de T∞i y Tse se inicia la segunda iteración.

2ª iteración:

- Cálculo de hi:
i  15,52  10-6 m2 / s
T  Ti 38,27  10 
Tfi  si   24,13 ºC Pri  0,7283
2 2 k  0,02543 W / mK
 i

uiL
Rei   231943  5  105  flujo laminar
i

Nui  0,064 Re1/i 2 Pri1/ 3  287,7

hiL
Nui   hi  12,2 W/m2K
ki

- Cálculo de he:

 e  13,71 10-6 m2 / s
Tse  Te 7,821  0 
Tfe    3,91 ºC Pre  0,7334
2 2 k  0,02392 W / mK
 e

u eL
Re e   875363  5  105  flujo turbulento
e

Nue  (0,037Re e4 / 5  871,3)Pre1/ 3  1107

h eL
Nue   he  44,13 W/m2K
ke
Por tanto:

Tsi  Tse Ti  38,29 º C

hi (Ti  Tsi )  k  hi (Ti  Tsi ) 
e Tse  7,821 º C

Estos valores de T∞i y Tse son prácticamente iguales que los obtenidos en la primera iteración, por
tanto, se asumen como correctos y ya no es necesario seguir iterando. La temperatura que se
pide en el enunciado es T∞i ≥ 38,29 ºC.

72
Problema 6.6
Por una larga tubería horizontal de diámetro 50 mm y espesor despreciable fluyen 0,3 kg/s de
agua líquida saturada a una temperatura uniforme de 80 ºC. La tubería se cubre con una capa de
material aislante de 5 mm de espesor y conductividad térmica 0,05 W/m∙K y se expone a aire que
se encuentra a 25 ºC, presión atmosférica y en reposo. Se mide la temperatura del aislante en
contacto con el aire, resultando un valor uniforme de 50 ºC. Calcular:
a) Coeficiente de transferencia de calor por convección para el agua.
b) Coeficiente de transferencia de calor por convección para el aire.
c) Calor por unidad de longitud, en W/m, transferido por radiación.

Solución apartado a:

Se trata de un caso de convección forzada, flujo interno.

  355,1  10 -6 Ns/m2

Se evalúan las propiedades necesarias: Tm  80 ºC Pr  2,27
k  0,06561 W / mK


umD 4m 
ReD    21522  2300  flujo turbulento
 D

Para calcular el coeficiente de transferencia de calor por convección se aplica, por ejemplo, la
correlación de Dittus-Boelter. El coeficiente n en este caso es 0,3 porque el fluido se enfría:

NuDi  0,023ReD4 / 5 Pr n  86,07

A partir del número de Nusselt se obtiene el coeficiente de transferencia de calor por convección:

hiD NuDik
NuDi   hi   1129 W / m2K
k D

Solución apartado b:

Se trata de un caso de convección libre.

  16,76  10-6 m2 / s

Pr  0,7251
T  Ts 
Se evalúan las propiedades necesarias: Tf   37,5 ºC k  0,0264 W / mK
2   23,12  10-6 m2 / s

  1/ Tf  1/ 310,5  0,0032 K1

g Ts  T D3
Número de Rayleigh: RaD  GrD Pr   387150


El número de Nusselt se obtiene a partir de la siguiente correlación:

2
 
 0,387Ra1/D 6 
NuDe  0,6  8 / 27   11,2
 1   0,559 / Pr 9 /16  
   

73
A partir del número de Nusselt se obtiene el coeficiente de transferencia de calor por convección:

heD kNuDe
NuDe   he   4,93 W / m2K
k D

Solución apartado c:

Este apartado se resolverá aplicando analogía eléctrica. El circuito eléctrico equivalente se indica
en la siguiente figura. La transferencia de calor por radiación no se puede representar mediante
resistencias porque no es lineal con la variación de temperatura, por lo que simplemente se indica
de modo ilustrativo para entender como se comporta.

Con los datos que da el enunciado, la única manera de calcular el calor transferido por radiación
es mediante la resta del calor total, q ' , menos el calor disipado por convección en la parte exterior
del cilindro, q 'conve :

'
q rad  q '  q conve
'

El calor total por unidad de longitud se calcula como:

Ti  Ts
q '   51,17 W / m
1 ln(De / Di )

Dihi 2k a

El calor por convección en la parte exterior por unidad de longitud se calcula como:

q 'conve  Dehe (Ts  Te )  23,22 W

Por tanto, el calor disipado por radiación por unidad de longitud resulta:

'
q rad  q '  q conve
'
 51,17  23,22  27,9 W

74
 TEMA 7

CONDENSACIÓN Y
EBULLICIÓN

75
Problema 7.1
25 tubos de diámetro 15 mm se disponen formando un haz cuadrado. Exteriormente a los tubos
circula vapor de agua a presión atmosférica. Asumir que la superficie de los tubos se mantiene a
una temperatura uniforme de 98ºC y utilizar las siguientes propiedades para el agua a presión
atmosférica: Tsat = 100ºC, ρl = 960 kg/m3, μl = 2,82∙10-4 kg/m∙s, kl = 0,66 W/m∙K, ρv = 0,6 kg/m3 y
hlv = 2257 kJ/kg. Calcular la masa de vapor de agua condensado por unidad de longitud de tubos.

Solución:

El coeficiente de transferencia de calor por convección viene dado por:

1/ 4 1/ 4
 g (  v )k l3hlv   9,81 960  (960  0,6)  0,663  2257000000 
h  0,729  l l   0,729   
 Nl (Tsat  Ts )D   5  2,82  104 (100  98)  0,015 

 14065,9 W / mK

El área de la superficie total por unidad de longitud de tubos es:

A
 ND  25    0,015  1,18 m
L

De modo que la transferencia de calor por unidad de longitud de tubos es:


Q A
q '   h (Tsat  Ts )  14065,9  1,18  (100  2)  33125,2 W / m
L L

Finalmente, el flujo másico de condensado resulta:


m q ' 33125,2
   0,015 kg / s  m
L hlv 2257000

Problema 7.2
En un condensador se condensa un flujo de 2,5 kg/s de R12 a 37,8º C. Se basa en un haz de
tubos cuadrado de 20 ˣ 20 tubos. Cada tubo tiene un diámetro de 10 mm. El R12 circula por el
exterior de los tubos, mientras que por el interior circula agua que mantiene la temperatura de la
superficie a una temperatura prácticamente uniforme de 32,2 ºC. Utilizar como propiedades del
R12 a la presión de operación las siguientes: Tsat = 37,8ºC, ρl = 1274 kg/m3, μl = 2,1∙10-4 kg/m∙s, kl
= 0,06 W/m∙K, ρv = 48,2 kg/m3 y hlv = 132003 J/kg. Calcular la longitud de tubo.

Solución:

El coeficiente de transferencia de calor por convección viene dado por:

1/ 4 1/ 4
 gl (l  v )k l3hlv   9,81 1274  (1274  48,2)  0,063  132003 
h  0,729    0,729   
 Nl (Tsat  Ts )D   20  2,1 104 (37,8  32,2)  0,01 

 851 W / mK

76
El calor se puede calcular como:

  mh
Q  lv  2,5  132003  330007,5 W

A su vez, el calor también viene dado por la siguiente expresión:

  hA(T  T )  hNDL(T  T )
Q sat s sat s

Igualando las expresiones anteriores se obtiene la longitud de tubo:

  hNDL(T  T )
Q  330007,5  851 400    0,01 L  (37,8  32,2)  L  5,51 m
sat s

Problema 7.3
Por el interior de un tubo de 25 mm de diámetro circula agua a 5 bar en condiciones de ebullición
de convección forzada, estando la superficie del tubo a una temperatura 10 ºC por encima de la
de saturación. Calcular la transferencia de calor por longitud de tubo.

Solución:

Se empleará la correlación de Levy, según la cual el flujo de calor por unidad de superficie viene
dada por:

q ''  283,2p4 / 3 T3

Siendo p la presión en MPa. Sustituyendo valores resulta:

q ''  283,2p4 / 3 T3  283,2  0,54 / 3  103  112388,2 W / m2

El flujo de calor por unidad de longitud de tubo resulta:

q '  q '' D  112388,2    0,025  8822,5 W / m

Problema 7.4
Un condensador se basa en un serpentín de 10 tubos horizontales. Por el interior de los tubos
circula lentamente el agua a condensar a 4,25 kPa. Cada tubo tiene un diámetro interno de 6 mm,
una longitud de 0,4 m y, para simplificar los cálculos, asumir que la temperatura de la superficie
interna se mantiene a una temperatura uniforme de 20ºC. Calcular la tasa de transferencia de
calor. Utilizar como propiedades del agua a 4,25 kPa las siguientes Tsat = 30 ºC, cl = 4183 J/kg∙K,
hlv = 2429000 J/kg, ρl = 995,6 kg/m3, ρv = 0,03 kg/m3, μl = 0,0008 kg/m∙s y kl =0,6 W/m∙K.

77
Solución:

El coeficiente de transferencia de calor se calculará mediante la correlación de Chato:

1/ 4
 g (  v )k l3hlv' 
h  0,555  l l 
 l (Tsat  Ts )D 

Dicha correlación es aplicable para números de Reynolds a la entrada menores de 35000. Con los
datos aportados por el enunciado no es posible calcular el número de Reynolds, pero como se
indica que el agua circula lentamente se asumirá que el número de Reynolds es inferior a ese
valor.

El parámetro viene dado por:

3
hlv'  hlv  c l (Tsat  Ts )  2447000 J / kg
8

Por tanto, sustituyendo valores el coeficiente de transferencia de calor por convección resulta:

1/ 4
 g (  v )k l3hlv' 
h  0,555  l l   6738 W/m2K
 l (Tsat  Ts )D 

Y, finalmente, la tasa de transferencia de calor resulta:

  hA(T  T )  hDNL(T  T )  6738    0,006  8  0,4  (30  20)  4062,2 W

Q sat s sat s

Problema 7.5
Agua a presión 84,5 kPa hierve en un cazo de acero inoxidable pulido mecánicamente. Se estima
que la base del cazo aporta al agua un flujo de calor por unidad de área de 100000 W/m2.
Determinar la temperatura de la base del cazo en contacto con el agua. Asumir como propiedades
del agua a 84,5 kPa: Tsat = 95 ºC, cl = 4212 J/kg∙K, hlv = 2270000 J/kg, ρl = 961 kg/m3, ρv = 0,5
kg/m3, σ = 0,06 kg/s2, μl = 0,0003 kg/m∙s y Prl = 1,8.

Solución:

Se aplicará la correlación de Rohsenow, la cual viene dada por:

m
c l (Ts  Tsat )  q '' σ 
 Csf   Prl
n

hlv  μh
l lv g(ρl  ρv ) 

Sustituyendo valores se obtiene la temperatura de base en contacto con el fluido:

78
 0,33
4212  (Ts  95)  100000 0,06 
 0,013   1,8  Ts  104,1 ºC
2270000  0,0003  2270000 9,81 (961  0,5) 

Problema 7.6
Una placa genera un flujo de calor uniforme de 50000 W/m2. La parte de abajo se encuentra
perfectamente aislada, mientras que la parte superior se refrigera mediante ebullición de un fluido
de propiedades Tsat = 55 ºC, cl = 1100 J/kg∙K, hlv = 85000 J/kg, ρl = 1600 kg/m3, ρv = 15 kg/m3, σ =
0,007 kg/s2, μl = 450∙10-6 kg/m∙s y Prl = 9. Utilizando la correlación de Rohsenow con las
constantes Csf = 0,005, n = 1,7 y m = 0,33 calcular la temperatura de la parte superior de la placa
en contacto con el fluido en estado estacionario.

Solución:

La correlación de Rohsenow viene dada por:

m
c l (Ts  Tsat )  q '' σ 
 Csf   Prl
n

hlv  μh
l lv g(ρl  ρv ) 

Sustituyendo valores se obtiene la temperatura de la superficie en contacto con el fluido:

0,33
1600  (Ts  55)  50000 0,007 
 0,005  6  91,7  Ts  65,66 ºC
85000  450  10  85000 9,81 (1600  15) 

Problema 7.7
Agua a presión atmosférica fluye con una velocidad media de 1,5 m/s por el interior de un tubo de
acero de 15 mm de diámetro cuya superficie se mantiene a 110 ºC. Calcular la transferencia de
calor por unidad de longitud de tubo. Utilizar como propiedades del agua a presión atmosférica:
Tsat = 100ºC, ρl = 958 kg/m3, μl = 2,82∙10-4 kg/m∙s, kl = 0,66 W/m∙K, ρv = 0,6 kg/m3 y hlv = 2257
kJ/kg.

Solución:

Se trata de un caso de convección forzada, flujo interno. Rohsenow y Griffith (1955) propusieron
que el mecanismo de cambio de fase de un líquido fluyendo en un tubo consiste en dos
fenómenos superpuestos: el de ebullición nucleada y el de ebullición estrictamente convectiva, tal
y como indica la siguiente expresión:

q ''  q ''en  q ''ec

El término de ebullición nucleada según la correlación de Rohsenow viene dado por:

m
c l (Ts  Tsat )  q ''en  
 Csf   Prl
n
 q ''en  131546 W / m2
hlv  lhlv g(l  v ) 

El término de ebullición convectiva viene dado por la correlación de Dittus-Boelter pero utilizando
un factor de 0,019 en lugar de 0,023:

79
 hecD
Nu   0,019ReD0,8 Prl0,4
kl

El número de Reynolds viene dado por:

uD uDl
ReD    770,14
l l

Con el número de Reynolds calculado, el número de Nusselt y h resultan:

hecD
Nu   0,019ReD0,8 Prl0,4 = 4,88  hec  8622 W / m  K
kl

Y el término de ebullición convectiva resulta:

q ''ec  hec (Ts  Tsat )  86220 W / m2

Por tanto, el calor como suma del término de ebullición nucleada más el de ebullición convectiva
resulta:

q ''  q ''en  q ''ec  131546  86220  217766 W / m2

Y el calor por unidad de longitud:

q '  q '' D  10261 W / m

80
 TEMA 8

RADIACIÓN TÉRMICA

81
Problema 8.1
Dos placas negras paralelas de 0,6 ˣ 0,6 m están separadas una distancia de 0,3 m. Una de las
placas se mantiene a 800 ºC y la otra a 400 ºC. Calcular el intercambio neto de radiación entre
ambas placas.

Solución:

El factor de forma se obtiene de la gráfica. Teniendo en cuenta que X/D e Y/D resultan 0,6/0,3 = 2,
la gráfica proporciona F aproximadamente 0,3.

La transferencia de calor por radiación resulta:

  AF (E  E )  AF (T 4  T 4 )  5,67  108  0,6  0,6  0,3  (10734  6734 )  6861 W

Q 12 b1 b2 12 1 2

Problema 8.2
Dos cilindros concéntricos de longitud 10 cm tienen diámetros 10 y 20 cm. Calcular el factor de
forma entre los extremos abiertos de los cilindros.

Solución:

Se utilizará la nomenclatura indicada en las siguientes figuras, y los extremos abiertos se

designarán como 3 y 4:

82
De acuerdo con los datos proporcionados en el enunciado, L/r2 = 20/10 = 2 y r1/r2 = D1/D2 = 10/20
= 0,5. Por tanto, a partir de las figuras se obtiene que F22 ≈ 0,3286 y F21 ≈ 0,4126.

Con la relación de reciprocidad se tiene:

D2 20
A1F12  A 2F21  F12  F21  0, 4126  0,8253
D1 10

Para la superficie 2 se tiene:

F21  F22  F23  F24  1

Por simetría F23 = F24, por tanto:

83
 1
F23  F24  (1  0,4126  0,3286)  0,1294
2

Aplicando nuevamente la relación de reciprocidad:

  0,2  0,2
A 2F23  A 3F32  F32  0,1294  0,6901
(0,22  0,12 ) / 4

Se cumple que F11 = F33 = F44 = 0, y para la superficie 3:

F31  F32  F34  1

Para la superficie 1:

F12  F13  F14  1

Y por simetría F13 = F14, de modo que:

1
F13  (1  0,8253)  0,0874
2

Utilizando la relación de reciprocidad se tiene:

  0,1 0,2
A1F13  A 3F31  F31  0,0874  0,233
(0,22  0,12 ) / 4

Por tanto:

F31  F32  F34  1  F34  1- 0,233  0,6901  0,0769

Problema 8.3
Dos planos paralelos muy grandes con emisividades 0,4 y 0,9 intercambian calor. Calcular el
porcentaje de reducción de transferencia de calor cuando se coloca entre ellos un apantallamiento
radiante de emisividad 0,05.

Solución:

La transferencia de calor sin apantallamiento viene dada por:

(T14  T24 )
q ''   0,323(T14  T24 )
1/ 1  / 2  1

El circuito de radiación es el indicado en la siguiente figura:

84
Las resistencias son las siguientes:

1  1 1  0,4
  1,5
1 0,4

1  2 1  0,9
  0,11
2 0,9

1  3 1  0,05
  19
3 0,05

La resistencia total viene dada por:

1,5  2  19  2  1  0,11  41,61

Y el calor total por unidad de área:

(T14  T24 )
q ''   0,024(T14  T24 )
41,61

Por tanto, la transferencia de calor se reduce en el siguiente porcentaje:

0,323  0,024
 100  92,6 % %.
0,323

85
 TEMA 9

INTERCAMBIADORES DE
CALOR

86
Problema 9.1
Por la carcasa de un intercambiador de calor de carcasa y tubos entra aceite a 160 ºC y sale a
100 ºC. Por los tubos circulan 5 kg/s de agua, que entra a 20 y sale a 50 ºC. Se trata de 30 tubos
de diámetro 25 mm que se disponen formando dos pasos, mientras que el flujo en la carcasa se
dispone de un único paso. Según el catálogo del fabricante, el coeficiente global de transferencia
de calor es de 300 W/m2ꞏK. Calcular:
a) Flujo de aceite.
b) Longitud de los tubos.

Solución apartado a:

En primer lugar se calculará el calor específico para el agua y para el aceite:

Tce  Tcs 160  100

Aceite (fluido caliente), Tc    130 º C  cpc  2350 J / kgK
2 2

Tfe  Tfs 20  50
Agua (fluido frío), T f    35 º C  c pf  4182 J / kgK
2 2

Despreciando pérdidas de calor hacia el exterior del intercambiador de calor, a partir de un

balance de energía se obtiene que el calor intercambiado entre el fluido caliente y el fluido frío
viene dado por:

 m
Q  c c pc (Tce  Tcs )  m
 f c pf (Tfs  Tfe )

A partir de la expresión anterior se obtiene el flujo de aceite (fluido caliente):

 c c pc (Tce  Tcs )  m
m  f c pf (Tfs  Tfe )  m
 f  2350  (160  100)  5  4182  (50  20)  m
 c  4,4 kg / s

Solución apartado b:

Este apartado se puede realizar utilizando dos métodos, el DTML y el ε–NUT. A continuación se
muestra el procedimiento para ambos métodos.

Método DTML:

  UAT
Q

siendo T  FTmlcontraflujo

El factor de corrección F se determina a partir de la siguiente gráfica:

87
En la gráfica anterior se entra con los siguientes datos y se obtiene que F es aproximadamente
0,96. Nótese que, según la figura, t representa el fluido que circula por el interior de los tubos (en
este caso el agua) y T el que circula exteriormente a los tubos (en este caso el aceite):

ts  te 50  20 
P   0,21
Te  t e 160  20 
 F  0,96
Te  Ts 160  100
R  2 
ts  te 50  20 

La diferencia de temperaturas media logarítmica para un intercambiador de calor a contraflujo es:

(Tce  Tfs )  (Tcs  Tfe ) (160  50)  (100  20)

Tmlcontraflujo    94,2 º C o K
 Tce  Tfs   160  50 
ln   ln  
 Tcs  Tfe   100  20 

  UAFT
Tal y como se indicó anteriormente, la transferencia de calor viene dada por Q mlcontraflujo .

Como se trata de 20 tubos, el área total es NπDL, siendo N el número de tubos. Por tanto resulta:

  UNDLFT
Q  627450  300  30    0,025  L  0,96  94,2  L  9,8 m
mlcontraflujo

Es decir, se necesitan 30 tubos de 9,8 m cada uno. Como cada tubo tiene dos pasos, serían 4,9 m
por cada paso.

88
Método ε–NUT:

Primeramente se calcularán las capacidades térmicas tanto para el fluido caliente como para el
fluido frío:

Cf  m
 f c pf  5  4182  20910 W/K

Cc  m
 c c pc  4,45  2350  10457 W/K

A continuación se calculará el cociente entre la capacidad térmica mínima, que en este caso se
corresponde con el fluido caliente, y la capacidad térmica máxima, que se corresponde con el
fluido frío:

Cmin 10457
  0,5
Cmax 20910

El calor máximo viene dado por:


Q max  Cmin (Tce  Tfe )  10457  (160  20)  1464050 W

El calor real viene dado por:

 m
Q  c c pc (Tce  Tcs )  m
 f c pf (Tfs  Tfe )  627450 W

El rendimiento del intercambiador de calor viene dado por:


Q 624750
   0,43

Qmax 1464050

De la siguiente gráfica se obtiene que el NUT es de aproximadamente 0,65:

A partir del NUT se obtiene el área según la siguiente expresión:

89
 UA CminNUT 10457  0,65
NUT   A   22,6 m2
Cmin U 300

Y a partir del área la longitud de cada tubo:

A  NDL  L  9,6 m

Como se puede observar, el resultado obtenido mediante ambos métodos es parecido pero no
igual debido a la limitada precisión que se puede obtener al utilizar las gráficas.

Problema 9.2
Un flujo de gases de escape calientes entra a un intercambiador de calor de tubos con aletas de
flujo cruzado a 300 ºC y sale a 100 ºC. Por los tubos fluye 1 kg/s de agua presurizada que entra a
35 ºC y sale a 125 ºC. Asumir el calor específico de los gases de escape como 1000 J/kg∙K y el
coeficiente global de transferencia de calor U = 100 W/m2∙K. Calcular:
a) Flujo de gases.
b) Área de intercambio de calor.

Solución apartado a:

Aplicando balance de energía:

 m
Q  c c pc (Tce  Tcs )  m
 f c pf (Tfs  Tfe )

El calor específico del fluido caliente (gases de escape) es dato del problema, cpc = 1000 kJ/kgK.
El calor específico del fluido frío (agua) se consulta en tablas evaluándolo a la temperatura
promedio del fluido frío, es decir:

Tfe  Tfs 35  125

Agua, T f    80 º C  c pf  4194 J / kgK
2 2

Por tanto:

 c c pc (Tce  Tcs )  m
m  f c pf (Tfs  Tfe )  m
 c  1000  (300  100)  1 4194  (125  35)  m
 c  1,89 kg / s

Solución apartado b:

Este apartado se puede realizar utilizando dos métodos, el DTML y el ε–NUT. A continuación se
muestra el procedimiento para ambos métodos.

Método DTML:

  UAT
Q

siendo T  FTmlcontraflujo

El factor de corrección F se determina a partir de la siguiente gráfica:

90
En la gráfica anterior se entra con los siguientes datos y se obtiene que F es aproximadamente
0,89. Nótese que, según la figura, t representa el fluido que circula por el interior de los tubos y T
el que circula exteriormente a los tubos:

t s  t e 125  35 
P   0,34 
Te  t e 300  35 
 F  0,89
Te  Ts 300  100
R   2,22
ts  te 125  35 

La diferencia de temperaturas media logarítmica para un intercambiador de calor a contraflujo es:

(Tce  Tfs )  (Tcs  Tfe )

Tmlcontraflujo   111 º C
 Tce  Tfs 
ln  
 Tcs  Tfe 

Por tanto:

  UAFT
Q  37460  100  A  0,89  111  A  38,2 m2
mlcontraflujo

91
Método ε–NUT:

Primeramente se calcularán las capacidades térmicas tanto para el fluido caliente como para el
fluido frío:

Cf  m
 f c pf  1 4194  4194 W/K

Cc  m
 c c pc  1,89  1000  1890 W/K

A continuación se calculará el cociente entre la capacidad térmica mínima, que en este caso se
corresponde con el fluido caliente, y la capacidad térmica máxima, que se corresponde con el
fluido frío:

Cmin 1890
  0,45
Cmax 4194

El calor máximo viene dado por:


Q max  Cmin (Tce  Tfe )  1890  (300  35)  500850 W

El calor real viene dado por:

 m
Q  c c pc (Tce  Tcs )  m
 f c pf (Tfs  Tfe )  377460 W

El rendimiento del intercambiador de calor viene dado por:


Q 377460
   0,75

Qmax 500850

De la siguiente gráfica se obtiene que el NUT es de aproximadamente 2:

A partir del NUT se obtiene el área según la siguiente expresión:

UA CminNUT 2  1890
NUT   A   37,8 m2
Cmin U 100

92
Como se puede observar, el área obtenida mediante ambos métodos es muy parecida. No
proporcionan el mismo resultado debido a la limitada precisión que se puede obtener al utilizar las
gráficas.

Problema 9.3
El radiador de un automóvil se puede ver como un intercambiador de calor en flujo cruzado con
ambos fluidos sin mezclar. Analizar el caso para el que el agua entra al radiador a 90 ºC y sale a
85 ºC con un flujo de 0,8 kg/s. El agua se enfría mediante aire que entra a 20 ºC con un flujo de
0,9 kg/s. El coeficiente global de transferencia de calor es 200 W/m2∙K. Para simplificar los
cálculos, asumir estado estacionario.
a) Calcular la temperatura de salida del aire.
b) Calcular el área de transferencia de calor.

Solución apartado a:

Aplicando balance de energía:

 m
Q  c c pc (Tce  Tcs )  m
 f c pf (Tfe  Tfs )

El calor específico de cada fluido se evalúa a la temperatura promedio del fluido. Del fluido
caliente (agua) resulta:

Tce  Tcs 90  85
Agua, T c    87,5 º C  c pc  4201 J / kgK
2 2

Del fluido frío (aire) se conoce la temperatura de entrada pero precisamente la temperatura de
salida es lo que se pide calcular en el apartado a. Se seguirá por tanto un procedimiento iterativo

1ª iteración:

Se asumirá que por ejemplo T f = 310 K, con lo cual resulta un calor específico de cpf = 1005
J/kgK.

Por tanto:

 c c pc (Tce  Tcs )  m
m  f c pf (Tfe  Tfs )  0,8  4201(90  85)  0,9  1005(Tfs  20)  Tfs  38,6 ºC

2ª iteración:

Se utilizará el valor de Tfs = 38,6 ºC obtenido en la primera iteración. Utilizando este valor, el calor
específico del agua resulta:

Tfe  Tfs 20  38,6

Aire, T f    29,3 º C  c pf  1007 J / kgK
2 2

Por tanto:

 c c pc (Tce  Tcs )  m
m  f c pf (Tfe  Tfs )  0,8  4201(90  85)  0,9  1007(Tfs  20)  Tfs  38,6 ºC

93
Como se puede observar, la temperatura de salida del aire presenta el mismo valor que en la
primera iteración (tomando un decimal no se percibe la diferencia), por lo que se asumirá este
valor como correcto y ya no es necesario realizar más iteraciones.

Solución apartado b:

Este apartado se puede realizar utilizando dos métodos, el DTML y el ε–NUT. A continuación se
muestra el procedimiento para ambos métodos.

Método DTML:

  UAT
Q

siendo T  FTmlcontraflujo

El factor de corrección F se determina a partir de la siguiente gráfica:

En la gráfica anterior se entra con los siguientes datos y se obtiene que F es aproximadamente
0,99. Nótese que, según la figura, t representa el fluido que circula por el interior de los tubos, es
decir, el agua y T el que circula exteriormente a los tubos, es decir, el aire:

t s  t e 85  90 
P   0,07 
Te  t e 20  85 
 F  0,99
T  Ts 20  38,6
R e   3,71
ts  te 85  90 

94
La diferencia de temperaturas media logarítmica para un intercambiador de calor a contraflujo es:

(Tce  Tfs )  (Tcs  Tfe )

Tmlcontraflujo   58 º C
 Tce  Tfs 
ln  
 Tcs  Tfe 

Por tanto:

  UAFT
Q  16806  200  A  0,98  58  A  1,46 m2
mlcontraflujo

Método ε–NUT:

Primeramente se calcularán las capacidades térmicas tanto para el fluido caliente como para el
fluido frío:

Cf  m
 f c pf  0,9  1007  905,9 W/K

Cc  m
 c c pc  0,8  4201  3361 W/K

A continuación se calculará el cociente entre la capacidad térmica mínima, que en este caso se
corresponde con el fluido frío, y la capacidad térmica máxima, que se corresponde con el fluido
caliente:

Cmin 905,9
  0,27
Cmax 3361

El calor máximo viene dado por:


Q max  Cmin (Tce  Tfe )  63415 W

El calor real viene dado por:

 m
Q  c c pc (Tce  Tcs )  m
 f c pf (Tfs  Tfe )  16806 W

El rendimiento del intercambiador de calor viene dado por:


Q 16806
   0,26

Q 63415
max

De la siguiente gráfica se obtiene que el NUT es de aproximadamente 0,3:

95
A partir del NUT se obtiene el área según la siguiente expresión:

UA C NUT
NUT   A  min  1,35 m2
Cmin U
Como se puede observar, el área obtenida mediante ambos métodos es muy parecida. No
proporcionan el mismo resultado debido a la limitada precisión que se puede obtener al utilizar las
gráficas.

Problema 9.4
Vapor saturado a 100 ºC se condensa en un intercambiador de calor de carcasa y tubos de un
paso por la carcasa y dos pasos por los tubos, con un área superficial de 0,5 m2 y coeficiente
global de transferencia de calor de 2000 W/m2∙K. Como fluido frío se utiliza un flujo de 0,5 kg/s de
agua que entra a 15 ºC. Determinar la temperatura de salida del agua fría.

Solución:

Se resolverá el problema utilizando al método ε–NUT porque con los datos que proporciona el
enunciado la resolución utilizando el método DTML sería muy tediosa.

Es necesario determinar el calor específico del fluido frío. Como solamente se conoce la
temperatura de entrada pero no la de salida se realizará un procedimiento iterativo.

- 1ª iteración´:

En esta primera iteración se asumirá, por ejemplo, que la temperatura promedio del fluido frío es
de 27 ºC. A esta temperatura se evalúa el calor específico, que resulta el indicado a continuación:

Tfe  Tfs
Agua, T f   27 º C  c pf  4183,4 J / kgK
2

Las capacidades térmicas tanto para el fluido caliente como para el fluido frío. Se indican a
continuación. Nótese que para el fluido caliente la capacidad térmica tiende a infinito por estar a
temperatura constante:

Cf  m
 f c pf  0,5  4183  2091,7 W/K

96
 Cc  

El cociente entre la capacidad térmica mínima resulta cero:

Cmin 2091,7
 0
Cmax 

El número de unidades de transferencia resulta:

UA 2000  0,5
NUT    0,48
Cmin 2091,7

Para un valor de Cmin/Cmax = 0 y NUT = 0,48, de la siguiente gráfica se obtiene la eficiencia, que
resulta aproximadamente 0,38:

Teniendo en cuenta el rendimiento, el calor transferido en esta primera iteración resulta:


Q
    Q


Q max  Cmin (Tce  Tfe )  0,38  2091,7  (100  15)  67561,9 W
Q max

Y, finalmente, la temperatura de salida del agua fría obtenida en esta primera iteración resulta:

 m
Q  f c pf (Tfs  Tfe )  67651,9  0,5  4183,4  (Tfs  15)  Tfs  47,3 ºC

- 2ª iteración:

En esta segunda iteración se utilizará el valor de la temperatura de salida del agua fría obtenida
en la primera iteración. El calor específico del agua fría resulta:

Tfe  Tfs 15  47,3

Agua, T f    31,1 º C  c pf  4183,3 J / kgK
2 2

97
Las capacidades térmicas resultan:

Cf  m
 f c pf  0,5  4183,3  2091,6 W/K

Cc  

El cociente entre la capacidad térmica mínima resulta cero:

Cmin 2091,6
 0
Cmax 

El número de unidades de transferencia resulta:

UA 2000  0,5
NUT    0,48
Cmin 2091,6

Para un valor de Cmin/Cmax = 0 y NUT = 0,48, de la gráfica anterior se obtiene que la eficiencia
resulta aproximadamente 0,38, al igual que en la primera iteración.

El calor transferido en esta segunda iteración resulta:


Q
    Q


Q max  Cmin (Tce  Tfe )  0,38  2091,6  (100  15)  67558,7 W
Q max

Y, finalmente, la temperatura de salida del agua fría obtenida en esta primera iteración resulta:

 m
Q  f c pf (Tfs  Tfe )  67558,7  0,5  4183,3  (Tfs  15)  Tfs  47,3 ºC

Como se puede observar, en la segunda iteración se ha obtenido el mismo valor de temperatura

de salida del agua fría que en la primera iteración, por lo cual se considera este resultado como
válido y ya no se continúa iterando.

Problema 9.5
Vapor saturado a 0,14 bar (Tsat = 52,5 ºC y hcf = 2376000 J/kg) se condensa en un intercambiador
de carcasa y tubos con un paso por la carcasa y dos pasos por los tubos que consisten en 130
tubos de bronce (k = 114 W/m∙K), cada uno con una longitud por paso de 2 m. Los tubos tienen
diámetros interior y exterior de 13,4 y 15,9 mm, respectivamente. El coeficiente de transferencia
de calor para la condensación en las superficies exteriores de los tubos es de 13500 W/m2∙K. El
agua de enfriamiento entra en los tubos a 20 ºC con una velocidad media de 1,25 m/s. Utilizar
como propiedades del agua de enfriamiento ν = 6,6∙10-7 m2/s, k = 0,6176 W/m∙K, Pr = 4,437, ρ =
995,7 kg/m3 y cp = 4183 J/kgꞏK.
Calcular:
a) Coeficiente global de transferencia de calor.
b) Temperatura de salida del agua de enfriamiento.
c) Flujo de condensación de vapor.

98
Solución apartado a:

El coeficiente global de transferencia de calor viene dado por la siguiente expresión:

1 1
U 
Rglobal A Rglobal 2reL

Contabilizando la convección en el interior de los tubos, la conducción a lo largo del espesor del
tubo y la convección exterior a los tubos, resulta:

1 ln(re / ri ) 1
Rglobal  Rconvi  Rcond  Rconve   
2rhL
i i 2 k bronce L 2 reheL

Por tanto:

1 1
U 
Rglobal 2reL 1 r ln(r e / ri ) 1
e
 re 
hi ri k bronce he

Cálculo de hi:

Se trata de un problema de convección, flujo interior. Es necesario calcular el número de Reynolds,

el cual resulta:

umDi 1,25  0,0134

ReD    25370
 6,602  107

Como el número de Reynolds resultó mayor de 2300 se trata de flujo turbulento. Por ejemplo,
empleando la correlación de Dittus-Boelter resulta:

NuD  0,023ReD4 / 5 Pr n  0,023ReD4 / 5 Pr 0,4  139,3

hD
NuD  i i
 hi  6421 W/m2K
k agua

99
Por tanto:

1 1 1
U    3693 W/m2K
Rglobal 2reL 1 e r ln(r / r ) 1 1 15,9 ln(15,9 /13,4) 0,0159 1
 e i
r   
hi ri k bronce e he 642113,4 114 2 13500

Solución apartado b:

Utilizando el método ε–NUT:

Cc  

Cf  m
 f c pf

El flujo másico de agua fría viene dado por:

Di2 0,01342
 f  f um A f N  f um
m N  995,7  1,25  130  22,81 kg / s
4 4

Por tanto:

Cf  m
 f c pf  22,81 4183  95404 W / K

A partir de las gráficas de ε-NUT se obtiene la eficiencia ε:

Cmin 
0 
Cmax 

   0,6

UA 3693DeLNNpasos 3693    0,0159  2  130  2
NUT    1 
Cmin 95404 95404 

El calor transferido se obtiene a partir de la siguiente expresión:


Q
    Q

max  Cmin (Tce  Tfe )  0,6  95404  (52,5  20)  1,86  10 W
6
Q

Q max

Y finalmente la temperatura de salida del agua fría resulta:

 m
Q  f c pf (Tfs  Tfe )  1,86  106  22,81 4183  (Tfs  20)  Tfs  39,5 º C

Solución apartado c:

El flujo de condensación de vapor se obtiene calculando el calor mediante la siguiente expresión:

 m
Q  c hcf  1,86  106  m
 c 2376000   c  0,783 kg/s
m

100