Ejercicios Resueltos de Calor

Colección de Problemas Propuestos y Resueltos de Transmisión de Calor
2 (Conductividad variable)
Considérese un muro compuesto por dos capas cuyas características son las siguientes:
Capa 1: espesor 0.4 m, conductividad: k 1 = 0.9 (1 + 0.006 T ) [W / m·K ]

Capa 2: espesor 0.05 m, conductividad: k 2 = 0.04 W / m·K
Y sometido a un flujo solar en la cara exterior de 300 W/m², esta cara se encuentra en
contacto con aire a 40°C (Coeficiente convectivo exterior 10 W/m²K). La cara interior se
encuentra en contacto con aire a 20°C (Coeficiente convectivo interior 5 W/m²K)
Calcular:
a. Flujo de calor por unidad de área que atraviesa el muro.

b. Temperatura en las dos superficies extremas y en la interfase entre las dos capas
Solución:
Capa 1: e1 = 0.4 m k 1 = 0.9 (1 + 0.006 T ) [W / m·K ]

Datos:
-
- Capa 2: e2 = 0.05 m k 2 = 0.04 W / m·K
- Condición de contorno exterior: q′sol
′ = 300 W / m² Text = 40°C hext = 10 W / m²K
- Condición de contorno interior: Tint = 20°C hint = 5 W / m²K
Incógnitas:
a. Flujo de calor por unidad de área que atraviesa el muro: q′′
b. Temperatura de las superficies: T1 , T2 , T3
Esquema:
q′sol
′
hext hint
Text T1 T2 T3 Tint
q′′
25
Desarrollo:
a. Flujo de calor por unidad de área que atraviesa el muro:
La ecuación diferencial en la capa 1 será la siguiente:
 k (T ) k (T ) q′′ = − k (T )
d  dT 
 = 0; = cte;
dT dT
dx  dx  dx dx
variables con la temperatura siguiendo una ley lineal del tipo: k (T ) = k 0 (1 + β T ) . Si

El flujo de calor por unidad de área debe ser constante. La conductividad es
integramos la ecuación anterior para toda la capa 1:
− ∫ q′′ dx = ∫ k 0 (1 + β T ) dT
e1 T2
0 T1
− q′′ e1 = k 0 (T2 − T1 ) + k 0 (
β 2
)
T2 − T12 = k 0 (T2 − T1 ) + k 0 (T2 − T1 )(T2 + T1 ) = k 0 (T2 − T1 )[1 + β Tmed ]
β
− q′′ e1 = k 0 [1 + β Tmed ] (T2 − T1 ) = k med (T2 − T1 )

2 2
Ahora impondremos las dos condiciones de contorno:
x=0 → hext (Text − T1 ) + q′sol

′ = q′′
Esta condición de contorno podemos expresarla como si fuera una condición de
contorno puramente convectiva contra una temperatura equivalente (Temperatura
x = 0 → hext  Text + sol  − T1  = hext (Tsol,aire − T1 ) = q′′

sol-aire) de 70°C
 q′′  
 hext  
En el otro contorno la condición será:

T − Tint
x = e1 → q′′ = 2
+
k2 1
e2 hint
Tenemos pues 3 ecuaciones con 3 incógnitas ( T1 , T2 , q′′ ):
 (T2 − T1 )
  T + T1 
− q′′ e1 = k 0 1 + β  2
  2 
(1)
q′′ = hext (Text − T1 ) + q′sol

′ (2)
T2 − Tint
q′′ = (3)
+
k2 1
e2 hint
Igualando la ecuación (2) con la (3) e introduciendo la (2) en la (1) tenemos 2
ecuaciones con 2 incógnitas:
hext (Text − T1 ) + q′sol

T2 − Tint
′ =
+
k2 1
e2 hint
− e1 (hext (Text − T1 ) + q′sol

′ ) = k 0 1 + β  2  (T2 − T1 )
  T + T1 
  2 
26
Si despejamos en la primera T2 y lo introducimos en la segunda tendremos una

ecuación cuadrática en T1:
T1 = 67.33°C T2 = 58.72°C q′′ = 26.7 W / m²
Finalmente podemos calcular la temperatura en la superficie 3:
T3 − Tint
q′′ = ; T3 = 25.34°C
1
hint
Si pintamos la distribución de temperaturas será la siguiente:
T1
T2
Text
T3
Tint
27
3 (Generación)
Una tubería de acero de 36cm de diámetro exterior, 34cm de diámetro interior y
conductividad térmica 40 kcal/h·m°C, transporta fueloil a 50°C a través de un local que se
encuentra a 10°C. Con objeto de mantener constante la temperatura del fueloil, se rodea la
tubería con una resistencia eléctrica asimilable a una capa de 1 cm de material de
conductividad térmica 200 cal/h·m°C, y una generación uniforme de calor G. Calcular:
A. Valor mínimo de G en kcal/h·m3 para que la pérdida de calor del fuel sea nula.
B. Distribución de temperatura en la tubería y en la resistencia.
Los coeficientes de película en el exterior e interior de la tubería son 15 y 45 kcal/h·m2°C

respectivamente.
Solución:
Datos:
- Capa 1 (tubería acero): D1 = 0.34 m D2 = 0.36 m k t = 40 kcal / h·m°C
- Capa 2 (resistencia eléctrica): D3 = 0.38 m k R = 200 kcal / h·m°C
- Condición de contorno exterior: Text = 10°C hext = 15 kcal / h·m²°C
- Condición de contorno interior: Tint = 50°C hint = 45 kcal / h·m²°C
Incógnitas:
A. Generación de calor volumétrica en la resistencia: G, para que las pérdidas sean
nulas
B. Distribución de temperaturas: T(r).
Esquema:
G
Text, hext
D3 D2 D1 Tint, hint
T1
T3
T2
28
Desarrollo:
A. G, para que la pérdidas sean nulas:

Para que las pérdidas de calor del fueloil sean nulas es necesario que el calor por
convección en el interior de la tubería sea cero, o lo que es lo mismo, que no exista
diferencia de temperaturas entre el fluido y la superficie interna del acero:
Tint = T1 = 50°C
Podemos decir también, que como la capa de acero no tiene generación interna el
flujo de calor por conducción (q) a través de ella debe ser constante y como en la
superficie interior es cero, debe ser cero en toda la capa cilíndrica, o lo que es lo
mismo la temperatura debe ser constante en toda la capa de acero, e igual a la del
fueloil:
Tint = T1 = T2 = 50°C
Para la segunda capa tenemos generación y por tanto la ecuación general de

transmisión de calor en este medio es:
1 d  dT 
k 2 r +G=0
r dr  dr 
Con condiciones de contorno:
r = r2 → − k2 =0
dT
dr r 2
r = r3 → − k2 = hext (T3 − Text )

dT
dr r 3
Integrando la ecuación diferencial una vez e imponiendo la primera condición de

contorno, tendremos:
= − r 2 + C1; 0=− r2 + C1; C1 =

dT G G 2 G 2
( )
k 2r r2
dr 2 2 2
= − r − r22
dT G 2
k 2r
dr 2
− r22 ln(r ) + C2

Integrando una segunda vez tendremos:
∫ ∫
 r22  G  r2 
= −  −  = − 
  
G
2 k2  r  2 k2  2 
dT r dr; T
Si imponemos la segunda condición de contorno:
=  r3 − 2  = hext (T3 − Text ) = hext  − ( )

G r2   G  r32  
− k2  −  + C2 − Text 
 2 k2  2  
dT
2 r3   
2
 
r2 ln r3
dr r
− r22 ln(r3 ) + Text

 r2  G  r32 
3
C2 =  r3 − 2  + 
  
G
2 hext  r3  2 k 2  2 
29
T (r ) = − ( ) 2 k  2 − r2 ln(r3 ) + 2 h
Luego:
 r2  G  r32   r2 
 − +    r3 − 2  + Text
2  r3 
G G
  2   
2 2
r2 ln r
2 k2 ext
Si ahora imponemos que la temperatura en la cara interior tiene que ser 50°C
( ) ( )
Tendremos una ecuación de la cual obtenemos el valor de la generación:
G  r22  G  r32  G  r22 
T2 = 50°C = −  −  +  −  +  
2 k 2  2
2   3  −
 3 r  + Text
 2 k2  2  2 hext  3 
2 2
r2 ln r r2 ln r r
G = 61598 kcal / h·m³
B. Distribución de temperaturas:
Y por tanto la distribución de temperaturas será:
T (r ) = −77 r 2 + 5 ln(r ) + 61.05 [°C]
T(°C)
T2
T3
T1
Text
r1 r2 r3 r (m)
30
4 (Cuestión)
La siguiente figura muestra la distribución de densidad de flujo de calor q’’ (W/m²) en el
espesor de un muro con tres capas. La conductividad de las tres capas es constante, siendo
la del material A, el doble (2k) a la del material C (k).
A. Calcular el valor de la generación volumétrica G en el material B.

B. Calcular que proporción existe entre dT/dx en el material A y el C.
C. Dibujar cualitativamente la distribución de temperatura en el muro en función de x.
q’’(W/m²)
2k k
A B C
50
0
-100 G
X
L L L
Solución:
Datos:
- Capa A: k A = 2 k , Capa C: k C = k
- Espesores: e A = eB = eC = L
- Distribución de flujo de calor por unidad de área: q′′(x )
Incógnitas:
A. Generación de calor volumétrica en el material B: G.
dTA dTC
B.
dx dx
C. Dibujar cualitativamente: T(r).
Desarrollo:
A. G en el material B:
Si realizamos un balance de energía en la capa B, tendremos (debemos suponer
−q′A′ + q′C′
que los flujos de calor son positivos en la dirección creciente de la coordenada x):
q′A′ A + G A L = q′C′ A; G = =
150
W / m³
L L
B. Proporción entre dT/dx en el material A y en el C:
31
q′′
q′A′ = −k A =− A = =
dTA dTA 100 50
;
dx dx kA 2k k
q′′
q′C′ = −k C =− C =−
dTC dTC 50
;
dx dx kC k
Luego las derivadas en ambos medios son de igual valor y signo contrario:
= −1 ; =− C
dTA dTC dTA dT
dx dx dx dx
C. Distribución de temperaturas:
La ecuación diferencial en el medio B es la siguiente:
d  dT 
k B  + G = 0; + =0
d2 T G
dx  dx  dx 2 K B
Luego la distribución de temperaturas tendrá forma de polinomio de 2º orden
T (x ) = −
(parábola):
x 2 + C1x + C2
G
2 KB
q′′(x ) = −k B
El flujo de calor será:
= G x − k BC1
dT
dx
q′′(0) = −100;
Si imponemos condiciones de contornor:
x=0 → − k BC1 = −100; C1 =
100
kB
x=L → q′′(L ) = 50; GL − 100 = 50; G=

150
L
q′′(x ) = x − 100; q′′(x ) = 0 → x =

La temperatura tendrá un máximo donde el flujo sea igual a 0:
L = L
150 100 2
L 150 3
2
T (°C)
L
3
B
A C
x
L L L
32
5 (Aletas)
Se separan aire y agua mediante una pared plana hecha de acero. Se propone aumentar la
razón de transferencia de calor entre estos 2 fluidos agregando aletas rectangulares rectas
de acero de 1,5 mm de espesor, 2,5 cm de longitud y espaciadas 1 cm entre los centros.
Calcular el porcentaje de aumento en la transferencia de calor al añadir aletas en:
A. Aire exterior
B. Lado del agua
C. Ambos lados de la pared plana
El coeficiente de película en aire es 9 kcal/h·m2ºC y en agua 200 kcal/h·m2ºC. La

conductividad del acero es 37 kcal/h·mºC.
Solución:
Datos:
- Lado del aire: ha = 9 kcal / h·m²°C
- Lado del agua (w): hw = 200 kcal / h·m²°C
- Conductividad del acero de la pared y las aletas: k = 37 kcal / h·m°C
- Dimensiones de aletas: L = 2.5 cm; S = 1 cm; δ = 1.5 mm
Incógnitas:
A. Porcentaje de aumento de la transferencia de calor con aletas en el aire
B. Porcentaje de aumento de la transferencia de calor con aletas en el agua
C. Porcentaje de aumento de la transferencia de calor con aletas en ambos lados
Esquema:
Elemento
S = 0.01 m repetitivo
δ = 0.0015 m
L = 0.025 m
33
Desarrollo:
= 8.612 ∆T [kcal / h·m²]

Calculemos en primer lugar la transferencia a través de la pared sin aletas:
∆T
q′sin
′ =
+
1 1
ha hw
Hemos supuesto que la resistencia asociada a la capa acero es despreciable frente
a las resistencias convectivas.
A. Aletas en el aire:
El flujo de calor en este caso será:
∆T
qaire =
+
1 1
ha A totalηs hw A sin
Donde la eficiencia global de la superficie aleteada se calcula como:
A a ηa + A p
ηs =
A total
Si calculamos las áreas para el elemento repetitivo tendremos los siguientes
Área de aleta: A a = (2 L + δ ) W [m²]

valores:
Área primaria: A p = (S − δ ) W [m²]

Área Total: A total = (S + 2 L ) W [m²]
Área del lado del agua: A sin = S W [m²]
Siendo W la longitud perpendicular al plano del dibujo
Al ser una aleta recta la eficiencia de aleta puede calcularse como:
tgh(m L C )
ηa = = 0.934 m= = 18.01, L c = L + δ / 2 = 0.02575 m
2 ha
kδ
Donde
mLC
= 0.4028 W ∆T [kcal / h]
A a ηa + A p ∆T
Por tanto ηs = = 0.9433 y qaire =
+
A total 1 1
ha A totalηs hw A sin
qsin = 8.612 S W ∆T = 0.08612 W ∆T [kcal / h]

El flujo de calor sin aletas para cada unidad repetitiva será:
y por tanto el porcentaje de aumento es: 368%
B. Aletas en el agua:
∆T
qagua =
+
1 1
hw A totalηs ha A sin
34
La obtención de la eficiencia global de la superficie aleteada se obtiene igual que
tgh(m L C )
antes pero usando el coeficiente de película del agua:
ηa = = 0.446 Donde m = = 84.9, L c = L + δ / 2 = 0.02575 m

2 hw
mLC kδ
A a ηa + A p
Por tanto ηs = = 0.5245 y
= 0.08873 W ∆T [kcal / h]
A total
∆T
qagua =
+
1 1
hw A totalηs ha A sin
Por tanto el porcentaje de aumento es: 3%
B. Ambos lados:
= 0.4712 W ∆T [kcal / h]
∆T
qambos =
+
1 1
hw A totalηs,w ha A totalηs,a
Con lo cual el aumento con respecto a la situación inicial es del: 447%.
35
6 (Cuestión)
La pared de un cilindro está compuesta por dos capas de materiales con conductividad kA y
kB. Ambos materiales están separados por una resistencia eléctrica muy delgada de muy alta
conductividad. Por el interior de la tubería circula un líquido a temperatura Ti y con un
coeficiente de película hi. En el exterior la temperatura y el coeficiente de película son
respectivamente Te y he.
A. Obtener la temperatura de la resistencia eléctrica cuando el calor disipado por ésta

es nulo.
B. Obtener la temperatura de la resistencia eléctrica cuando el calor disipado por ésta
es q’’c (W/m²).
Tc, q’’c
Ti , hi
Te, he
A B
Solución:
Datos:
- Resistencia eléctrica muy delgada de alta conductividad que genera: q′c′ [W / m²]
- Capa A: k A , Capa B: k B
- Condición de contorno exterior: Te , he

- Condición de contorno interior: Ti , hi
Incógnitas:
A. Tc cuando q′c′ = 0
B. Tc cuando q′c′ ≠ 0
Desarrollo:
A:
Utilizando la analogía eléctrica de conducción, podemos expresar el flujo de calor

desde la superficie intermedia hacia el interior y hacia el exterior:
Tc − Ti Tc − Ti
ln(rc / ri ) ln(rc / ri )
qi = =
+ +
1 1
hi A i 2π k BL hi 2π riL 2π k BL
Tc − Te Tc − Te
ln(re / rc ) ln(re / rc )
qe = =
+ +
1 1
he A e 2π k A L he 2π reL 2π k A L
36
la suma de ambos debe ser igual a cero: qi + qe = 0

Ambos han sido expresados como flujos salientes de la superficie intermedia, luego
Si despejamos la temperatura de la interfase de la expresión anterior tendremos:
ln(rc / ri ) ln(re / rc )
+
Ti Te
+ +
1 1
Tc =
hi ri kB he re kA
+
1 1
+ +
1 1
hi ri kB he re kA
B:
Para el caso de que exista una generación de energía superficial el balance de

energía en esa superficie sería el siguiente:
qi + qe = q′c′ 2π rcL
Y por tanto al despejar la temperatura tendríamos:
q′c′ rc + +
Ti Te
+ +
1 1
Tc =
hi ri kB he re kA
+
1 1
+ +
1 1
hi ri kB he re kA
37
2.2. Problemas resueltos de convección
7 (Convección forzada flujo externo)

Una superficie plana horizontal de ancho w = 1 m, se mantiene a una temperatura
uniforme de 230°C, mediante el uso de resistencias eléctricas controladas
independientemente. Cada resistencia tiene una longitud de 50 mm. Si sobre la superficie
circula aire atmosférico a 25°C, con una velocidad de 60 m/s, determinar la resistencia que
presenta un mayor consumo y el valor del mismo.
Solución:
Datos:
- Flujo externo de aire a presión atmosférica: u∞ = 60 m / s T∞ = 25°C
- Temperatura superficial: Ts = 230°C
- Longitudes: w = 1 m L i = 0.05 m
Incógnitas:
La resistencia i-esima presenta la mayor potencia eléctrica, calcular cual es “i” y y
cuanto vale esta potencia qi .
Esquema:
T∞ = 25°C
u∞ = 60 m / s
Ts = 230°C
L i = 50 mm
q1 q2 q3 q4 q5 q6
Hipótesis:
- Régimen permanente
- Efectos de radiación despreciables
- Superficie inferior de la placa adiabática
38
Desarrollo:
La potencia eléctrica consumida por cada una de las resistencias será cedida por
convección al aire, debido a que la cara inferior de la placa se encuentra
perfectamente aislada. Por tanto, buscar la placa con máxima potencia eléctrica es
lo mismo que buscar la placa con flujo de calor por convección máximo.
Veamos primero donde se produce la transición a régimen turbulento:
Las propiedades en las correlaciones de convección forzada flujo externo se evalúan

en la mayoría de los casos a la temperatura media de película:
T∞ + Ts
Tmp = = 127.5 °C
2
µ = 22.95·10 −6 N·s / m²

ρ = 0.881 kg / m³
Propiedades del aire a 127.5°C (Tabla 4.4.1) 
k = 32.88·10 W / m·K
−3
Pr = 0.706

ρ u∞ x cr Re cr µ
Re cr = = 500000; x cr = = 0.217
µ ρ u∞
La transición se produce por tanto en el 5º elemento calentador.
El flujo de calor transferido en cada uno de los elementos será el siguiente:
qi = hi w L i (Ts − T∞ )
Por tanto, el flujo de calor será máximo allí donde el coeficiente de película
promedio sea máximo. Si recordamos como varía el coeficiente de película sobre
una placa plana flujo externo, concluimos que sólo existen tres posibilidades:
h(x )
h
39
1. Calentador nº 1: Corresponde al mayor coeficiente de convección local en

régimen laminar
2. Calentador nº 5: Se produce la transición al turbulento y aparece el mayor
coeficiente de convección local en régimen turbulento.
3. Calentador nº 6: Al ser todo el calentador turbulento puede ocurrir que el
promedio sea mayor que el anterior.
Calentador nº 1:
En este calentador la corriente es laminar y a temperatura superficial constante,

por tanto podemos usar la correlación Polhausen (nº 5, tabla 6.1):
 ρu L 
Nu1 = = 0.664 Re10.5 Pr 1 / 3 = 0.664 ∞ 1  Pr 1 / 3 = 200.64
0.5
h1 L 1
k  µ 
h1 = = 131.94 W / m²K
Nu1 k
q1 = h1 w L 1 (Ts − T∞ ) = 1352 W
L1
Calentador nº 5:
( )
El calor cedido en este elemento calefactor lo calcularemos por diferencia:
q5 = q1−5 − q1− 4 = h1−5 w L 1−5 (Ts − T∞ ) − h1− 4 w L 1− 4 (Ts − T∞ ) = h1−5L 1−5 − h1− 4L 1− 4 w (Ts − T∞ )
 ρu L 
Nu1− 4 = = 0.664 Re10−.54 Pr 1 / 3 = 0.664 ∞ 1− 4  Pr 1 / 3 = 401.29
0.5
h1− 4 L 1− 4
k  µ 
h1− 4 = = 65.97 W / m²K

Nu1− 4 k
L 1− 4
Al final del elemento quinto ya hemos entrado en zona turbulenta y por tanto
debemos usar otra correlación (nº 9, tabla 6.1):
( )
  ρ u∞L 1−5  
= = 0.037 Re1−5 − 871 Pr 
= 0.037  − 871 Pr 1 / 3 = 561.12
0.8

h1−5 L 1−5
 µ 
Nu1−5
 
0.8 1/3
 
k
h1−5 = = 73.80 W / m²K

Nu1−5 k
( )
q5 = q1−5 − q1− 4 = h1−5L 1−5 − h1− 4L 1− 4 w (Ts − T∞ ) = 1077 W
L 1 −5
Calentador nº 6:
De forma similar obtenemos la potencia disipada por el sexto elemento restando al
calor total entre el primer y el sexto elemento, el calor entre el primer y el quinto
( )
elemento:
q6 = q1−6 − q1−5 = h1−6L 1−6 − h1−5L 1−5 w (Ts − T∞ )
40
( )
  ρu L  
Nu1−6 = = 0.037 Re10−.86 − 871 Pr 1 / 3 =  0.037 ∞ 1−6  − 871 Pr 1 / 3 = 771.02
0.8
h1−6 L 1−6
  µ  
 
k
h1−6 = = 84.50 W / m²K

Nu1−6 k
( )
q6 = q1−6 − q1−5 = h1−6L 1−6 − h1−5L 1−5 w (Ts − T∞ ) = 1414 W
L 1− 6
Por tanto el elemento con máxima pérdida de calor es sexto elemento:
q6 = 1414 W > q1 = 1352 W > q5 = 1077 W
41
8 (Convección forzada flujo interno)

Se desea calentar 3 kg/s de agua desde 10°C hasta 66°C, manteniendo la temperatura de la
superficie interna de la tubería a 82°C. Si el diámetro interior de la tubería es de 5 cm,
determinar:
a. Longitud de tubería necesaria para alcanzar la temperatura requerida
b. Coeficiente de transferencia de calor en la superficie.
Solución:
Datos:
- Caudal de agua: m & = 3 kg / s
- Condiciones de entrada y salida del agua: Tm,ent = 10°C Tm,sal = 66°C
- Temperatura de la superficie interior del conducto: Tsup = 82°C
- Diámetro interior del conducto: Di = 0.05 m
Incógnitas:
a. Longitud de la tubería: L
b. Coeficiente de transferencia de calor en la superficie: h
Esquema:
3 kg/s 3 kg/s
10°C 66°C
T (°C)
82°C
66°C
10°C
L x(m)
Hipótesis:
Desarrollo:
a. Longitud de la tubería: L
Realizando un balance de energía sobre el volumen de agua podemos calcular el

calor ganado por esta:
42
& cp (Tm,sal − Tm,ent ) = 701.232 kW

q=m
entre la entrada y la salida 38°C cp = 4.174 kJ / kg·K (Tabla 4.5 de la colección de

Donde el calor específico del agua líquida se ha evaluado a la temperatura media
tablas, gráficas y ecuaciones de transmisión de calor)
La ecuación de transferencia para un conducto con temperatura superficial

constante dice:
q = h A ∆Tlm = h π Di L ∆Tlm
∆Tlm = DTLM =
∆Tent − ∆Tsal
=
(Tsup − Tm,ent ) − (Tsup − Tm,sal ) = 37.232°C
 ∆T   Tsup − Tm,ent 
ln ent  ln 
 ∆Tsal   Tsup − Tm,sal 
 
Calculo del coeficiente de película:
Las propiedades en las correlaciones de convección forzada flujo interno se evalúan

en la mayoría de los casos a la temperatura media de masas:
Tm,ent + Tm,sal
Tm,med = = 38 °C
2
µ = 678.6·10 −6 N·s / m²

ρ = 993 kg / m³
Propiedades del agua a 38°C (Tabla 4.5) 
k = 627.7·10 W / m·K
−3
Pr = 4.521

ReD = = 1.126·105
4m&
π Di µ
El régimen es claramente turbulento (mayor que 2300), Realizamos la hipótesis de

flujo completamente desarrollado L/D>10, que comprobaremos posteriormente.
Utilizando la correlación de Dittus-Boelter (correlación 27, Tabla 6.6):
NuD = = 0.023 ReD0.8 Pr 0.4 ; h = 5805 W / m²K

h Di
k
Volviendo a la ecuación de transferencia despejamos la longitud de tubería

necesaria:
L = = 20.65 m
q
h π Di ∆Tlm
L/D = 413, por tanto la hipótesis de flujo completamente desarrollado es válida.
Se propone utilizar alguna otra correlación válida para este caso y comparar los
resultados.
43
9 (Convección forzada, banco de tubos)

A menudo se dispone de agua presurizada a temperaturas elevadas, la cual se puede usar
para calefacción de locales o aplicaciones en procesos industriales. En tales casos es normal
usar un haz de tubos en el que el agua se hace pasar por éstos, mientras que también pasa
aire en flujo cruzado sobre ellos. Considérese una disposición de los tubos cruzada con un
diámetro exterior de los tubos de 16.4 mm y don los espaciados longitudinales y
transversales valen SL = 34.3 mm y ST = 31.3 mm respectivamente. Hay siete filas de tubos
en la dirección del flujo de aire y ocho tubos en cada una de las filas. En condiciones de
operación típicas la temperatura superficial de los tubos es de 70°C, mientras que la
temperatura del flujo de aire a contracorrientes es de 15°C y su velocidad 6 m/s. Determine
el coeficiente de convección del lado del aire y la transferencia de calor para el haz de
tubos.
Solución:
Datos:
- Condiciones del aire a la entrada: u∞ = 6 m / s T∞ = Tent = 15°C
- Temperatura superficial de los tubos: Tsup = 70°C
- Geometría del haz de tubos: De = 0.0164 m SL = 0.0343 m S T = 0.0313 m
- Nº de filas: 7
- Nº de tubos: 8
Incógnitas:
a. Coeficiente de película medio del lado del aire: h
b. Calor total transferido en el haz de tubos: q
Esquema:
Tsup = 70°C
6 m/s
15°C
Hipótesis:
44
Desarrollo:
Para obtener e coeficiente de película medio usaremos la correlación nº 17 de la

tabla 6.5:
NuD = C1 C2 ReDm,max Pr 0.36 (Pr/ Prs )

1/4
Donde todas las propiedades del fluido se evalúan a la temperatura del fluido
menos Prs que se evalúa a la temperatura de la superficie. Teóricamente
deberíamos evaluar las propiedades a la temperatura media del fluido entre la
entrada y la salida, comenzaremos por evaluarlas a la temperatura de entrada y
más tarde recalcularemos las mismas si la temperatura de salida es
sustancialmente diferente a la de entrada:
µ = 18.02·10 −6 N·s / m²

ρ = 1.225 kg / m³
Propiedades del aire seco a 15°C (Tabla 4.5) 
k = 24.76·10 W / m·K
−3
Pr = 0.7323

A la temperatura de la superficie 70°C: Prs = 0.7177
El Reynolds está basado en la máxima velocidad alcanzada por el aire:
SL
SD
D
ST
u∞
A1
A2
A1 = (S T − Dext ) w = (14.9 mm) w

La máxima velocidad se producirá en el lugar de sección de paso mínima:
S 
S = S +  T  ; SD = 37.7 mm
2
 2 
2 2
A 2 = 2 (SD − Dext ) w = (42.6 mm) w

D L
u∞ S T = v max (S T − Dext ); v max =

Por tanto la velocidad máxima se produce en la sección A1:
= 12.6 m / s
ST
S T − Dext
ρ v maxDext
ReD,max = = 14047
µ
45
Este valor está dentro de los límites de la correlación usada y si miramos las
constantes de la correlación valdrán:
C2 = 0.95 (Tabla 6.5.4)
= 0.915 < 2; C1 = 0.35(S T / SL ) = 0.3438; m = 0.6

ST 1 /5
(Tabla 6.5.3)
SL
NuD = 0.3438·0.95·ReD0.,6max Pr 0.36 (Pr/ Prs ) = 90.376

1/4
NuD = h = = 136.45 W / m²K

h Dext NuD k
;
k Dext
Flujo de calor:
& cp (Tsal − Tent )

Un balance de energía sobre la corriente de aire diría que:
q=m
El flujo de calor sobre un intercambiador a temperatura superficial constante puede

establecerse también como:
q = h A ext ∆Tlm
∆Tlm = DTLM =
∆Tent − ∆Tsal
=
(Tsup − Tent ) − (Tsup − Tsal )
=
Tsal − Tent
 ∆T   T − Tent   ∆T 
ln ent  ln sup  ln ent 
 ∆Tsal   Tsup − Tsal   ∆Tsal 
 
q = h A ext
Tsal − Tent & cp (Tsal − Tent )
=m
 ∆Tent 
ln 
 ∆Tsal 
 ∆T  h A ext ∆Tent
ln ent  = =e
h A ext
 ∆Tsal ∆Tsal
& cp
m

;
& cp
m
A ext = π Dext w Nt = 2.88·w m²

Área exterior de transferencia de calor:
& = ρ u∞ (8S T + Dext / 2) w = 1.9·w kg / s

Caudal de aire:
m
Tsal = 25.21°C
Luego la temperatura de salida valdrá:
& cp (Tsal − Tent ) = 19.53·w kW

Y por tanto: q = m
Deberíamos rehacer el problema calculando las propiedades a la temperatura media

T + Tent
del fluido: Tmed = sal = 20.1°C , pero el cambio es pequeño y el coeficiente de
2
película saldrá aproximadamente el mismo.
46
10 (Convección forzada + libre, conducto circular)

Por el interior de una tubería de 1” de diámetro y 100 m de longitud, circula agua
procedente de una caldera a una velocidad de 1.5 m/s. Calcular el espesor de aislamiento
necesario (Conductividad del aislante: k = 0.040 W/m·K), si la caída máxima de
temperatura permitida en el agua es de 0.5°C. La temperatura de salida del agua de la
caldera es de 90°C y el ambiente exterior se encuentra a 10°C.
Solución:
Datos:
- Condiciones del agua a la entrada: um,ent = 1.5 m / s Tm,ent = 90°C
- Dimensiones de la tubería: D = 1" = 0.0254 m L = 100 m
- Conductividad del aislante: k a = 0.040 W / m·K
- Máxima caída de temperatura en el agua: ∆Tm = Tm,ent − Tm,sal = 0.5°C
- Ambiente exterior: T∞ = 10°C
Incógnitas:
Espesor de aislante: e a
Esquema:
1.5 m/s
89.5°C
90°C
100 m
Sección media
Del conducto
qconv,int
Hipótesis:
- La tubería es de un espesor muy pequeño ( Dint = Dext )

47
Desarrollo:
Conocida la velocidad a la entrada calculamos el caudal másico de agua:
& = ρ um,ent π Dint = 0.7338 kg / s

2
m
Agua a presión atmosférica y 90°C, Tabla 4.5: ρ = 965.5 kg / m³
4
Si realizamos un balance de energía sobre el volumen de agua que circula por el

interior del conducto, obtenemos la potencia perdida por el agua en su
& cp (Tm,ent − Tm,sal ) = m

enfriamiento:
q=m & cp ∆Tm = 1.542 kW
Ese calor que pierde el agua es igual al calor que hacia el exterior atraviesa el
conducto, utilizando la analogía eléctrica para conducción en la sección media del
conducto (supondremos que la temperatura en la sección media es la media de
masas):
Tent + Tsal Tent + Tsal

− T∞ − T∞
ln(Dext / Dint ) ln(Dext / Dint )
q= 2 = 2 (1)
+ + + +
1 1 1 1
hint A int 2 π ka L hext A ext hint πDintL 2 π ka L hext πDextL
Las incógnitas en la ecuación anterior son: Dext , hint , hext

Podemos plantear dos ecuaciones más para cerrar el problema que son:
→ Nuint = f (Re, Pr )
Nuext = f (Gr, Pr )
hint
hext →
Con estas tres ecuaciones el problema queda cerrado (3 ecuaciones con 3

incógnitas).
Procedamos en primera instancia al cálculo del coeficiente de película interior:
Coeficiente de película interior:

Re = = 1.169 · 105
4m&
π Dint µ
Propiedades del agua a la temparatura media de masas 89.75°C (aprox. 90°C,
Tabla 4.5)
µ = 314.6·10 −6 N·s / m²

ρ = 965.5 kg / m³

k = 0.6755 W / m·K
Pr = 1.958

El régimen es turbulento ( ReD > 2300 ) con el flujo completamente desarrollado

( L / Dint = 3937 > 10 ), usaremos la correlación (27) de la tabla 6.6:
48
NuD = 0.023 ReD0.8 ⋅ Pr 0.3 ; ⇒hnti = 8479 W / m²K

hint Dint
k
Observamos que la resistencia asociada al coeficiente de película interior, se puede

considerar despreciable frente a la del aislante.
R cv,i = = 1.478 × 10 −5 K / W
1
hint A int
Para calcular el coeficiente de película externo es necesario conocer previamente la

temperatura de la superficie exterior del conducto y su diámetro exterior. Por ello
se hace necesario establecer un proceso iterativo de resolución:
1. Suponer un coeficiente de película exterior inicial hext ≈ 5 W / m²K

2. Calcular el diámetro exterior utilizando la ecuación (1)
3. Calcular la temperatura superficial exterior, planteando una ecuación de
transferencia que contenga dicha incógnita
Tent + Tsal
− Tsup
ln(Dext / Dint )
q= 2 (2)
+
1
hint A int 2 π ka L
′
4. Calcular el coeficiente de película exterior hext
′ ≈ hext se terminó el proceso iterativo en caso contrario volver al punto 2.
5. Sí hext
Dext = 0.075 m
El diámetro exterior calculado resolviendo de forma iterativa la ecuación (1) es:
Calculemos ahora la temperatura superficial exterior media despejando de la

ecuación (2): Tsup = 23.3°C .
Calcular el coeficiente de película exterior:
Podemos utilizar la correlación de Morgan (40, tabla 6.9) donde las propiedades
deben evaluarse a la temperatura media de película (23.3+10)/2 = 16.6°C
υ = 14.71·10 −6 kg/ms

k = 0.02476 W/mK
Propiedades del aire a 15°C 
Pr = 0.7323

β = = 0.0035
1
 16 .6 + 273.15
g β ∆T D3ext
RaD = Pr = 6.5198 x 105 , con este valor del número de Rayleigh la
υ2
correlación de Morgan toma el siguiente valor:
NuD = = 0.48 Ra1D/ 4 → hext = 4.503 W / m²K

hextDext
k
49
Es ligeramente inferior al que habíamos estimado, si volvemos a calcular el

diámetro exterior usando este coeficiente de película y la ecuación (1):
Dext = 0.073 m e a = 2.38 cm
50
2.3. Problemas resueltos de radiación

11 (Ecuaciones de intercambio radiante)
El proceso de fabricación de una superficie curva de absorción solar de área A2 = 15 m²,
consiste en la fijación de un recubrimiento especial sobre dicha superficie. Para fijar dicho
recubrimiento se cura mediante la exposición a un calentador infrarrojo de ancho W = 1 m.
El absorbedor y el calentador son cada uno de longitud L = 10 m y se encuentran
separados H = 1 m.
El calentador se encuentra a T1 = 1000 K y tiene una emisividad ε1 = 0.9, mientras que el

absorbedor está a T2 = 600 K y tiene una emisividad ε2 = 0.5. Todo el sistema se encuentra
en un local de grandes dimensiones cuyas paredes pueden considerarse a 300 K. ¿Cuál es la
transferencia neta de calor sobre la superficie de absorción?.
A2, T2, ε2
Superficie absorbedora
Paredes local
T3
A1, T1, ε1
Calentador
Solución:
Datos:
- Calentador (superficie 1): T1 = 1000 K ε1 = 0.9 A 1 = W L = 10 m²
- Superficie absorbedora (superficie 2): T2 = 600 K ε 2 = 0.5 A 2 = 15 m²
- Paredes del Local (superficie 3): T3 = 300 K
- Dimensiones: W = 1 m H =1m L = 10 m
Incógnitas:
Flujo de calor por radiación sobre la superficie absorbedora: q2
51
Esquema:
Superficie A’2
Hipótesis:
- Los flujos de calor por convección son despreciables
Las paredes del local pueden suponerse como una superficie negra ( ε 3 = 1 ) que
- Las superficies son difusas y grises
-
cierra el recinto a T3 = 300 K.
Desarrollo:
El problema se constituye como un recinto formado por tres superficies que

intercambian calor por radiación entre ellas, de las cuales conocemos sus
temperaturas. Podemos expresar el flujo de calor sobre 2 como:
M20 − J2 σ T24 − J2
q2 =
(1 − ε2 ) (1 − ε2 )
=
A 2 ε2 A 2 ε2
La única incógnita en la ecuación anterior es la radiosidad sobre la superficie 2, J2.

Para calcular todas las radiosidades plantearemos un sistema de ecuaciones donde
las radiosidades sobre cada una de las superficies sean las únicas incógnitas:
Ji = εi Mi0 + (1 − εi )∑ Fij J j
N
j=1
Si aplicamos la ecuación anterior a cada una de las superficies tendremos:
J1 = ε1 M10 + (1 − ε1 ) [F12 J2 + F13 J3 ]

J2 = ε2 M20 + (1 − ε2 ) [F21 J1 + F22 J2 + F23 J3 ]
J3 = M30 = σ T34
52
Debemos observar que la superficie 1 es plana y por tanto F11 = 0, y que sin
embargo la superficie 2 es convexa, F22 > 0.
Factores de forma:
Si observamos la geometría del problema veremos que los factores de forma F12 =
F12’. Debido a que toda la radiación que sale de 1 y llega a 2’, es la misma que
saliendo de 1 llega a 2.
Para calcular el factor de forma F12’ utilizaremos la gráfica 7.6 de la colección de
= 1/1 = 1. F12' = F12 = 0.4 .

tablas, gráficas y ecuaciones de transmisión de calor. Donde Y/L = 10/1 = 10 y X/L
Aplicando reciprocidad: A1F12 = A 2F21; F21 = F12 = 0.4 = 0.267

A1 10 m²
A2 15 m²
Aplicando adición: F13 = 1 − F12 = 0.6
Y por reciprocidad: F31 = F13 = 0.6 = 0.2727 ( A 3 ≈ 2 H L + 2 W L = 22 m² )

A1 10 m²
A3 22 m²
Por simetría podemos ver que F31 = F32' = F32 = 0.2727
Por último si aplicamos de nuevo reciprocidad y adición:
F23 = 3 F32 = 0.2727 = 0.4
A 22 m²
A2 15 m²
F22 = 1 − F21 − F23 = 1 − 0.267 − 0.4 = 0.333
Sustituyendo en las ecuaciones de radiación tendremos:
J1 = 51569 W / m²
J2 = 12778 W / m²
J3 = 459 W / m²
Y volviendo a la ecuación del flujo de calor sobre 2:
q2 = −81445 W = −81.445 kW
Este flujo es negativo debido a que la superficie 2 absorbe más radiación de la que
emite.
53
12 (Analogía Eléctrica + Fuente pequeña)

Un elemento cilíndrico de calentamiento de 10 mm de diámetro y gran longitud, se utiliza
en un horno tal como se indica en la figura. La temperatura del elemento es T1 = 1500 K y
difusa, con ε2 = 0,6 mantenida a una temperatura de 500 K. El resto de las superficies están
se puede considerar como un cuerpo negro. La superficie inferior A2 es una superficie gris
constituidas por un material refractario adiabático, con una emisividad de 0,9. La longitud
del horno en dirección normal es muy grande comparada con el ancho (w = 1 m) y la
altura (h = 0,87 m). Despreciando el intercambio por convección, determinar:
1. La potencia por unidad de longitud en W/m, que ha de suministrar el elemento de

calentamiento para mantener las condiciones de operación.
2. La temperatura de las paredes del horno.
60°
h = 0,87m
w=1m
Solución:
Datos:
- Calentador (superficie 1): T1 = 1500 K ε1 = 1 D1 = 0.01 m
A1 = π D1 L = 0.0314 L m²
- Superficie inferior (superficie 2): T2 = 500 K ε2 = 0.6 A 2 = w L = L m²
- Superficie refractaria (superficie 3): ε 3 = 0.9
- Dimensiones: w = 1 m h = 0.87 m α = 60°
Incógnitas:
1. Flujo de calor por radiación por unidad de longitud sobre la superficie 1: q1

(W/m).
2. Temperatura de la superficie 3: T3
54
Hipótesis:
- Los flujos de calor por convección son despreciables
La superficie refractaria es adiabática ( qcd,3 = 0 ) y puede considerarse el flujo de

- Las superficies son difusas y grises
-
calor por convección despreciable ( qcv,3 = 0 ), por tanto el flujo de calor por
radiación debe ser también nulo ( qrd,3 = 0 ). A este tipo de superficie se le suele
llamar “rerradiantes”.
Desarrollo:
El problema se constituye como un recinto formado por tres superficies que

intercambian calor por radiación entre ellas, de las cuales conocemos las
temperaturas de dos de ellas y el flujo de calor por radiación sobre la tercera. Este
tipo de recintos con 3 superficies donde una de ellas es rerradiante pueden
resolverse utilizando la analogía eléctrica.
1 − ε1 1 1 − ε2
ε1 A 1 A1F1,2 ε2 A 2
M20
M10
1 1
A1F1,3 A 2F2,3
q3 = 0
Podemos obtener el flujo de calor sobre la superficie 1 calculado la resistencia
( )
equivalente del circuito:
σ T14 − T24
q1 =
R eq
Cálculo de factores de forma:
Si consideramos la superficie 1, muy pequeña (un punto) frente a las superficies 2

y 3, y recordemos que el factor de forma no es más que la fracción de radiación
que saliendo de una superficie alcanza otra tendremos que:
60°
F12 = = 1 / 6 = 0.1667
360°
Y si aplicamos reciprocidad y adición:
F13 = 1 − F12 = 5 / 6 = 0.8333 F21 = F12 = 0.0052 F23 = 1 − F21 = 0.9948

A1
A2
55
Las resistencias por tanto valdrán:
R eq = 33.193 m K / W
190.948 0.6667
0
M20 M10 M20
M10
38.199 1.0052
q3 = 0
q1 =
(
σ T14 − T24 )
= 8541 W / m = 8.541 kW / m
R eq
Para calcular la temperatura de la superficie 3 sabemos que:
q3 =
A 3 ε3 0
1 − ε3
(
M3 − J3 = 0; ) J3 = M30 = σ T34
q3 = A 3 ∑ F (J − J j ) = 0; A 3F31 (J3 − J1 ) + A 3F32 (J3 − J2 ) = 0

3
3j 3
j=1
Para calcular la radiosidad sobre la superficie 3 y por tanto su temperatura

debemos conocer primero la radiosidad sobre las otras superficies:
J1 = M10 = σ T14
q2 = −8541 W / m =
1 − ε2
(
A 2 ε2 0
M2 − J2 ; ) J2 = 9238 W / m²
A1F13 J1 + A 2F23 J2
J3 = J3 = 10693 W / m² = σ T34
A1F13 + A 2F23
;
Por tanto:
T3 = 656 K
Comentarios:
Si consideramos que una superficie es muy pequeña comparada con las restantes y
que además se encuentra a una temperatura mucho mayor que el resto (bombillas,
resistencias de calentamiento de tamaño pequeño etc...), podemos suponer que la
energía que emiten es mucho mayor que la que absorben y por tanto que esta es
independiente de las superficies que la rodean y función únicamente de su
temperatura y propiedades radiantes:
q1 = A 1ε1M10 = A1σT14 = 9018 W / m
Error cometido: 5.6 %
56
2.4. Problemas resueltos de mecanismos combinados

13 (Mecanismos combinados, radiación + conducción)
Un conducto circular con los gases de escape de una caldera cruza verticalmente una
habitación. La temperatura superficial del metal que constituye la tubería es de 200°C. Una
pared de 2 m de ancho, con baja conductividad térmica (aislada) y negra a efectos radiantes
se encuentra situada próxima a la tubería (ver figura 1). El efecto radiante de la tubería
produce grietas por dilataciones en la pared anteriormente mencionada. Para subsanar este
problema se debe mantener la pared a una temperatura no superior a 30°C y se plantean
dos soluciones posibles:
1. Aislar exteriormente la tubería con una capa de aislante de conductividad 0.05

W/m·K negro en su cara exterior. Calcular el espesor de aislante necesario para
mantener la pared a 30°C.
2. Se coloca una placa metálica (muy alta conductividad) a modo de pantalla de
considerar negra y la enfrentada a la tubería tiene una emisividad ε. La pared

radiación tal como muestra la figura 2. La cara enfrentada a la pared se puede
desnuda de la tubería puede considerarse también negra. Determinar ε para que la

pared se siga manteniendo a 30°C.
El resto de la habitación puede suponerse un cuerpo negro a 20°C, y el aire se encuentra a

la misma temperatura.
1m 1m 1m 1m
30°C 30°C 0,25 m

0,5 m
0,25 m
0,5 m 0,5 m
200°C
200°C
Taire = 20°C Taire =20°C

20°C 20°C
Figura 1 Figura 2
Notas:
- Considerar todos los coeficientes de película iguales a 10 W/m²K (tubería y pared),
menos los de ambas caras de la placa metálica en el apartado 2 que serán iguales a 5
W/m²K
- Para el cálculo de factores de forma considerar, en el apartado 1, despreciable el espesor
de aislamiento
57
Solución:
Datos:
ε1 = 1 k 1 = 0.05 W / m·K
Apartado 1:
- Superficie exterior del aislamiento (1):
- Temperatura superficial del conducto: Tsup = 200°C
- Pared adyacente (2): ε 2 = 1 T2 = 30°C qCD,2 = 0 (aislada)
- Aire y resto de la habitación (3): Taire = 20°C Trm,hab = T3 = 20°C ε3 = 1
- Todos los coeficientes de película: h1 = h2 = 10 W / m²K
ε1 = 1 T1 = 200°C
Apartado 2:
- Superficie exterior de la tubería (1):
- Pared adyacente (2): ε 2 = 1 T2 = 30°C qCD,2 = 0 (aislada)
- Aire y resto de la habitación (3): Taire = 20°C Trm,hab = T3 = 20°C ε3 = 1
- Cara enfrentada al conducto de la placa metálica (4): h4 = 5 W / m²K
- Cara enfrentada a la pared de la placa metálica (5): ε5 = 1 h5 = 5 W / m²K
- Resto de coeficientes de película: h1 = h2 = 10 W / m²K
Los datos geométricos se tomarán de las figuras del enunciado
Incógnitas:
- Apartado 2: Emisividad de la cara enfrentada a la tubería de la placa metálica, ε 4 ,

- Apartado 1: Espesor de aislamiento, e1, para mantener la pared adyacente a 30°C
para mantener la pared adyacente a 30°C
Hipótesis:
- Conducción unidimensional en el aislamiento exterior de la tubería para el
apartado 1
- Superficies isotermas, grises y difusas para el intercambio radiante
- Suponer despreciable el espesor de aislamiento al calcular los factores de forma
del recinto
Esquema:
2
1 qconv,2
qrad,1
qcond,1 qrad,2
qconv,1
58
Desarrollo:
Apartado 1:
Balance de energía sobre la superficie exterior del aislamiento (1) y sobre la pared
adyacente (2) de acuerdo con los sentidos de los flujos dibujados en el esquema:
qrad,1 + qconv,1 + qcond,1 = 0 (1)

qrad,2 + qconv,2 = 0 (2)
Las ecuaciones de transferencia de cada uno de los flujos de calor son las
siguientes:
qconv,2 = h2 A 2 (T2 − Taire ) = 10·2 L (30 − 20) = 200 L [W]
qconv,1 = h1 A1 (T1 − Taire ) = 10 A1 (T1 − 20) = 10 π (D + 2e1 )(T1 − 20)
qcond,1 =
(T − Tsup )
=
2π 0.05 L
(T1 − 200)
 D + 2e1   0.5 + 2e1 
1
ln  ln 
 D   0.5 
2πk 1L 2πk 1L
[
qrad,1 = A1 M10 − (F1,2 J2 + F1,3 J3 ) ]
qrad,2 = A [M − (F
2
0
2 J + F2,3
2,1 1 J )]
3
Como todas las superficies son negras las radiosidades son de muy fácil cálculo:
Ji = Mi0 = σ Ti4
[ ( )]
Luego:
qrad,1 = A1σ T14 − F1,2 T24 + F1,3 T34

qrad,2 = A σ[T − (F
2
4
2
4
2,1 1T + F2,3 T )]
4
3

Tabla 7.4: S1 = −1 m S2 = 1 m L = 0.75 m r = 0.25 m
0.25  −1 
F2,1 =  − tan−1 = 0.23 F2,3 = 1 − F2,1 = 1 − 0.23 = 0.77
0.75 
tan−1
1
2  0.75
F1,2 = F2,1 = 0.23 = 0.29 F1,3 = 1 − F1,2 = 1 − 0.29 = 0.71

A1 2L
A2 π 0.5 L
En la ecuación (2) la única incógnita es T1, si despejamos:
T1 = 374.4 k = 101.4°C
Volviendo a la ecuación 1 la única incógnita es el espesor de aislamiento, aunque la

ecuación no es lineal. Si resolvemos esta ecuación de forma iterativa obtendremos:
59
e1 = 0.0033 m = 3.3 mm
Apartado 2:
Esquema:
2
qrad,5 qconv,2
qconv,5 qrad,2
5
4 qconv,4
qrad,4
1
qrad,1 qconv,1
qcond,1
Si tomamos el balance de energía sobre la superficie 2:
qrad,2 + qconv,2 = 0
El flujo de calor por convección es el mismo del apartado anterior pero el de
[ ]
radiación es diferente al aparecer nuevas superficies en el recinto.
qrad,2 = A 2 σT24 − (F2,1 J1 + F2,3 J3 + F2,4 J4 + F2,5 J5 )
La superficie 2 “no ve” a las superficies 1 y 4 y por tanto los factores de forma:
F2,1 = 0 F2,4 = 0
[ (
qrad,2 = A 2 σ T24 − + F2,3 T34 + F2,5 T54 )]
[(8 + 8) ] [ ]
Tabla 7.4: w1 = w2 = 2 m L = 0.25 m W1 = w1 / L = W2 = w2 / L = 8
+4 − (8 − 8) + 4
F2,5 = F5,2 = = 0.88 F2,3 = 1 − F2,5 = 1 − 0.88 = 0.12
2 1/2 2 1/2
2·8
Tabla 7.4: S1 = −1 m S2 = 1 m L = 0.5 m r = 0.25 m

0.25  −1
F4,1 = − tan−1 = 0.28 F4,3 = 1 − F4,1 = 1 − 0.28 = 0.72
2  0.5 
tan−1
1
0.5
En la ecuación de balance sobre la superficie la única incógnita es T5, si

despejamos:
60
T5 = 320.7 k = 47.7°C
La temperatura de la superficie 5 debe ser igual a la de la superficie 4 puesto que

ambos superficies son constituyen las dos caras de una placa metálica de muy alta
conductividad.
Si planteamos a continuación un balance de energía sobre toda la placa, superficie

4 más superficie 5:
qrad,4 + qconv,4 + qrad,5 + qconv,5 = 0
siguientes:
qconv,4 = qconv,5 = h A 4 (T4 − Taire ) = 10·2 L (47.7 − 20) = 554 L [W]

[
qrad,4 = A 4 ε 4M04 − ε 4 (F4,1J1 + F4,3 J3 )]
qrad,5 = A 5 [M − (F
0
5 J + F5,3 J3 )
5,2 2 ]
La única incógnita en la ecuación de balance anterior es la emisividad de la
superficie 4, que si sustituimos:
ε 4 = 0.82
61
14 (Mecanismos combinados, Analogía eléctrica + h rad.)

El aire de un local acondicionado se encuentra a una temperatura de 20°C, en él se
encuentra situado un radiador (superficie 1) de emisividad superficial 0.9, y una pared
conectada con el exterior (superficie 2), cuya cara interior y exterior tienen un emisividad de
0.8. La cara interior del muro (superficie 2) tiene un coeficiente de película de hint = 3
W/m²K y en su cara exterior intercambia calor por convección con el aire ambiente a 10°C
(el coeficiente de película exterior puede considerarse de hext = 10 W/m²K) y por radiación
con los alrededores que también pueden considerarse a una temperatura de 10°C. Las
restantes superficies del local, excepto el radiador y el muro exterior deben considerarse
rerradiantes.
Calcular la temperatura superficial del radiador para que las pérdidas de calor por
conducción en el muro conectado con el exterior no superen el 80% del calor total emitido
radiantemente por el radiador.
La composición del muro exterior es:
Espesor Conductividad
Capa (cm) (W/mK)
Enlucido de cemento 2 1.0
Ladrillo macizo 25 0.7
Enlucido de cemento 2 1.0
1.5 m
1 1.5 m
Text = 10°C
2 Tint = 20°C
3m
3m
Solución:
Datos:
- Temperatura del aire del local: Tint = 20°C
62
- Temperatura del aire exterior: Text = 10°C

- Superficie 1 (radiador): ε1 = 0.9 A1 = 3 m · 1.5 m = 4.5 m²
Cara interior (superficie 2): ε 2 = 0.8 h2 = 3 W / m²·K A 2 = 9 m²

- Pared al exterior:
Cara exterior: ε ext = 0.8 hext = 10 W / m²·K A ext = A 2

- Temperatura media radiante del exterior: Tmr,ext = 10°C
- Las restantes superficies del local son rerradiantes.
e1 = 0.02 m; k 1 = 1.0 W / m·K

- Composición del muro exterior: e2 = 0.25 m; k 2 = 0.7 W / m·K
e3 = 0.02 m; k 3 = 1.0 W / m·K
Incógnitas:
Temperatura de la superficie 1 (radiador), para que el flujo de calor por conducción
en el muro sea el 80% del calor emitido radiantemente por 1.
T1 ? → qcond = 0.8 qrad,1
Hipótesis:
- Conducción unidimensional en el muro exterior
Esquema:
qrad,ext qrad,2
qconv,ext qcond qconv,2
Desarrollo:
Balance de energía sobre las superficies interior (2) y exterior (s,ext) de acuerdo
con los sentidos de los flujos dibujados en el esquema:
qconv,2 = qcond + qrad,2 (1)

qcond = qconv,ext + qrad,ext (2)
siguientes:
qconv,2 = h2 A 2 (Tint − T2 )
63
qcond = A 2
(T − Ts,ext )
= A2
(T2 − Ts,ext )
∑
2
+ +
ej e1 e2 e3
k k k2 k3
qconv,ext = hext A 2 (Ts,ext − Text )

capas j 1
qrad,ext = hrad A 2 (Ts,ext − Trm,ext );

Ts,ext + Trm,ext
donde hrad = 4 σ Tm3 ε ext ; Tm =
2
La temperatura media para el cálculo del h radiante debe estar cercana a 10°C
comencemos por suponer este valor y retomaremos esta hipótesis al final del
problema para su revisión: hrad = 4.1189 W / m²·K
Para calcular el flujo de calor por radiación en 2 debemos resolver el intercambio

radiante en el local interior, por ser un recinto con 3 superficies con una de ellas
rerradiantes podemos utilizar la analogía eléctrica (la superficie 3 está constituida
por el resto de las superficies del recinto que no son ni 1 ni 2).
1 − ε1 1 − ε2
( )
1
ε1 A 1 ε2 A 2
σ T24 − T14
A1F1,2
qrad,2 =
R eq
M20
M10
1 1
A 1F1,3 A 2F2,3
Cálculo de los factores de forma:
A 2F2,1 + A 2F2,A = A 2F2,1A

A
F2,1 = F2,A F2,1 =
F2,1A
2
W / D = 1; L / D = 1 → F2,1A = 0.2
De la Gráfica 7.7
F2,1 = = 0.1
F2,1A
2 1
= 2 F2,1 = 0.2
A
F1,2
A1
F1,3 = 1 − F1,2 = 0.8
F2,3 = 1 − F2,1 = 0.9
2
Luego:
0.0247 1.111 0.0278
M20
M10
0.2778 0.1235
64
R eq = 0.0247 + + 0.0278 = 0.3473 K / W

1
+
1 1
1.111 (0.2778 + 0.1235)
qrad,2 =
(
σ T24 − T14 )
[W]
0.3473
Si sustituimos las ecuaciones de transferencia en las dos ecuaciones de balance (1)

y (2), tendremos un sistema de dos ecuaciones con 3 incógnitas (las tres
temperaturas superficiales: T1 , T2 y Ts,ext ).
h2 A 2 (Tint − T2 ) = A 2
(T (
− Ts,ext ) σ T24 − T14
+
)
∑
2
(1)
ej R eq
(T − Ts,ext )
capas k j
= hext A 2 (Ts,ext − Text ) + hradA 2 (Ts,ext − Trm,ext )

∑
2
A2 (2)
ej
capas kj
Para cerrar el sistema necesitamos otra ecuación que la obtenemos de la condición
( )
impuesta por el problema:
σ T14 − T24
qcond = 0.8 qrad,1 ; qcond = 0.8 (3)
0.3473
( )
Sustituyendo:
27 (20 − T2 ) = 9
(T2 − Ts,ext )
+
σ T24 − T14
(1)
0.3971 0.3473
(T 2 − Ts,ext )
= 10 (Ts,ext − 10 ) + 4.1189 (Ts,ext − 10 ) (2)
0.3971
Si introducimos las ecuaciones (3) en la (1):
 (T − Ts,ext )
27 (20 − T2 ) = qcond − qcond = −0.25 9 2 
1
0.8  0.3971 
Con esta ecuación y la ecuación (2) tenemos un sistema lineal de dos ecuaciones
con dos incógnitas que resolvemos:
T2 = 22.167 °C; Ts,ext = 11.842 °C
En este punto deberíamos recalcular el coeficiente de película radiante con la nueva

temperatura media, si lo hacemos: hrad = 4.1592 W / m²·K observamos que el
cambio es muy pequeño y que no merece la pena rehacer todos los números.
Volviendo a la ecuación (3) obtenemos la temperatura superficial del radiador:
T1 = 311.354 K = 38.204 °C
65
Comentarios:
Al plantear el balance de energía sobre la superficie (2) podríamos haber expresado

el flujo de calor por conducción como:
qcond = A 2
(T2 − Text )
∑ +
ej 1
capas k j hext + hrad
Eliminando la temperatura de a superficie exterior como incógnita. Esto es posible

debido a que en este problema la temperatura exterior del aire y la temperatura
radiante media exterior coinciden y a que el intercambio radiante puede expresarse
usando un coeficiente de película radiante.
66
15 (Mecanismos combinados, Analogía eléctrica + Bal. aire)

Una fila de elementos de calentamiento cilíndricos regularmente espaciados (superficie 1)
se usa para curar un recubrimiento superficial que se aplica a una de las caras a una lámina
de metal (superficie 2) colocado por debajo de los elementos. La lámina tiene un espesor de
1 cm y una conductividad k = 30 W/m·K. Un segundo panel (superficie 3), cuya superficie
negros y se mantienen a T1 = 600 K, mientras que la lámina tiene una emisividad ε2 = 0.5 y
superior está bien aislada, se coloca por encima de los elementos. Los elementos son
se mantiene a T2 = 400 K. La cavidad se llena con un gas no participativo teniendo lugar

una transferencia de calor convectiva en las superficies 1 y 2, con h1 = 10 W/m2K y
h2 = 2 W/m2K . (La convección en el panel aislado se puede despreciar).
1. Evalúe la temperatura media del gas, Tg .
2. ¿Qué potencia eléctrica por unidad de longitud axial debe suministrársele a cada
elemento para mantener la temperatura establecida?
3. ¿Cuál es el coeficiente de transferencia de calor por convección en la cara inferior

de la lámina de metal (superficie 4) (ε4 = 0.5) si la temperatura del aire en contacto
directo con él y la temperatura media radiante de los alrededores es Ta = 300 K?.
3
s = 50 mm
1
D = 25 mm
Gas, Tg 2
4
Aire, Ta
Solución:
Datos:
- Superficie 1: T1 = 600 K ε1 = 1 h1 = 10 W / m²K
- Superficie 2: T2 = 400 K ε 2 = 0.5 h2 = 2 W / m²K k 2 = 30 W / m·K e2 = 1 cm
- Superficie 3: Aislada → qcond = 0; Convección despreciable → qconv = 0
- Superficie 4: ε 4 = 0.5
- Aire en contacto con la superficie 4: Ta = Trm,alr = 300 K
Los datos geométricos se encuentran en la figura del enunciado
Incógnitas:
- Apartado 1: Temperatura del gas en contacto con las superficies 1, 2 y 3: Tg
q1′
- Apartado 2: Potencia eléctrica por unidad de longitud [W/m] para cada cilindro:
67
- Apartado 3: Coeficiente de transferencia de calor en 4: h4
Hipótesis:
- Conducción unidimensional entre las superficies 2 y 4
- Al existir un número elevado de elementos de calentamiento, estudiaremos el
comportamiento de un elemento intermedio como representativo del problema.
Esquema:
3
qCV1
1
qCV2 2
Desarrollo:
Apartado 1:
Tomamos un volumen de control que contenga solamente un elemento de

calentamiento, y realizamos un balance de energía sobre el gas de este volumen
(régimen permanente)
qCV1 + qCV2 = 0
h1 A1 (T1 − Tg ) + h2 A 2 (T2 − Tg ) = 0
Despejando la temperatura del gas de la ecuación anterior tenemos:
h1A1 T1 + h2 A 2 T2 10·π·0.025·600 + 2·0.05·400

Tg = = = 577.411 K
h1A1 + h2 A 2 10·π·0.025 + 2·0.05
Tg = 577.411 K
qRD1
Apartado 2:
qCV1
Balance de energía sobre la superficie 1: q1

q1′ = q′CV1 + q′RD1
q1′ : Potencia eléctrica consumida por cada elemento por unidad de longitud axial
q′CV1 = h1 πD1 (T1 − Tg ) = 10·π·0.025·(600 − 577 .411) = 17.741 W/m
Para calcular el flujo de calor por radiación sobre la superficie 1 utilizamos la

analogía eléctrica, ya que la superficie 3 puede considerarse rerradiante (aislada y
68
con convección despreciable). En el cálculo de las resistencias utilizaremos las áreas
( )
partidas por L.
1 − ε1 1− ε2
σ T14 − T24
1
q′RD1 = ε1 A1 A2 F2,1 ε 2 A2
R eq
M20
M10
1 1
A3 F3,1 A2 F2 ,3
Factor de forma entre una placa y un conjunto infinito de cilindros equidistantes

(Tabla 7.4):
 D  D −1  s − D 
+   tan  
1/2
= 1 − 1 −    = 0.6576
2 1/2

2 2
  s   s  D 
F2,1 2
F2,2 =0 F2,3 = 1 − 0.6576 = 0.3424 F3,1 = F2,1 = 0.6576
Luego la analogía eléctrica queda como:
0 30.4136 20
M20
M10
30.4136 58.4112
R eq = 42 .6561 K/W
q′RD1 =
(
σ T14 − T24
=
) (
5.67 e − 8 600 4 − 400 4 )
= 138.24 W/m
42.6561 42.6561
Volviendo a la ecuación de balance sobre 1:
q1′ = 17.741 + 138.24 = 155.981 W/m qRD2
qCV2
Apartado 3:
qCD
Balance de energía sobre 2 y 4:
qCV4
qCV2 + qRD2 + qCD = 0
qRD4
qCV 4 + qRD4 = qCD
Del balance sobre la superficie 2 podemos calcular la temperatura de la superficie 4

al ser la única incógnita:
69
h2 A 2 (T2 − Tg ) +
(
σ T24 − T14) (T − T4 ) = 0
+ 2
R eq e
(400 − T4 ) = 0 ;
kA 4
2·0.05(400 − 577.411) − 138.24 +

0.01
30·0.05
T4 = 398.96 K ; q′CD = 155.981 W/m
( )
Y del balance sobre la superficie 4 despejamos el coeficiente de película:
h4 A 4 (T4 − Ta ) + ε 4 A 4 σ T44 − Ta4 = 155 .981 W/m
( )
h4 0.05(398.96 − 300 ) + 0.5·0.05·5.68e − 8 398.96 4 − 300 4 = 155.981 W/m
h4 = 26.587 W/m²K
70
16 (Mecanismos combinados)
Un conjunto de placas como las de la figura se encuentran en el interior de un túnel de
calentamiento muy largo, la parte alta del túnel es una superficie calefactora a 80°C
(superficie 3) y puede considerarse negra a efectos radiantes. La superficie inferior
(superficie 4) es también negra y se encuentra a 20°C. El aire fluye sobre las placas a 20°C.
1. Calcule el coeficiente de película medio (supuesto el mismo en ambas caras) sobre

las placas si la cara superior (superficie 1) se encuentra a 27°C (negra) y la inferior
(superficie 2) a 25°C ( ε2 = 0,8 )
2. Calcule la velocidad del aire sobre la cara superior de la placa si ésta se encuentra en
convección forzada
3. Las placas están compuestas de dos capas de materiales diferentes de espesor 1mm,
la capa superior tiene una conductividad de 1W/m·K. Calcule la conductividad del
material inferior.
Nota: Considerar que las capas límites que se desarrollan sobre las superficies de las placas
no se interfieren entre sí
T3 = 80°C
......
T2 = 25°C
Taire = 20°C
H = 20 cm
T1 = 27°C e1 = 1mm
e2 = 1mm
...
w = 20 cm
T4 = 20°C
71
Solución:
Datos:
- Superficie 1: T1 = 27°C e1 = 1 mm k 1 = 1 W / m·K ε1 = 1
- Superficie 2: T2 = 25°C e2 = 1 mm ε 2 = 0.8
- Superficie 3: T3 = 80°C ε3 = 1
- Superficie 4: T4 = 20°C ε4 = 1
- Aire que fluye sobre las placas: Taire = 20°C

Dimensiones en la figura del enunciado
Incógnitas:
- Apartado 1: Coeficiente de película medio sobre las placas uno y dos h = h1 = h2 .
- Apartado 2: Velocidad del aire antes de llegar a las placas u∞ .
- Apartado 3: Conductividad del material inferior de las placas k 2 .
Hipótesis:
- Conducción unidimensional en las placas
- Al existir un número elevado de placas, estudiaremos el comportamiento de una
placa intermedia como representativo del problema.
Esquema:
qconv,1 qrad,1
qconv,2 qrad,2
Desarrollo:
Apartado 1:
Realizaremos un balance de energía sobre toda la placa (superficie superior más

superficie inferior:
qrad,1 + qconv,1 + qrad,2 + qconv,2 = 0 (1)
siguientes:
qconv,1 = h1 A1 (T1 − Taire ) = h w L (T1 − Taire ) = 1.4 h L
72
qconv,2 = h2 A 2 (T2 − Taire ) = h w L (T2 − Taire ) = h L




[
qrad,1 = A1  ε1M10 − ε1 ∑ F1, j J j  = A1 ε1σT14 − ε1 (F1,2 J2 + F1,3 J3 ) ]
 
∑F [
J  = A 2 ε 2 σT24 − ε 2 (F2,1 J1 + F2,4 J4 ) ]
j
 
qrad,2 = A 2  ε 2M20 − ε 2
 
 
2, j j
j
Planteamiento del intercambio radiante en el recinto:

Plantearemos una ecuación de radiosidad por cada superficie
J1 = ε1M10 = σT14 = 460 .189 W/m²

J2 = ε 2M20 + (1 − ε2 )∑ F2, j J j = ε2 σT24 + (1 − ε2 )(F2,1J1 + F2,4 J4 )
j
J3 = ε 3M = σT = 881 .901 W/m²

0 4
J4 = ε 4M04 = σT44 = 418 .738 W/m²

3 3
W1 = w1 / H = 0.2 / 0.2 = 1 W2 = w2 / H = 0.2 / 0.2 = 1

Placas paralelas con las líneas medias en la misma perpendicular (Tabla 7.4)
F1,2 =
[(W
1 + W2 ) + 4
2
]
1/2
[
− (W2 − W1 ) + 4
2
]
1/2
=
[(1 + 1) 2
+4 ]
1/2
[
− (1 − 1) + 4
2
]
1/2
= 0.4142
2W1 2·1
F2,1 = F1,2 = 0.4142

F1,3 = 1 − F1,2 = 0.5858
F2,4 = 1 − F2,1 = 0.5858
( )
Volviendo a la ecuación de la radiosidad sobre 2:
J2 = ε2 σT24 + (1 − ε2 )(F2,1J1 + F2,4 J4 ) = ε2 σT24 + (1 − ε2 ) 0.4142 σT14 + 0.5858 σT44 = 445.618 W/m²
[ ] [ ( )]
Así los flujos de calor radiantes valdrán:
qrad,1 == A1 ε1σT14 − ε1 (F1,2 J2 + F1,3 J3 ) = wL σT14 − F1,2 J2 + F1,3 σT34 = −48.2 L
[ ] [ (
qrad,2 = A 2 ε 2 σT24 − ε 2 (F2,1 J1 + F2,4 J4 ) = wL σT24 − F2,1σT14 + F1,3 σT44 = 1.94 L )]
Sustituyendo en la ecuación de balance (1):
−48.2 L + 1.4 h L + 1.94 L + h L = 0
Despejando de aquí el coeficiente de película obtenemos:
h = 19.27 W/m²K
73
Apartado 2:
Debido al bajo valor del coeficiente de película supondremos que este se encuentra
en convección forzada flujo externo (placa plana, flujo paralelo), régimen laminar.
Por tanto la correlación apropiada es la correlación de Polhausen (correlación 5 de
la tabla 6.1), está hipótesis deberá ser revisada si una vez calculado el nº de
Reynolds, éste es mayor que 500000 y por lo tanto debe cambiarse a una
correlación turbulenta.
Las propiedades se evalúan a la temperatura media de película:

+ T1
Tmp = aire = 23.5°C
T
2
µ = 18.42·10 −6 kg/m·s

ρ = 1.19 kg / m³
Propiedades del gas (aire) a 23.5°C (Tabla 4.4.1) 
k = 25.4·10 W/mK
−3
Pr = 0.73

NuL = = 151.7323 = 0.664 Re1L / 2 Pr 1 / 3

hw
k
Despejando Reynolds: ReL = 6.4408·10 4 < 5·105
Luego la hipótesis de régimen laminar era correcta.
ρ u∞ w
Despejemos la velocidad de la expresión del Reynolds: ReL =
µ
u∞ = 4.985 m/s
Apartado 3:
qconv,1 qrad,1
qcond
qconv,2 qrad,2
Realizando un balance sobre cualquiera de las dos caras.
qcond = −qconv,1 − qrad,1 = qconv,2 − qrad,2 = 21.22 L
Utilizando la ecuación de transferencia de conducción.

T − T2
qcond = 21.22 L = w L 1
+
e1 e2
k1 k 2
k 2 = 0.056 W/m·K
74
En una caldera para calentamiento de agua los tubos de humos tienen un diámetro de 7 cm
y una temperatura superficial igual a 385 K. A través de estos conductos circulan los gases
de escape a una temperatura de masa de 900 K. Para mejorar la transferencia de calor desde
el gas al agua, se coloca una pared delgada de alta conductividad en el plano medio del
conducto. Las propiedades físicas de los gases pueden aproximarse usando las del aire.
Tubo, Ts, εs
&g
Gases, Tm,g, m
Partición, Tp, εp
4. Sin la partición y con un flujo de gases igual a 0.05 kg/s. ¿Cuál es la transferencia de
calor por unidad de longitud del conducto?.
5. Para el mismo flujo de gases que en el caso anterior y colocando la partición (la
emisividad de la superficie interior del conducto y de la partición es 0.5).
Determinar la temperatura de la partición y la transferencia de calor por unidad de
longitud del conducto.
6. Explicar físicamente a que se debe el aumento de la transferencia de calor.

c. ¿Será así para otros valores de la emisividad? Razónelo. ¿Qué valor de la
emisividad hace máxima la transferencia?
d. Si el caudal variara, ¿aumentaría siempre la transferencia de calor por
convección?. Cuando así fuere ¿en qué proporción?
Solución:
Datos:
- Caudal másico de gas: & g = 0.05 kg / s
m
- Diámetro interno: Di = 0.07 m
- Superficie interna del conducto: Ts = 385 K ε s = 0.5
- Temperatura de masa del fluido: Tm,g = 900 K
- Partición de conducto: εp = 0.5
Incógnitas:
exterior sin partición: q′sin

- Apartado 1: Flujo de calor por unidad de longitud perdido por el conducto hacia el
75
- Apartado 2: Temperatura de la partición: Tp , y el flujo de calor por unidad de

longitud perdido por el conducto hacia el exterior con partición: q′con
- Apartado 3: Explicación física del aumento de transferencia de calor
Hipótesis:
- Temperatura superficial uniforme y constante en todos los apartados
- En el apartado 1 no existe intercambio de calor por radiación.
- Superficies isotermas, grises y difusas para el intercambio radiante del 2º
apartado
- Flujo completamente desarrollado
Esquema:
q′con = q′cond q′conv

q′sin = q′cond
q′conv
q′rad
Desarrollo:
Apartado 1:
Desde la pared de la tubería al fluido solo hay convección ya que al estar toda la
superficie interior del conducto a una temperatura uniforme el intercambio de calor
por radiación es nulo. El calor transferido por convección desde el fluido a la pared
será igual al calor perdido por conducción a través de la pared del conducto.
q′sin = q′conv = A i hi / L (Tm,g − Ts ) = π Di hi (Tm,g − Ts )

Pérdidas de calor por unidad de longitud:
Cálculo del coeficiente de película, h:

µ = 39.3·10 −6 kg / ms

ν = 100·10 m / s
−6
Propiedades del gas (aire) a 900 K (626.85°C) 
2
k = 62.3·10 W / mK
−3
(Tabla 4.4.1)
Pr = 0.706

Nº de Reynolds para conducto circular: ReD = = 23141 .4

4m &g
π Di µ
El flujo es claramente turbulento y supondremos que el flujo es completamente

desarrollado. Utilizaremos la correlación de Dittus-Boelter (nº 27, tabla 6.6) con un
exponente n = 0.3 puesto que la superficie está más fría que el gas.
NuD = 0.023 ReD4 / 5 Pr 0.3 = 64.2511
76
hi = = 57.184 W / m2K
NuDk
D
Luego: q′sin = π Di hi (Tm,g − Ts ) = 6476.3 W / m
Apartado 2:
& g = 0.025 kg / s ). Las temperaturas superficiales y de masa siguen siendo las

El caudal másico por cada una de las dos mitades del conducto se reduce a la mitad
(m
mismas del apartado anterior.
En este segundo apartado aparece un nuevo flujo de calor por radiación debido a
que insertamos una partición que se encuentra a una temperatura diferente, Tp, de
la superficie del conducto. Para determinar esta temperatura será necesario
plantear un balance de energía sobre la partición.
Cálculo del nuevo coeficiente de película, h (las propiedades del gas son las mismas
del apartado anterior):
Tomando la mitad del conducto tendremos un conducto no circular con el siguiente

diámetro equivalente
πDi2
π
De = = = Di = 0.043 m
4
4A
πDi π+2
8
+ Di
P
2
Si calculamos Reynolds en este nuevo conducto no circular, tendremos:

ρ umDe
ReD = = = 14215
& g De
m
µ AT µ
Podemos utilizar la misma correlación del apartado anterior pero usando el

diámetro equivalente:
NuD = 0.023 ReD4 / 5 Pr 0.3 = 43.508 hi = = 63.036 W / m2K
NuDk
D
Realicemos un balance de energía sobre la partición (debemos sumar los flujos por
las dos caras, que son iguales ya que el problema es simétrico):
qconv + qrad = 0 (1)
q′conv = 2 D hi (Tp − Tm,g )
( )
Si aplicamos la analogía eléctrica para un recinto con dos superficies:
1 − εp 1 1 − εs
σ Tp4 − Ts4 εp A p εs A s
q′rad = 2
A pFp,s
1 − εp 1 − εs
+ +
1
εp A p A pFp,s εsA s
Cálculo de los factores de forma: Fp,s = 1
77
q′rad = 2
(
σ Tp4 − Ts4 )
(T )=0
14.286 + 14.286 + 9.094
8.825 (Tp − 900 ) + σ

4
p − 385 4
Sí sustituimos en la ecuación (1) tendremos:
18.833
Tp = 780.7 K = 507.55°C
(T ) = 7139.1 + 1052.3 = 8191.4 W / m

Pérdidas de calor por unidad de longitud:
q′con = q′conv + q′rad = A ihi (Tm,g − Ts ) + σ

4
p − Ts4
18.833
Apartado 3:
El aumento de la transferencia de calor se debe a dos razones: Primero al aumento

del coeficiente de película, manteniéndose la diferencia de temperaturas constante,
con lo cual aumenta el flujo de calor por convección, y segundo a la aparición de un
flujo de calor por radiación positivo que no existía en el apartado 1.
a. El flujo de calor aumentará para todos los valores de la emisividad, ya que la

convección será siempre mayor y el flujo de calor por radiación será siempre
positivo, ya que la temperatura de la partición estará siempre entre 385 y 900 K.
Para la emisividad igual a 1 se hará mínima la resistencia equivalente del circuito
radiante y el flujo de calor por radiación será máximo.
b. El flujo de calor por convección aumentará para todo valor de caudal ya que
siempre aumenta el coeficiente de película, a no ser que el caudal sea tal que al
pasar a la segunda configuración el Reynolds pase a ser laminar y sea necesario
utilizar una correlación de régimen laminar para el segundo caso, si no es así y
siempre se utiliza una correlación de turbulento el flujo disminuirá en la misma
proporción que el h
 π 
& g,1 / 2 Di
 π + 2  = π Re
m
ReD,2 = =
& g,2De
m
A Tµ πDi2 π+2
µ
D,1
 π 
NuD,2 = 0.023 ReD4 ,/25 Pr 0.3 =  
4 /5
 π + 2
NuD,1
 π 
  NuD,1k
4 /5
NuD,2k  π + 2 
= = = 1.103 hi,1
 π 
hi,2
 Di
De
π + 2
78
Durante la transfusión de sangre a un paciente se utiliza un dispositivo compuesto de dos
partes: un conducto circular de caucho (k = 0.1 W/m K) de 10 mm de diámetro interior,
14 mm de diámetro exterior y longitud 5 m, y posteriormente un calentador de sangre.
3. Un caudal de sangre de 200 ml/min entra en el conducto a 2°C procedente de un

banco de sangre y atraviesa los 5 m de longitud en una habitación, cuyo aire y
paredes pueden considerarse a 25°C. Suponiendo que la emisividad superficial
exterior del caucho es 0.9, calcular la temperatura de la sangre al final del conducto
(Ts).
4. Posteriormente ese mismo caudal de sangre entra en un calentador de sangre. El

calentador está compuesto de una matriz sólida que se mantiene a una temperatura
constante, 40°C en nuestro caso, en la cual se encuentra embebido el conducto.
Calcular la longitud de conducto necesaria para que la sangre salga a 37°C del
calentador. El conducto embebido es metálico con una conductividad muy alta y
entre éste y la matriz existe una resistencia de contacto de 0.001 m² K/W.
Matriz
a 40 °C
37°C
25°C
25°C
5m
2°C Ts
RC = 0.001 m² K/W
14 mm
10 mm
10 mm
Nota: Suponer que el flujo está completamente desarrollado y que la sangre tiene las
mismas propiedades que el agua
79
Solución:
Datos:
L =5m Dint = 0.01 m Dext = 0.014 m k = 0.1 W / m·K

Apartado 1:
- Conducto de caucho:
- Caudal de sangre a la entrada: Vent = 200 ml / min Tent = 2°C
&
- Aire y resto de la habitación: Taire = 25°C Trm,hab = T3 = 25°C εhab = 1
- Superficie exterior del conducto: ε ext = 0.9
Tmatriz = 40°C R t,cont = 0.001 m²K / W

Apartado 2:
- Matriz sólida:
Incógnitas:
- Apartado 1: Temperatura de salida de la sangre al final del conducto, Tsal .
- Apartado 2: Longitud del conducto embebido en la matriz, Lmatriz, para que la
temperatura de salida de la matriz sea 37°C.
Hipótesis:
- Conducción unidimensional en la pared del conducto
Esquema:
Sección media qrad,ext

Del conducto
Tent Tsal
qconv,ext
&
m qconv,int
Desarrollo:
Apartado 1:
El caudal volumétrico a la entrada lo convertimos en caudal másico:
m& = ρV
& ent = 0.0034 kg / s
Agua a presión atmosférica y 2°C, Tabla 4.5: ρ = 1005 kg / m³
Si realizamos un balance de energía sobre el volumen de sangre que circula por el
q=m & cp (Tsal − Tent )

interior del conducto:
80
Ese calor que gana la sangre es igual al calor que desde el exterior atraviesa el
conducto, utilizando la analogía eléctrica para conducción en la sección media del
conducto (supondremos que la temperatura en la sección media es la media de
masas):
Tent + Tsal
Thab −
ln(Dext / Dint )
q= 2
+ +
(hext + hr,ext ) A ext
1 1
hint A int 2 πk L
Igualando esta ecuación con la anterior tendremos:
Tent + Tsal
Thab −
& cp (Tsal − Tent ) =
ln(Dext / Dint )
2
(hext + hr,ext ) A ext

m (1)
+ +
1 1
hint A int 2 πk L
Las incógnitas en esta ecuación son: Tsal , hint , hext , hr,ext

Podemos plantear tres ecuaciones más para cerrar el problema que son:
→ Nuint = f (Re, Pr )
Nuext = f (Gr, Pr )
hint
hext →
 Thab + Tsup 
= 4 σ ε ext  
3

 
hr ,ext
2
Con estas cuatro ecuaciones el problema queda cerrado pero al ser fuertemente no
lineal estableceremos un proceso iterativo para resolverlo.
Tsal = 10°C .
1. Suponer la temperatura de salida, entre la de entrada y la de la habitación,
2. Calcular el coeficiente de película interior.

3. Calcular el coeficiente de película exterior.
4. Calcular el coeficiente de película radiante exterior.
5. Calcular la temperatura de salida de la ecuación (1)
6. Si esta temperatura es muy diferente de la inicial volver a 2.
Coeficiente de película interior:

Re = = 287.451
4m&
π Dint µ
Propiedades de la sangre (agua) a la tª media de masas (aprox. 5°C, Tabla 4.5)
µ = 1506·10 −6 kg/ms

ρ = 1004 kg/m³

k = 0.5748 W/mK
Pr = 11.0

El régimen es laminar con el flujo completamente desarrollado, dadas estas

condiciones podemos usar los valores para temperatura superficial constante o flujo
de calor constante. Nuestro caso está más cercano al problema de temperatura
81
superficial constante puesto que la temperatura de la habitación es constante a lo

largo del tubo.
Usaremos la correlación (23) de la tabla 6.6:

NuD = int int = 3.66 hint = 210.377 W / m²K .
h D
k
Calcular el coeficiente de película exterior:
Si suponemos que la resistencia controlante será la asociada al coeficiente de

película exterior podemos estimar que la temperatura superficial media del exterior
del conducto será muy parecida a la temperatura media de masas (10+2)/2 = 6°C.
Podemos utilizar la correlación de Morgan (40, tabla 6.9) donde las propiedades
deben evaluarse a la temperatura media de película (25+6)/2 = 15.5°C
υ = 15·10 −6 kg/ms

k = 0.025 W/mK
Propiedades del aire a 15°C 
Pr = 0732

β = = 0.0035
1
 15 .5 + 273.15
g β ∆T D3ext
RaD = Pr = 5824 , con este valor del número de Rayleigh la correlación de
υ2
Morgan toma el siguiente valor:
NuD = = 0.85 RaD0.188 → hext = 7.746 W / m²K

hextDext
k
Calcular el coeficiente de película radiante exterior:

 Thab + Tsup 
= 4 σ ε ext   = 4.91 W / m²K
3

 
hr ,ext
2
Volvemos a la ecuación (1) para calcular la temperatura de salida de la sangre:
Tsal + 2°C
25°C −
0.034 kg / s 4201 J / kg·K (Tsal − 2°C ) =
(0.0303 + 0.1071 + 0.359) K / W
2
Tsal = 5.03°C ≈ 5 °C
Ahora deberíamos volver a recalcular los coeficientes de película pero estos no van
a cambiar apreciablemente.
Apartado 2:
El flujo de calor absorbido por la sangre en el proceso de calentamiento en el

interior del calentador de sangre será:
& cp (37°C − 5°C ) = 0.0034 kg / s 4180 J / kg·K 32°C = 454.784 W

q=m
82
Propiedades de la sangre (agua) a la tª media de masas (aprox. 20°C, Tabla 4.5)
µ = 1001·10 −6 kg/ms

ρ = 999.5 kg/m³

k = 0.5998 W/mK
cp = 4180 J / kg·K

Este calor es igual al transferido desde las paredes del calentador
 1 
q  + R t,cont 
Tmatriz − 21°C
q= A int = π Dint L matriz =  int 
h
Tmatriz − 21°C
+ t,cont
1 R
hint A int A int
Si calculamos el número de Reynolds veremos que el problema sigue siendo

laminar y que por tanto:
NuD = int int = 3.66 hint = 219.527 W / m²K
h D
k
Despejando de la ecuación anterior:
L matriz = 4.233 m
83

Ejercicios Resueltos de Calor

Cargado por

Copyright:

Formatos disponibles

Ejercicios Resueltos de Calor

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ejercicios Resueltos de Calor

Cargado por

Copyright:

Formatos disponibles

Colección de Problemas Propuestos y Resueltos de Transmisión de Calor

Capa 1: espesor 0.4 m, conductividad: k 1 = 0.9 (1 + 0.006 T ) [W / m·K ]

a. Flujo de calor por unidad de área que atraviesa el muro.

Capa 1: e1 = 0.4 m k 1 = 0.9 (1 + 0.006 T ) [W / m·K ]

a. Flujo de calor por unidad de área que atraviesa el muro:

La ecuación diferencial en la capa 1 será la siguiente:

variables con la temperatura siguiendo una ley lineal del tipo: k (T ) = k 0 (1 + β T ) . Si

integramos la ecuación anterior para toda la capa 1:

− q′′ e1 = k 0 [1 + β Tmed ] (T2 − T1 ) = k med (T2 − T1 )

Ahora impondremos las dos condiciones de contorno:

x=0 → hext (Text − T1 ) + q′sol

x = 0 → hext  Text + sol  − T1  = hext (Tsol,aire − T1 ) = q′′

En el otro contorno la condición será:

Tenemos pues 3 ecuaciones con 3 incógnitas ( T1 , T2 , q′′ ):

q′′ = hext (Text − T1 ) + q′sol

hext (Text − T1 ) + q′sol

− e1 (hext (Text − T1 ) + q′sol

Si despejamos en la primera T2 y lo introducimos en la segunda tendremos una

T1 = 67.33°C T2 = 58.72°C q′′ = 26.7 W / m²

Finalmente podemos calcular la temperatura en la superficie 3:

Si pintamos la distribución de temperaturas será la siguiente:

Los coeficientes de película en el exterior e interior de la tubería son 15 y 45 kcal/h·m2°C

A. G, para que la pérdidas sean nulas:

Para la segunda capa tenemos generación y por tanto la ecuación general de

Con condiciones de contorno:

r = r3 → − k2 = hext (T3 − Text )

Integrando la ecuación diferencial una vez e imponiendo la primera condición de

= − r 2 + C1; 0=− r2 + C1; C1 =

− r22 ln(r ) + C2

Si imponemos la segunda condición de contorno:

=  r3 − 2  = hext (T3 − Text ) = hext  − ( )

− r22 ln(r3 ) + Text

G = 61598 kcal / h·m³

Y por tanto la distribución de temperaturas será:

T (r ) = −77 r 2 + 5 ln(r ) + 61.05 [°C]

A. Calcular el valor de la generación volumétrica G en el material B.

B. Proporción entre dT/dx en el material A y en el C:

Luego la distribución de temperaturas tendrá forma de polinomio de 2º orden

x=L → q′′(L ) = 50; GL − 100 = 50; G=

q′′(x ) = x − 100; q′′(x ) = 0 → x =

El coeficiente de película en aire es 9 kcal/h·m2ºC y en agua 200 kcal/h·m2ºC. La

= 8.612 ∆T [kcal / h·m²]

El flujo de calor en este caso será:

Donde la eficiencia global de la superficie aleteada se calcula como:

Área de aleta: A a = (2 L + δ ) W [m²]

Área primaria: A p = (S − δ ) W [m²]

Al ser una aleta recta la eficiencia de aleta puede calcularse como:

qsin = 8.612 S W ∆T = 0.08612 W ∆T [kcal / h]

y por tanto el porcentaje de aumento es: 368%

El flujo de calor en este caso será:

La obtención de la eficiencia global de la superficie aleteada se obtiene igual que

ηa = = 0.446 Donde m = = 84.9, L c = L + δ / 2 = 0.02575 m

El flujo de calor en este caso será:

Con lo cual el aumento con respecto a la situación inicial es del: 447%.

A. Obtener la temperatura de la resistencia eléctrica cuando el calor disipado por ésta

- Condición de contorno exterior: Te , he

Utilizando la analogía eléctrica de conducción, podemos expresar el flujo de calor

la suma de ambos debe ser igual a cero: qi + qe = 0

Si despejamos la temperatura de la interfase de la expresión anterior tendremos: