Guía de Ejercicios Resueltos - 2015 AJ
Guía de Ejercicios Resueltos - 2015 AJ
Guía de Ejercicios Resueltos - 2015 AJ
(a) Si una pared de 20 cm de espesor y sus superficies externas son mantenidas isotérmicamente
a 40 y 60 ºC respectivamente, calcule el flujo de calor que se conduce a través de la pared.
(b) Considere que 40 y 60ºC son las temperaturas ambientales a ambos lados de la pared y que
los coeficientes de transferencia de calor son 6 y 10 W/m2 ·ºC respectivamente ¿Cuál es ahora
el flujo de calor? ¿Cuáles son las temperaturas de las paredes?
SOLUCIÓN:
𝜕2𝑇 𝑇 (0) = 𝑇1
=0 ⇒ ( )
𝑇 𝑥 = 𝐶1 𝑥 + 𝐶2 sujeto a {
𝜕𝑥 2
𝑇(𝐿) = 𝑇2
𝑇(0) = 𝑇1 = 𝐶1 ∙ 0 + 𝐶2 ⇒ 𝐶2 = 𝑇1
𝑇2 − 𝑇1
𝑇(𝐿) = 𝑇2 = 𝐶1 𝐿 + 𝐶2 ⇒ 𝑇2 = 𝐶1 𝐿 + 𝑇1 ⇒ 𝐶1 =
𝐿
𝑇2 − 𝑇1
𝑇 (𝑥) = 𝑥 + 𝑇1
𝐿
𝜕𝑇 𝜕 𝑇2 − 𝑇1 𝑇2 − 𝑇1 𝑇1 − 𝑇2
𝑄̇ = −𝑘𝐴 ⇒ 𝑄̇ = −𝑘𝐴 ( 𝑥 + 𝑇1 ) = −𝑘𝐴 = 𝐿
𝜕𝑥 𝜕𝑥 𝐿 𝐿
𝑘𝐴
(a) Se supone que la conductividad térmica varía linealmente
W W
𝑇 = 20 ℃ ⇒ 𝑘 = 5,65 𝑇1 = 60 ℃ ⇒ 𝑘1 = 5,72
m∙K } ⇒ { m∙K
W W
𝑇 = 200 ℃ ⇒ 𝑘 = 5,97 𝑇2 = 40 ℃ ⇒ 𝑘 2 = 5,69
m ∙K m∙ K
La conductividad térmica 𝑘en la pared es la conductividad promedio entre𝑇1 y 𝑇2 definida por:
1 𝑇2 𝑘 + 𝑘 2 5,72 + 5,69 W
𝑘̅ = ∫ 𝑘 (𝑇)𝑑𝑇 = 1 = = 5,705
∆𝑇 𝑇1 2 2 m ∙K
Conducción
material
𝑇1 − 𝑇2 𝑄̇ 𝑇1 − 𝑇2 60 − 40 W
𝑄̇ = 𝐿 ⇒ = 𝑄̇ ′′ = 𝐿 = 0,2 = 570,5
𝐴 m2
𝑘𝐴 𝑘 5,705
𝑄̇ ′′ 66,29
𝑇1 = 𝑇∞,1 − = 60 − = 53,4 ℃
ℎ1 10
𝑄̇ ′′ 66,29
𝑇2 = 𝑇∞,2 + = 40 + = 51,0 ℃
ℎ2 6
PROBLEMA 2
Un cable de un material muy conductor tiene una resistencia de 4 Ω por metro de longitud. El cable
tiene 4 mm de diámetro y por el circula una corriente de 3 A. Si el coeficiente de transmisión por
convección al ambiente es de 10 W/m2 ºC y la temperatura del ambiente es de 25ºC determine la
temperatura de la superficie del cable. Si se toman en cuenta los efectos de radiación, ¿Cómo se
modifica el resultado anterior? Suponga que para el cable 𝜀 = 0,92.
SOLUCIÓN:
El calor por convección viene dado por
𝑄̇𝑐𝑜𝑛𝑣 = ℎ𝐴 (𝑇𝑆 − 𝑇∞ ) = ℎ(𝜋𝐷𝐿 )(𝑇𝑆 − 𝑇∞ )
El calor por radiación viene dado por
𝑄̇𝑟𝑎𝑑 = 𝜀𝜎𝐴 (𝑇𝑆4 − 𝑇∞4 ) = 𝜀𝜎 (𝜋𝐷𝐿 )(𝑇𝑆4 − 𝑇∞4 )
El calor debido a la Ley de Joule del calor generado por una corriente eléctrica viene dado
por:
𝑄̇ = 𝐼 2 𝑅𝑒
Despreciando la radiación, la ecuación de transferencia de calor es
𝑄̇ = 𝑄̇𝑐𝑜𝑛𝑣 ⇒ 𝐼 2 𝑅𝑒 = ℎ (𝜋𝐷𝐿 )(𝑇𝑆 − 𝑇∞ )
Despejando 𝑇𝑆
𝐼 2 (𝑅𝑒 /𝐿 ) 𝐼 2 (𝑅𝑒 /𝐿) 32 ∙ 4
𝑇𝑆 = + 𝑇∞ = + 𝑇∞ = + 25 = 311,5 ℃
ℎ (𝜋𝐷) ℎ (𝜋𝐷) 10 ∙ 𝜋 ∙ 0,004
Considerando la radiación, la ecuación de transferencia de calor es
𝑄̇ = 𝑄̇𝑐𝑜𝑛𝑣 + 𝑄̇𝑟𝑎𝑑 ⇒ 𝐼 2 𝑅𝑒 = ℎ(𝜋𝐷𝐿 )(𝑇𝑆 − 𝑇∞ ) + 𝜀𝜎 (𝜋𝐷𝐿 )(𝑇𝑆4 − 𝑇∞4 )
Arreglando se obtiene que
𝐼 2 (𝑅𝑒 /𝐿 )
= ℎ(𝑇𝑆 − 𝑇∞ ) + 𝜀𝜎 (𝑇𝑆4 − 𝑇∞4 )
𝜋𝐷
Sustituyendo los datos (temperaturas en grados Kelvin) y la constante de Stefan-Boltzman:
W
𝜎 = 5,67 ∙ 10−8
m2
∙ K4
32 ∙ 4
= 10(𝑇𝑆 − 298) + 0,92 ∙ 5,67 ∙ 10−8 (𝑇𝑆4 − 2984 )
𝜋 ∙ 0,004
Resolviendo la temperatura 𝑇𝑆 con un método numérico:
𝑇𝑆 = 436,4 K = 163 ℃
PROBLEMA 3
Se tiene una tubería metálica (𝑘 𝑡 = 53 W/m ∙ K) de diámetro interno 8 cm y diámetro externo 9 cm.
Dentro de la tubería entra vapor saturado a 𝑇𝑠𝑎𝑡 = 205 ℃ ( ℎ𝑓𝑔 = 1912 kJ/kg ) con un flujo másico
de 0,1 kg/s y un coeficiente convectivo de 20000 W/m2 ·K. Por fuera de la tubería circula aire a 𝑇∞ =
25 ℃ y un coeficiente convectivo de 1000 W/m2 ·K .
(a) Determine la longitud de la tubería para que la calidad en la salida de la tubería sea de 0,3 y
calcule la temperatura superficial del tubo.
SOLUCIÓN:
(a) Primero se plantean las resistencias térmicas desde la temperatura del fluido hasta la
temperatura de la corriente libre
𝐷
1 ln ( 𝑡𝑒 ) 1 1 ln (9/8) 1 K∙m
𝐷𝑡𝑖
𝑅 ′𝑡 = + + = + + = 0,00409
ℎ 𝑖 ∙ 𝜋𝐷𝑡𝑖 2𝜋𝑘 𝑡 ℎ 𝑒 ∙ 𝜋𝐷𝑡𝑒 20000 ∙ 𝜋 ∙ 0,08 2𝜋 ∙ 53 1000 ∙ 𝜋 ∙ 0,09 W
∆𝑇 𝑇𝑖 − 𝑇∞ 205 − 25 W
𝑄̇ ′ = ′ = ′ = = 44016
𝑅𝑡 𝑅𝑡 0,00409 m
ℎ 𝑒 = ℎ𝑓 + 𝜒𝑒 ℎ𝑓𝑔
{ ⇒ ℎ 𝑒 − ℎ 𝑠 = ∆ℎ = ( 𝜒𝑒 − 𝜒𝑠 )ℎ𝑓𝑔 = ℎ𝑓𝑔 ∆𝜒
ℎ 𝑠 = ℎ𝑓 + 𝜒𝑠 ℎ𝑓𝑔
kg J
𝑄̇ = 𝑚̇ (ℎ 𝑒 − ℎ 𝑠 ) = 𝑚̇ℎ𝑓𝑔 ∆𝜒 = 0,1 ∙ 1912000 ∙ (1 − 0,3) = 133840 W
s kg
𝑄̇ 133840 W
𝐿= = = 3,04 m
𝑄̇ ′ 44016 W/m
(b) Al añadir aislante a la tubería entonces se agregan dos resistencias térmicas más, la
resistencia de conducción del aislante y la resistencia de contacto:
h ∙ ft 2 ∙ ℉ K ∙ m2
𝑅′′
𝑡𝑐 = 0,001 = 1,76 ∙ 10−4
BTU W
El calor transferido es
∆𝑇 𝑇𝑖 − 𝑇∞ 205 − 25
𝑄̇ = = = = 35698 W = 35,698 kW
𝑅𝑡 𝑅𝑡 0,00504
𝑄̇ 𝑄̇ 35,698
𝑄̇ = 𝑚̇ ℎ𝑓𝑔 ∆𝜒 ⇒ ∆𝜒 = ⇒ 𝜒𝑠 = 𝜒𝑒 − =1− = 0,813
𝑚̇ℎ𝑓𝑔 𝑚̇ ℎ𝑓𝑔 0,1 ∙ 1912
𝐷
𝑇𝑖 − 𝑇𝑠𝑒 𝑇𝑖 − 𝑇𝑠𝑒 1 ln ( 𝑡𝑒)
𝐷𝑡𝑖
𝑄̇ = = ⇒ 𝑇𝑠𝑒 = 𝑇𝑖 − 𝑄̇ ( + )
𝑅𝑡 𝐷 ℎ 𝑖 ∙ 𝜋𝐷𝑡𝑖𝐿 2𝜋𝑘 𝑡 𝐿
1 ln ( 𝑡𝑒 )
𝐷𝑡𝑖
+
ℎ 𝑖 ∙ 𝜋𝐷𝑡𝑖 𝐿 2𝜋𝑘 𝑡𝐿
1 ln(9/8)
𝑇𝑠𝑒 = 205 − 35698 ( + ) = 185,3 ℃
20000 ∙ 𝜋 ∙ 0,08 ∙ 3,04 2𝜋 ∙ 53 ∙ 3,04
La pared de un horno de inducción es como la mostrada en la figura, y está compuesta por dos capa
de refractarios de conductividad térmica k1 = 1,2 W/m·K y k2 = 0,30 W/m·K respectivamente. La
estructura del horno es de un chapa de acero (k3 = 60 W/m·K) de 10 mm de espesor. Si los gases
en el interior del horno están a 1000 °C y producen un coeficiente de convección de 100 W/m2 ·K.
El aire en el exterior se encuentra a 40 °C y produce un coeficiente de convección en la pared del
horno de 15 W/m2 ·K. Considerando la radiación debido a las resistencias inductivas que están a
1000 °C y la pared interior tiene una emisividad de 0,8, determine:
(a) La pérdida de calor por unidad de área que ocurra a través de las paredes del horno.
L1 = L2 = 120 mm
L3 = 10 mm
k1 = 1,2 W/m°C
k2 = 0,3 W/m°C
k3 = 60 W/m°C
Ti = 1000 °C
Te = 40 °C
hi = 100 W/m2 K
he = 15 W/m2 K
ε = 0,8
Ta = 1000°C (1273 K)
SOLUCIÓN:
Radiación
𝑇1 − 𝑇𝑒
𝑄̇ = 𝐿 𝐿 𝐿
̇ ̇
1 = 𝑄𝐶 + 𝑄𝑅
1
+ 2 + 3 +
𝑘1 𝐴 𝑘 2 𝐴 𝑘 3 𝐴 ℎ 𝑒 𝐴
𝑄̇𝐶 = ℎ 𝑖 𝐴(𝑇𝑖 − 𝑇1 )
𝑇1 − 𝑇𝑒 4)
( ) ( 4
𝐿1 𝐿2 𝐿3 1 = ℎ 𝑖 𝐴 𝑇𝑖 − 𝑇1 + 𝜀𝜎𝐴 𝑇𝑎 − 𝑇1
+ + +
𝑘1 𝐴 𝑘 2 𝐴 𝑘 3 𝐴 ℎ 𝑒𝐴
𝑇1 − 𝑇𝑒 4)
( ) ( 4
𝐿1 𝐿2 𝐿3 1 = ℎ 𝑖 𝑇𝑖 − 𝑇1 + 𝜀𝜎 𝑇𝑎 − 𝑇1
+ + +
𝑘1 𝑘 2 𝑘 3 ℎ 𝑒
𝑇𝑒 = 40 + 273 = 313 𝐾
𝑇1 − 313 −8 ( 4 4)
( )
0,12 0,12 0,01 1 = 100 1273 − 𝑇1 + 0,8 ∙ 5,67 ∙ 10 1273 − 𝑇1
+ + +
1,2 0,3 60 15
Conocida esta temperatura, se calcula la transferencia de calor por unidad de área de pared del horno
𝑄̇ 𝑇1 − 𝑇𝑒 1269 − 313 W
= 𝑄̇ ′′ = 𝐿 𝐿 𝐿 1 = 0,12 0,12 0,01 1 = 1686,6 2
𝐴 1
+ 2+ 3+ + + + m
𝑘1 𝑘 2 𝑘 3 ℎ 𝑒 1,2 0,3 60 15
𝑇4 − 𝑇𝑒 1 1686,6
𝑄̇ ′′ = 1 ⇒ 𝑇4 = 𝑇𝑒 + 𝑄̇ ′′ ( ) = 40 + = 152 ℃
ℎ𝑒 15
ℎ𝑒
PROBLEMA 5
(a) Determine la fracción del calor producido que se transfiere a cada uno de los dos fluidos.
(b) Determine las temperaturas superficiales del cilindro y la temperatura máxima en el interior
del mismo.
SOLUCIÓN:
Suponiendo que el proceso se lleva a cabo en estado estacionario, la conductividad térmica constante
y la conducción de calor es unidireccional en la dirección radial, se tiene que el balance de energía
diferencial viene dado por:
1𝜕 𝜕𝑇 𝑞̇
(𝑟 ) + = 0
𝑟 𝜕𝑟 𝜕𝑟 𝑘
Donde 𝑞̇ es la generación de calor por unidad de volumen del material. Integrando la ecuación una
vez se tiene que
1𝜕 𝜕𝑇 𝑞̇ 𝜕 𝜕𝑇 𝑞̇ 𝜕𝑇 𝑞̇ 𝑞̇
(𝑟 ) + = 0 ⇒ (𝑟 ) = − 𝑟 ⇒ 𝑟 = ∫ − 𝑟𝑑𝑟 + 𝐶1 = − 𝑟 2 + 𝐶1
𝑟 𝜕𝑟 𝜕𝑟 𝑘 𝜕𝑟 𝜕𝑟 𝑘 𝜕𝑟 𝑘 2𝑘
Integrando otra vez
𝜕𝑇 𝑞̇ 𝐶1 𝑞̇ 𝐶1 𝑞̇
= − 𝑟+ ⇒ 𝑇(𝑟) = ∫ ( − 𝑟 + ) 𝑑𝑟 + 𝐶2 = − 𝑟 2 + 𝐶1 ln 𝑟 + 𝐶2
𝜕𝑟 2𝑘 𝑟 2𝑘 𝑟 4𝑘
El radio 𝑅0 es el radio tal que la temperatura en el sólido es máxima (𝑇0 ). La condición para encontrar
la temperatura máxima es
𝜕𝑇 𝑞̇ 𝐶1 𝑞̇ 𝐶1 𝑞̇ 2
| =0 ⇒ (− 𝑟 + )| =0 ⇒ − 𝑅0 + =0 ⇒ 𝐶1 = 𝑅
𝜕𝑟 𝑟 = 𝑅 2𝑘 𝑟 𝑟=𝑅 2𝑘 𝑅0 2𝑘 0
0 0
𝑞̇ 2 𝑞̇ 2 𝑞̇ 𝑞̇
𝑇0 = (− 𝑟 + 𝑅 0 ln 𝑟 + 𝐶2) | = − 𝑅20 + 𝑅20 ln 𝑅0 + 𝐶2
4𝑘 2𝑘 𝑟=𝑅 4𝑘 2𝑘
0
𝑞̇ 2 𝑞̇ 2 𝑞̇ 𝑞̇ 𝑞̇ 𝑞̇ 2 𝑟
𝑇(𝑟) − 𝑇0 = (− 𝑟 + 𝑅0 ln 𝑟 + 𝐶2 ) − (− 𝑅20 + 𝑅 20 ln 𝑅0 + 𝐶2 ) = ( 𝑅20 − 𝑟 2 ) + 𝑅0 ln ( )
4𝑘 2𝑘 4𝑘 2𝑘 4𝑘 2𝑘 𝑅0
Arreglando la ecuación anterior, el perfil de temperaturas se convierte en (aun con las constantes 𝑅0
y 𝑇0 desconocidas):
𝑞̇ 𝑅 20 𝑟 2 𝑟
𝑇(𝑟) − 𝑇0 = (1 − ( ) + 2 ln ( ))
4𝑘 𝑅0 𝑅0
𝑞̇ 𝑅 20 𝑅𝑒 2 𝑅
si 𝑟 = 𝑅𝑒 ⇒ 𝑇 = 𝑇𝑠𝑒 𝑇𝑠𝑒 − 𝑇0 = ( 1− ( ) + 2 ln ( 𝑒 ))
4𝑘 𝑅0 𝑅0
{ ⇒ 2 2
si 𝑟 = 𝑅 𝑖 ⇒ 𝑇 = 𝑇𝑠𝑖 𝑞̇ 𝑅 𝑅𝑖 𝑅𝑖
𝑇𝑠𝑖 − 𝑇0 = 0 (1 − ( ) + 2 ln ( ))
{ 4𝑘 𝑅0 𝑅0
𝑞̇ 𝑅20 𝑅𝑖 2 𝑅 2 𝑅
𝑇𝑠𝑒 − 𝑇𝑠𝑖 = (( ) − ( 𝑒 ) + 2 ln ( 𝑒))
4𝑘 𝑅0 𝑅0 𝑅𝑖
𝑄̇ 𝑖 = 𝑞̇ 𝜋(𝑅 2𝑒 − 𝑅 20 )𝐿 = ℎ 𝑒 𝐴𝑒 (𝑇𝑠𝑒 − 𝑇∞ )
𝑄̇ = ∭ 𝑞̇ 𝑑𝑉 = 𝑞̇ 𝑉 = 𝑄̇ 𝑖 + 𝑄̇𝑒 ⇒ {
𝑉 𝑄̇ 𝑒 = 𝑞̇ 𝜋(𝑅 20 − 𝑅 2𝑖 )𝐿 = ℎ 𝑖 𝐴𝑖 (𝑇𝑠𝑖 − 𝑇𝑖 )
𝑞̇ 𝜋
𝑇𝑠𝑒 − 𝑇∞ = (𝑅 2 − 𝑅20 )𝐿
ℎ 𝑒 𝐴𝑒 𝑒 𝑞̇ 𝜋 𝑞̇ 𝜋
⇒ 𝑇𝑠𝑒 − 𝑇𝑠𝑖 = (𝑅2𝑒 − 𝑅20 )𝐿 − (𝑅 2 − 𝑅 2𝑖 )𝐿
𝑞̇ 𝜋 ℎ 𝑒 𝐴𝑒 ℎ 𝑖 𝐴𝑖 0
𝑇𝑠𝑖 − 𝑇𝑖 = (𝑅 2 − 𝑅 2𝑖 )𝐿
{ ℎ𝑖 𝐴𝑖 0
𝑞̇ 𝜋 𝑞̇ 𝜋 𝑞̇ 𝑅2 𝑅 2 𝑅 2 𝑅
( 𝑅2𝑒 − 𝑅 20 )𝐿 − (𝑅 20 − 𝑅2𝑖 )𝐿 = 0 (( 𝑖 ) − ( 𝑒 ) + 2 ln ( 𝑒))
ℎ 𝑒 𝐴𝑒 ℎ 𝑖 𝐴𝑖 4𝑘 𝑅0 𝑅0 𝑅𝑖
𝜋 𝜋 1 𝑅
(𝑅2𝑒 − 𝑅20 )𝐿 − (𝑅02 − 𝑅 2𝑖 )𝐿 = (𝑅2𝑖 − 𝑅 2𝑒 + 2𝑅20 ln ( 𝑒 ))
ℎ 𝑒 𝐴𝑒 ℎ 𝑖 𝐴𝑖 4𝑘 𝑅𝑖
𝑅 2𝑒 𝑅2 1 1 𝑅
4𝜋𝑘𝐿 ( + 𝑖 − 𝑅 20 ( + )) = 𝑅 2𝑖 − 𝑅 2𝑒 + 2𝑅20 ln ( 𝑒 )
ℎ 𝑒 𝐴𝑒 ℎ 𝑖 𝐴 𝑖 ℎ 𝑖 𝐴𝑖 ℎ 𝑒 𝐴𝑒 𝑅𝑖
𝑅 2𝑒 𝑅 2𝑖 1 1 𝑅
4𝜋𝑘𝐿 ( + − 𝑅 20 ( + )) = 𝑅2𝑖 − 𝑅2𝑒 + 2𝑅 20 ln ( 𝑒 )
ℎ 𝑒 ∙ 2𝜋𝑅𝑒 𝐿 ℎ 𝑖 ∙ 2𝜋𝑅𝑖 𝐿 ℎ 𝑖 ∙ 2𝜋𝑅𝑖 𝐿 ℎ 𝑒 ∙ 2𝜋𝑅𝑒𝐿 𝑅𝑖
𝑅𝑒 𝑅𝑖 1 1 𝑅
4𝑘 ( + − 𝑅20 ( + )) = 𝑅 2𝑖 − 𝑅 2𝑒 + 2𝑅20 ln ( 𝑒 )
2ℎ 𝑒 2ℎ 𝑖 ℎ 𝑖 ∙ 2𝑅𝑖 ℎ 𝑒 ∙ 2𝑅 𝑒 𝑅𝑖
Sustituyendo valores y resolviendo para 𝑅0
0,04 0,02 1 1 4
4∙ 8 ∙( + − 𝑅 20 ( + )) = 0,022 − 0,042 + 2𝑅20 ln ( )
2 ∙ 240 2 ∙ 60 60 ∙ 2 ∙ 0,02 240 ∙ 2 ∙ 0,04 2
0,0046
𝑅 20 = ⇒ 𝑅0 = 0,0237 m = 2,37 cm
ln (2) + 7,5
𝑞̇ 𝜋 𝑞̇ 𝜋 𝑞̇ (𝑅2𝑒 − 𝑅20 )
𝑇𝑠𝑒 = 𝑇∞ + (𝑅2𝑒 − 𝑅 20 )𝐿 = 𝑇∞ + (𝑅2𝑒 − 𝑅20 )𝐿 = 𝑇∞ +
ℎ 𝑒 𝐴𝑒 ℎ 𝑒 ∙ 2𝜋𝑅𝑒 𝐿 ℎ 𝑒 ∙ 2𝑅 𝑒
𝑞̇ 𝜋 2 2 𝑞̇ 𝜋 2 2
𝑞̇ (𝑅 20 − 𝑅 2𝑖 )
𝑇𝑠𝑖 = 𝑇𝑖 + (𝑅 − 𝑅𝑖 )𝐿 = 𝑇𝑖 + (𝑅 − 𝑅𝑖 )𝐿 = 𝑇𝑖 +
ℎ 𝑖 𝐴𝑖 0 ℎ 𝑖 ∙ 2𝜋𝑅 𝑖 𝐿 0 ℎ 𝑖 ∙ 2𝑅 𝑖
𝑞̇ 𝑅20 𝑅 2 𝑅 𝑞̇ 𝑅 20 𝑅 2 𝑅
𝑇𝑠𝑒 − 𝑇0 = (1 − ( 𝑒 ) + 2 ln ( 𝑒 )) ⇒ 𝑇0 = 𝑇𝑠𝑒 − (1 − ( 𝑒) + 2 ln ( 𝑒 ))
4𝑘 𝑅0 𝑅0 4𝑘 𝑅0 𝑅0
Dos paredes paralelas, a temperaturas T1 y T2 , están separadas entre sí una distancia L. Entre ambas
se coloca una barra cilíndrica de radio R y conductividad térmica k. El ambiente se encuentra a
temperatura T∞ y por el mecanismo de convección recibe calor de la barra con coeficiente de
convección h.
SOLUCIÓN:
𝑑𝑞𝑐𝑜𝑛𝑑
𝑞𝑐𝑜𝑛𝑑 = 𝑞𝑐𝑜𝑛𝑑 + 𝑑𝑥 + 𝑑𝑞𝑐𝑜𝑛𝑣
𝑑𝑥
𝑑𝑇
𝑞𝑐𝑜𝑛𝑑 = −𝑘𝐴 𝑡
𝑑𝑥
𝑑𝑞𝑐𝑜𝑛𝑣 = ℎ (𝑇 − 𝑇∞ )𝑑𝐴 𝑠 = ℎ(𝑇 − 𝑇∞ )𝑃𝑑𝑥
𝑑 𝑑𝑇 𝑑 𝑑𝑇 ℎ𝑃
0= (−𝑘𝐴 𝑡 ) 𝑑𝑥 + ℎ(𝑇 − 𝑇∞ )𝑃𝑑𝑥 ⇒ (𝐴 )− (𝑇 − 𝑇∞ ) = 0
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑡 𝑑𝑥 𝑘
𝑑 2 𝑇 ℎ𝑃
− (𝑇 − 𝑇∞ ) = 0
𝑑𝑥 2 𝑘𝐴 𝑡
𝑇 − 𝑇∞ 𝑑 2𝜃 ℎ𝑃 𝑑 2𝜃
𝜃≡ ⇒ − 𝜃=0 ⇒ − 𝑚2 𝜃 = 0
𝑇1 − 𝑇∞ 𝑑𝑥 2 ⏟𝑘𝐴 𝑡 𝑑𝑥 2
𝑚2
𝑑 2𝜃
− 𝑚2 𝜃 = 0
𝑑𝑥 2
𝜃 = 1 si 𝑥 = 0
𝑇2 − 𝑇∞
𝜃 = 𝜃2 = si 𝑥 = 𝐿
{ 𝑇1 − 𝑇∞
𝑇2 − 𝑇∞ 𝜃2 − cosh(𝑚𝐿)
𝑥=𝐿 ⇒ 𝜃2 = = 𝐶1 cosh(𝑚 ∙ 𝐿) + 𝐶2 senh (𝑚 ∙ 𝐿) ⇒ 𝐶2 =
𝑇1 − 𝑇∞ senh (𝑚𝐿)
𝜃2 − cosh(𝑚𝐿)
𝜃 = cosh(𝑚𝑥) + senh(𝑚𝑥)
senh (𝑚𝐿)
cosh(𝑚𝐿) 1 𝑇2 − 𝑇∞ cosh(𝑚𝐿) 1
ℎ𝑃 (𝑇1 − 𝑇∞ ) ( − +( )( − )) =
𝑚 senh(𝑚𝐿) 𝑚 senh (𝑚𝐿) 𝑇1 − 𝑇∞ 𝑚senh (𝑚𝐿) 𝑚 senh (𝑚𝐿)
cosh(𝑚𝐿) − 1 𝑇 − 𝑇∞
𝑄𝑐𝑜𝑛𝑣 = ℎ𝑃 (𝑇1 − 𝑇∞ ) (1 + 2 )
𝑚 senh (𝑚𝐿) 𝑇1 − 𝑇∞
Método 2. Mediante balance de energía y calculando los calores por conducción encontrarndo
indirectamente el de convección
𝑑𝑇 𝑑𝑇 𝑑𝑇 𝑑𝑇
𝑄𝑐𝑜𝑛𝑣 = −𝑘𝐴 𝑡 | − (−𝑘𝐴 𝑡 | ) = 𝑘𝐴 𝑡 ( | − | )=
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑥 =0 𝑥=𝐿 𝑥 =𝐿 𝑥=0
𝑑 𝑇 − 𝑇∞ 𝑑 𝑇 − 𝑇∞
= 𝑘𝐴 𝑡 (𝑇1 − 𝑇∞ ) ( ( )| − ( )| )
𝑑𝑥 𝑇1 − 𝑇∞ 𝑑𝑥 𝑇1 − 𝑇∞
𝑥=𝐿 𝑥 =0
𝑑𝜃 𝑑𝜃
= 𝑘𝐴 𝑡 (𝑇1 − 𝑇∞ ) ( | − | )
𝑑𝑥 𝑑𝑥
𝑥=𝐿 𝑥 =0
Calculando la derivada
1 𝑇2 − 𝑇∞ cosh(𝑚𝐿) cosh(𝑚𝐿) 𝑇2 − 𝑇∞ 1
𝑄𝑐𝑜𝑛𝑣 = 𝑘𝐴 𝑡 (𝑇1 − 𝑇∞ )𝑚 (− +( ) − (− +( ) ))
(
senh 𝑚𝐿 ) (
𝑇1 − 𝑇∞ senh 𝑚𝐿 ) (
senh 𝑚𝐿 ) 𝑇1 − 𝑇∞ senh(𝑚𝐿)
1 𝑇2 − 𝑇∞ cosh(𝑚𝐿) cosh(𝑚𝐿) 𝑇2 − 𝑇∞ 1
𝑄𝑐𝑜𝑛𝑣 = 𝑘𝐴 𝑡 (𝑇1 − 𝑇∞ )𝑚 (− +( ) + −( ) )
(
senh 𝑚𝐿 ) (
𝑇1 − 𝑇∞ senh 𝑚𝐿 ) (
senh 𝑚𝐿 ) 𝑇1 − 𝑇∞ senh (𝑚𝐿)
cosh(𝑚𝐿) − 1 𝑇 − 𝑇∞
𝑄𝑐𝑜𝑛𝑣 = 𝑘𝐴 𝑡 (𝑇1 − 𝑇∞ )𝑚 (1 + 2 )
senh (𝑚𝐿) 𝑇1 − 𝑇∞
𝑑𝜃 cosh(𝑚(𝐿 − 𝑥 0 )) 𝑇2 − 𝑇∞ cosh(𝑚𝑥0 )
| =0 ⇒ −𝑚 +𝑚( ) =0
𝑑𝑥 𝑥 = 𝑥 senh (𝑚𝐿) 𝑇1 − 𝑇∞ senh (𝑚𝐿)
0
𝑇2 − 𝑇∞
− cosh(𝑚(𝐿 − 𝑥 0 )) + ( ) cosh(𝑚𝑥0 ) = 0
𝑇1 − 𝑇∞
PROBLEMA 7
Conforme se colocan más y más componentes en un solo circuito integrado (chip), la cantidad de
calor que se disipa continúa en aumento. Sin embargo, este incremento está limitado por la
temperatura máxima permisible de operación del chip, que es alrededor de 75°C. Para maximizar la
disipación de calor se propone que un arreglo de 4×4 aletas rectas circulares de cobre se una
metalúrgicamente a la superficie externa de un clip cuadrado que tiene 12,7 mm de lado.
SOLUCIÓN:
𝑇𝐶 − 𝑇∞,𝑖 𝑇𝐶 − 𝑇∞,𝑖
𝑄̇𝑖 = = ′′
𝑅′′
𝑡,𝑐 𝐿 1 𝑅 𝑡,𝑐 𝐿 1
+ 𝑏 + + 𝑏 +
𝐴 𝑘 𝑏𝐴 ℎ 𝑖 𝐴 𝑤 2 𝑘 𝑏𝑤 2 ℎ 𝑖 𝑤 2
75 − 20
𝑄̇𝑖 = = 0,29 W
1∙ 10−4 0,005 1
+ +
0,01272 1 ∙ 0,01272 40 ∙ 0,01272
𝑇𝐶 − 𝑇∞,0 𝑇𝐶 − 𝑇∞,0
𝑄̇ 0 = 1 = 1
𝑁 1 𝑁 𝜋
+ 1 + ℎ 0 (𝑤 2 − 𝐷𝑃2 𝑁)
𝑅𝑡𝑓 𝑅𝑡𝑓 4
𝜋 2
ℎ0 (𝐴 − 𝐷𝑃 𝑁)
4
𝐷𝑃 1,5
𝐿 𝐶 = 𝐿𝑃 + = 15 + = 15,375 mm = 0,015375 m
4 4
4ℎ 0 4 ∙ 1000
𝑚=√ =√ = 81,65 m−1
𝑘𝑓𝐷𝑃 400 ∙ 0,0015
1 1 K
𝑅𝑡𝑓 = = −5 = 20,39
𝜂𝑓 𝐴𝑓 ℎ0 0,677 ∙ 7,245 ∙ 10 ∙ 1000 W
75 − 20
𝑄̇0 = 1 = 50,47 W
16 𝜋
+ 1000 ∙ (0,01272 − ∙ 0,00152 ∙ 16)
20,39 4
Una tarjeta de circuitos disipa una potencia de 5 W. Suponiendo que la superficie inferior disipa una
cantidad de calor muy pequeña, se puede asumir condición adiabática para simplificar el problema.
El sistema está compuesto de la manera siguiente (de abajo hacia arriba): La tarjeta de 1 mm de
espesor (kT = 12 W/m·K), un plástico de 1 mm de espesor (kA = 1,8 W/m·K) que fija la parte siguiente,
un cuerpo de aletas rectas de perfil rectangular (N = 15 aletas, kf =237 W/m·K). Por el área aleteada
se expone a un fluido a una temperatura de 37 °C y un coeficiente convectivo de 45 W/m2 ·K.
SOLUCIÓN:
Conducción Conducción
tarjeta adhesivo Aletas
Area sin
aletas
Al considerar la tarjeta como una pared plana en estado estacionario, se tiene que su ecuación
unidimensional de conducción viene dada por
𝜕2𝑇
𝑘𝑇 + 𝑞̇ = 0
𝜕𝑥 2
𝑇(𝐿 𝑇 ) = 𝑇𝑇
𝜕𝑇
| =0
{ 𝜕𝑥 𝑥 = 0
𝜕2𝑇 𝑑𝑇 𝑞̇ 𝑞̇ 2
𝑘𝑇 = −𝑞̇ ⇒ = − 𝑥 + 𝐶1 ⇒ 𝑇(𝑥) = − 𝑥 + 𝐶1 𝑥 + 𝐶2
𝜕𝑥 2 𝑑𝑥 𝑘𝑇 2𝑘 𝑇
𝜕𝑇 𝑞̇
| =0 ⇒ 0=− ∙ 0 + 𝐶1 ⇒ 𝐶1 = 0
𝜕𝑥 𝑘𝑇
𝑥 =0
Y la condición isotérmica
𝑞̇ 2 𝑞̇ 2
𝑥 = 𝐿𝑇 ⇒ 𝑇 = 𝑇𝑇 ⇒ 𝑇𝑇 = − 𝐿 + 𝐶2 ⇒ 𝐶2 = 𝑇𝑇 + 𝐿
2𝑘 𝑇 𝑇 2𝑘 𝑇 𝑇
𝑞̇ 2 𝑞̇ 2 𝑞̇
𝑇(𝑥) = − 𝑥 + 𝑇𝑇 + 𝐿 𝑇 = 𝑇𝑇 + (𝐿2 − 𝑥 2 )
2𝑘 𝑇 2𝑘 𝑇 2𝑘 𝑇 𝑇
Si 𝑥 = 0 ⇒ 𝑇 = 𝑇𝑚𝑎𝑥 :
𝑞̇ 𝑞̇ 𝐿2
𝑇𝑚𝑎𝑥 = 𝑇(0) = 𝑇𝑇 + (𝐿2𝑇 − 02 ) = 𝑇𝑇 + 𝑇
2𝑘 𝑇 2𝑘 𝑇
𝑄̇
𝑄̇ = ∭ 𝑞̇ 𝑑𝑉 = 𝑞̇ 𝑉 = 𝑞̇ 𝐿 𝑇 𝑤𝑏 ⇒ 𝑞̇ =
𝑉 𝐿 𝑇 𝑤𝑏
𝑄̇ 𝐿2 𝑇𝑚𝑎𝑥 − 𝑇𝑇
𝑇𝑚𝑎𝑥 = 𝑇𝑇 + ( ) 𝑇 ⇒ 𝑄̇ =
𝐿 𝑇 𝑤𝑏 2𝑘 𝑇 𝐿𝑇
2𝑘 𝑇 𝑤𝑏
𝐿𝑇 0,001 K
𝑅𝑡,𝑡𝑎𝑟𝑗𝑒𝑡𝑎 = = = 0,0185
2𝑘 𝑇 𝑤𝑏 2 ∙ 12 ∙ 0,05 ∙ 0,045 W
𝐿𝐴 𝐿𝐴 0,001 K
𝑅 𝑡,𝑎𝑑ℎ𝑒𝑠𝑖𝑣𝑜 = = = = 0,2469
𝑘𝐴 𝐴 𝑘𝐴 𝑤𝑏 1,8 ∙ 0,05 ∙ 0,045 W
2ℎ 2 ∙ 45
𝑚=√ =√ = 43,57 m−1
𝑘𝑓 𝑡 237 ∙ 0,0002
𝑡 0,0002
𝐿𝐶 = 𝐿 + = 0,02 + = 0,0201 m
2 2
1 1 K
𝑅 𝑡𝑓 = = = 13,748
𝜂𝑓 𝐴𝑓ℎ 0,8042 ∙ 0,00201 ∙ 45 W
1 1 K
𝑅𝑡,𝑐𝑜𝑛𝑣 = = = 10,582
ℎ𝐴 𝑒 45 ∙ 0,0021 W
1 1 K
𝑅𝑡 = 𝑅𝑡,𝑡𝑎𝑟𝑗𝑒𝑡𝑎 + 𝑅𝑡,𝑎𝑑ℎ𝑒𝑠𝑖𝑣𝑜 + 𝑁 1 = 0,0185 + 0,2469 + 15 = 1,1089
1 W
+ +
𝑅𝑡𝑓 𝑅𝑡,𝑐𝑜𝑛𝑣 13,748 10,582
∆𝑇 𝑇𝑚𝑎𝑥 − 𝑇∞
𝑄̇ = = ⇒ 𝑇𝑚𝑎𝑥 = 𝑇∞ + 𝑄̇𝑅𝑡
𝑅𝑡 𝑅𝑡
∆𝑇 𝑇𝑚𝑎𝑥 − 𝑇𝑇
𝑄̇ = = ⇒ 𝑇𝑇 = 𝑇𝑚𝑎𝑥 − 𝑄̇𝑅 𝑡,𝑡𝑎𝑟𝑗𝑒𝑡𝑎
𝑅 𝑡,𝑡𝑎𝑟𝑗𝑒𝑡𝑎 𝑅𝑡,𝑡𝑎𝑟𝑗𝑒𝑡𝑎
𝑞̇ = 𝑞̇ 0 (1 − (𝑟/𝑅) 2 )
SOLUCIÓN:
Transferencia de calor
𝑄̇ = ∭ 𝑞̇ 𝑑𝑉 = ∭ 𝑞̇ 𝑟 2 sen 𝜙 𝑑𝑟𝑑𝜙𝑑𝜃
𝑉 𝑉
4
Como el volumen de una esfera es 𝑉 = 𝜋𝑟 3 ⇒ 𝑑𝑉 = 4𝜋𝑟 2 𝑑𝑟
3
𝑅 𝑅
𝑟 2
𝑄̇ = ∭ 𝑞̇ 𝑑𝑉 = ∫ 𝑞̇ ∙ 4𝜋𝑟 2 𝑑𝑟 = ∫ 𝑞̇ 0 (1 − ( ) ) ∙ 4𝜋𝑟 2 𝑑𝑟 =
𝑉 0 0 𝑅
𝑅
𝑟4 𝑟3 𝑟5 𝑅 𝑅3 𝑅5 8𝜋
= 4𝜋𝑞̇ 0 ∫ (𝑟 2 − ) 𝑑𝑟 = 4𝜋𝑞̇ 0 ( − ) | = 4𝜋𝑞̇ 0 ( − 2) = 𝑞̇ 𝑅3
0 𝑅 2 3 5𝑅 2 3 5𝑅 15 0
0
8𝜋 8𝜋
𝑄̇ = 𝑞̇ 0 𝑅3 = ∙ 1 ∙ 106 ∙ 0,33 = 45239 W
15 15
(b) Se emplea un gradiente de temperaturas desde la temperatura externa del núcleo (𝑇𝑠𝑖 ) hasta
la temperatura ambiental (𝑇∞), entonces el calor transferido puede ser expresado como
Donde 𝐴 𝑜 es el área de la superficie externa quitándole las áreas de las bases de las aletas
𝜋
𝐴 𝑜 = 4𝜋𝑅22 − 𝑁𝑓 ∙ 𝐷 2
4
Aletas:
tanh(𝑚𝐿𝑐 )
𝜂𝑓 =
𝑚𝐿𝑐
Donde la longitud corregida 𝐿𝑐 es la longitud tal que el área total de la aleta se conserva (la punta de
la aleta se considera adiabática)
𝜋
𝜋 𝜋𝐷𝐿𝑓 + 𝐷 2 𝐷
Area aleta = 𝜋𝐷𝐿𝑓 + 𝐷 2 = 𝜋𝐷𝐿𝑐 ⇒ 𝐿𝑐 = 4 = 𝐿𝑓 +
4 𝜋𝐷 4
𝐷 0,008
𝐿𝑐 = 𝐿𝑓 + = 0,06 + = 0,062 m
4 4
𝑄̇ 𝑄̇
𝑇𝑠𝑒 = 𝑇∞ + = 𝑇∞ + 𝜋 =
ℎ(𝑁𝑓 𝜂𝑓 𝐴𝑓 + 𝐴 𝑜 ) ℎ (𝑁𝑓𝜂𝑓𝐴𝑓 + 4𝜋𝑅 22 − 𝑁𝑓 ∙ 𝐷 2 )
4
45239
= 25 + 𝜋 = 204,96 ℃ ≈ 205 ℃
100 ∙ (600 ∙ 0,58954 ∙ 0,001508 + 4𝜋 ∙ 0,42 − 600 ∙ ∙ 0,0082 )
4
𝑄̇ 45239
𝑇𝑠𝑖 = 𝑇𝑠𝑒 + 1 1 = 205 + 1 1 = 1569 ℃
− −
𝑅1 𝑅 2 0,3 0,4
4𝜋𝑘 𝑚 4𝜋 ∙ 0,002
1 𝜕 2 𝜕𝑇
𝑘 (𝑟 ) + 𝑞̇ = 0
𝑟 2 𝜕𝑟 𝜕𝑟
Arreglando e integrando una vez
𝜕 2 𝜕𝑇 𝑞̇ 𝜕𝑇 𝑞̇
(𝑟 ) = − 𝑟2 ⇒ 𝑟2 = ∫ ( − 𝑟 2 ) 𝑑𝑟 + 𝐶1
𝜕𝑟 𝜕𝑟 𝑘 𝜕𝑟 𝑘
𝜕𝑇 𝑞̇ 0 𝑟 2 𝑞̇ 0 𝑟4
𝑟2 = − ∫ (1 − ( ) ) 𝑟 2 𝑑𝑟 + 𝐶1 = − ∫ (𝑟 2 − 2 ) 𝑑𝑟 + 𝐶1
𝜕𝑟 𝑘 𝑅 𝑘 𝑅
𝜕𝑇 𝑞̇ 0 𝑟 3 𝑟 5 𝜕𝑇 𝑞̇ 0 𝑟 𝑟3 𝐶1
𝑟2 = − ( − 2 ) + 𝐶1 ⇒ = − ( − 2) + 2
𝜕𝑟 𝑘 3 5𝑅 𝜕𝑟 𝑘 3 5𝑅 𝑟
𝑞̇ 0 𝑟 𝑟3 𝐶 𝑞̇ 𝑟 2 𝑟4 𝐶
𝑇 = ∫ (− ( − 2 ) + 12 ) 𝑑𝑟 + 𝐶2 = − 0 ( − 2
) − 1 + 𝐶2
𝑘 3 5𝑅 𝑟 𝑘 6 20𝑅 𝑟
Como la temperatura tiene que ser finita en r = 0 entonces necesariamente 𝐶1 = 0, entonces el perfil
de temperaturas viene dado por:
𝑞̇ 0 𝑟 2 𝑟4
𝑇=− ( − ) + 𝐶2
𝑘 6 20𝑅2
𝑇 = 𝑇𝑠𝑖 si 𝑟 = 𝑅
𝑞̇ 0 𝑅 2 𝑅4
𝑇𝑠𝑖 = − ( − ) + 𝐶2
𝑘 6 20𝑅2
𝑇 = 𝑇𝑜 si 𝑟 = 0
𝑇𝑜 = 0 + 𝐶2
𝑞̇ 0 𝑅 2 𝑅4 𝑞̇ 0 𝑅2 1 1
𝑇𝑠𝑖 = − ( − ) + 𝑇𝑜 ⇒ 𝑇𝑜 = 𝑇𝑠𝑖 + ( − )
𝑘 6 20𝑅 2 𝑘 6 20
1 ∙ 106 ∙ 0,32 1 1
𝑇𝑜 = 1569 + ( − ) = 1832 ℃
40 6 20
PROBLEMA 10
Considere que se desea estudiar la transferencia de calor en un cable protegido por un aislante
plástico. En el cable circula 250 amperios de corriente eléctrica y se transfiere energía a los
alrededores por mecanismos de radiación y convección. El cable tiene un diámetro de 30 mm, tiene
una conductividad térmica de 200 W/m·K y una resistividad de 0,005 /m. El aislante tiene un
diámetro de 31 mm, una conductividad térmica de 0,15 W/m·K y una emisividad promedio de 0,9.
El aire de la corriente libre se encuentra a una temperatura de 25 °C con coeficiente convectivo de 25
W/m2 ·K. La superficie de los alrededores tiene una temperatura de 35 °C.
T0
TS2
R1
R2 TS1
(a) Obtenga las expresiones analíticas para los perfiles de temperaturas y de transferencia de
calor en el cable y el aislante en función del radio. Establezca claramente todas las
suposiciones empleadas para hallar dichas expresiones.
(b) Calcule la temperatura exterior del aislante (TS2 ) y la temperatura de la interface cable-
aislante (TS1 ) suponiendo que no hay resistencia térmica de contacto en esta.
(c) Determine la temperatura en el centro del cable T0 .
(d) Demuestre que la temperatura promedio en el cable es igual al promedio aritmético entre la
temperatura del centro y de la superficie del cable.
SOLUCIÓN:
(a) Perfiles de temperatura en el cable y aislante.
Suposiciones:
Proceso en estado estacionario.
Propiedades térmicas del cable y el aislante constantes.
Conducción térmica unidimensional.
Generación de calor en el cable uniformemente distribuida.
Resistencia térmica de contacto en la interface despreciable.
donde 𝑞̇ es la generación de calor en el cable por unidad de volumen del mismo. Esta ecuación está
sujeto a las condiciones de borde:
𝑇(𝑅1 ) = 𝑇𝑆1
{𝜕𝑇
| =0
𝜕𝑟
𝑟=0
Integrando la ecuación una vez se tiene que
1𝜕 𝜕𝑇 𝑞̇ 1𝜕 𝜕𝑇 𝑞̇ 𝜕 𝜕𝑇 𝑞̇
(𝑟 ) + =0⇒ (𝑟 ) = − ⇒ (𝑟 ) = − 𝑟
𝑟 𝜕𝑟 𝜕𝑟 𝑘𝐶 𝑟 𝜕𝑟 𝜕𝑟 𝑘 𝐶 𝜕𝑟 𝜕𝑟 𝑘𝐶
𝜕𝑇 𝑞̇ 𝑞̇ 2
𝑟 = ∫ − 𝑟𝑑𝑟 + 𝐶1 = − 𝑟 + 𝐶1
𝜕𝑟 𝑘𝐶 2𝑘 𝐶
𝜕𝑇 𝑞̇ 𝐶1
| =0 ⇒ lim ( − 𝑟+ )=0 ⇒ 𝐶1 = 0
𝜕𝑟 𝑟→0 2𝑘 𝐶 𝑟
𝑟=0
Para que se cumpla la igualdad anterior la constante 𝐶1 debe ser nula.
𝑞̇ 2 𝑞̇ 2
𝑇(𝑅1 ) = 𝑇𝑆1 ⇒ 𝑇(𝑅1 ) = 𝑇𝑆1 = − 𝑅1 + 0 ∙ ln 𝑟 + 𝐶2 ⇒ 𝐶2 = 𝑇𝑆1 + 𝑅
4𝑘 𝐶 4𝑘 𝐶 1
Por lo que
𝑞̇ 2 𝑞̇ 2 𝑞̇ 2 𝑞̇
𝑇(𝑟) = − 𝑟 + 𝐶2 = − 𝑟 + 𝑇𝑆1 + 𝑅1 = 𝑇𝑆1 + ( 𝑅2 − 𝑟 2 )
4𝑘 𝐶 4𝑘 𝐶 4𝑘 𝐶 4𝑘 𝐶 1
También es útil expresar el perfil de temperatura del cable en función de la temperatura del centro
del cable:
𝑞̇ 𝑞̇ 2
𝑇(0) = 𝑇0 = 𝑇𝑆1 + (𝑅12 − 02 ) = 𝑇𝑆1 + 𝑅
4𝑘 𝐶 4𝑘 𝐶 1
Y se dividen las dos ecuaciones:
𝑞̇
𝑇(𝑟) − 𝑇𝑆1 = (𝑅 2 − 𝑟 2 )
4𝑘 𝐶 1 𝑇(𝑟) − 𝑇𝑆1 𝑅12 − 𝑟 2 𝑟 2
⇒ = = 1 − ( )
𝑞̇ 2 𝑇0 − 𝑇𝑆1 𝑅12 𝑅1
𝑇0 − 𝑇𝑆1 = 𝑅
{ 4𝑘 𝐶 1
𝑇(𝑅1 ) = 𝑇𝑆1
{
𝑇(𝑅2 ) = 𝑇𝑆2
𝑅2 𝑇𝑆2 − 𝑇𝑆1
𝑇𝑆2 − 𝑇𝑆1 = 𝐶1 ln 𝑅2 − 𝐶1 ln 𝑅1 = 𝐶1 ln ( ) ⇒ 𝐶1 = 𝑅
𝑅1 ln ( 2 )
𝑅1
𝑇𝑆2 − 𝑇𝑆1
𝐶2 = 𝑇𝑆1 − 𝐶1 ln 𝑅1 = 𝑇𝑆1 − 𝑅 ln 𝑅1
ln ( 2 )
𝑅1
Así se obtiene el perfil de temperaturas en el aislante
𝑇𝑆2 − 𝑇𝑆1 𝑇𝑆2 − 𝑇𝑆1 𝑇𝑆2 − 𝑇𝑆1
𝑇(𝑟) = 𝑅 2 ln 𝑟 + 𝑇𝑆1 − 𝑅2 ln 𝑅1 = 𝑇𝑆1 + 𝑅
(ln 𝑟 − ln 𝑅1 )
ln ( ) ln ( ) ln ( 2 )
𝑅1 𝑅1 𝑅1
𝑇𝑆2 − 𝑇𝑆1 𝑟
𝑇(𝑟) = 𝑇𝑆1 + 𝑅 ln ( )
ln ( 2 ) 𝑅1
𝑅1
𝜕𝑇 𝑞̇ 𝑄̇(𝑟)
𝑄̇(𝑟) = −𝑘 𝐶𝐴 = −𝑘 𝐶2𝜋𝑟𝐿 ( − 𝑟) ⇒ = 𝜋𝑞̇ 𝑟 2
𝜕𝑟 2𝑘 𝐶 𝐿
Para el aislante:
𝜋𝑞̇ 𝑟 2 si 0 ≤ 𝑟 < 𝑅1
𝑄̇(𝑟)
= 2𝜋𝑘𝐴
𝐿 ( )
𝑅2 𝑇𝑆1 − 𝑇𝑆2 si 𝑅1 ≤ 𝑟 ≤ 𝑅2
ln ( )
{ 𝑅1
(b) Realizando un balance de energía tomando como volumen de control la interface cable-
aislante, se obtiene por continuidad de la función de calor la siguiente expresión:
𝑄̇𝐸 𝑅𝐸
𝑄̇𝐸 = 𝐼2 𝑅𝐸 ⇒ = 𝐼2 ( ) = 2502 ∙ 0,005 = 312,5 W/m
𝐿 𝐿
Y la generación de calor por unidad de volumen es
Para encontrar las temperaturas 𝑇𝑆1 y 𝑇𝑆2 se realiza un balance de energía en el cable
𝑄̇ = 𝑄̇ 𝑐𝑜𝑛𝑣 + 𝑄̇𝑟𝑎𝑑
𝑄̇𝑐𝑜𝑛𝑣 = ℎ𝐴 (𝑇𝑆2 − 𝑇∞ )
Así
𝑄̇ 𝑄̇ 𝑐𝑜𝑛𝑣 𝑄̇𝑟𝑎𝑑
= +
𝐿 𝐿 𝐿
𝜋𝑞̇ 𝑅12 = ℎ(2𝜋𝑅2)(𝑇𝑆2 − 𝑇∞ ) + 𝜀𝜎(2𝜋𝑅 2 )(𝑇𝑆2
4 4 )
− 𝑇𝑎𝑙𝑟
𝑞̇ 𝑅12 4 4 )
= ℎ(𝑇𝑆2 − 𝑇∞ ) + 𝜀𝜎(𝑇𝑆2 − 𝑇𝑎𝑙𝑟
2𝑅 2
442097 ∙ 0,0152 4
= 25 ∙ (𝑇𝑆2 − 298) + 0,9 ∙ 5,67 ∙ 10 −8 (𝑇𝑆2 − 3084 )
2 ∙ 0,0155
Una aproximación inicial a la solución de la ecuación anterior se obtiene despreciando el término de
radiación:
442097 ∙ 0,0152 442097 ∙ 0,0152
= 25 ∙ (𝑇𝑆2 − 298) ⇒ 𝑇𝑆2 ~298 + = 426 K
2 ∙ 0,0155 2 ∙ 0,0155 ∙ 25
Usando como aproximación inicial, se resuelve la ecuación de grado 4:
𝑄̇ 2𝜋𝑘𝐴
= 𝜋𝑞̇ 𝑅12 = (
𝑅2 𝑇𝑆1 − 𝑇𝑆2
)
𝐿 ln ( )
𝑅1
1
𝑇̅ = ∭ 𝑇(𝑟) 𝑑𝑉
𝑉 𝑉
𝑇̅ − 𝑇𝑆1
𝜃̅ =
𝑇0 − 𝑇𝑆1
Se cumple que
1
𝜃̅ = ∭ 𝜃 (𝑟)𝑑𝑉
𝑉 𝑉
𝑇(𝑟) − 𝑇𝑆1 𝑟 2
𝜃(𝑟) = = 1 −( )
𝑇0 − 𝑇𝑆1 𝑅1
𝑉 = 𝜋𝑟 2 𝐿 ⇒ 𝑑𝑉 = 2𝜋𝑟𝐿𝑑𝑟
𝑅1
1 1 2𝜋𝐿 𝑅1 𝑟 2
𝜃̅ = ∭ 𝜃 (𝑟)𝑑𝑉 = ∫ 𝜃 (𝑟) 2𝜋𝑟𝐿𝑑𝑟 = ∫ ( 1 − ( ) ) 𝑟𝑑𝑟
𝑉 𝑉 𝜋𝑅12 𝐿 0 𝜋𝑅12 𝐿 0 𝑅1
2 𝑅1 𝑟 2 2 𝑅1 𝑟3 2 𝑟 2 𝑟 4 𝑅1 2 𝑅2 𝑅4
= 2 ∫ (1 − ( ) ) 𝑟𝑑𝑟 = 2 ∫ (𝑟 − 2 ) 𝑑𝑟 = 2 ( − 2 ) | = 2 ( 1 − 12 )
𝑅1 0 𝑅1 𝑅1 0 𝑅1 𝑅1 2 4𝑅1 𝑅1 2 4𝑅1
0
2 𝑅12 𝑅12 1
= 2( − )=
𝑅1 2 4 2
Por lo tanto
𝑇̅ − 𝑇𝑆1 1 1 𝑇0 + 𝑇𝑆1
𝜃̅ = = ⇒ 𝑇̅ = 𝑇𝑆1 + (𝑇0 − 𝑇𝑆1 ) =
𝑇0 − 𝑇𝑆1 2 2 2
Así, la temperatura promedio en este caso es el promedio aritmético de la temperatura del centro y
superficie del cable.
𝑇0 + 𝑇𝑆1 133,1 + 132,9
𝑇̅ = = = 133 ℃
2 2
PROBLEMA 11
La densidad promedio del agua es de 969 kg/m3 y su calor específico promedio es de 4200 J/kgºC.
SOLUCIÓN:
𝑇𝑖 − 𝑇𝑆 𝑄̇ 𝑇𝑖 − 𝑇𝑆
𝑄̇ = ⇒ =
1 ln(𝐷𝑡𝑒 /𝐷𝑡𝑖 ) 1 𝐿 1 ln(𝐷𝑡𝑒 /𝐷𝑡𝑖 ) 1
+ + + +(
ℎ 𝑖(𝜋𝐷𝑡𝑖 𝐿) 2𝜋𝑘 𝑡 𝐿 𝑆𝑘 𝑆 ℎ 𝑖(𝜋𝐷𝑡𝑖 ) 2𝜋𝑘 𝑡 𝑆/𝐿)𝑘𝑆
2𝜋𝐿
𝑆= con 𝑧 > 𝑟
cosh−1 (𝑧/𝑟)
𝑆 2𝜋 2𝜋
= = = 1,7923 m/m
𝐿 cosh (𝑧/𝑟) cosh−1 ( 1 )
−1
0,12/2
𝑄̇ 90 − 20
= = 187,3 W/m
𝐿 1 ln(12/10) 1
+ +
4500 ∙ 𝜋 ∙ 0,1 2𝜋 ∙ 25 1,7923 ∙ 1,5
(b) La longitud de la tubería se encuentra a partir del balance de energía del agua
𝜋 2
𝑚̇ 𝐶𝑃 (𝑇𝑒𝑛𝑡 − 𝑇𝑠𝑎𝑙 ) (1,2 ∙ 969 ∙ 4 ∙ 0,1 ) ∙ 4200 ∙ (90 − 79)
𝑄̇ = 𝑚̇𝐶𝑃 (𝑇𝑒𝑛𝑡 − 𝑇𝑠𝑎𝑙 ) ⇒ 𝐿= = = 2253 m
𝑄̇/𝐿 187,3
PROBLEMA 12
SOLUCIÓN:
𝑄̇ 𝑇𝑖 − 𝑇∞
=
𝐿 1 ln (𝐷𝑡𝑒/𝐷𝑡𝑖) 1 1
+ +( +
ℎ 𝑖 (𝜋𝐷𝑡𝑖 ) 2𝜋𝑘𝐴 𝑆/𝐿)𝑘 𝐵 ℎ 𝑒 (4𝑧)
El diámetro externo es 𝐷𝑡𝑒 = 0,2 m por lo tanto el interno es de 𝐷𝑡𝑒 = 0,2 − 2 ∙ 0,02 = 0,16 m
𝑆 2𝜋
= = 13,02 m/m
𝐿 ln ( 0,54 ∙ 0,3 )
0,2/2
Sustituyendo los valores
𝑄̇ 50 − 25
= = 535,4 W/m
𝐿 1 ln(20/16) 1 1
+ + +
875 ∙ 𝜋 ∙ 0,16 2𝜋 ∙ 25 13,02 ∙ 15 22 ∙ 4 ∙ 0,3
(b) Mediante un balance de energía del agua se busca el cambio de temperatura del agua. El flujo
másico de agua viene dado por (densidad del agua a 50 ºC es 988 kg/m3 y el calor específico
es 4195 J/kgºC)
𝜋
𝑚̇ = 𝜌𝑢𝐴 𝑡 = 988 ∙ 1 ∙ ∙ 0,162 = 19,86 kg/s
4
𝑄̇
𝑄̇ = 𝑚̇𝐶𝑃 (𝑇𝑒𝑛𝑡 − 𝑇𝑠𝑎𝑙 ) = 𝑚̇ 𝐶𝑃 ∆𝑇 ⇒ ∆𝑇 =
𝑚̇𝐶𝑃
Suponiendo que la tasa de transferencia de calor no cambia mucho con la longitud (lo cual es falso
ya que la temperatura interna 𝑇𝑖 varía a lo largo de la longitud, pero aquí se supone que la tasa de
transferencia de calor es la misma calculada en la parte (a)):
W
𝑄̇ 535,4 ∙ 30 m
m
∆𝑇 = = = 0,19 K
𝑚̇𝐶𝑃 19,86 kg ∙ 4195 J
s kg ∙ ℃
Por lo que la variación de temperatura del agua es de 0,2 K, verificando que es válida la aproximación
de tasa de transferencia de calor constante a lo largo de la tubería.
PROBLEMA 13
Agua a una temperatura de 130 °C y una calidad de 99% entra con un caudal de 1,1 kg/s a una tubería
de un material A (kA = 10 W/m·K) de 200 m de longitud y de 4,0 cm y 4,2 cm sus respectivos
diámetros. Está recubierta por un material B (kB = 10 W/m·K) de forma cuadrada de arista 6,0 cm.
La superficie posee 12 aletas longitudinales (a lo largo de toda la estructura) de espesor 2 mm y
longitud 4,0 cm con una conductividad térmica de 15 W/m·K.
SOLUCIÓN:
Area sin
aletas
𝐷
1 ln ( 𝑡𝑒 ) 1 1
𝐷𝑡𝑖
𝑅𝑡 = + + + 𝑁 1
ℎ 𝑖 𝐴𝑖 2𝜋𝑘𝐴 𝐿 𝑆𝑘 𝐵 + 1
𝑅𝑡𝑓
ℎ𝑒 𝐴𝑒
2𝜋 ∙ 200
𝑆= = 2897,9 m
0,54 ∙ 0,06
ln ( )
0,042/2
Luego se calcula la resistencia debido a cada aleta:
𝑡 = 0,002 m
𝑤 = 𝐿 = 200 m
𝐿𝑎 = 0,04 m
W
ℎ 𝑒 = 40 2
m ∙K
W
𝑘𝑓 = 15
m∙K
2ℎ 𝑒 2 ∙ 40
𝑚=√ =√ = 51,64 m−1
𝑘𝑓 𝑡 15 ∙ 0,002
𝑡 0,002
𝐿 𝐶 = 𝐿𝑎 + = 0,04 + = 0,041 m
2 2
1 1 K
𝑅𝑡𝑓 = = = 3,322 ∙ 10−3
𝜂𝑓 𝐴𝑓ℎ𝑒 0,4588 ∙ 16,4 ∙ 40 W
4,2
1 ln ( ) 1 1 K
4,0
𝑅𝑡 = + + + 12 1 = 2,336 ∙ 10−4
5000 ∙ 25,1 2𝜋 ∙ 10 ∙ 200 2897,9 ∙ 10 + W
3,322 ∙ 10−3 1
40 ∙ 43,2
∆𝑇 𝑇𝑖 − 𝑇∞ 130 − 20
𝑄̇ = = = = 470882 W = 470,9 W
𝑅𝑡 𝑅𝑡 2,336 ∙ 10 −4
ℎ 𝑒𝑛𝑡 = ℎ𝑓 + 𝜒𝑒 ℎ𝑓𝑔
{ ⇒ ℎ 𝑒𝑛𝑡 − ℎ 𝑠𝑎𝑙 = ∆ℎ = ( 𝜒𝑒 − 𝜒𝑠 )ℎ𝑓𝑔 = ℎ𝑓𝑔∆𝜒
ℎ 𝑠𝑎𝑙 = ℎ𝑓 + 𝜒𝑠 ℎ𝑓𝑔
𝑄̇ 470,9
𝑄̇ = 𝑚̇ ℎ𝑓𝑔 (𝜒𝑒 − 𝜒𝑠 ) ⇒ 𝜒𝑠 = 𝜒𝑒 − = 0,99 − = 0,793 (79,3%)
𝑚̇ ℎ𝑓𝑔 1,1 ∙ 2174,2
PROBLEMA 14
Una caja tiene dimensiones externas 40cm x 50 cm x 60 cm y sus paredes tienen 3 cm de espesor y
está construida de un material de conductividad 0,22 W/m·K. En su interior hay hielo a 0 ºC y el
ambiente externo está a 30 ºC. El coeficiente de convección interno de la caja es de 5 W/m2 ∙ K y el
exterior 8 W/m2 ∙ K . Determine cuanto tiempo tarda en derretirse 500 g de hielo.
La entalpía de fusión del hielo a 0 °C es 335,2 kJ/kg y la densidad del hielo es 916,8 kg/m3 .
SOLUCIÓN:
1
𝑅𝑡 =
1 1 1 1 1 1 1
2( ) + 2( )+2( )+ 4 ( )+4( )+4( ) + 8( )
𝑅𝑡,𝑃𝑆 𝑅 𝑡,𝑃𝐿1 𝑅𝑡,𝑃𝐿2 𝑅𝑡,𝐴60 𝑅𝑡,𝐴40 𝑅𝑡,𝐴50 𝑅 𝑡,𝐸
Pared superior e inferior 40x50 cm:
𝐴 𝑃𝑆 = (0,4 − 2 ∙ 0,03)(0,4 − 2 ∙ 0,03) = 0,1496 m2 𝐴 0,1496
𝑆𝑃𝑆 = = = 4,99 m
𝐿 0,03
1 1 1 K
𝑅𝑡,𝑃𝑆 = + + = 3,084
ℎ 𝑖 𝐴𝑃𝑆 𝑆𝑃𝑆𝑘 ℎ 𝑒 𝐴𝑃𝑆 W
Arista de 60 cm:
𝐴𝐴60 = 2 ∙ 0,03 ∙ (0,6 − 2 ∙ 0,03) = 0,0162 m2 𝑆𝐴60 = 0,54𝐷 = 0,54 ∙ 0,54 = 0,2916 m
1 1 K
𝑅𝑡,𝐴60 = + = 23,30
𝑆𝐴60 𝑘 ℎ 𝑒 𝐴𝐴60 W
Arista de 40 cm:
𝐴𝐴40 = 2 ∙ 0,03 ∙ (0,4 − 2 ∙ 0,03) = 0,0102 m2 𝑆𝐴40 = 0,54𝐷 = 0,54 ∙ 0,36 = 0,1836 m
1 1 K
𝑅𝑡,𝐴40 = + = 37,01
𝑆𝐴40 𝑘 ℎ 𝑒 𝐴𝐴40 W
Arista de 50 cm:
𝐴𝐴50 = 2 ∙ 0,03 ∙ (0,5 − 2 ∙ 0,03) = 0,0132 m2 𝑆𝐴50 = 0,54𝐷 = 0,54 ∙ 0,44 = 0,2376 m
1 1 K
𝑅𝑡,𝐴50 = + = 28,60
𝑆𝐴50 𝑘 ℎ 𝑒 𝐴𝐴50 W
Esquina:
1 1 K
𝑅 𝑡,𝐸 = + = 1056
𝑆𝐸 𝑘 ℎ 𝑒 𝐴𝐸 W
masa de líquido
𝜒=
masa de líquido + hielo
ℎ1 = ℎ 𝑠 + 𝜒1 (ℎ𝑓 − ℎ 𝑠 )
{ ⇒ ∆ℎ = ℎ 2 − ℎ1 = (𝜒2 − 𝜒1 )ℎ𝑓𝑠
ℎ 2 = ℎ 𝑠 + 𝜒2 (ℎ𝑓 − ℎ 𝑠 )
kJ
𝑄 = 𝑚∆ℎ = 𝑚(𝜒2 − 𝜒1 )ℎ𝑓𝑠 = 74,06 kg ∙ (0,00675 − 0) ∙ 335,2 = 167,6 kJ = 167600 J
kg
𝑄 167600 J
𝑡= = = 1925 𝑠 ~ 32 min
𝑄̇ 87,05 J/s
PROBLEMA 15
SOLUCIÓN:
Primero se calcula el volumen y área superficial de la pieza (que está en contacto con el fluido) que
tiene de arista 𝐿, altura 𝑧 y diámetro de orificio 𝐷:
3√3 𝜋 3√3 𝜋 𝐷 2
𝑉 = 𝐿2 ∙ ∙ 𝑧 − 𝐷 2 𝑧 = 𝑧𝐿2 ( − ( ) )
2 4 2 4 𝐿
3 √3 𝜋 𝑧 𝜋 𝐷 2 𝐷 𝑧
𝐴 𝑠 = 6𝐿𝑧 + 2 (𝐿2 ∙ ) − 2 ( 𝐷 2 ) + 𝜋𝐷𝑧 = 𝐿2 (6 + 3√3 − ( ) + 𝜋 ( ) ( ))
2 4 𝐿 2 𝐿 𝐿 𝐿
𝑉 = 21,633 cm3
𝐴 𝑠 = 82,765 cm2
𝑉
𝐿𝐶 = = 0,26138 cm
𝐴𝑠
𝑇 − 𝑇∞ ℎ𝐴 𝑠 𝑇 − 𝑇∞ ℎ𝐴 𝑠
ln ( )=− 𝑡 ⇒ 𝜃≡ = exp (− 𝑡)
𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃𝑉 𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃 𝑉
𝑇 − 𝑇∞ 30 − 20
𝜃= = = 0,04348
𝑇𝑖 − 𝑇∞ 250 − 20
(c) El calor retirado de las doscientas piezas la absorbe el aceite. Entonces la masa de aceite que
experimenta un cambio en su temperatura de 0,5 K es
𝑁𝑃 𝑄 200 ∙ 20941
𝑚 𝑎 𝐶𝑃𝑎∆𝑇 = 𝑁𝑃 𝑄 ⇒ 𝑚𝑎 = = = 2327 kg
𝐶𝑃𝑎∆𝑇 3600 ∙ 0,5
(d) El valor máximo del coeficiente convectivo para que el problema pueda aproximarse como
resistencia interna despreciable viene dado por
𝜕𝑇
𝜌𝐶𝑃 𝑉 = −ℎ𝐴 𝑠 (𝑇 − 𝑇∞ )
𝜕𝑡
Si se desea resolver este modelo matemático se despeja 𝑇∞ de la primera ecuación y se sustituye en
la segunda
𝜕𝑇 𝑁𝑃 𝜌𝐶𝑃 𝑉 𝑁𝑃 𝜌𝐶𝑃 𝑉
𝜌𝐶𝑃 𝑉 = −ℎ𝐴 𝑠 (𝑇 − (𝑇∞𝑖 + 𝑇𝑖 ) + 𝑇)
𝜕𝑡 𝑚 𝑎 𝐶𝑃𝑎 𝑚𝑎 𝐶𝑃𝑎
𝑑𝑇 ℎ𝐴 𝑠
=− 𝑑𝑡
𝑁𝑃 𝜌𝐶𝑃 𝑉 𝑁 𝜌𝐶 𝑉 𝜌𝐶𝑃 𝑉
(1 + ) 𝑇 − (𝑇∞𝑖 + 𝑃 𝑃 𝑇𝑖 )
𝑚 𝑎 𝐶𝑃𝑎 𝑚 𝑎 𝐶𝑃𝑎
𝑇 𝑡
𝑑𝑇 ℎ𝐴 𝑠
∫ =∫ − 𝑑𝑡
𝑁𝑃 𝜌𝐶𝑃 𝑉 𝑁 𝜌𝐶 𝑉 𝜌𝐶𝑃 𝑉
𝑇𝑖 (1 + ) 𝑇 − (𝑇∞𝑖 + 𝑃 𝑃 𝑇𝑖 ) 0
𝑚 𝑎𝐶𝑃𝑎 𝑚 𝑎 𝐶𝑃𝑎
𝑁 𝜌𝐶 𝑉 𝑁 𝜌𝐶 𝑉
1 (1 + 𝑃 𝑃 ) 𝑇 − (𝑇∞𝑖 + 𝑃 𝑃 𝑇𝑖 ) ℎ𝐴 𝑠
𝑚 𝑎 𝐶𝑃𝑎 𝑚𝑎 𝐶𝑃𝑎
ln ( )=− 𝑡
𝑁 𝜌𝐶 𝑉 𝑁 𝜌𝐶 𝑉 𝑁 𝜌𝐶 𝑉 𝜌𝐶𝑃 𝑉
1+ 𝑃 𝑃 (1 + 𝑃 𝑃 ) 𝑇𝑖 − (𝑇∞𝑖 + 𝑃 𝑃 𝑇𝑖 )
𝑚 𝑎 𝐶𝑃𝑎 𝑚 𝑎 𝐶𝑃𝑎 𝑚𝑎 𝐶𝑃𝑎
PROBLEMA 16
(a) El tiempo que se requiere para que el centro de la cara superior alcance los 30 ºC.
(b) La temperatura del centro geométrico en ese mismo momento.
(c) La temperatura media en ese mismo momento.
(d) Temperatura en una de las esquinas en ese mismo momento.
(e) Calor transferido hasta ese instante.
(f) El tiempo necesario para que el centro de la cara superior alcance los 30 ºC si el coeficiente
de convección se reduce a 20 W/m2 ºC.
8 cm
16 cm
8 cm
SOLUCIÓN:
𝑘 4 m2
𝛼= = = 8,889 ∙ 10−7
𝜌𝐶𝑃 3000 ∙ 1500 s
𝜃 = 𝜃𝑥 𝜃𝑦 𝜃𝑧
𝛼𝑡 𝜁 𝑥 2 𝛼𝑡 )
𝜁1𝑦 𝑦 2 𝛼𝑡 ) 𝜁1𝑧 𝑧
2
𝜃 = (𝐶1𝑥 exp (−𝜁1𝑥 2
) cos ( 1𝑥 ))(𝐶1𝑦 exp (−𝜁1𝑦 2 cos ( ))(𝐶1𝑧 exp ( −𝜁1𝑧 2 cos ( ))
𝐿𝑥 𝐿𝑥 𝐿𝑦 𝐿𝑦 𝐿𝑧 𝐿𝑧
Arreglando términos
2
𝜁1𝑥 2 𝜁1𝑦 𝜁1𝑧 2 𝜁1𝑥 𝑥 𝜁1𝑦 𝑦 𝜁 𝑧
𝜃 = 𝐶1𝑥 𝐶1𝑦 𝐶1𝑧 exp (− (( ) + ( ) + ( ) ) 𝛼𝑡) cos ( ) cos ( ) cos ( 1𝑧 )
𝐿𝑥 𝐿𝑦 𝐿𝑧 𝐿𝑥 𝐿𝑦 𝐿𝑧
𝜃
− ln ( 𝜁1𝑦 𝑦 )
𝜁 𝑥 𝜁 𝑧
𝐶1𝑥 𝐶1𝑦 𝐶1𝑧 cos ( 1𝑥 ) cos ( ) cos ( 1𝑧 )
𝐿𝑥 𝐿𝑦 𝐿𝑧
𝑡= 2 2 2
𝜁 𝜁 𝜁
(( 1𝑥 ) + ( 1𝑦 ) + ( 1𝑧 ) ) 𝛼
𝐿𝑥 𝐿𝑦 𝐿𝑧
𝑇(0,0, 𝐿𝑧 ) − 𝑇∞ 30 − 20
𝜃(0,0, 𝐿𝑧 ) = = = 0,035714
𝑇𝑖 − 𝑇∞ 300 − 20
Sustituyendo valores
0,035714
− ln ( )
1,247882 ∙ 1,26497 ∙ cos(0) ∙ cos(0) ∙ cos(1,45050 ∙ 1)
𝑡= 2 2 = 816 s
1,34955 1,45050 −7
(2 ∙ ( ) +( ) ) ∙ 8,889 ∙ 10
0,04 0,08
𝑡 = [0 horas][13 min][36 s]
2 2
1,34955 1,45050
= 1,247882 ∙ 1,26497 ∙ exp (− (2 ∙ ( ) +( ) ) ∙ 8,889 ∙ 10 −7 ∙ 816) = 0,29761
0,04 0,08
𝑇𝑜 − 𝑇∞
𝜃𝑜 = ⇒ 𝑇𝑜 = 20 + 0,29761 ∙ (300 − 20) = 103,33 ℃
𝑇𝑖 − 𝑇∞
Temperatura media:
1 1 𝐿𝑥 1 𝐿𝑦 1 𝐿𝑧
𝜃̅ = ∭ 𝜃𝑑𝑉 = ( ∫ 𝜃𝑥 𝑑𝑥) ( ∫ 𝜃𝑦 𝑑𝑦) ( ∫ 𝜃𝑧 𝑑𝑧) = 𝜃̅𝑥 𝜃𝑦̅ 𝜃𝑧̅
𝑉 𝑉 𝐿𝑥 0 𝐿𝑦 0 𝐿𝑧 0
2 2
1 𝐿𝑥 1 𝐿𝑥 𝜁
−( 1𝑥) 𝛼𝑡 𝜁1𝑥 𝑥 𝜁
−( 1𝑥) 𝛼𝑡
1
𝜃̅𝑥 = ∫ 𝜃𝑥 𝑑𝑥 = ∫ 𝐶1𝑥 𝑒 𝐿𝑥 cos ( ) 𝑑𝑥 = 𝐶1𝑥 𝑒 𝐿𝑥 ∫ cos(𝜁1𝑥 𝑥) 𝑑𝑥 =
𝐿𝑥 0 𝐿𝑥 0 𝐿𝑥 0
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,00172 ∙ (300 − 20) = 20,48 ℃
𝑇𝑖 − 𝑇∞
Calor transferido:
𝑄 = 𝜌𝑉𝐶𝑃 (𝑇𝑖 − 𝑇̅) = 3000 ∙ 0,082 ∙ 0,16 ∙ 1500 ∙ (300 − 49,81) = 1152876 J
Como se reduce mucho el valor del coeficiente convectivo, puede ser que el problema se aproxime
como resistencia interna despreciable por obtener un número de Biot bajo:
ℎ𝐿𝐶 20 ∙ 0,016
𝐵𝑖𝐶 = = = 0,08 < 0,1 ⇒ es de RID
𝑘 4
El balance de energía suponiendo que la temperatura no depende de los efectos espaciales es
𝜕𝑇 𝑑𝑇 ℎ𝐴 𝑆
𝜌𝐶𝑃 𝑉 = −ℎ𝐴 𝑆 (𝑇 − 𝑇∞ ) ⇒ =− 𝑑𝑡
𝜕𝑡 𝑇 − 𝑇∞ 𝜌𝐶𝑃 𝑉
𝑇 𝑡
𝑑𝑇 ℎ𝐴 𝑆 𝑇 − 𝑇∞ ℎ𝐴 𝑆 𝜌𝐶𝑃 𝑉 𝑇 − 𝑇∞
∫ =∫ − 𝑑𝑡 ⇒ ln ( )= − 𝑡 ⇒ 𝑡 =− ln ( )
𝑇𝑖 𝑇 − 𝑇∞ 0 𝜌𝐶𝑃 𝑉 𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃 𝑉 ℎ𝐴 𝑆 𝑇𝑖 − 𝑇∞
t = [3 horas][19 min][56 s]
PROBLEMA 17
SOLUCIÓN:
𝑘 40 m2
𝛼= = = 1,111 ∙ 10−5
𝜌𝐶𝑃 2000 ∙ 1800 s
0,18
𝐿 𝑥 = 𝐿 𝑦 = 𝐿𝑧 = = 0,09 m
2
𝛼𝑡 𝜁 𝑥 2 𝛼𝑡 )
𝜁1𝑦 𝑦 2 𝛼𝑡 ) 𝜁1𝑧 𝑧
2
𝜃 = (𝐶1𝑥 exp (−𝜁1𝑥 2
) cos ( 1𝑥 ))(𝐶1𝑦 exp (−𝜁1𝑦 2 cos ( ))(𝐶1𝑧 exp ( −𝜁1𝑧 2 cos ( ))
𝐿𝑥 𝐿𝑥 𝐿𝑦 𝐿𝑦 𝐿𝑧 𝐿𝑧
Arreglando términos
2
𝜁1𝑥 2 𝜁1𝑦 𝜁1𝑧 2 𝜁1𝑥 𝑥 𝜁1𝑦 𝑦 𝜁 𝑧
𝜃 = 𝐶1𝑥 𝐶1𝑦 𝐶1𝑧 exp (− (( ) + ( ) + ( ) ) 𝛼𝑡) cos ( ) cos ( ) cos ( 1𝑧 )
𝐿𝑥 𝐿𝑦 𝐿𝑧 𝐿𝑥 𝐿𝑦 𝐿𝑧
𝜃
− ln ( 𝜁1𝑦 𝑦 )
𝜁1𝑥 𝑥 𝜁1𝑧 𝑧
𝐶1𝑥 𝐶1𝑦 𝐶1𝑧 cos ( ) cos ( ) cos ( )
𝐿𝑥 𝐿𝑦 𝐿𝑧
𝑡=
𝜁 2 𝜁 2 𝜁 2
(( 1𝑥 ) + ( 1𝑦 ) + ( 1𝑧 ) ) 𝛼
𝐿𝑥 𝐿𝑦 𝐿𝑧
𝑇(0,0,0) − 𝑇∞ 𝑇𝑜 − 𝑇∞ 30 − 20
𝜃 (0,0,0) = = = = 0,035714
𝑇𝑖 − 𝑇∞ 𝑇𝑖 − 𝑇∞ 300 − 20
Sustituyendo valores
0,035714
− ln ( 3 )
1,06419 ∙ cos(0) ∙ cos(0) ∙ cos(0)
𝑡= = 2192,9 s
0,62444 2
(3 ∙ ( ) ) ∙ 1,111 ∙ 10−5
0,09
𝑡 = [0 horas][36 min][33 s]
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,01907 ∙ (300 − 20) = 25,34 ℃
𝑇𝑖 − 𝑇∞
Evaluando la solución analítica multidimensional en la coordenada del centro de una de las caras
(0,0, 𝐿𝑧 ):
𝜃 = 𝜃𝑜 cos(0) cos(0) cos(𝜁1𝑧 ) = 0,035714 ∙ cos(0,62444) = 0,02897
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,02897 ∙ (300 − 20) = 28,11 ℃
𝑇𝑖 − 𝑇∞
Temperatura media:
3
1 sen 𝜁1𝑥 sen 𝜁1𝑦 sen 𝜁1𝑧 sen(0,62444)
𝜃̅ = ∭ 𝜃𝑑𝑉 = 𝜃𝑜 ∙ ∙ ∙ = 0,035714 ∙ ( ) = 0,02931
𝑉 𝑉 𝜁1𝑥 𝜁1𝑦 𝜁1𝑧 0,62444
𝑇̅ − 𝑇∞
𝜃̅ = ⇒ 𝑇̅ = 20 + 0,02931 ∙ (300 − 20) = 28,21 ℃
𝑇𝑖 − 𝑇∞
𝑄 = 𝜌𝑉𝐶𝑃 (𝑇𝑖 − 𝑇̅) = 2000 ∙ 0,183 ∙ 1800 ∙ (300 − 28,21) = 5706285 J = 5,7 MJ
0,18 𝐿𝑧 = 0,18 m
𝐿 𝑥 = 𝐿𝑦 = = 0,09 m
2
Se calcula el nuevo número de Biot, autovalor y constante para la coordenada 𝑧 (las otras coordenadas
siguen manteniendo los mismos valores anteriores):
𝜁1𝑧 tan 𝜁1𝑧 = 𝐵𝑖𝑧 ⇒ 𝜁1𝑧 tan 𝜁1𝑧 = 0,9 ⇒ 𝜁1𝑧 = 0,82740
𝑇(0,0, 𝐿𝑧 /2) − 𝑇∞ 30 − 20
𝜃 (0,0, 𝐿𝑧 /2) = = = 0,035714
𝑇𝑖 − 𝑇∞ 300 − 20
Sustituyendo valores
0,035714
− ln ( )
1,064192 ∙ 1,11067 ∙ cos(0) ∙ cos(0) ∙ cos(0,82740 ∙ 1/2)
𝑡= = 2662,6 s ≈ 2663 s
0,62444 2 0,82740 2 −5
(2 ∙ ( ) +( )
) ∙ 1,111 ∙ 10
0,09 0,18
𝑡 = [0 horas][44 min][23 s]
La temperatura en el centro de la cara inferior o cara adiabática (cuya coordenada coincide con el
origen de coordenadas) se obtiene a partir de:
1 𝜃(0,0, 𝐿𝑧 /2) 0,035714
𝜃 (0,0, 𝐿𝑧 /2) = 𝜃𝑜 cos(0) cos(0) cos (𝜁1𝑧 ∙ ) ⇒ 𝜃𝑜 = =
2 cos(𝜁1𝑧 ∙ 1/2) cos(0,82740 ∙ 1/2)
= 0,03900
𝑇𝑜 − 𝑇∞
𝜃𝑜 = ⇒ 𝑇𝑜 = 20 + 0,039 ∙ (300 − 20) = 30,92 ℃
𝑇𝑖 − 𝑇∞
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,02640 ∙ (300 − 20) = 27,39 ℃
𝑇𝑖 − 𝑇∞
1
𝜃 = 𝜃𝑜 cos(𝜁1𝑥 ∙ 0) cos(𝜁1𝑦 ∙ 1) cos (𝜁1𝑧 ∙ ) = 0,039 ∙ cos(0,62444) ∙ cos(0,82740) = 0,02897
2
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,02897 ∙ (300 − 20) = 28,11 ℃
𝑇𝑖 − 𝑇∞
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,01737 ∙ (300 − 20) = 24,86 ℃
𝑇𝑖 − 𝑇∞
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 20 + 0,02567 ∙ (300 − 20) = 27,19 ℃
𝑇𝑖 − 𝑇∞
La temperatura media es
2
sen 𝜁1𝑥 sen 𝜁1𝑦 sen 𝜁1𝑧 sen(0,62444) sen(0,82740)
𝜃̅ = 𝜃𝑜 ∙ ∙ ∙ = 0,039 ∙ ( ) ∙ = 0,03042
𝜁1𝑥 𝜁1𝑦 𝜁1𝑧 0,62444 0,82740
𝑇̅ − 𝑇∞
𝜃̅ = ⇒ 𝑇̅ = 20 + 0,03042 ∙ (300 − 20) = 28,52 ℃
𝑇𝑖 − 𝑇∞
𝑄 = 𝜌𝑉𝐶𝑃 (𝑇𝑖 − 𝑇̅) = 2000 ∙ 0,183 ∙ 1800 ∙ (300 − 28,52) = 5699777 J = 5,7 MJ
PROBLEMA 18
(a) Determine el tiempo necesario para que el centro geométrico llegue a 50 ºC.
(b) Determine en este mismo tiempo cual es la temperatura en el centro de la cara inferior y en
el centro de la tapa superior.
(c) Determine la temperatura en cualquier punto del borde de ambas tapas.
(d) ¿Cuál es la temperatura media en ese momento?
(e) ¿Cuál es la cantidad de energía que ha disipado al ambiente?
SOLUCIÓN:
𝑘 2,5 m2
𝛼= = = 5,556 ∙ 10−7
𝜌𝐶𝑃 3000 ∙ 1500 s
ℎ 𝑥 𝐿 50 ∙ 0,1
𝐵𝑖𝑥 = = = 3,2
𝑘 2,5
𝜁1𝑥 tan 𝜁1𝑥 = 𝐵𝑖𝑥 ⇒ 𝜁1𝑥 tan 𝜁1𝑥 = 3,2 ⇒ 𝜁1𝑥 = 1,20944
ℎ 𝑟 𝑅 80 ∙ 0,05
𝐵𝑖𝑟 = = =1
𝑘 2,5
𝐽1 (𝜁1𝑟 ) 𝐽1 (𝜁1𝑟 )
𝜁1𝑟 = 𝐵𝑖𝑟 ⇒ 𝜁1𝑟 =1 ⇒ 𝜁1𝑟 = 1,25578
𝐽0 (𝜁1𝑟 ) 𝐽0 (𝜁1𝑟 )
2 𝐽1 (𝜁1𝑟 ) 2 𝐽1 (1,25578)
𝐶1𝑟 = ∙ 2 = ∙ 2 = 1,20709
𝜁1𝑟 𝐽0 (𝜁1𝑟 ) + 𝐽1 (𝜁1𝑟 ) 1,25578 𝐽0 (1,25578) + 𝐽12 (1,25578)
2
𝛼𝑡 𝜁 𝑥 2 𝛼𝑡 ) ( 𝜁1𝑟 𝑟))
2
𝜃 = ( 𝐶1𝑥 exp (−𝜁1𝑥 ) cos ( 1𝑥 )) (𝐶1𝑟 exp (−𝜁1𝑟 𝐽
𝐿2 𝐿 𝑅2 0 𝑅
Arreglando términos
𝜃
− ln ( 𝜁1𝑥 𝑥 𝜁 𝑟
)
𝐶1𝑥 𝐶1𝑟 cos ( ) 𝐽0 ( 1𝑟 )
𝑡= 𝐿 𝑅
𝜁1𝑥 2 𝜁1𝑟 2
(( ) + ( ) ) 𝛼
𝐿 𝑅
𝑇 − 𝑇∞ 50 − 25
𝜃= = = 0,25
𝑇𝑖 − 𝑇∞ 125 − 25
0,25
− ln ( 1
)
1,21471 ∙ 1,20709 ∙ cos (1,20944 ∙ ) ∙ 𝐽0 (0)
𝑡= 2 = 3645,5 ≈ 3646 s
2 2
1,20944 1,25578 −7
(( ) +( ) ) ∙ 5,556 ∙ 10
0,1 0,05
𝜁1𝑥 𝑥 𝜁 𝑟 𝜃
𝜃 = 𝜃𝑜 cos ( ) 𝐽0 ( 1𝑟 ) ⇒ 𝜃𝑜 =
𝐿 𝑅 𝜁1𝑥 𝑥 𝜁 𝑟
cos ( ) 𝐽0 ( 1𝑟 )
𝐿 𝑅
Evaluando la ecuación anterior en la coordenada del centro geométrico cuya temperatura es conocida
(𝑟 = 0, 𝑥 = 𝐿/2):
0,25
𝜃𝑜 = 1 = 0,30389
cos (1,20944 ∙ ) ∙ 𝐽0 (0)
2
𝑇𝑜 − 𝑇∞
𝜃𝑜 = ⇒ 𝑇𝑜 = 25 + 0,30389 ∙ (125 − 25) = 55,39 ℃
𝑇𝑖 − 𝑇∞
𝜁1𝑥 ∙ 𝐿 𝜁 ∙0
𝜃 = 𝜃𝑜 cos ( ) 𝐽0 ( 1𝑟 ) = 𝜃𝑜 cos(𝜁1𝑥 ) 𝐽0 (0) = 0,30389 ∙ cos(1,20944) ∙ 1 = 0,10744
𝐿 𝑅
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 25 + 0,10744 ∙ (125 − 25) = 35,74 ℃
𝑇𝑖 − 𝑇∞
𝜁1𝑥 ∙ 𝐿 𝜁 ∙𝑅
𝜃 = 𝜃𝑜 cos ( ) 𝐽0 ( 1𝑟 ) = 𝜃𝑜 cos(𝜁1𝑥 ) 𝐽0 (𝜁1𝑟 ) = 0,30389 ∙ cos(1,20944) ∙ 𝐽0 (1,25578)
𝐿 𝑅
= 0,069078
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 25 + 0,069078 ∙ (125 − 25) = 31,91 ℃
𝑇𝑖 − 𝑇∞
𝜁1𝑥 ∙ 0 𝜁 ∙𝑅
𝜃 = 𝜃𝑜 cos ( ) 𝐽0 ( 1𝑟 ) = 𝜃𝑜 cos(0) 𝐽0 (𝜁1𝑟 ) = 0,30389 ∙ 1 ∙ 𝐽0 (1,25578) = 0,19539
𝐿 𝑅
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 25 + 0,19539 ∙ (125 − 25) = 44,54 ℃
𝑇𝑖 − 𝑇∞
Temperatura media:
𝐿 𝑅
1 2𝜋
𝜃̅ = ∭ 𝜃𝑑𝑉 = (∫ 𝜃 𝑑𝑥 )(∫ 𝜃𝑟 ∙ 𝑟𝑑𝑟)
𝑉 𝑉 𝜋𝑅2 𝐿 0 𝑥 0
𝐿 𝑅
2𝜋 𝜁1𝑥 𝑥 𝜁 𝑟
= 2
(∫ 𝜃𝑜𝑥 cos ( ) 𝑑𝑥) (∫ 𝜃𝑜𝑟 𝐽0 ( 1𝑟 ) 𝑟𝑑𝑟) =
𝜋𝑅 𝐿 0 𝐿 0 𝑅
𝐿 𝑅 1 1
2 𝜁1𝑥 𝑥 𝜁1𝑟 𝑟
= 𝜃
⏟ 𝜃 ∫ cos ( ) 𝑑𝑥 ∫ 𝐽 ( ) 𝑟𝑑𝑟 = 2𝜃 ∫ cos(𝜁 𝑥 ) 𝑑𝑥 ∫ 𝐽0 (𝜁1𝑟 𝑟)𝑟𝑑𝑟 =
𝑅 2 𝐿 𝑜𝑥 𝑜𝑟 0 𝐿 0
0
𝑅 𝑜
0
1𝑥
0
𝜃𝑜
𝑇̅ − 𝑇∞
𝜃̅ = ⇒ 𝑇̅ = 25 + 0,19165 ∙ (125 − 25) = 44,17 ℃
𝑇𝑖 − 𝑇∞
Calor transferido:
𝑄 = 𝜌𝑉𝐶𝑃(𝑇𝑖 − 𝑇̅) = 3000 ∙ (𝜋 ∙ 0,052 ∙ 0,1) ∙ 1500 ∙ (125 − 44,17) = 285677 J = 0,286 MJ
PROBLEMA 19
Considere una esfera de 20 cm de diámetro que tiene una temperatura inicial de 275 ºC que se sumerge
en un fluido en reposo a temperatura constante de 𝑇∞ = 25 ℃. Suponga que la esfera cae a la
velocidad terminal o en estado estacionario. En la región de Newton, el factor de fricción viene dado
por:
Las propiedades del sólido y las del fluido se reportan en la siguiente tabla:
SOLUCIÓN:
1
0 = 𝜌𝑠 𝑉𝑠 𝑔 − 𝜌𝑓 𝑉𝑓 𝑔 − 𝐶𝑓𝜌𝑓 𝑢2∞ 𝐴𝐹
2
𝜋
𝜌𝑠 𝑉𝑠 𝑔 − 𝜌𝑓 𝑉𝑓 𝑔 2(𝜌𝑠 − 𝜌𝑓 )𝑉𝑔 2(𝜌𝑠 − 𝜌𝑓 ) (6 𝐷 3 ) 𝑔 4𝐷 𝜌𝑠
𝑢2∞ = 1 = = 𝜋 = ( − 1) 𝑔
𝐶𝑓𝜌𝑓 𝐴𝐹 𝐶𝑓𝜌𝑓 𝐴𝐹 𝐶𝑓𝜌𝑓 ( 𝐷 2 ) 3𝐶𝑓 𝜌𝑓
2 4
Sustituyendo los valores
4𝐷 𝜌𝑠 4 ∙ 0,2 1138
𝑢∞ = √ ( − 1) 𝑔 = √ ( − 1) ∙ 9,8 = 0,9 m/s
3𝐶𝑓 𝜌𝑓 3 ∙ 0,445 1000
̅̅̅̅̅̅
𝑁𝑢 𝐷 ∙ 𝑘𝑓 692,8 ∙ 0,5 W
ℎ= = = 1732 2
𝐷 0,2 m ∙℃
ℎ𝑅 1732 ∙ 0,1
𝐵𝑖𝑅 = = = 4,33
𝑘𝑠 40
𝑘𝑠 0,5 m2
𝛼= = = 1,406 ∙ 10−5
𝜌𝑠 𝐶𝑃𝑠 1138 ∙ 2500 s
𝜁1 𝜁1
1− = 𝐵𝑖 ⇒ 1− = 4,33 ⇒ 𝜁1 = 2,4980
tan 𝜁1 tan 𝜁1
𝑇𝑜 − 𝑇∞ 50 − 25
𝜃𝑜 = = = 0,1
𝑇𝑖 − 𝑇∞ 275 − 25
𝜃 0,1
𝛼𝑡 − ln ( 𝑜 ) − ln ( )
𝐶1 1,7450
𝜃𝑜 = 𝐶1 exp (−𝜁12 2 ) ⇒ 𝑡= = = 326 s
𝑅 𝜁1 2 −5 2,4980 2
𝛼( ) 1,406 ∙ 10 ∙ ( )
𝑅 0,1
𝑧 180
𝑡= = = 200 s
𝑢∞ 0,9
𝑇𝑜 − 𝑇∞
𝜃𝑜 = ⇒ 𝑇𝑜 = 25 + 0,30183 ∙ (275 − 25) = 100,46 ℃
𝑇𝑖 − 𝑇∞
𝛼𝑡 𝜁𝑅
𝜃𝑅 = 𝜃 | = 𝐶1 exp (−𝜁12 2
) cos ( 1 ) = 𝜃𝑜 cos 𝜁1 = 0,30183 ∙ cos(2,4980) = 0,07251
𝑟=𝑅 𝑅 𝑅
𝑇𝑅 − 𝑇∞
𝜃𝑅 = ⇒ 𝑇𝑅 = 25 + 0,07251 ∙ (275 − 25) = 43,13 ℃
𝑇𝑖 − 𝑇∞
𝑟
1 1 𝑅
𝛼𝑡 sen (𝜁1 𝑅)
𝜃̅ = ∭ 𝜃𝑑𝑉 = 4 2
∫ 𝐶1 exp ( −𝜁1 2 ) 𝑟 ∙ 4𝜋𝑟 2 𝑑𝑟 =
𝑉 𝑉 3
𝜋𝑅 0 𝑅 𝜁1
3 𝑅
𝑇̅ − 𝑇∞
𝜃̅ = ⇒ 𝑇̅ = 25 + 0,15094 ∙ (275 − 25) = 62,73 ℃
𝑇𝑖 − 𝑇∞
4 4
𝑄 = 𝜌𝑠 𝑉𝐶𝑃𝑠 (𝑇𝑖 − 𝑇̅) = 𝜌𝑠 ( 𝜋𝑅 3 ) 𝐶𝑃𝑠 (𝑇𝑖 − 𝑇̅) = 1138 ∙ ( 𝜋 ∙ 0,13 ) ∙ 2500 ∙ (275 − 62,73)
3 3
= 2529589 J = 2529,6 kJ
̅̅̅̅̅̅
𝑁𝑢 𝐷 ∙ 𝑘𝑓 83,56 ∙ 0,25 W
ℎ= = = 104,5 2
𝐷 0,2 m ∙℃
Como se reduce mucho el valor del coeficiente convectivo, puede ser que el problema se aproxime
como resistencia interna despreciable por obtener un número de Biot bajo:
4 3
𝑉 3 𝜋𝑅 𝑅 0,1 ℎ𝐿𝐶 104,5 ∙ 0,1/3
𝐿𝐶 = = 2 = = ⇒ 𝐵𝑖𝐶 = = = 0,09 < 0,1 ⇒ es de RID
𝐴 𝑆 4𝜋𝑅 3 3 𝑘 40
𝜕𝑇 𝑑𝑇 ℎ𝐴 𝑆
𝜌𝐶𝑃 𝑉 = −ℎ𝐴 𝑆 (𝑇 − 𝑇∞ ) ⇒ =− 𝑑𝑡
𝜕𝑡 𝑇 − 𝑇∞ 𝜌𝐶𝑃 𝑉
𝑇 𝑡
𝑑𝑇 ℎ𝐴 𝑆 𝑇 − 𝑇∞ ℎ𝐴 𝑆 𝜌𝐶𝑃 𝑉 𝑇 − 𝑇∞
∫ =∫ − 𝑑𝑡 ⇒ ln ( )= −− 𝑡 ⇒ 𝑡=− ln ( )
𝑇𝑖 𝑇 − 𝑇∞ 0 𝜌𝐶𝑃 𝑉 𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃 𝑉 ℎ𝐴 𝑆 𝑇𝑖 − 𝑇∞
La velocidad que tiene la esfera se halla a partir de este valor de número de Reynolds:
m2
𝜈 = 0,1 St = 1 ∙ 10 −5
s
𝑢∞ 𝐷 1000 ∙ 1 ∙ 10−5
𝑅𝑒𝐷 = = 1000 ⇒ 𝑢∞ = = 0,05 m/s
𝜈 0,2
¿A qué profundidad debe enterrarse una tubería por la que circula agua si el piso que inicialmente
está a 15 ºC se enfría repentinamente en la superficie a −5 ℃ ? Suponga que el piso tiene
conductividad térmica 2,5 W/mºC, densidad 2500 kg/m3 y calor específico 3950 J/kgºC. La finalidad
es la de evitar que el agua llegue a congelarse antes de 24 horas. La temperatura de fusión del agua a
presión atmosférica es de 0 ºC.
SOLUCIÓN:
𝑇 − 𝑇𝑠 𝑥 2 𝜂 −𝜇2
𝜃≡ = erf ( ) donde erf (𝜂) = ∫ 𝑒 𝑑𝜇
𝑇𝑖 − 𝑇𝑠 2√𝛼𝑡 √𝜋 0
𝑇 − 𝑇𝑠 0 − (−5) 𝑥
𝜃= = = 0,25 ⇒ = erf −1 (𝜃) = erf −1 (0,25)
𝑇𝑖 − 𝑇𝑠 15 − (−5) 2√𝛼𝑡
Utilizando la Función Gaussiana de Error inversa erf −1 (𝜂) se halla el valor del argumento para
encontrar la distancia:
√𝜋 𝜋 7𝜋 2 5 127𝜋 3 7 4369𝜋 4 9
erf −1 (𝜂) = (𝜂 + 𝜂 3 + 𝜂 + 𝜂 + 𝜂 )
2 12 480 40320 5806080
√𝜋 𝜋 7𝜋 2 127𝜋 3 4369𝜋 4
erf −1 (0,25) = (0,25 + 0,253 + 0,255 + 0,257 + 0,259 ) = 0,22531
2 12 480 40320 5806080
La distancia a la cual debe enterrarse la tubería tal que la temperatura del suelo es la temperatura de
fusión del agua en 24 horas es:
𝑥
= erf −1 (𝜃) ⇒ 𝑥 = 2√𝛼𝑡 ∙ erf −1 (𝜃) = 2√1,406 ∙ 10−5 ∙ 24 ∙ 3600 ∙ 0,22531 = 0,497 m
2√𝛼𝑡
PROBLEMA 21
Una viga cilíndrica semi-infinita de aluminio de 5 cm de diámetro está inicialmente a una temperatura
uniforme de 200°C. Este cilindro es sometido de forma rápida, a una condición de contorno
convectiva a 70 ºC cuyo coeficiente convectivo es de 5250 W/m2 K. La difusividad térmica del
cilindro es 8,4.10-5 m2 /s y la conductividad es 215 W/mK.
(a) Calcúlese la temperatura en la superficie del cilindro a 10 cm de la base, 10 segundos después
de la exposición al ambiente.
(b) Si el coeficiente convectivo se reduce 10 veces, calcule la temperatura en la superficie del
cilindro a 10 cm de la base, 10 segundos después de la exposición al ambiente.
La solución analítica de un sólido semi-infinito en un medio convectivo es:
𝑇 − 𝑇∞ 𝑥 ℎ𝑥 ℎ 2 𝛼𝑡 𝑥 ℎ √𝛼𝑡
𝜃≡ = erf ( ) + exp ( + 2 )(1 − erf ( + ))
𝑇𝑖 − 𝑇∞ 2√𝛼𝑡 𝑘 𝑘 2√𝛼𝑡 𝑘
SOLUCIÓN:
(a) Se resuelve este problema por superposición de las soluciones de un medio semi-infinito y la
solución de un cilindro infinito
ℎ𝑅 5250 ∙ 0,025
𝐵𝑖 = = = 0,6105
𝑘 215
Con este valor del número de Biot, se encuentran el valor del 𝐶1 𝑦 𝜁1:
𝐽1 (𝜁𝑛 ) 𝐽1 (𝜁𝑛 )
𝜁𝑛 = 𝐵𝑖 ⇒ 𝜁𝑛 = 0,6105 ⇒ 𝜁1 = 1,0260
𝐽0 (𝜁𝑛 ) 𝐽0 (𝜁𝑛 )
2 𝐽1 (𝜁𝑛) 2 𝐽1 (1,0260)
𝐶𝑛 = ∙ 2 = ∙ 2 = 1,1366
𝜁𝑛 𝐽0 (𝜁𝑛 ) + 𝐽1 (𝜁𝑛) 1,0260 𝐽0 (1,0260) + 𝐽12 (1,0260)
2
𝛼𝑡 8,4 ∙ 10−5 ∙ 10
𝐹𝑜 = = = 1,344 > 0,2
𝑅2 0,0252
𝛼𝑡 𝜁𝑟
𝜃 = 𝐶1 exp (−𝜁12 2
) 𝐽0 ( 1 ) = 1,1366 ∙ exp(−1,02602 ∙ 1,344) 𝐽0 (1,0260 ∙ 1) = 0,2081
𝑅 𝑅
La solución semi-infinita es
𝑥 ℎ𝑥 ℎ 2 𝛼𝑡 𝑥 ℎ √𝛼𝑡
𝜃 = erf ( ) + exp ( + 2 ) (1 − erf ( + ))
2√𝛼𝑡 𝑘 𝑘 2√𝛼𝑡 𝑘
𝑥 0,1 𝑥
= = 1,7252 ⇒ erf ( ) = erf(1,7252) = 0,9853
2√𝛼𝑡 2√8,4 ∙ 10−5 ∙ 10 2√𝛼𝑡
ℎ √𝛼𝑡 5250√8,4 ∙ 10−5 ∙ 10
= = 0,7077
𝑘 215
𝑥 ℎ√𝛼𝑡 𝑥 ℎ √𝛼𝑡
+ = 1,7252 + 0,7077 = 2,4392 ⇒ erf ( + ) = erf(2,4392) = 0,9994
2 √𝛼𝑡 𝑘 2√𝛼𝑡 𝑘
ℎ𝑥 ℎ 2 𝛼𝑡 5250 ∙ 0,1 52502 ∙ 8,4 ∙ 10 −5 ∙ 10
+ 2 = + = 2,9427
𝑘 𝑘 215 2152
𝑇 − 𝑇∞
𝜃≡ = 𝜃cilindro infinito 𝜃semi−infinito = 0,2081 ∙ 0,9967 = 0,2074
𝑇𝑖 − 𝑇∞
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 70 + 0,2074 ∙ (200 − 70) = 97,0 ℃
𝑇𝑖 − 𝑇∞
(b) Como el coeficiente de convección se reduce mucho, el problema puede convertirse en resistencia
interna despreciable. Calculando la longitud característica y el número de Biot:
𝑇 − 𝑇∞ ℎ𝐴 𝑆 525 ∙ 10
𝜃≡ = exp (− 𝑡) = exp (− ) = 0,8487
𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃𝑉 (215/8,4 ∙ 10 −5 ) ∙ 0,0125
𝑇 − 𝑇∞
𝜃= ⇒ 𝑇 = 70 + 0,8487 ∙ (200 − 70) = 180,3 ℃
𝑇𝑖 − 𝑇∞
PROBLEMA 22
Demuestre que el perfil de temperaturas en una caja es el producto de los perfiles unidimensionales
SOLUCIÓN:
𝜕𝜃
= 𝛼∇2 𝜃(𝑡, 𝑥, 𝑦, 𝑧)
𝜕𝑡
𝜕2 𝜕2 𝜕2
∇2 = + +
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2
Sustituyendo
𝜕
(𝜗 (𝑡,𝑥) 𝜗𝑦 (𝑡, 𝑦)𝜗𝑧 (𝑡, 𝑧) ) = 𝛼∇2 (𝜗𝑥 (𝑡, 𝑥)𝜗𝑦 (𝑡, 𝑦) 𝜗𝑧 (𝑡, 𝑧))
𝜕𝑡 𝑥
𝜕𝜗𝑥 (𝑡, 𝑥) 𝜕𝜗𝑦 (𝑡, 𝑦) 𝜕𝜗𝑧 (𝑡,𝑧)
𝜗𝑦 (𝑡, 𝑦) 𝜗𝑧 (𝑡, 𝑧) + 𝜗𝑥 (𝑡, 𝑥) 𝜗𝑧 (𝑡, 𝑧) + 𝜗𝑥 (𝑡, 𝑥)𝜗𝑦 (𝑡, 𝑦) =
𝜕𝑡 𝜕𝑡 𝜕𝑡
Para que la ecuación anterior sea consistente matemáticamente, cada una de las funciones debe ser
idénticamente cero ya que cada término es independiente entre sí, implicando las soluciones
unidimensionales en el volumen 𝑉 = {(𝑥, 𝑦, 𝑧): 0 < 𝑥 < 𝑎, 0 < 𝑦 < 𝑏, 0 < 𝑧 < 𝑐}
Y la condición de convección en las paredes para 𝜃 (𝑡, 𝑥, 𝑦, 𝑧) implica la condición de borde para
cada una de las soluciones unidimensionales, cancelando términos ambos lados de la igualdad:
La temperatura media
1 1
𝜃̅(𝑡) = ∭ 𝜃 (𝑡, 𝑥, 𝑦, 𝑧) 𝑑𝑥𝑑𝑦𝑑𝑧 = ∭ 𝜗 (𝑡, 𝑥) 𝜗𝑦 (𝑡, 𝑦)𝜗𝑧 (𝑡, 𝑧) 𝑑𝑥𝑑𝑦𝑑𝑧 =
𝑉 𝑉 𝑎𝑏𝑐 𝑉 𝑥
1 𝑎 1 𝑏 1 𝑐
̅̅̅
= ( ∫ 𝜗𝑥 (𝑡, 𝑥)𝑑𝑥) ( ∫ 𝜗𝑦 (𝑡, 𝑦) 𝑑𝑦) ( ∫ 𝜗𝑧 (𝑡, 𝑧)𝑑𝑧) = 𝜗 ̅̅̅ ̅̅̅
𝑥 (𝑡)𝜗𝑦(𝑡)𝜗𝑧 (𝑡)
𝑎 0 𝑏 0 𝑐 0
𝑇̅(𝑡) − 𝑇∞
𝑄 = 𝜌𝑉𝐶𝑃(𝑇𝑖 − 𝑇̅(𝑡)) = 𝜌𝑉𝐶𝑃 ((𝑇𝑖 − 𝑇∞ ) − (𝑇̅(𝑡) − 𝑇∞ )) = 𝜌𝑉𝐶𝑃 (𝑇𝑖 − 𝑇∞ ) (1 − ( ))
𝑇𝑖 − 𝑇∞
𝑄
𝑄 = 𝜌𝑉𝐶𝑃(𝑇𝑖 − 𝑇∞ ) (1 − 𝜃̅ (𝑡)) ⇒ Φ= = 1 − 𝜃̅ (𝑡) ⇒ 𝜃̅ (𝑡) = 1 − Φ
𝑄𝑚𝑎𝑥
̅̅̅
Como 𝜗 ̅̅̅ ̅̅̅
𝑥 (𝑡) = 1 − Φ𝑥 , 𝜗𝑦 (𝑡) = 1 − Φ𝑦 , 𝜗𝑧 (𝑡) = 1 − Φ𝑧 entonces
1 − Φ = 𝜃̅ (𝑡) = 𝜗
̅̅̅ ̅̅̅ ̅̅̅
𝑥 (𝑡)𝜗𝑦 (𝑡)𝜗𝑧(𝑡) = (1 − Φ𝑥 )(1 − Φ𝑦 )(1 − Φ𝑧 )
PROBLEMA 23
Una tubería de acero comercial (de conductividad térmica 𝑘 𝑡 = 51,3 W/m ∙ K) de diámetro interno
3,5 cm y diámetro externo 4,0 cm, está empotrada en concreto (de conductividad térmica 𝑘 𝐶 =
1,40 W/m ∙ K) de corte transversal en forma de cuadrado de arista 6,0 cm. La tubería transporta agua
a una velocidad promedio de 1,5 m/s y una temperatura promedio de 100 °C. Si en el exterior se
encuentra aire a 10 °C circulando transversalmente sobre la columna a una velocidad promedio de 10
m/s, calcule el flujo de calor por metro de tubería que pierde el agua al ambiente.
SOLUCIÓN:
Conducción bidimensional:
2𝜋𝐿
𝑆= con 𝑤 > 𝐷𝑜
ln(1,08𝑤/𝐷𝑜 )
𝑆 2𝜋 2𝜋
= = = 13,024 m/m
𝐿 ln(1,08𝑤/𝐷𝑜 ) ln(1,08 ∙ 6/4)
Flujo interno:
Calculando el número de Reynolds
𝜌𝑢𝐷𝑖 957,9 ∙ 1,5 ∙ 0,035
𝑅𝑒𝐷 = = = 180250 (turbulento)
𝜇 0,000279
𝜇𝐶𝑃 0,000279 ∙ 4217
𝑃𝑟 = = = 1,730
𝑘 0,680
Se calcula el factor de fricción en la tubería (tubería lisa)
Luego el número de Nusselt válido para 3000 ≤ 𝑅𝑒𝐷 ≤ 5 ∙ 106 , 0,5 ≤ 𝑃𝑟 ≤ 2000 , 𝐿/𝐷𝑖 ≥ 10
𝑓 0,01596
(𝑅𝑒𝐷 − 1000) 𝑃𝑟 ∙ (180250 − 1000) ∙ 1,73
̅̅̅̅̅̅
𝑁𝑢 8 8
𝐷= = = 494,9
𝑓 1/2 2/3 0,01596 1/2 2/3
1 + 12,7 ( ) (𝑃𝑟 − 1) 1 + 12,7 ∙ ( (
) ∙ 1,73 − 1 )
8 8
ℎ 𝑖𝐷𝑖 ̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 494,9 ∙ 0,680 W
̅̅̅̅̅̅
𝑁𝑢 𝐷= ⇒ ℎ𝑖 = = = 9615 2
𝑘 𝐷𝑖 0,035 m ∙℃
Flujo externo:
̅̅̅̅̅̅
𝑁𝑢 𝑚
𝐷 = 𝐶𝑅𝑒𝐷 𝑃𝑟
1/3
Aire
𝐶 = 0,102 , 𝑚 = 0,675
Propiedades a 𝑇𝑓
1 ≤ 𝑅𝑒𝐷 ≤ 106 , 0,7 ≤ 𝑃𝑟 ≤ 500
Se calcula el número de Reynolds del aire, número de Prandtl y Nusselt. Primero se buscan las
propiedades del aire a la temperatura de la película (𝑇𝑓 ). Como no se tiene el valor de la temperatura
en la superficie, entonces como primera estimación se considera que 𝑇𝑓 ~𝑇∞ = 10 °C = 283 K. Las
propiedades del aire a esta temperatura son (interpolando linealmente):
𝜌 = 1,2407 kg/m3
𝐶𝑃 = 1007 J/kg ∙ ℃
aire a 283 K
𝑘 = 0,02494 W/m ∙ ℃
{𝜇 = 176,1 ∙ 10−7 Pa ∙ s
Calculando los números de Reynolds, Prandtl y Nusselt
𝜌𝑉∞ 𝐷 1,2407 ∙ 10 ∙ 0,06
𝑅𝑒𝐷 = = = 42273
𝜇 176,1 ∙ 10−7
ℎ𝑜𝐷 ̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 120,7 ∙ 0,02494 W
̅̅̅̅̅̅
𝑁𝑢 𝐷= ⇒ ℎ𝑜 = = = 50,2 2
𝑘 𝐷 0,06 m ∙℃
Calor transferido:
La resistencia total de transferencia de calor es
1 1 ln (𝐷𝑜 /𝐷𝑖 ) 1 1 1 ln( 𝐷𝑜 /𝐷𝑖 ) 1 1
= + + + = + + 𝑆 +
𝑈𝐴 ℎ 𝑖 𝐴𝑖 2𝜋𝑘 𝑡 𝐿 𝑆𝑘 𝐶 ℎ 𝑜 𝐴𝑜 ℎ 𝑖 (𝜋𝐷𝑖 𝐿) 2𝜋𝑘 𝑡 𝐿 𝐿 ( ) 𝑘 𝐶 ℎ 𝑜 (4𝑤𝐿)
𝐿
Escribiendo la expresión anterior por unidad de longitud
𝐿 1 ln(𝐷𝑜 /𝐷𝑖 ) 1 1
= + + 𝑆 + =
(
𝑈𝐴 ℎ 𝑖 𝜋𝐷𝑖 ) 2𝜋𝑘 𝑡 ( ) 𝑘 𝐶 ℎ 𝑜 4𝑤)
(
𝐿
1 ln(4/3,5) 1 1 K∙m
= + + + = 0,1392
9615 ∙ 𝜋 ∙ 0,035 2𝜋 ∙ 51,3 13,024 ∙ 1,4 50,2 ∙ 4 ∙ 0,06 W
El calor transferido es
𝑄̇ 𝑇𝑖 − 𝑇∞ 100 − 10
𝑄̇ = 𝑈𝐴 (𝑇𝑖 − 𝑇∞ ) ⇒ = = = 646,5 W/m
𝐿 𝐿/𝑈𝐴 0,1392
Temperatura de la superficie externa:
La temperatura de la superficie externa de la columna de concreto viene dada por
𝑄̇ 𝑄̇/𝐿 646,5
𝑄̇ = ℎ 𝑜 𝐴 𝑜 (𝑇𝑠 − 𝑇∞ ) ⇒ 𝑇𝑠 = 𝑇∞ + = 𝑇∞ + = 10 + = 63,7 °C
ℎ 𝑜 𝐴𝑜 ℎ 𝑜 (4𝑤) 50,2 ∙ 4 ∙ 0,06
Repetir los cálculos con una segunda iteración:
Ahora que ya se tiene el valor de la temperatura de la superficie externa de la columna de concreto,
es posible encontrar la temperatura de la película:
𝜌 = 1,1281 kg/m3
𝑇𝑠 + 𝑇∞ 63,7 + 10 𝐶𝑃 = 1007 J/kg ∙ ℃
aire a 𝑇𝑓 = = = 36,8 ℃ ≈ 310 K
2 2 𝑘 = 0,02704 W/m ∙ ℃
{𝜇 = 189,3 ∙ 10−7 Pa ∙ s
Calculando los números de Reynolds, Prandtl y Nusselt y luego el coeficiente convectivo externo
𝜌𝑉∞ 𝐷 1,1281 ∙ 10 ∙ 0,06
𝑅𝑒𝐷 = = = 35756
𝜇 189,3 ∙ 10−7
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 107,5 ∙ 0,02704 W
ℎ𝑜 = = = 48,5 2
𝐷 0,06 m ∙℃
Calculando nuevamente el valor de la resistencia total de transferencia de calor
𝐿 1 ln(4/3,5) 1 1 K∙m
= + + + = 0,1421
𝑈𝐴 9615 ∙ 𝜋 ∙ 0,035 2𝜋 ∙ 51,3 13,024 ∙ 1,4 48,5 ∙ 4 ∙ 0,06 W
El calor transferido es
𝑄̇ 𝑇𝑖 − 𝑇∞ 100 − 10
𝑄̇ = 𝑈𝐴 (𝑇𝑖 − 𝑇∞ ) ⇒ = = = 633,3 W/m
𝐿 𝐿/𝑈𝐴 0,1421
La temperatura de la superficie externa de la columna de concreto viene dada por
𝑄̇ /𝐿 646,5
𝑇𝑠 = 𝑇∞ + = 10 + = 63,7 °C
ℎ 𝑜 (4𝑤) 50,2 ∙ 4 ∙ 0,06
PROBLEMA 24
Calcule:
SOLUCIÓN:
T = 20 ºC
h = 42 W/m2 ºC
Ti = 80 ºC
To = ?
ui = 1,2 m/s
L = 80 m
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 394,7 ∙ 0,5 W
ℎ𝑖 = = = 1973,7 2
𝐷𝑖 0,1 m ∙℃
El flujo másico es
𝜋 𝜋 kg
𝑚̇ = 𝜌𝑢𝐴 𝑡 = 𝜌𝑢 ( 𝐷𝑖2 ) = 800 ∙ 1,2 ∙ ∙ 0,12 = 7,54
4 4 s
Luego se halla el coeficiente global:
1 1 ln(𝐷𝑜 /𝐷𝑖 ) 1 K
= ( + + ) = 8,8518 ∙ 10−4
80 1973,7 ∙ 𝜋 ∙ 0,1 2𝜋 ∙ 50 42 ∙ 𝜋 ∙ (0,1 + 2 ∙ 0,005) W
1 W
𝑈𝐴 = −4 = 1129,7
8,8518 ∙ 10 K
𝑈𝐴 1129,7 W
𝑈𝐴 = 𝑈𝑜 𝐴 𝑜 = 𝑈𝑖 𝐴 𝑖 ⇒ 𝑈𝑜 = = = 40,86 2
𝐴𝑜 𝜋 ∙ 0,11 ∙ 80 m ∙℃
𝑄̇ = 𝑚̇𝐶𝑃(𝑇𝑖 − 𝑇𝑜 )
{ ⇒ 𝑄̇ = 𝑚̇𝐶𝑃 (𝑇𝑖 − 𝑇𝑜 ) = 𝑈𝐴∆𝑇𝑚𝑙
𝑄̇ = 𝑈𝐴∆𝑇𝑚𝑙
(𝑇𝑖 − 𝑇∞ ) − (𝑇𝑜 − 𝑇∞ ) 1
𝑄̇ = 𝑚̇𝐶𝑃(𝑇𝑖 − 𝑇𝑜 ) = 𝑈𝐴 𝑇 − 𝑇∞ ⇒ 𝑚̇ 𝐶𝑃 = 𝑈𝐴 𝑇 − 𝑇∞
ln ( 𝑖 ) ln ( 𝑖 )
𝑇𝑜 − 𝑇∞ 𝑇𝑜 − 𝑇∞
𝑇𝑖 − 𝑇∞ 𝑈𝐴 𝑇𝑜 − 𝑇∞ 𝑈𝐴
ln ( )= ⇒ = exp (− )
𝑇𝑜 − 𝑇∞ 𝑚̇ 𝐶𝑃 𝑇𝑖 − 𝑇∞ 𝑚̇ 𝐶𝑃
𝑈𝐴 1129,7
𝑇𝑜 = 𝑇∞ + (𝑇𝑖 − 𝑇∞ ) exp (− ) = 20 + (80 − 20) exp (− ) = 77,08 ℃
𝑚̇𝐶𝑃 7,54 ∙ 3000
Por lo tanto como la variación de temperatura del fluido es pequeña, es adecuada la aproximación
con propiedades constantes.
𝑄̇ 66118
𝑄̇ = ℎ 𝑜 𝐴𝑜 (𝑇𝑠𝑜 − 𝑇∞ ) ⇒ 𝑇𝑠𝑜 = 𝑇∞ + = 20 + = 76,94 ℃
ℎ 𝑜 𝐴𝑜 42 ∙ 𝜋 ∙ (0,1 + 2 ∙ 0,005) ∙ 80
T = 20 ºC
w= ?
h = 42 W/m2 ºC
Ti = 80 ºC
ui = 1,2 m/s To = 79 ºC
L = 80 m
𝑄̇ 22620 W
𝑄̇ = 𝑈𝐴∆𝑇𝑚𝑙 ⇒ 𝑈𝐴 = = = 380,17
∆𝑇𝑚𝑙 59,50 K
Como el sistema tiene una resistencia más a la transferencia de calor, debido al aislante añadido
entonces
𝐷𝐴 − 𝐷𝑜 0,1331 − 0,11
𝐷𝐴 = 𝐷𝑜 + 2𝑤 ⇒ 𝑤= = = 0,01155 m ≈ 12 mm
2 2
PROBLEMA 25
Aceite de motor limpio (sus propiedades se muestran en la tabla) con un flujo másico de 1 kg/s y
temperatura de entrada de 120 ºC entra por una tubería de acero comercial (la conductividad es de
𝑘 𝑡 = 51,9 W/m ∙ K) de diámetro externo 5 cm y 5 mm de espesor. La tubería tiene una longitud de
50 m. Esta se recubre con un aislante de conductividad 0,2 W/m ∙ K. Alrededor de la tubería circula
aire a 10 m/s y una temperatura de 27 ºC.
(a) ¿Cuál debe ser el espesor mínimo de aislamiento para que la temperatura de la salida del
aceite no sea menor de 118 ºC? ¿Cuál es la temperatura de la superficie externa del aislante?
(b) ¿Cuál sería la temperatura del aceite de salida si la tubería no se recubre con el aislante?
Aire a 1 atm
27 ºC 1,1614 1007 0,0263 184,6 ∙ 10−7
Nota: La viscosidad del aceite no es lineal, puede interpolarse con la siguiente fórmula:
𝑇 𝑛 ln (𝜇 2 /𝜇1 )
𝜇 = 𝜇1 ( ) donde 𝑛= con 𝑇1 < 𝑇 < 𝑇2
𝑇1 ln (𝑇2 /𝑇1 )
SOLUCIÓN:
La temperatura media entre entrada y salida para evaluar las propiedades del aceite es
𝑇𝑖 + 𝑇𝑜 120 + 118
𝑇̅ = = = 119 ℃ = 392 K
2 2
Interpolando linealmente las propiedades excepto la viscosidad:
𝜌 = 829,5 kg/m3
𝐶𝑃 = 2303 J/kg ∙ ℃
𝑘 = 0,1348 W/m ∙ ℃
𝜌𝑢𝐷𝑖 𝜋 4 𝐷𝑖 4 𝑚̇ 4 1
𝑅𝑒𝐷 = = 𝜌𝑢 ( 𝐷𝑖2 ) ∙ 2 ∙ = ∙ = ∙ = 3032
𝜇 ⏟ 4 𝜋𝐷𝑖 𝜇 𝜋 𝜇𝐷𝑖 𝜋 0,0105 ∙ 0,04
𝑚̇
𝑓 0,0454
(𝑅𝑒𝐷 − 1000) 𝑃𝑟 ∙ (3032 − 1000) ∙ 179,4
̅̅̅̅̅̅
𝑁𝑢 = 8 = 8 = 67,86
𝐷
𝑓 1/2 2/3 0,0454 1/2 2/3
1 + 12,7 ( ) (𝑃𝑟 − 1) 1 + 12,7 ∙ ( ) ∙ (179,4 − 1)
8 8
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 67,86 ∙ 0,1348 W
ℎ𝑖 = = = 228,7 2
𝐷𝑖 0,04 m ∙℃
Como se conoce la temperatura de salida del aceite es posible encontrar el valor del calor que se
escapa al ambiente
𝑄̇ 4606 W
𝑄̇ = 𝑈𝐴∆𝑇𝑚𝑙 ⇒ 𝑈𝐴 = = = 50,07
∆𝑇𝑚𝑙 92 K
Diámetro de aislante requerido:
En la que el diámetro del aislante es desconocido, y por lo tanto también el coeficiente convectivo
externo ℎ 𝑜 es desconocido. Se requiere un método de tanteo para resolver el diámetro de aislante 𝐷𝐴
o formular una ecuación que permita encontrar dicho valor. Utilizando la correlación de Churchill
para el coeficiente convectivo externo
−1/4 4/5
ℎ 𝑜 𝐷𝐴 1/2 0,4 2/3 𝑅𝑒𝐷 5/8
̅̅̅̅̅̅
𝑁𝑢𝐷 ≡ 1/3
= 0,3 + 0,62𝑅𝑒𝐷 𝑃𝑟 (1 + ( ) ) (1 + ( ) )
𝑘 𝑃𝑟 282000
𝜌𝑢∞ 𝐷𝐴 1,1614 ∙ 10 ∙ 𝐷𝐴
𝑅𝑒𝐷 = = = 629144𝐷𝐴
𝜇 184,6 ∙ 10−7
El producto ℎ 𝑜 𝐷𝐴 es una función exclusiva sólo de 𝐷𝐴. Esto permite resolver la ecuación:
𝐷𝐴 = 0,156 m
El espesor de aislante es
𝐷𝐴 − 𝐷𝑜 0,156 − 0,05
𝐷𝐴 = 𝐷𝑜 + 2𝑤 ⇒ 𝑤= = = 0,053 m ≈ 5,3 cm
2 2
𝜌𝑢∞𝐷𝐴
𝑅𝑒𝐷 = = 629144 ∙ 0,156 = 98231
𝜇
−1/4 4/5
0,4 2/3 98231 5/8
̅̅̅̅̅̅ 1/2 1/3
𝑁𝑢𝐷 = 0,3 + 0,62 ∙ 98231 ∙ 0,707 (1 + ( ) ) (1 + ( ) ) = 212,4
0,707 282000
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 212,4 ∙ 0,0263 W
ℎ𝑜 = = = 35,8 2
𝐷𝐴 0,156 m ∙℃
𝑄̇ 4606
𝑄̇ = ℎ 𝑜 𝐴𝐴 (𝑇𝑠𝑜 − 𝑇∞ ) ⇒ 𝑇𝑠𝑜 = 𝑇∞ + = 27 + = 32,25 ℃
ℎ 𝑜 ∙ 𝜋𝐷𝐴 𝐿 35,8 ∙ 𝜋 ∙ 0,156 ∙ 50
1 1 1 ln(𝐷𝑜 /𝐷𝑖 ) 1
= ( + + )
𝑈𝐴 𝜋𝐿 ℎ 𝑖 𝐷𝑖 2𝑘 𝑡 ℎ 𝑜 𝐷𝑜
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 103,3 ∙ 0,0263 W
ℎ𝑜 = = = 54,4 2
𝐷𝑜 0,05 m ∙℃
El coeficiente convectivo interno no se conoce su valor ya que la temperatura de salida del aceite es
desconocida. Sin embargo puede tomarse ℎ 𝑖 como aproximación inicial, el mismo valor calculado
anteriormente con 𝑇𝑜 = 118 ℃ y todas las propiedades del aceite evaluadas entre 118 ºC y 120 ºC:
1 W
𝑈𝐴 = = 327,9
0,00305 K
𝑇𝑜 − 𝑇∞ 𝑈𝐴 𝑈𝐴
= exp ( − ) ⇒ 𝑇𝑜 = 𝑇∞ + (𝑇𝑖 − 𝑇∞ ) exp (− )
𝑇𝑖 − 𝑇∞ 𝑚̇ 𝐶𝑃 𝑚̇𝐶𝑃
𝑇𝑖 + 𝑇𝑜 120 + 107,66
𝑇̅ = = = 114 ℃ = 387 K
2 2
Interpolando propiedades:
𝜌 = 832,2 kg/m3
ln( 0,011/0,0141)
387 ln( 390/380)
𝐶𝑃 = 2281 J/kg ∙ ℃ 𝜇 = 0,0141 ∙ ( ) = 0,0118 Pa ∙ s
380
𝑘 = 0,1353 W/m ∙ ℃
Calculando el número de Prandtl y el número de Reynolds
4 𝑚̇ 4 1
𝑅𝑒𝐷 = ∙ = ∙ = 2698
𝜋 𝜇𝐷𝑖 𝜋 0,0118 ∙ 0,04
𝑓 0,0472
(𝑅𝑒𝐷 − 1000) 𝑃𝑟 ∙ (2698 − 1000) ∙ 201,6
̅̅̅̅̅̅
𝑁𝑢 = 8 = 8 = 59,9
𝐷
𝑓 1/2 2/3 0,0472 1/2 2/3
1 + 12,7 ( ) (𝑃𝑟 − 1) 1 + 12,7 ∙ ( ) ∙ (201,6 − 1)
8 8
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 59,9 ∙ 0,1353 W
ℎ𝑖 = = = 202,6 2
𝐷𝑖 0,04 m ∙℃
𝑈𝐴 318,5
𝑇𝑜 = 𝑇∞ + (𝑇𝑖 − 𝑇∞ ) exp (− ) = 27 + (120 − 27) exp (− ) = 107,88 ℃
𝑚̇𝐶𝑃 1 ∙ 2281
La nueva temperatura de salida del aceite es 𝑇𝑜 = 107,88 ℃ (no varió casi nada con la iteración
anterior 𝑇𝑜 = 107,66 ℃ por lo que no es necesario hacer más iteraciones) y la temperatura de salida
se fija en:
𝑇𝑜 = 108 ℃
PROBLEMA 26
Aire
u = 8 m/s
30 cm
T = 27 ºC
8 cm
12 cm
Flujo interno:
El flujo másico es
𝜋 𝜋 kg
𝑚̇ = 𝜌𝑢𝐴 𝑡 = 𝜌𝑢 ( 𝐷𝑖2 ) = 1100 ∙ 0,5 ∙ ∙ 0,112 = 5,227
4 4 s
Se calcula el número de Reynolds:
𝑓 0,01653
(𝑅𝑒𝐷 − 1000)𝑃𝑟 ∙ (151250 − 1000) ∙ 2
̅̅̅̅̅̅
𝑁𝑢 8 8
𝐷= = = 463,6
𝑓 1/2 2/3 0,01653 1/2 2/3
1 + 12,7 ( ) ( 𝑃𝑟 − 1) 1 + 12,7 ∙ ( ) ∙ (2 − 1)
8 8
̅̅̅̅̅̅
𝑁𝑢 𝐷 ∙𝑘 463,6 ∙ 0,8 W
ℎ𝑖 = = = 3372 2
𝐷𝑖 0,11 m ∙℃
Flujo externo:
̅̅̅̅̅̅
𝑁𝑢 𝐷 ∙𝑘 285,16 ∙ 0,0263 W
ℎ𝑜 = = = 25 2
𝐷𝐶 0,3 m ∙℃
Conducción en 2 dimensiones:
2𝜋𝐿
𝑆= con
𝐷 2 + 𝑑 2 − 4𝑧 2
cosh−1 ( )
2𝐷𝑑
D
z
𝑆 2𝜋
d = = 14,57 m/m
𝐿 302 + 122 − 4 ∙ 82
cosh−1 ( )
2 ∙ 30 ∙ 12
Calor transferido:
𝑄̇ 209080 W
𝑄̇ = 𝑈𝐴∆𝑇𝑚𝑙 ⇒ 𝑈𝐴 = = = 3613,5
∆𝑇𝑚𝑙 57,86 K
Longitud de la tubería:
1 ln(𝐷𝑜 /𝐷𝑖 ) 1 1
𝐿 = 𝑈𝐴 ( + + 𝑆 + )=
ℎ 𝑖 ∙ 𝜋𝐷𝑖 2𝜋𝑘 𝑡 ( ) 𝑘 𝐶 ℎ 𝑜 ∙ 𝜋𝐷𝐶
𝐿
1 ln(12/11) 1 1
= 3613,5 ∙ ( + + + ) = 405,7 m
3372 ∙ 𝜋 ∙ 0,11 2𝜋 ∙ 40 14,57 ∙ 5 25 ∙ 𝜋 ∙ 0,3
Caída de presión:
1 𝐿 𝜇 −0,14 1 405,7
∆𝑃 = 𝑓 𝜌𝑢2 ( ) = ∙ 0,01653 ∙ ∙ 1100 ∙ 0,52 ∙ 1 = 8383 Pa = 8,4 kPa
2 𝐷𝑖 𝜇𝑆 2 0,11
Se calcula el nuevo número de Reynolds con la nueva velocidad del aire y el diámetro de aislante:
̅̅̅̅̅̅
𝑁𝑢 𝐷∙𝑘 418,2 ∙ 0,0263 W
ℎ𝑜 = = = 32,3 2
𝐷𝐴 0,34 m ∙℃
1 1 ln(12/11) 1 ln(34/30) 1 K
= ( + + + + ) = 2,7192 ∙ 10−4
405,7 3372 ∙ 𝜋 ∙ 0,11 2𝜋 ∙ 40 14,57 ∙ 5 2𝜋 ∙ 0,3 32,3 ∙ 𝜋 ∙ 0,34 W
1 W
𝑈𝐴 = −4 = 3678
2,7192 ∙ 10 K
𝑄̇ = 𝑚̇𝐶𝑃(𝑇𝑖 − 𝑇𝑜 )
{ ⇒ 𝑄̇ = 𝑚̇𝐶𝑃 (𝑇𝑖 − 𝑇𝑜 ) = 𝑈𝐴∆𝑇𝑚𝑙
𝑄̇ = 𝑈𝐴∆𝑇𝑚𝑙
(𝑇𝑖 − 𝑇∞ ) − (𝑇𝑜 − 𝑇∞ ) 1
𝑄̇ = 𝑚̇𝐶𝑃(𝑇𝑖 − 𝑇𝑜 ) = 𝑈𝐴 𝑇 − 𝑇∞ ⇒ 𝑚̇ 𝐶𝑃 = 𝑈𝐴 𝑇𝑖 − 𝑇∞
ln ( 𝑖 ) ln ( )
𝑇𝑜 − 𝑇∞ 𝑇𝑜 − 𝑇∞
𝑇𝑖 − 𝑇∞ 𝑈𝐴 𝑇𝑜 − 𝑇∞ 𝑈𝐴
ln ( )= ⇒ = exp (− )
𝑇𝑜 − 𝑇∞ 𝑚̇ 𝐶𝑃 𝑇𝑖 − 𝑇∞ 𝑚̇ 𝐶𝑃
𝑈𝐴 3678
𝑇𝑜 = 𝑇∞ + (𝑇𝑖 − 𝑇∞ ) exp (− ) = 27 + (90 − 27) exp (− ) = 79,84 ℃
𝑚̇ 𝐶𝑃 5,227 ∙ 4000
PROBLEMA 27
Considere un intercambiador de calor basado en banco de tubos. En los tubos circula vapor saturado
a 200 °C para calentar una corriente de aire que está a 20 °C. El banco de tubos está formado por 100
tubos de cobre, de diámetro interno = 2 cm, diámetro externo 2,5 cm, longitud = 1 m en arreglo
cuadrado con 5 tubos en transversal y 20 tubos en la dirección del flujo de aire, cuya distancia entre
centro y centro de los tubos es del doble del diámetro exterior. Suponga que el coeficiente de
transferencia de calor vapor-tubería es muy elevado.
(a) Calcule el caudal de aire que puede producir este equipo si se requiere que salga a 50 °C.
(b) Calcule el flujo de vapor, si se supone que sale líquido saturado.
𝑆𝐿
𝑁𝑇 = 5
𝑇∞,𝑖 = 20 ℃ 𝑇∞,𝑜 = 50 ℃
𝑇𝑖 = 200 ℃ 𝑇𝑠
𝑁𝐿 = 20
Datos adicionales:
SOLUCIÓN:
Suponiendo que el equipo está aislado, sólo hay intercambio de calor entre el vapor y el aire. Debido
a que el agua entra como vapor saturado y sale como líquido saturado su calor es igual a la entalpía
de vaporización. En este sentido, el balance de energía viene dado por:
(1) Suponer 𝑣∞ .
(2) Con 𝑣∞ se calcula el número de Reynolds y el coeficiente convectivo externo ℎ 𝑜 en el banco
de tubos.
(3) Con el coeficiente convectivo externo, se calcula el calor mediante las ecuaciones de
transferencia de calor:
Donde 𝑇𝑖 es la temperatura del vapor que permanece constante debido a que es un cambio de
fase 𝑇𝑖 = 200 ℃, ∆𝑇𝑚𝑙 es la temperatura media logarítmica y 𝑈𝐴 el coeficiente global
definido por las resistencias térmicas:
1 ln(𝐷𝑜 /𝐷𝑖 ) 1
= +
𝑈𝐴 2𝜋𝑘 𝑡 𝐿𝑁 ℎ 𝑜 𝐴 𝑜
se obtienen los valores de los flujos 𝑚̇𝐴 y 𝑚̇ 𝑉. Por otro lado, el flujo másico de aire puede
encontrarse también a partir del área transversal, la densidad y la velocidad mediante:
𝑚̇𝐴
𝑚̇𝐴 = 𝜌𝑣∞ (𝑆𝑇 𝑁𝑇𝐿) ⇒ 𝑣∞ =
𝜌(𝑆 𝑇𝑁𝑇 𝐿)
(5) Se compara la velocidad del aire calculada en (4) con la supuesta en (1).
Al principio se supone que la temperatura de la superficie de los tubos 𝑇𝑠 es muy cercana a la del
vapor ya que el coeficiente interno es muy alto (ℎ 𝑖 = ∞) y el cobre es un material muy conductor.
Calculando la temperatura de película, que es la temperatura promedio entre entrada y superficie de
los tubos, y también entre salida y superficie de los tubos, suponiendo al principio un 𝑇𝑠 ≈ 200 ℃:
Del balance de masa de aire, entre entrada e intersticios de los tubos del banco, se obtiene la velocidad
máxima que alcanza en el banco de tubos
𝑆 𝑇 /𝐷𝑜
𝑚̇𝐴 = 𝜌𝑣∞ (𝑆𝑇 𝑁𝑇𝐿) = 𝜌𝑣𝑚𝑎𝑥 (𝑆𝑇 − 𝐷𝑜 )𝑁𝑇 𝐿 ⇒ 𝑣𝑚𝑎𝑥 = 𝑣
𝑆𝑇 /𝐷𝑜 − 1 ∞
𝑆𝑇 /𝐷𝑜 2
𝑣𝑚𝑎𝑥 = 𝑣∞ = ∙ 10 = 20 m/s
𝑆𝑇 /𝐷𝑜 − 1 2 −1
𝜇𝐶𝑃
𝑃𝑟 = = 0,693
𝑘
Correlación de coeficiente convectivo (Propiedades a 𝑇𝑓 )
𝑚
̅̅̅̅𝐷 = 1,13 ∙ 𝐶1 ∙ 𝐶2 ∙ (𝑅𝑒𝑚𝑎𝑥 ) 𝑃𝑟 1/3
𝑁𝑢
Donde las constantes 𝐶1 , 𝑚 de la [Tabla 7.5 del Incropera], 𝐶2 de la [Tabla 7.6 5 del Incropera] .
Válida para 2000 < 𝑅𝑒𝑚𝑎𝑥 < 40000, 𝑃𝑟 > 0,7
𝑆𝑇 𝑆𝐿 𝐶1 = 0,229
= =2 𝑁𝐿 = 25 ≥ 10 ⇒ { 𝑚 = 0,632
𝐷𝑜 𝐷𝑜 𝐶2 = 1
𝑈𝐴 = 1290,5 W/K
𝑄̇ 212295
𝑄̇ = 𝑚̇𝐴 𝐶𝑃𝐴(𝑇∞,𝑜 − 𝑇∞,𝑖 ) ⇒ 𝑚̇𝐴 = = = 7,0273 kg/s
𝐶𝑃𝐴(𝑇∞,𝑜 − 𝑇∞,𝑖 ) 1007 ∙ 30
A partir del flujo de aire anterior, la densidad del aire a la entrada y el área transversal del equipo, se
obtiene 𝑣∞ y se compara con el valor supuesto inicialmente:
𝑚̇𝐴 7,0273
𝑣∞ = = = 23,73 m/s
𝜌(𝑆 𝑇𝑁𝑇𝐿) 1,1847 ∙ (2 ∙ 0,025) ∙ 5 ∙ 1
𝑇∞,𝑖 − 𝑇∞,𝑜 20 − 50
𝑄̇ = ℎ 𝑜 𝐴𝑜 ∙ ⇒ 212295 = 164,5 ∙ (𝜋 ∙ 0,025 ∙ 1 ∙ 100) ∙
𝑇 − 𝑇∞,𝑜 𝑇 − 50
ln ( 𝑠 ) ln ( 𝑠
𝑇𝑠 − 20
)
𝑇𝑠 − 𝑇∞,𝑖
𝑇𝑠 = 199,77 ℃ ≈ 200 ℃
∆𝑇𝑚𝑙
𝑄̇ = 𝑚̇𝐴 𝐶𝑃𝐴 (𝑇∞,𝑜 − 𝑇∞,𝑖 ) ⇒ = 𝜌𝑣∞ (𝑆𝑇𝑁𝑇 𝐿)𝐶𝑃𝐴 (𝑇∞,𝑜 − 𝑇∞,𝑖 )
1 ln(𝐷𝑜 /𝐷𝑖) 1
( + )
𝜋𝐿𝑁 2𝑘 𝑡 ℎ 𝑜 𝐷𝑜
Como todos los parámetros se mantienen constantes menos el 𝑣∞ , entonces el coeficiente exterior
es proporcional a (𝑣∞ )0,632 según la correlación para el coeficiente de convección:
𝑣∞ 0,632
ℎ 𝑜 = 157,2 ∙ ( )
10
164,5
= 𝑣∞ ∙ 1,1847 ∙ (2 ∙ 0,025) ∙ 5 ∙ 1 ∙ 1007 ∙ 30
1 ln(2,5/2) 1
( + )
𝜋 ∙ 1 ∙ 100 2 ∙ 413 𝑣 0,632
157,2 ∙ ( ∞ ) ∙ 0,025
10
𝑣∞ = 91,69 m/s
Resultados
W 𝑄̇ = 𝑈𝐴∆𝑇𝑚𝑙 = 820433 W
ℎ 𝑜 = 637,8 2
m ∙℃
Flujo másico de aire Flujo másico de vapor
𝑄̇ 𝑄̇ 820433
𝑚̇𝐴 = = 27,16 kg/s 𝑚̇𝑉 = = = 0,4205 kg/s
𝐶𝑃𝐴 (𝑇∞,𝑜 − 𝑇∞,𝑖 ) ∆𝐻𝑓𝑔 1951000
PROBLEMA 28
SOLUCIÓN:
La condición para resistencia interna despreciable (RID) es que el número de Biot es menor que 0,1
Calculando la longitud característica de la placa (con la cara de abajo adiabática, por esta superficie
no hay transferencia de calor):
𝑉 volumen 2 ∙ 1 ∙ 0,02
𝐿𝐶 ≡ = = = 0,01887 m
𝐴 𝑠 área expuesta 1 ∙ 2 + 2 ∙ (0,02 ∙ +0,02 ∙ 1)
Por lo tanto, para el siguiente rango de coeficientes convectivos el problema sería de resistencia
interna despreciable:
Los coeficientes convectivos del aire varían entre 1 y 500 W/m2 ∙ K por lo que puede suponerse que
este problema es de resistencia interna despreciable y se verifica luego de calcular ℎ.
El balance de energía considerando RID, es decir que 𝑇(𝑥, 𝑦, 𝑧, 𝑡) ≈ 𝑇(𝑡) viene dado por
𝜕𝑇
𝜌𝐶𝑃 𝑉 = −ℎ𝐴 𝑠 (𝑇 − 𝑇∞ )
𝜕𝑡
Separando variables e integrando entre la condición inicial y una arbitraria
𝑇
𝑑𝑇 𝑑𝑇 ℎ𝐴 𝑠 𝑡 𝑇 − 𝑇∞ ℎ𝐴 𝑠 𝑡
𝜌𝐶𝑃𝑉 = −ℎ𝐴 𝑠 (𝑇 − 𝑇∞ ) ⇒ ∫ =− ∫ 𝑑𝑡 ⇒ ln ( )=−
𝑑𝑡 𝑇𝑖 𝑇 − 𝑇∞ 𝜌𝐶𝑃 𝑉 0 𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃𝑉
Despejando el tiempo
𝜌𝐶𝑃 𝑉 𝑇 − 𝑇∞ 𝜌𝐶 𝐿 𝑇 − 𝑇∞
𝑡 =− ln ( ) = − 𝑃 𝐶 ln ( )
ℎ𝐴 𝑠 𝑇𝑖 − 𝑇∞ ℎ 𝑇𝑖 − 𝑇∞
̅̅̅̅̅̅ 1/4
𝑁𝑢𝐿 = 0,27𝑅𝑎 𝐿
𝜌 = 0,9628 kg/m3
𝐶𝑃 = 1010 J/kg ∙ ℃
𝑘 = 0,03099 W/m ∙ ℃
aire a 𝑇̅𝑓 = 363 K
𝜇 = 213,9 ∙ 10 −7 Pa ∙ s
1 1
𝛽= ̅ = = 0,002755 K −1
{ 𝑇𝑓 363
El número de Nusselt es
̅̅̅̅̅̅ 1/4
𝑁𝑢𝐿 = 0,27𝑅𝑎𝐿 = 0,27 ∙ (1,695 ∙ 108 ) 1/4 = 30,81
̅̅̅̅̅̅
𝑁𝑢𝐿 ∙ 𝑘 30,81 ∙ 0,03099 W
ℎ= = = 2,9 2
𝐿 1/3 m ∙K
Vapor saturado
120 ºC
Datos adicionales: la entalpía de vaporización del agua a 120 ºC es ℎ𝑓𝑔 = 2202 kJ/kg y el calor
específico del agua líquida fría es 𝐶𝑃𝑐 = 4180 J/kgºC en el condensador y 𝐶𝑃𝑐 = 4188 J/kgºC en el
intercambiador. El calor específico del agua líquida caliente es de 𝐶𝑃ℎ = 4203 J/kgºC en el
intercambiador.
SOLUCIÓN:
Cálculos en el condensador:
Como es un condensador, entonces la temperatura de la corriente de salida del lado caliente es 120
ºC por ser líquido saturado (el cambio de fase ocurre a temperatura constante), así:
W
𝐶𝑚𝑖𝑛 = 𝐶𝑐 = 𝑚̇ 𝑐𝐶𝑃𝑐 = 120 ∙ 4180 = 501600
K
W W
𝑈𝐴 = 3250 2 ∙ 80 m2 = 260000
m ∙K K
El número de unidades de transferencia NTU es
𝑈𝐴 260000
NTU = = = 0,51834
𝐶𝑚𝑖𝑛 501600
Luego el flujo másico de vapor 𝑚̇ℎ se calcula por el balance de energía del vapor:
𝑄̇ 19883,6 kJ/s kg
𝑄̇ = 𝑚̇ ℎ ℎ𝑓𝑔 ⇒ 𝑚̇ℎ = = = 9,03
ℎ𝑓𝑔 2202 kJ/kg s
W W
𝐶ℎ = 𝑚̇ ℎ 𝐶𝑃ℎ = 9,03 ∙ 4203 = 37953 𝐶𝑚𝑖𝑛 = min {𝐶ℎ ,𝐶𝑐 } = 37953
K K
W W
𝐶𝑐 = 𝑚̇ 𝑐 𝐶𝑃𝑐 = 120 ∙ 4188 = 1753944 𝐶𝑚𝑎𝑥 = max {𝐶ℎ ,𝐶𝑐} = 17539344
K K
𝐶𝑚𝑖𝑛 37953
𝐶𝑟 = = = 0,02164
𝐶𝑚𝑎𝑥 17539344
W W
𝑈𝐴 = 1800 2 ∙ 60 m2 = 108000
m ∙K K
𝑈𝐴 108000
NTU = = = 2,8456
𝐶𝑚𝑖𝑛 0,02164
1 − exp(−NTU(1 − 𝐶𝑟 ))
𝜀= = 0,94166
1 − 𝐶𝑟 exp(−NTU(1 − 𝐶𝑟))
𝑄̇ 2085727
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 120 − = 65,04 ℃
𝐶ℎ 37953
𝑄̇ 2085727
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 61,64 + = 61,76 ℃
𝐶𝑐 1753944
PROBLEMA 30
En el esquema adjunto se presentan dos maneras de operar una batería de tres intercambiadores en
serie. Cada uno de ellos tiene 100 m2 de área externa y están construidos con tubos de 1,25 cm de
diámetro externo y 1 cm de diámetro interno. La conductividad del material es de 250 W/m-ºC. El
calor específico del fluido de enfriamiento es 𝐶𝑃𝑐 = 3650 J/kg ∙ K y el del fluido caliente
𝐶𝑃ℎ = 3250 J/kg ∙ K. El coeficiente interno 2250 W/m2 ∙ ℃ y el externo 1950 W/m2 ∙ ℃.
Complete los cuadros adjuntos y determine cuál de las dos opciones es mejor (% de velocidad
adicional de transferencia de calor).
Operación en paralelo:
W W
𝐶ℎ = 𝑚̇ ℎ 𝐶𝑃ℎ = 10 ∙ 3250 = 32500 𝐶𝑚𝑖𝑛 = min {𝐶ℎ ,𝐶𝑐 } = 32500
K K
W W
𝐶𝑐 = 𝑚̇ 𝑐 𝐶𝑃𝑐 = 30 ∙ 3650 = 109500 𝐶𝑚𝑎𝑥 = max {𝐶ℎ ,𝐶𝑐 } = 109500
K K
𝐶𝑚𝑖𝑛 32500
𝐶𝑟 = = = 0,2968
𝐶𝑚𝑎𝑥 109500
𝐴 𝑜 = 𝜋𝐷𝑜 𝐿𝑁 𝐴 𝑜 𝐷𝑜 𝐷𝑖 1
} ⇒ = ⇒ 𝐴𝑖 = 𝐴𝑜 = ∙ 100 = 80 m2
𝐴 𝑖 = 𝜋𝐷𝑖 𝐿𝑁 𝐴 𝑖 𝐷𝑖 𝐷𝑜 1,25
𝐴𝑜 100
𝐴 𝑜 = 𝜋𝐷𝑜 𝐿𝑁 ⇒ 𝐿𝑁 = = = 2546,5 m
𝜋𝐷𝑜 𝜋 ∙ 0,0125
1 1 ln(𝐷𝑜 /𝐷𝑖 ) 1 1 ln(1,25/1) 1
= + + = + +
𝑈𝐴 ℎ 𝑖𝐴 𝑖 2𝜋𝑘 𝑡 𝐿𝑁 ℎ 𝑜 𝐴 𝑜 2250 ∙ 80 2𝜋 ∙ 250 ∙ 2456,5 1950 ∙ 100
W
𝑈𝐴 = 93114
K
𝑈𝐴 93114
NTU = = = 2,865
𝐶𝑚𝑖𝑛 32500
Método NTU-eficiencia:
NTU
𝜀= si 𝐶𝑟 = 1
1 + NTU
𝑄̇ 1515991
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 80 − = 33,35 ℃
𝐶ℎ 32500
𝑄̇ 1515991
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 18 + = 31,84 ℃
𝐶𝑐 109500
𝑄̇ 1177583
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 80 − = 43,77 ℃
𝐶ℎ 32500
𝑄̇ 1177583
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 31,84 + = 42,59 ℃
𝐶𝑐 109500
𝑄̇ 914730
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 80 − = 51,85 ℃
𝐶ℎ 32500
𝑄̇ 914730
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 42,59 + = 50,94 ℃
𝐶𝑐 109500
𝑄̇ 1818245
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 80 − = 24,05 ℃
𝐶ℎ 32500
𝑄̇ 1818245
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 18 + = 34,60 ℃
𝐶𝑐 109500
𝑄̇ 1331278
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 80 − = 39,04 ℃
𝐶ℎ 32500
𝑄̇ 1331278
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 34,60 + = 46,76 ℃
𝐶𝑐 109500
Cálculos I-3: 𝑇𝑐𝑖 = 46,76 ℃
𝑄̇ 974733
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 80 − = 50,01 ℃
𝐶ℎ 32500
𝑄̇ 974733
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 46,76 + = 55,66 ℃
𝐶𝑐 109500
Resumen de valores
Tho1 Tho2
Tho1 = ?
Tho = ?
SOLUCIÓN:
𝐶𝑚𝑖𝑛 32000
𝐶𝑟 = = = 0,5333
𝐶𝑚𝑎𝑥 60000
𝑄 2160000
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 200 − = 132,5 ℉
𝑚̇ℎ 𝐶𝑃ℎ 32000
𝑄 2160000
𝜀= = = 0,5114
𝑄𝑚𝑎𝑥 4224000
1 𝜀−1 1 0,5114 − 1
NTU = ln ( )= ∙ ln ( ) = 0,8522
𝐶𝑟 − 1 𝜀𝐶𝑟 − 1 0,5333 − 1 0,5114 ∙ 0,5333 − 1
𝑈𝐴
NTU = ⇒ 𝑈𝐴 = NTU ∙ 𝐶𝑚𝑖𝑛 = 0,8522 ∙ 32000 = 27270 BTU/h ∙ ℉
𝐶𝑚𝑖𝑛
𝑄
𝜀= ⇒ 𝑄 = 𝜀𝑄𝑚𝑎𝑥 = 0,5114 ∙ 3072000 = 1570909 BTU/h
𝑄𝑚𝑎𝑥
𝑄 1570909
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 200 − = 150,9 ℉
𝑚̇ℎ 𝐶𝑃ℎ 32000
𝑄 1570909
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 104 + = 130,2 ℉
𝑚̇𝑐 𝐶𝑃𝑐 60000
𝑚̇ℎ3 = 𝑚̇ ℎ1 + 𝑚̇ℎ2
{ ⇒ (𝑚̇ ℎ1 + 𝑚̇ℎ2 )𝐻3 = 𝑚̇ ℎ1 𝐻1 + 𝑚̇ ℎ2 𝐻2
𝑚̇ℎ3 𝐻3 = 𝑚̇ ℎ1 𝐻1 + 𝑚̇ ℎ2𝐻2
𝐶𝑚𝑖𝑛 60000
𝐶𝑟 = = = 0,9375
𝐶𝑚𝑎𝑥 64000
𝑈𝐴 27270
NTU = = = 0,4545
𝐶𝑚𝑖𝑛 60000
Además 𝑇𝑐𝑖2 = 𝑇𝑐𝑜1 y 𝑇ℎ𝑜2 = 𝑇ℎ𝑖1, por lo que se resuelven cuatro ecuaciones con cuatro icógnitas
Sustituyendo valores
Resolviendo el sistema de ecuaciones lineal se obtienen 𝑇ℎ𝑖1 = 𝑇ℎ𝑜2 = 169,9 ℉, 𝑇ℎ𝑜1 = 141,6 ℉,
𝑇𝑐𝑖2 = 𝑇𝑐𝑜1 = 98,1 ℉, 𝑇𝑐𝑜2 = 130,3 ℉. Calculando los calores transferidos
En la figura se muestra un esquema para condensar y subenfriar 7200 kg/h de un vapor saturado que
condensa en el condensador a 80 °C y luego se subenfría en un intercambiador en contracorriente
(IC) hasta 25 °C. El calor de condensación del vapor es de 525 kJ/kg. El calor específico del
condensado es de 3250 J/kg ∙ ℃. El agua para enfriar y condensar está disponible a 18 °C y no debe
salir del condensador a más de 52 °C. El agua para condensar debe pasar a 1,2 m/s a través de los
tubos de 1,6 cm de diámetro interno y 2 mm de espesor. El coeficiente de transferencia de calor del
vapor que condensa es de 5250 W/m2 ∙ ℃. La conductividad de la tubería es de 50 W/m ∙ K.
Determine:
Vapor saturado
7200 kg/h
80 °C
Agua
Condensador
52 °C
CP = 3250 J/kg-°C
IC
Condensado
25 °C
Agua
18 °C
El calor específico promedio del agua entre 18 °C y 52 °C es 𝐶𝑃𝐶 = 4181 J/kg ∙ ℃ (tabla de
propiedades del agua líquida). En el condensador, el vapor sale como líquido saturado a la misma
temperatura del vapor 80 °C (ya que en un cambio de fase no cambia la temperatura). Haciendo un
balance de energía en el sistema completo suponiendo que ambos intercambiadores son adiabáticos
(se toma el sistema como una caja negra, donde sólo se consideran las entradas y salidas del mismo
y omitiendo las corrientes intermedias):
Vapor saturado
7200 kg/h
80 °C
Agua
Condensador
52 °C
CP h = 3250 J/kg-°C
IC
Líquido saturado Condensado
7200 kg/h 25 °C
80 °C
Agua
18 °C
CP c = 4181 J/kg-°C
kg kg
𝑚̇ ℎ = 7200 =2
h s
Agua
Tco = ?
Líquido saturado
mh = 2 kg/s
Thi = 80 °C CP = 3250 J/kg-°C
IC
Condensado
Tho = 25 °C
Agua
Tci = 18 °C
mc = 9,9012 kg/s
CP c = 4181 J/kg-°C
𝑚̇ ℎ 𝐶𝑃ℎ 2 ∙ 3250
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + (𝑇ℎ𝑖 − 𝑇ℎ𝑜 ) = 18 + (80 − 25) = 26,64 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 9,9012 ∙ 4181
𝜌 = 993 kg/m3
𝐶𝑃 = 4178 J/kg ∙ ℃
agua a 310 K 𝑘 = 0,628 W/m ∙ ℃
𝜇 = 0,000695 Pa ∙ s
{
ℎ 𝑖 𝐷𝑖 ̅̅̅̅̅̅
𝑁𝑢 𝐷 ∙𝑘 150,8 ∙ 0,628 W
̅̅̅̅̅̅
𝑁𝑢 𝐷= ⇒ ℎ𝑖 = = = 5918 2
𝑘 𝐷𝑖 0,016 m ∙K
Vapor saturado
7200 kg/h
Thi = 80 °C
Agua Agua
Condensador
Tco = 52 °C Tci = 26,64 °C
Líquido saturado
Tho = 80 °C
𝑄 1050000 W
𝑄 = 𝑈𝐴∆𝑇𝑚𝑙 ⇒ 𝑈𝐴 = = = 26700
∆𝑇𝑚𝑙 39,33 K
𝑈𝐴 1 ln(𝐷𝑜 /𝐷𝑖 ) 1
𝐿𝑁 = ( + + )
𝜋 ℎ 𝑖 𝐷𝑖 2𝑘 𝑡 ℎ 𝑜 𝐷𝑜
26700 1 ln(2/1,6) 1
𝐿𝑁 = ( + + ) = 189,7 m
𝜋 5918 ∙ 0,016 2 ∙ 50 5250 ∙ 0,02
El área externa es
𝐴 𝑜 = 𝜋𝐷𝑜 𝐿𝑁 = 𝜋 ∙ 0,02 ∙ 151,7 = 11,92 m2
𝑚̇𝑐 9,9012
𝑁= = = 41,32 ≈ 42 tubos
𝑚̇𝑡 0,2396
189,7 189,7
𝐿𝑁 = 189,7 m ⇒ 𝐿= = = 4,52 m/tubo
𝑁 42
PROBLEMA 33
Un aceite (de calor específico promedio de 2500 J/kg ∙ ℃), se enfría en un intercambiador de calor
tubos y coraza 1-2 (un paso por la coraza y dos pasos por los tubos). El aceite entra al intercambiador
(circula por la coraza) con un flujo másico de 10 kg/s a una temperatura de 180 ºC y sale a una
temperatura de 90,8 ºC. Se utiliza agua como fluido de enfriamiento (calor específico promedio de
4200 J/kg ∙ ℃) el cual entra al intercambiador a 30 ºC con un flujo másico de 10 kg/s (entra por los
tubos).
(a) Ahora se requiere enfriar más el aceite, para esto se decide cambiar la temperatura de entrada
del agua de enfriamiento a 5 ºC. Utilizando sólo estos datos, determine la nueva temperatura
de salida del aceite.
(b) ¿Cuántos intercambiadores de tubos y coraza arreglados en serie se necesitan para enfriar el
aceite hasta una temperatura menor o igual a 35 ºC si se utiliza agua de enfriamiento a 5 ºC?
SOLUCIÓN:
𝐶𝑚𝑖𝑛 25000
𝐶𝑟 = = = 0,59524
𝐶𝑚𝑎𝑥 42000
𝑄 2230000
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 30 + = 83,10 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 42000
𝑄 2230000
𝜀= = = 0,5947
𝑄𝑚𝑎𝑥 3750000
2 1,519 − 1
− (1 + 0,59524) ln ( )
0,5947 1,519 + 1
𝐸= = 1,519 ⇒ NTU = − = 1,3571
(1 + 0,595242 ) 1/2 (1 + 0,595242 )1/2
𝑈𝐴
NTU = ⇒ 𝑈𝐴 = NTU ∙ 𝐶𝑚𝑖𝑛 = 1,3571 ∙ 25000 = 33929 W/K
𝐶𝑚𝑖𝑛
(b) Como los flujos no cambian y los calores específicos son funciones débiles de la temperatura,
entonces el valor de 𝑈𝐴 y 𝐶𝑟 son constantes y por lo tanto NTU y 𝜀 son los mismos que los calculados
anteriormente. Calculando de nuevo el calor máximo y con el balance de energía obteniendo las
temperaturas de salidas se obtiene:
𝑄
𝜀= ⇒ 𝑄 = 𝜀𝑄𝑚𝑎𝑥 = 0,5947 ∙ 4375000 = 2601667 W
𝑄𝑚𝑎𝑥
𝑄 2601667
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 180 − = 75,93 ℃
𝑚̇ℎ 𝐶𝑃ℎ 25000
𝑄 2601667
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 5+ = 66,94 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 42000
(c) Añadiendo otro intercambiador tubos y coraza idéntico y con un servicio de agua de enfriamiento
de 10 kg/s disponible a 5 ºC, calculando de nuevo la temperatura de salida del aceite:
𝑄
𝜀= ⇒ 𝑄 = 𝜀𝑄𝑚𝑎𝑥 = 0,5947 ∙ 1773333 = 1054542 W
𝑄𝑚𝑎𝑥
𝑄 1054542
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 180 − = 33,75 ℃ < 35 ℃
𝑚̇ℎ 𝐶𝑃ℎ 25000
𝑄 1054542
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 5+ = 30,11 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 42000
Entonces con un arreglo de dos intercambiadores en serie del lado del aceite es posible bajar la
temperatura de 180 ºC hasta una temperatura menor o igual a 35 ºC.
PROBLEMA 34
Se desea calentar una corriente de agua de proceso desde una temperatura de 34 ºC hasta una
temperatura de 45 ºC con 5 kg/s de una corriente de agua cruda disponible a 147 ºC cuya temperatura
de salida es de 50 ºC utilizando un intercambiador de calor de tubos y coraza con configuración
cocorriente-contracorriente con 4 pasos por tubo y un paso por la coraza.
Tubos
Catálogo BWG 18 – 0,75" Número de tubos 176
tubos/coraza
Diámetro externo [in] 0,75 Número de pasos por tubo 4 pasos/tubo
Diámetro interno [in] 0,652 Triangular
Conductividad térmica 51,9 Arreglo de los tubos
[W/m ∙ ℃]
Espaciamiento centro-centro 1"
Coraza
Diámetro de la coraza 17,25" Número de pasos de corazas 1 paso/coraza
Deflectores
Espaciamiento entre 17,25" Corte de deflectores (%) 25 %
deflectores
(a) Determine la longitud requerida para los tubos del intercambiador de calor y calcule la caída
de presión en el lado de la coraza y el lado de los tubos y exprésela en [bar].
(b) Hallar las temperaturas de salida si el flujo de agua de proceso se reduce un 40% y el flujo
de agua cruda se reduce un 50%, y se mantienen los mismos factores de obstrucción de
diseño.
SOLUCIÓN:
Primero se calcula el valor de 𝑈𝐴 requerido para la transferencia de calor. Para esto se determina la
diferencia de temperaturas media logarítmica y el factor de corrección.
T1 = 34 ºC
t2 = 50 ºC
t1 = 147 ºC
T2 = 45 ºC
𝑇1 − 𝑇2 34 − 45 𝑡2 − 𝑡1 50 − 147
𝑅= = = 0,1134 𝑆= = = 0,8584
𝑡2 − 𝑡1 50 − 147 𝑇1 − 𝑡1 34 − 147
El factor de configuración térmica (configuración cocorriente-contracorriente 1-2) es:
1−𝑆
√1 + 𝑅2 ln ( )
𝐹= 1 − 𝑅𝑆 =
2 − 𝑆(𝑅 + 1 − √1 + 𝑅2 )
(𝑅 − 1) ln ( )
2 − 𝑆(𝑅 + 1 + √1 + 𝑅2 )
1 − 0,8584
√1 + 0,11342 ∙ ln ( )
1 − 0,1134 ∙ 0,8584
= = 0,8913
2 − 0,8584 (0,1134 + 1 − √1 + 0,11342 )
(0,1134 − 1) ∙ ln ( )
2 − 0,8584 (0,1134 + 1 + √1 + 0,11342 )
La diferencia de temperaturas media logarítmica es
A partir del balance de energía se busca el calor transferido y el flujo másico de la corriente fría:
𝑄 2061250
𝑚̇𝑐 = = = 44,7543 kg/s
(
𝐶𝑃𝑐 𝑇2 − 𝑇1 ) 4187 ∙ (45 − 34)
𝑄 2061250 W
𝑄 = 𝑈𝐴𝐹∆𝑇𝑚𝑙 ⇒ 𝑈𝐴 = = = 49811
𝐹∆𝑇𝑚𝑙 0,8913 ∙ 46,43 K
Cálculos del lado de los tubos:
El flujo másico por los tubos es 𝑚̇ 𝑡 = 𝑚̇ℎ = 5 kg/s. El diámetro interior de los tubos es 0,652" =
0,01656 m. Por lo tanto el área de flujo en el lado de los tubos es
𝜋 2 𝑁 𝜋 176
𝐴𝑓𝑡 = 𝐷𝑖 = ∙ 0,016562 ∙ = 0,009477 m2
4 𝑁𝑃𝑡 4 4
𝜌 𝑚̇ 𝑡 𝑚̇ 𝐷 5 ∙ 0,01656
𝑅𝑒𝑡 = ( )𝐷𝑖 = 𝑡 𝑖 = = 22992
𝜇 𝜌𝐴𝑓𝑡 𝜇𝐴𝑓𝑡 0,38 ∙ 10−3 ∙ 0,009477
0,4 si es calentamiento
𝑁𝑢𝑡 = 0,023𝑅𝑒𝑡0,8 𝑃𝑟 𝑛
̅̅̅̅̅ con 𝑛 = {
0,3 si es enfriamiento
̅̅̅̅̅𝑡 = 0,023 ∙ (22992) 0,8 ∙ (2,45) 0,3 = 92,84
𝑁𝑢
ℎ 𝑖 𝐷𝑖 ̅̅̅̅̅
𝑁𝑢𝑡 ∙ 𝑘 92,84 ∙ 0,66 W
̅̅̅̅̅𝑡 =
𝑁𝑢 ⇒ ℎ𝑖 = = = 3700 2
𝑘 𝐷𝑖 0,01656 m ∙K
P Do
𝑃 − 𝐷𝑜 𝐷𝑜
B 𝐴𝑓𝑠 = 𝐷𝑠 𝐵 ∙ = 𝐷𝑠 𝐵 (1 − )
𝑃 𝑃
Ds
𝐷𝑜 3
𝐴𝑓𝑠 = 𝐷𝑠 𝐵 (1 − ) = 0,43815 ∙ 0,43815 ∙ (1 − ) = 0,048 m2
𝑃 4
𝑚̇ 𝑠 𝐷𝐻 44,7543 ∙ 0,01829
𝑅𝑒𝑠 = = = 25457
𝜇𝐴𝑓𝑠 0,67 ∙ 10−3 ∙ 0,048
El número de Nusselt se encuentra a partir de la gráfica del coeficiente de transferencia de calor del
lado de la coraza con deflectores segmentados 25% (Kern, 1950):
𝜇 0,14
̅̅̅̅̅𝑠 = 0,36𝑅𝑒𝑠0,55 𝑃𝑟 1/3 (
𝑁𝑢 ) = 0,36 ∙ (25457) 0,55 ∙ (4,38) 1/3 ∙ (1) 0,14 = 156,1
𝜇𝑤
𝑈𝐴 1 𝑅𝑖 ln(𝐷𝑜 /𝐷𝑖 ) 𝑅𝑜 1
𝐿= ( + + + + )=
𝜋𝑁 ℎ 𝑖𝐷𝑖 𝐷𝑖 2𝑘 𝑡 𝐷𝑜 ℎ 𝑜 𝐷𝑜
Y las pérdidas de presión de regreso se consideran como 4 cabezales de velocidad por paso por los
tubos (Kern, 1950), así
2 2
𝜌𝑢2 1 𝑚̇𝑡 1 5
∆𝑃𝑐 = 4𝑁𝑃𝑡 = 2𝑁𝑃𝑡 ( ) = 2 ∙ 4 ∙ ( ) = 2334 Pa
2 𝜌 𝐴𝑓𝑡 954 0,009477
𝐿 6,62
𝑁𝐵 + 1 = = = 15,1~15 espacios entre deflectores
𝐵 0,43815
Los nuevos flujos son 𝑚̇′𝑡 = 0,5𝑚̇𝑡 y 𝑚̇ ′𝑠 = 0,6𝑚̇ 𝑠 . Como varían los flujos, hay que encontrar los
nuevos valores de los coeficientes convectivos. Sabiendo que en el lado de los tubos, las propiedades
y otros parámetros se mantienen constantes, entonces el coeficiente convectivo es proporcional a 𝑚̇ 0,8
𝑡
por lo que se emplea esta regla de escalamiento
0,8 0,8
ℎ ′𝑖 𝑚̇′ 𝑚̇′ W
= ( 𝑡) ⇒ ℎ ′𝑖 = ℎ 𝑖 ( 𝑡) = 3700 ∙ (0,5) 0,8 = 2125
ℎ𝑖 𝑚̇𝑡 𝑚̇𝑡 m2 ∙ K
𝑈𝐴 = 41399 W/K
Como el intercambiador está diseñado y se necesita encontrar las temperaturas de salida, entonces se
aplica el método de NTU-eficiencia para encontrar las temperaturas de salida.
𝐶𝑚𝑖𝑛 10625
𝐶𝑟 = = = 0,0945
𝐶𝑚𝑎𝑥 112432
𝑈𝐴 41399
NTU = = = 3,896
𝐶𝑚𝑖𝑛 10625
2
𝐸−1 − (1 + 𝐶𝑟 )
NTU = −(1 + 𝐶𝑟2 ) −1/2 ln ( ) donde 𝐸 = 𝜀
𝐸+1 (1 + 𝐶𝑟2 )1/2
Y se resuelve para 𝜀
2
𝐸−1 3,896 − (1 + 0,0945)
= exp (− ) ⇒ 𝐸 = 1,0422 ⇒ 𝐸 = 𝜀 ⇒ 𝜀 = 0,9346
𝐸+1 √1 + 0,09452 √1 + 0,09452
𝑄
𝜀= ⇒ 𝑄 = 𝜀𝑄𝑚𝑎𝑥 = 0,9346 ∙ 1200625 = 1122141 W
𝑄𝑚𝑎𝑥
𝑄 1122141
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 147 − = 41,39 ℃
𝑚̇ℎ 𝐶𝑃ℎ 10625
𝑄 1122141
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 34 + = 43,98 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 112432
PROBLEMA 35
Un intercambiador de calor de flujo cruzado está operando actualmente enfriando un flujo másico de
aceite de 8 kg/s con una temperatura de 120 ºC usando un flujo másico de 27 kg/s de aire atmosférico
disponible a 25 ºC. Las características del intercambiador de calor y las propiedades del aceite y del
aire (supóngalas constantes) se muestran en las siguientes tablas:
(a) Calcule el calor transferido y la temperatura de salida del aceite y del aire. Puede considerarse
que el factor de obstrucción en el exterior de los tubos es cero y el factor de obstrucción
generado por el aceite es de 0,0001 m2 ∙ K/W.
Por un cambio en el proceso, se requiere que el aceite se siga enfriando hasta 50 ºC utilizando una
corriente de 6 kg/s de agua de enfriamiento disponible a 20 ºC. Para esto, se debe instalar un
intercambiador de doble tubo 3 × 2 construido con tubos catálogo 40 de 6 m por tramo. El tubo
externo tiene un diámetro nominal de 3" (un diámetro interno de 3,068 in) y tubo interno con diámetro
nominal de 2" (un diámetro interno de 2,067 in y un diámetro externo de 2,375 in). El aceite pasa por
los tubos internos y el agua de enfriamiento por la región anular y una configuración contracorriente.
Los factores de obstrucción en ambos lados son iguales a cero.
(b) Determine el número de tramos de doble tubo necesario para el nuevo intercambiador de
calor y calcule la temperatura de salida del agua de enfriamiento. Puede considerar las
propiedades de los fluidos constantes.
Agua de enfriamiento
L = 12 m 6 kg/s
20 ºC
Intercambiador
de doble tubo
Aceite Aceite
8 kg/s 50 ºC
120 ºC
Aire
27 kg/s
25 ºC
SOLUCIÓN:
Coeficiente convectivo del lado de los tubos del intercambiador de flujo cruzado:
El flujo másico por los tubos es 𝑚̇ 𝑡 = 𝑚̇ ℎ = 8 kg/s . El área de flujo del lado de los tubos es
𝜋 𝑁 𝜋 144
𝐴𝑓𝑡 = 𝐷𝑖2 = ∙ 0,022 ∙ = 0,04524 m2
4 𝑁𝑃𝑡 4 1
𝜌 𝑚̇𝑡 𝑚̇ 𝐷 8 ∙ 0,02
𝑅𝑒𝑡 = ( ) 𝐷𝑖 = 𝑡 𝑖 = = 393 (flujo laminar)
𝜇 𝜌𝐴𝑓𝑡 𝜇𝐴𝑓𝑡 0,009 ∙ 0,04524
𝐷𝑖 0,02
( ) 𝑅𝑒𝑡 𝑃𝑟 = ( ) ∙ 393 ∙ 128,6 = 84,23 > 10
𝐿 12
0,0668(𝐷𝑖 /𝐿)𝑅𝑒𝑡 𝑃𝑟 0,0668 ∙ 84,23
̅̅̅̅̅
𝑁𝑢𝑡 = 3,66 + = 3,66 + = 6,841
2/3 1 + 0,04(84,23) 2/3
1 + 0,04((𝐷𝑖 /𝐿)𝑅𝑒𝑡𝑃𝑟)
El número de Reynolds en el banco de tubos está basado en la velocidad máxima que tiene el fluido
dentro del mismo, esto es cuando el área es la mínima. Mediante un balance de masa del fluido se
concluye que
𝑆𝑇 𝑆 𝑇 /𝐷𝑜
𝑚̇ = 𝜌𝑣∞ (𝑆 𝑇𝑁𝑇 𝐿) = 𝜌𝑣𝑚𝑎𝑥 (𝑆𝑇 − 𝐷𝑜 )𝑁𝑇𝐿 ⇒ 𝑣𝑚𝑎𝑥 = 𝑣∞ = 𝑣
𝑆𝑇 − 𝐷𝑜 𝑆𝑇 /𝐷𝑜 − 1 ∞
𝑚̇𝑜 27
𝑣∞ = = = 4,513 m/s
𝜌 (𝑆𝑇 𝑁𝑇𝐿) 1,108 ∙ (1,25 ∙ 0,03) ∙ 12 ∙ 12
1,25
𝑣𝑚𝑎𝑥 = ∙ 4,513 = 22,56 m/s
1,25 − 1
Para un arreglo cuadrado con 𝑆 𝑇/𝐷𝑜 = 𝑆𝐿 /𝐷𝑜 = 1,25 y con más de 10 hileras de tubos a lo largo
(𝑁𝐿 ≥ 10) se tiene que las constantes de la correlación de banco de tubos son 𝐶 = 0,386, 𝑛 = 0,592
válida para 2000 < 𝑅𝑒𝑚𝑎𝑥 < 40000, 𝑃𝑟 > 0,7
̅̅̅̅𝑜 = 1,13𝐶𝑅𝑒𝑚𝑎𝑥
𝑁𝑢 𝑛 𝑃𝑟 1/3 = 1,13 ∙ 0,348 ∙ (39022) 0,592 ∙ (0,705) 1/3 = 182,9
1 1 1 0,0001 ln(3/2) 1 K
= ∙( + + + ) = 2,3056 ∙ 10−4
𝑈𝐴 𝜋 ∙ 12 ∙ 144 47,9 ∙ 0,02 0,02 2 ∙ 51,9 167,6 ∙ 0,03 W
𝑈𝐴 = 4337 W/K
𝐶𝑚𝑖𝑛 16000
𝐶𝑟 = = = 0,5879
𝐶𝑚𝑎𝑥 27216
𝑈𝐴 4337
NTU = = = 0,2711
𝐶𝑚𝑖𝑛 16000
Como el fluido con 𝐶𝑚𝑖𝑛 no se mezcla y el fluido con 𝐶𝑚𝑎𝑥 se mezcla, entonces
𝑄
𝜀= ⇒ 𝑄 = 𝜀𝑄𝑚𝑎𝑥 = 0,2216 ∙ 1520000 = 336832 W
𝑄𝑚𝑎𝑥
𝑄 336832
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 25 + = 37,38 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 27216
𝑄 336832
𝑇ℎ𝑜 = 𝑇ℎ𝑖 − = 120 − = 98,95 ℃
𝑚̇ ℎ 𝐶𝑃ℎ 16000
Haciendo un balance de energía en este intercambiador se obtiene la temperatura de salida del agua
de enfriamiento
𝑄 783200
𝑇𝑐𝑜 = 𝑇𝑐𝑖 + = 20 + = 51,23 ℃
𝑚̇ 𝑐 𝐶𝑃𝑐 6 ∙ 4179
𝑄 783200
𝑄 = 𝑈𝐴∆𝑇𝑚𝑙 ⇒ 𝑈𝐴 = = = 20515 W/K
∆𝑇𝑚𝑙 38,18
Coeficiente convectivo del lado del tubo interno del intercambiador de doble tubo:
𝑚̇𝑡 𝐷𝑖 4 𝑚̇𝑡 4 8
𝑅𝑒𝑖 = 𝜋 2 = 𝜋 ∙ 𝜇𝐷 = 𝜋 ∙ 0,009 ∙ 0,0525 = 21557 (turbulento)
𝜇 ( 𝐷𝑖 ) 𝑖
4
Coeficiente convectivo del lado de la región anular del intercambiador de doble tubo:
𝑚̇ 𝑜 𝐷𝐻 4 𝑚̇𝑜 4 6
𝑅𝑒𝑜 = 𝜋 = ∙ = ∙ = 64625
𝜇 (𝐷𝑜2 − 𝐷𝑖 ) 𝜋 𝜇 (𝐷𝑜 + 𝐷𝑖 ) 𝜋 0,000855 ∙ (0,06033 + 0,07793)
2
4
Y el número de Prandtl del fluido del lado de la región anular
̅̅̅̅̅̅
𝑁𝑢𝑜 = 0,023 ∙ (64625) 0,8 ∙ (5,829) 0,4 = 328,3
ℎ 𝑜 𝐷𝐻 ̅̅̅̅̅̅
𝑁𝑢𝑜 ∙ 𝑘 328,3 ∙ 0,613 W
̅̅̅̅̅̅
𝑁𝑢𝑜 = ⇒ ℎ𝑜 = = = 11435 2
𝑘 𝐷𝐻 0,07793 − 0,06033 m ∙K
Con esto se calcula el número de tramos de doble tubo necesarios para satisfacer el requerimiento:
𝑈𝐴 1 𝑅𝑖 ln(𝐷𝑜 /𝐷𝑖) 𝑅𝑜 1
𝑁= ( + + + + )=
𝜋𝐿 ℎ 𝑖 𝐷𝑖 𝐷𝑖 2𝑘 𝑡 𝐷𝑜 ℎ 𝑜 𝐷𝑜
20515 1 ln(0,06033/0,0525) 1
= ( +0 + +0 + ) = 29,9 ≈ 30 tramos
𝜋 ∙ 6 771,5 ∙ 0,0525 2 ∙ 51,9 11435 ∙ 0,06033
PROBLEMA 36
SOLUCIÓN:
Así la emisividad es
𝐹(0→1,5) = 0,2732
Sugerencia: en la calculadora HPP49G / 50G se edita F0𝜆(𝜆, 𝑇) con 𝜆 en [m] y 𝑇 en grados [K]
PROBLEMA 37
A partir de los balances de energía sobre cada una de las paredes de un recinto constituido por dos
superficies obtenga la expresión que permita calcular la velocidad de transferencia de calor neta entre
ambas superficies en función de las áreas, emisividades y temperaturas de cada superficie.
SOLUCIÓN:
Suposiciones: Notación:
𝑄̇ : Transferencia neta de calor
Transferencia de calor en estado estacionario 𝑄̇12 : Transferencia efectiva que emite la
Intercambio de radiación en el vacío superficie 1 y recibe la superficie 2
Las dos superficies forman un recinto cerrado 𝜀𝑖 : Emisividad de la superficie 𝑖
Las superficies son cuerpos grises 𝛼 ≈ 𝜀 𝛼𝑖 : Absortividad de la superficie 𝑖
𝜌𝑖 : Reflectividad de la superficie 𝑖
𝐸𝑏𝑖 : Potencia emisiva de cuerpo negro
1
2 𝐽𝑖 : Radiosidad (energía total emitida por la
superficie 𝑖)
𝐺𝑖 : Irradiación (energía total recibida por la
superficie 𝑖 de otras superficies)
𝐹𝑖𝑗 : Factor de vista
Balance de energía en la superficie 1:
𝑄̇
=⏟
𝜀1 𝐸𝑏1 + 𝜌1 𝐺1 − 𝐺1
𝐴1
𝐽1
Radiosidad
Calor retirado para
Emisión mantener la superficie
2 isotérmicamente a
𝑄̇
= 𝜀⏟2 𝐸𝑏2 + 𝜌2 𝐺2 − 𝐺2
𝐴2
𝐽2
Por otro lado suponiendo que la transmisividad𝜏 = 0 y que el cuerpo se aproxima como cuerpo gris
𝛼 = 𝜀 (cuerpo con emisividad independiente de la longitud de onda)
𝑄̇ = 𝜀1 𝐴1 (𝐸𝑏1 − 𝐺1 ) = 𝜀2 𝐴 2 (𝐸𝑏2 − 𝐺2 )
𝐸𝑏1 − 𝐺1 𝐸𝑏2 − 𝐺2
𝑄̇ = 1 = 1
𝜀1 𝐴1 𝜀2 𝐴2
𝐽1 − 𝜀1 𝐸𝑏1
𝐽1 = 𝜀1 𝐸𝑏1 + (1 − 𝜀1 )𝐺1 ⇒ 𝐺1 =
1 − 𝜀1
𝐽1 − 𝜀1 𝐸𝑏1
𝐸𝑏1 − 𝐸𝑏1 − 𝐽1
1 − 𝜀1
𝑄̇ = 1 =
1 − 𝜀1
𝜀1 𝐴1 𝜀1 𝐴1
Por lo tanto
SOLUCIÓN:
Como la placa 1 "observa completamente" la placa 2 ya que son placas infinitas, entonces el factor
de vista 𝐹12 es 1 (porque toda la radiación de la placa 1 llega a la placa 2). Además las placas "no se
miran a sí mismas" por lo que 𝐹11 = 𝐹22 = 0 .
(a) El intercambio de radiación entre las dos superficies viene representado mediante el siguiente
diagrama de resistencias térmicas
ℎ1 𝑇1 + ℎ 2 𝑇2 12 ∙ 150 + 12 ∙ 40
ℎ 2 𝐴 2(𝑇∞ − 𝑇2 ) = ℎ1 𝐴1 (𝑇1 − 𝑇∞ ) ⇒ 𝑇∞ = = = 95 ℃
ℎ1 + ℎ 2 12 + 12
𝜎(𝑇14 − 𝑇24 )
𝑄̇ = 𝑄̇𝑟𝑎𝑑 + 𝑄̇𝑐𝑜𝑛𝑣,1 = + ℎ1 𝐴1 (𝑇1 − 𝑇∞ )
1 − 𝜀1 1 1 − 𝜀2
+ +
𝜀1 𝐴1 𝐹12 𝐴1 𝜀2 𝐴2
(a) Si la esfera interior se encuentra a 177 °C y la exterior a 47 °C, determine la velocidad neta
de transferencia de calor entre ambas superficies esféricas.
(b) En qué porcentaje disminuye la transferencia de calor si se intercala una superficie esférica
muy delgada de 30 cm de diámetro (considere la misma emisividad de las dos anteriores por
cada lado) ¿La velocidad de transferencia sería la misma si en lugar de 30 cm de diámetro
tuviese 35 cm?
SOLUCIÓN:
Como toda la radiación de la esfera interna (superficie 1) llega a la esfera externa (superficie 2)
entonces en factor de vista 𝐹12 es 1.
2 2
1 1
S
Así la transferencia de calor con esta nueva superficie 𝑄̇ 𝑆viene dada por
𝐸𝑏1 − 𝐸𝑏2
𝑄̇𝑆 =
1 − 𝜀1 1 1 − 𝜀𝑆 1 − 𝜀𝑆 1 1 − 𝜀2
+ + + + +
𝜀1 𝐴1 𝐹1𝑆 𝐴1 𝜀𝑆 𝐴𝑆 𝜀𝑆 𝐴𝑆 𝐹𝑆2 𝐴𝑆 𝜀2 𝐴 2
𝑄̇𝑆 1 1
= =
𝑄̇ 1
(2
1 − 𝜀𝑆
+ 1) 2
1 − 𝜀𝑆
+1
𝐴 𝜀𝑆 𝐴1 𝜀𝑆
1 + 1 1 −𝑆𝜀 𝐴 1 − 𝜀2 1 + ∙
1 𝐴 𝑆 1 − 𝜀1 𝐴 1 − 𝜀2
( +1 + 1∙ ) + 1+ 1 ∙
𝐴1 𝜀1 𝐴2 𝜀2 𝜀1 𝐴2 𝜀2
Sustituyendo valores
𝑄̇𝑆 1
= = 0,657 (66 %)
𝑄̇ 2
1 − 0,75
2∙ +1
0,2 0,75
1+( ) ∙
0,3 1 − 0,75 0,2 2 1 − 0,75
+1 +( ) ∙
0,75 0,4 0,75
Por lo tanto, al intercalar una esfera con la misma emisividad de las anteriores de un diámetro de 30
cm se reduce un 34% la transferencia de calor. Como la fórmula anterior depende del diámetro de la
esfera intercalada 𝐷𝑆 entonces la velocidad de transferencia no es la misma al agregar una esfera de
diámetro 35 cm en lugar de una de 30 cm, de hecho mientras 𝐷𝑆 aumente reduce menos la
transferencia de calor.
PROBLEMA 40
Un horno tiene forma cúbica con aristas de 50 cm. La temperatura de la base es de 300ºC y la de la
tapa superior y las laterales 100ºC. La emisividad de la base es 0,8, la de la tapa superior 0,75 y la de
las paredes laterales 0,6, determine:
(a) La velocidad neta de transmisión de calor entre la base y las otras cinco caras.
(b) La velocidad neta de transferencia de calor entre la base y las otras superficies si las cuatro
caras laterales se encuentran aisladas exteriormente. Halle la temperatura de las caras
laterales en este caso. Compare las dos situaciones.
SOLUCIÓN:
La transferencia de calor entre una superficie real 𝑖 con las otras 𝑛 − 1 superficies viene dado por
𝑛
𝐸𝑏𝑖 − 𝐽𝑖
𝑄̇𝑖 = = ∑𝐹𝑖𝑗 𝐴 𝑖 (𝐽𝑖 − 𝐽𝑗 )
1 − 𝜀𝑖
𝑗=1
𝜀𝑖 𝐴 𝑖
2 6
4 3 5
1
(a) Formulando la expresión anterior para las superficies 1, 2 y 3 (las superficies 4, 5 y 6 son
iguales a la superficie 3):
𝐸𝑏1 − 𝐽1
𝑄̇1 = = 𝐹12 𝐴1 (𝐽1 − 𝐽2 ) + 4𝐹13 𝐴1 (𝐽1 − 𝐽3 )
1 − 𝜀1
𝜀1 𝐴1
𝐸𝑏2 − 𝐽2
𝑄̇2 = = 𝐹21 𝐴 2 (𝐽2 − 𝐽1 ) + 4𝐹23 𝐴 2 (𝐽2 − 𝐽3 )
1 − 𝜀2
𝜀2 𝐴2
𝐸𝑏3 − 𝐽3
𝑄̇3 = = 𝐹31 𝐴 3 (𝐽3 − 𝐽1 ) + 𝐹32 𝐴 3 (𝐽3 − 𝐽2 )
1 − 𝜀3
𝜀3 𝐴3
Arreglando en un sistema lineal de ecuaciones:
𝜀1 𝜀1
(𝐹12 + 4𝐹13 + ) 𝐽 − 𝐹12 𝐽2 − 4𝐹13 𝐽3 = 𝐸
1 − 𝜀1 1 1 − 𝜀1 𝑏1
𝜀2 𝜀2
−𝐹21 𝐽1 + ( 𝐹21 + 4𝐹23 + ) 𝐽2 − 4𝐹23 𝐽3 = 𝐸
1 − 𝜀2 1 − 𝜀2 𝑏2
𝜀3 𝜀3
−𝐹31 𝐽1 − 𝐹32 𝐽2 + (𝐹31 + 𝐹32 + ) 𝐽3 = 𝐸
{ 1 − 𝜀3 1 − 𝜀3 𝑏3
El factor de vista 𝐹12 se busca a partir de la gráfica de factor de vista entre dos placas paralelas. Como
es un cubo, en la siguiente figura se pone 𝑋 = 𝑌 = 𝐿. 𝑋/𝐿 = 1 , 𝑌/𝐿 = 1 . En esta figura se lee el
factor de vista entre las superficies 1 y 2 de 𝐹12 = 0,2
Como las placas no se ven a sí mismas 𝐹11 = 𝐹22 = 𝐹33 = 0 y por la regla de recinto cerrado:
𝐹11 + 𝐹12 + 4𝐹13 = 1 ⇒ 4𝐹13 = 1 − 𝐹11 − 𝐹12 = 1 − 0 − 0,2 = 0,8 ⇒ 𝐹13 = 0,2
𝐹21 + 𝐹22 + 4𝐹23 = 1 ⇒ 4𝐹23 = 1 − 𝐹21 − 𝐹22 = 1 − 0,2 − 0 = 0,8 ⇒ 𝐹23 = 0,2
𝐴2
𝐴 2 𝐹23 = 𝐴 3 𝐹32 ⇒ 𝐹32 = 𝐹 = 𝐹23 = 0,2
𝐴 3 23
𝐴1
𝐴1 𝐹13 = 𝐴 3 𝐹31 ⇒ 𝐹31 = 𝐹 = 𝐹13 = 0,2
𝐴 3 13
(b) Formulando la expresión anterior para la superficie 3 (las superficies 4, 5 y 6 son iguales a la
superficie 3). Además la superficie de las caras laterales son aisladas por lo que la
transferencia neta de calor en estas es cero 𝑄̇ 3 = 0
Sólo se modifica los dos últimos valores del sistema de ecuaciones anteriores
2
3
1
SOLUCIÓN:
Método analítico:
𝑅 2
𝑋 − √𝑋2 − 4 ( 2 )
𝑅1
𝐹12 =
2
Donde
𝑟1
𝑅1 =
𝐿
𝑟2
𝑅2 =
𝐿
1 + 𝑅22
𝑋=
𝑅12
𝐸𝑏1 − 𝐽1
𝑄̇1 = = 𝐹12 𝐴1 (𝐽1 − 𝐽2 ) + 𝐹13 𝐴1 (𝐽1 − 𝐽3 )
1 − 𝜀1
𝜀1 𝐴1
𝐸𝑏2 − 𝐽2
𝑄̇ 2 = = 𝐹21 𝐴 2 (𝐽2 − 𝐽1 ) + 𝐹23 𝐴 2 (𝐽2 − 𝐽3 )
1 − 𝜀2
𝜀2 𝐴2
𝐸𝑏3 − 𝐽3
𝑄̇ 3 = = 𝐹31 𝐴 3 (𝐽3 − 𝐽1 ) + 𝐹32 𝐴 3 (𝐽3 − 𝐽2 )
1 − 𝜀3
𝜀3 𝐴3
Arreglando en un sistema lineal de ecuaciones:
𝑄̇1 300
=𝜋 = 9549,3 W/m2
𝐴1 ∙ 0,22
4
Sólo se modifica los dos últimos valores del sistema de ecuaciones anteriores