Unidad I: Conducción de

calor unidimensional en
estado estacionario
Operaciones de Transferencia II

1
 Introducción
En esta sección se analizarán aplicaciones de la ley de Fourier de la conducción de
calor relacionados con el cálculo del flujo de calor en algunos sistemas
unidimensionales simples. En el caso de sistemas cilíndricos y esféricos son
unidimensionales cuando la temperatura en el cuerpo es una función solamente de
la distancia radial y es independiente del ángulo azimutal o de la distancia axial.
En algunos problemas bidimensionales, el efecto de una segunda coordenada
espacial puede ser tan pequeño como para justificar el análisis unidimensional y el
problema de flujo de calor multidimensional puede ser aproximado con la
descripción de una coordenada espacial. En estos casos las ecuaciones diferenciales
son simplificadas, lo cual conduce a una solución mucho más fácil como resultado de
esta simplificación.
2
 Pared plana
Considere una pared plana donde se aplica directamente la ley de
Fourier: 𝒅𝑻
calor 𝑸 = −𝒌𝑨
𝒅𝒙
Perfil de
temperatura
𝑻𝟏 Por integración (considerando k constante):
𝑻𝟐 −𝒌𝑨
𝑸= 𝑻𝟐 − 𝑻𝟏
∆𝒙

∆𝒙 𝑨 : Área transversal a la dirección de flujo de calor

3
 Pared plana
Si k presenta alguna dependencia con la temperatura, por ejemplo:

calor
𝒌 = 𝒌𝟎 𝟏 + 𝜷𝑻
Perfil de Por integración:
temperatura 𝒙𝟏 𝑻𝟏
𝑻𝟏
𝑸 𝒅𝒙 = − 𝒌𝟎 𝑨 𝟏 + 𝜷𝑻 𝒅𝑻
𝑻𝟐 𝒙𝟎 𝑻𝟎

𝒌𝟎 𝑨 𝜷 𝟐
𝑸=− 𝑻𝟐 − 𝑻𝟏 + 𝑻𝟐 − 𝑻𝟐𝟏
∆𝒙 𝟐
∆𝒙
4
 Pared plana compuesta
Si la pared está compuesta de varias capas de material (multicapas):
En estado estacionario:
𝑸𝑨 = 𝑸𝑩 = 𝑸𝑪 = 𝑸
A Perfil de
B
temperatura
𝑸 Por integración (con k constante):
C
−𝒌𝑨 𝑨 −𝒌𝑩 𝑨
𝑸= 𝑻𝟐 − 𝑻𝟏 = 𝑻𝟑 − 𝑻𝟐
∆𝒙𝑨 ∆𝒙𝑩
∆𝒙𝑨 ∆𝒙𝑩 ∆𝒙𝑪 −𝒌𝑪 𝑨
= 𝑻𝟒 − 𝑻𝟑
∆𝒙𝑪
𝟏 𝟐 𝟑 4
𝑨 : Área transversal a la dirección de flujo de calor 5
 Pared plana compuesta
Análogo a un circuito eléctrico:
𝑸
𝑹𝑨 𝑹𝑩 𝑹𝑪
A Perfil de
B
temperatura
𝑸 𝑻𝟏 ∆𝒙𝑨 𝑻𝟐 ∆𝒙𝑩 ∆𝒙𝑪
𝑻𝟑 𝑻𝟒
C 𝒌𝑩 𝑨 𝒌𝑪 𝑨
𝒌𝑨 𝑨

∆𝒙𝑨 ∆𝒙𝑩 ∆𝒙𝑪 𝑻𝟏 − 𝑻𝟒

Resolviendo 𝑸=
∆𝒙𝑨 ∆𝒙𝑩 ∆𝒙𝑪
triple igualdad: + +
𝟏 𝟐 𝟑 4 𝒌𝑨 𝑨 𝒌𝑩 𝑨 𝒌𝑪 𝑨
6
 Conceptual
La velocidad de transferencia de calor puede considerarse compuesta por un
gradiente térmico y una resistencia al flujo que surge de la combinación de
conductividad térmica, espesor de material y área. La ecuación de Fourier se puede
escribir de forma análoga a ley de Ohm (teoría de circuito eléctrico):
𝒈𝒓𝒂𝒅𝒊𝒆𝒏𝒕𝒆 𝒕é𝒓𝒎𝒊𝒄𝒐
𝑭𝒍𝒖𝒋𝒐 𝒅𝒆 𝒄𝒂𝒍𝒐𝒓 =
𝒓𝒆𝒔𝒊𝒔𝒕𝒆𝒏𝒄𝒊𝒂 𝒕é𝒓𝒎𝒊𝒄𝒂
La analogía eléctrica puede utilizarse para resolver problemas más complejos que
implican resistencias térmicas tanto en serie como en paralelo. La ecuación
unidimensional de flujo de calor para este tipo de problema puede escribirse:
∆𝑻 º𝑪 º𝑭 ∙ 𝒉
𝑸= 𝑹𝒊 ≡ ≡
𝑹𝒊 𝑾 𝑩𝒕𝒖
7
 Resistencias en serie y paralelo
Imagine un sistema de cierta complejidad como se indica en la Figura.

A B F 𝑸 𝑹𝑩
𝑹𝑬 𝑹𝑭
𝑸 E
C
𝑻𝟏 𝑹𝑨 𝑹𝑪
G 𝑹𝑫
D
𝑻𝟒 𝑹𝑮 𝑻𝟓
𝑻𝟐 𝑻𝟑
𝟏 𝟓
𝟐 𝟑 𝟒
∆𝑻 º𝑪 ∙ 𝒎𝟐 º𝑭 ∙ 𝒉 ∙ 𝒇𝒕𝟐
𝒄𝒐𝒏: 𝑹 ≡ ≡
La resistencia se define según: 𝑹 = 𝑸 𝑨 𝑾 𝑩𝒕𝒖
Se necesita resistencia en términos del área transversal debido a que no es constante como en
caso anterior 8
 Sistemas Radiales
Considere un cilindro largo de radio interior ri, radio exterior r0 y longitud
L, el cual se somete a un diferencial de temperatura Ti - T0.

Para un cilindro de longitud muy grande

en comparación con el diámetro, se
𝑸
puede suponer que el calor fluye sólo en
𝒓𝟎 la dirección radial, de modo que la única
𝒓 coordenada espacial necesaria para
𝒓𝒊
especificar el sistema es “r”.
𝒅𝒓
𝑳
9
 Sistemas Radiales
El área de flujo de calor en el sistema cilíndrico es: 𝑨𝒓 = 𝟐𝝅𝒓𝑳 Área del manto !!!

De la Ley de Fourier:
𝑸 𝒅𝑻 𝒅𝑻
𝑸 = −𝒌𝑨 = −𝟐𝝅𝒌𝒓𝑳
𝒅𝒓 𝒅𝒓
𝒓𝟎
Integrando (k cte):
𝒓𝒊 𝒓
𝒓𝒊 𝑻𝒊
𝒅𝒓 𝒅𝒓
𝑸 = −𝟐𝝅𝒌𝑳 𝒅𝑻
𝒓𝟎 𝒓 𝑻𝟎
𝑳

10
 Sistemas Radiales
𝒓𝒊
Resolviendo: 𝑸 𝒍𝒏 𝒓 = −𝟐𝝅𝒌𝑳 𝑻𝒊 − 𝑻𝟎
𝟎
𝑻𝒊 − 𝑻𝟎
Ordenando: 𝑸 =
𝑸
𝒓𝟎
𝒍𝒏
𝒓𝒊
𝒓𝟎 𝟐𝝅𝒌𝑳
𝒓𝒊 𝒓 𝑸 𝒓𝟎
𝒍𝒏
𝒅𝒓 𝒓𝒊
𝑹=
𝑻𝒊 𝟐𝝅𝒌𝑳
𝑳 𝑻𝟎

11
 Sistemas Radiales
En el caso de pared compuesta…. 𝑹𝑨 𝑹𝑩 𝑹𝑪
𝑻𝟒 𝑸 𝑻𝟏 𝑻𝟐 𝑻𝟑 𝑻𝟒
𝑻𝟑
𝑸
𝑻𝟐 𝒓 𝒓𝟑 𝒓𝟒
𝒍𝒏 𝒓𝟐 𝒍𝒏 𝒓 𝒍𝒏
𝒓𝟑
𝒓𝟐 𝟏 𝟐

A 𝒓𝟏 𝟐𝝅𝒌𝑨 𝑳 𝟐𝝅𝒌𝑩 𝑳 𝟐𝝅𝒌𝑪 𝑳

B
C Condiciones de borde:
𝑻𝟏 𝒓𝟑
𝑻 = 𝑻𝒊 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝒊
𝒓𝟒
𝒅𝒐𝒏𝒅𝒆 𝒊 = 𝟏, 𝟐, 𝟑 𝒚 𝟒

12
 Sistemas Radiales
Flujo de calor para las 3 secciones:

𝑸 𝒓𝟐
𝑻𝟏 − 𝑻𝟐 = 𝒍𝒏
𝟐𝝅𝒌𝑨 𝑳 𝒓𝟏 Sumando:
𝑸 𝒓𝟑
𝑻𝟐 − 𝑻𝟑 = 𝒍𝒏 𝟐𝝅𝑳 𝑻𝟏 − 𝑻𝟒
𝟐𝝅𝒌𝑩 𝑳 𝒓𝟐 𝑸=
𝒓𝟐 𝒓𝟑 𝒓𝟒
𝑸 𝒓𝟒 𝒍𝒏 𝒍𝒏 𝒍𝒏
𝒓𝟏 𝒓𝟐 𝒓𝟑
𝑻𝟑 − 𝑻𝟒 = 𝒍𝒏 + +
𝟐𝝅𝒌𝑪 𝑳 𝒓𝟑 𝒌𝑨 𝒌𝑩 𝒌𝑪

13
 Esferas
En sistemas esféricos es posible el tratamiento unidimensional cuando la
temperatura cambia con el radio.
El área de flujo de calor en 𝑨 = 𝟒𝝅𝒓𝟐
el sistema esférico es:
𝑸
𝟐
𝒅𝑻
De la Ley de Fourier: 𝑸 = −𝟒𝝅𝒌𝒓
𝒅𝒓
𝒅𝒓 Integrando (con k constante):
𝒓𝒊 𝒓𝒊
𝒅𝒓 𝑻𝒊
𝒓𝟎 𝑸 𝟐
= −𝟒𝝅𝒌 𝒅𝑻
𝒓𝟎 𝒓 𝑻𝟎
14
Esferas
𝟒𝝅𝒌 𝑻𝒊 − 𝑻𝟎
Resolviendo integral: 𝑸=
𝟏 𝟏
𝑸 −
𝒓𝒊 𝒓𝟎

𝒅𝒓
Análogo a un circuito eléctrico:
𝑸
𝒓𝒊
𝒓𝟎
𝑻𝒊 𝑻𝟎
𝟏 𝟏
−
𝒓𝒊 𝒓𝟎
𝑹=
𝟒𝝅𝒌
15
 Condiciones de Borde: Convección
El flujo de calor por convección puede ser determinado según:
𝑻𝒘 : temperatura de la superficie del sólido
𝑸 = 𝒉 𝑨 𝑻𝒘 − 𝑻∞
𝑻∞ : temperatura del fluido

Planteando la analogía de resistencia eléctrica:

𝑻𝒘 − 𝑻∞ Entonces la resistencia se define según:

𝑸= 𝟏
𝟏 𝑹=
𝒉𝑨 𝒉𝑨
16
 Coeficiente global de transferencia de calor
Considere una pared plana expuesta a fluido caliente A a un lado y un
fluido frío B al otro lado.
Q El calor transferido será.
𝑻𝑨 𝒌𝑨
𝑸 = 𝒉𝟏 𝑨 𝑻𝑨 − 𝑻𝟏 = 𝑻𝟏 − 𝑻𝟐 = 𝒉𝟐 𝑨 𝑻𝟐 − 𝑻𝑩
Fluido B

𝒉𝟏 ∆𝒙

𝑻𝟏 𝑻𝟐 Sumando estas ecuaciones:

Fluido A

𝒉𝟐 𝑻𝑨 − 𝑻𝑩
𝑸=
𝑻𝑩 𝟏 ∆𝒙 𝟏
+ +
∆𝒙
𝒉𝟏 𝑨 𝒌𝑨 𝒉𝟐 𝑨
17
 Coeficiente global de transferencia de calor
Por analogía con un 𝟏 ∆𝒙 𝟏
circuito eléctrico: 𝒉𝟏 𝑨 𝒉𝟐 𝑨
𝒌𝑨
𝑸 𝑻𝑨 𝑻𝟏 𝑻𝟐 𝑻𝑩
Q
𝑻𝑨
Fluido B

𝒉𝟏 El coeficiente global de transferencia de global considera

la sumatoria de las resistencias térmicas.
𝑻𝟏 𝑻𝟐
Fluido A

−𝟏
𝒉𝟐 𝟏 ∆𝒙 𝟏
𝑸 = 𝑼𝑨 𝑻𝑨 − 𝑻𝑩 𝑼= + +
𝑻𝑩 𝒉𝟏 𝒌 𝒉𝟐
∆𝒙
18
 Coeficiente global de transferencia de calor
Para el caso de un cilindro hueco que experimenta transferencia de calor
convección por ambas superficies.
Se asume que la temperatura de
Fluido B , 𝑻𝑩
ambos fluidos es constante!!
𝒓𝟎 1.- Convección interna:

𝒓𝒊 𝑸 = 𝒉𝒊 𝑨𝒊 𝑻𝑨 − 𝑻𝒊 = 𝒉𝒊 𝟐𝝅𝒓𝒊 𝑳 𝑻𝑨 − 𝑻𝒊
𝑳 Ordenando:
Condiciones de borde: 𝑸 𝑸
𝑻 = 𝑻𝒊 𝑨 = 𝑨𝒊 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝒊 𝑻𝑨 − 𝑻𝒊 = =
𝟐𝝅𝒓𝒊 𝑳𝒉𝒊 𝒉𝒊 𝑨𝒊
𝑻 = 𝑻𝟎 𝑨 = 𝑨𝟎 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝟎
19
 Coeficiente global de transferencia de calor
Fluido B , 𝑻𝑩 2.- Conducción:
𝑻𝒊 − 𝑻𝟎
𝒓𝟎 𝑸 =
𝒓𝟎
𝒍𝒏
𝒓𝒊
𝒓𝒊
𝟐𝝅𝒌𝑳
𝑳
Condiciones de borde: Ordenando:
𝑻 = 𝑻𝒊 𝑨 = 𝑨𝒊 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝒊
𝒍𝒏 𝒓𝟎 𝒓𝒊
𝑻 = 𝑻𝟎 𝑨 = 𝑨𝟎 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝟎 𝑻𝒊 − 𝑻𝟎 = 𝑸
𝟐𝝅𝒌𝑳

20
 Coeficiente global de transferencia de calor
Fluido B , 𝑻𝑩 2.- Convección externa:

𝒓𝟎 𝑸 = 𝒉𝟎 𝑨𝟎 𝑻𝟎 − 𝑻𝑩 = 𝒉𝟎 𝟐𝝅𝒓𝟎 𝑳 𝑻𝟎 − 𝑻𝑩

𝒓𝒊
𝑳 Ordenando:
Condiciones de borde: 𝑸 𝑸
𝑻 = 𝑻𝒊 𝑨 = 𝑨𝒊 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝒊 𝑻𝟎 − 𝑻𝑩 = =
𝟐𝝅𝒓𝟎 𝑳𝒉𝟎 𝒉𝟎 𝑨𝟎
𝑻 = 𝑻𝟎 𝑨 = 𝑨𝟎 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝟎

21
 Coeficiente global de transferencia de calor
∆𝑻
Recordando que: 𝑸=
𝑹𝒊
Fluido B , 𝑻𝑩 El calor es transferido en la dirección radial:
𝒓𝟎 𝟏 𝒓𝟎 𝟏
𝒍𝒏
𝒓𝒊
𝒉𝒊 𝑨𝒊 𝒉𝟎 𝑨𝟎
𝒓𝒊 𝑸 𝑻𝑨 𝟐𝝅𝒌𝑳 𝑻𝑩
𝑳 𝑻𝒊 𝑻𝟎
Condiciones de borde:
𝑻 = 𝑻𝒊 𝑨 = 𝑨𝒊 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝒊
𝑻 = 𝑻𝟎 𝑨 = 𝑨𝟎 𝒑𝒂𝒓𝒂 𝒓 = 𝒓𝟎
22
 Coeficiente global de transferencia de calor
Sumando las tres ecuaciones:
𝑸 𝑻𝑨 − 𝑻𝑩
𝑻𝑨 − 𝑻𝒊 = 𝑸=
𝒉𝒊 𝑨𝒊 𝟏 𝒍𝒏 𝒓𝟎 𝒓𝒊 𝟏
+ +
𝒉𝒊 𝑨𝒊 𝟐𝝅𝒌𝑳 𝒉𝟎 𝑨𝟎
𝒍𝒏 𝒓𝟎 𝒓𝒊
𝑻𝒊 − 𝑻𝟎 = 𝑸
𝟐𝝅𝒌𝑳
𝑸 𝟐𝝅𝑳 𝑻𝑨 − 𝑻𝑩
𝑻𝟎 − 𝑻𝑩 = 𝑸=
𝒉𝟎 𝑨𝟎 𝟏 𝒍𝒏 𝒓𝟎 𝒓𝒊 𝟏
+ +
𝒉𝒊 𝒓𝒊 𝒌 𝒉𝟎 𝒓𝟎

23
 Coeficiente global de transferencia de calor
El flujo de calor se plantea en término del 𝑸 = 𝑼𝑨 ∆𝑻
coeficiente global de transferencia de calor:
El coeficiente de transferencia puede ser definido por ejemplo en
función de la superficie externa: 𝑨𝟎 = 𝟐𝝅𝒓𝟎 𝑳
𝑸 = 𝑼𝟎 𝟐𝝅𝒓𝟎 𝑳 𝑻𝑨 − 𝑻𝑩 −𝟏
𝒓𝟎
Igualando con expresión anterior: 𝟏 𝒍𝒏 𝟏
𝒓𝒊
𝑼𝟎 = 𝒓−𝟏
𝟎 + +
𝟐𝝅𝑳 𝑻𝑨 − 𝑻𝑩 𝒉𝒊 𝒓𝒊 𝒌 𝒉𝟎 𝒓𝟎
𝑸=
𝟏 𝒍𝒏 𝒓𝟎 𝒓𝒊 𝟏
+ +
𝒉𝒊 𝒓𝒊 𝒌 𝒉𝟎 𝒓𝟎 24
 Espesor crítico de aislación
Considere una capa de aislante ubicado alrededor de una tubería circular
(de espesor despreciable). La temperatura interna de la tubería (Ti) es
conocida y constante.
𝑻∞ : temperatura ambiente
𝑻∞ 𝑸

𝑻 𝑸 𝑻𝒊 𝑻 𝑻∞
𝒓𝟎
𝑻𝒊 𝒓𝟎
𝒍𝒏 𝟏
𝒓𝒊 𝒓𝒊
𝟐𝝅𝒓𝟎 𝑳 𝒉
𝟐𝝅𝒌𝑳
aislante Conducción !! Convección !!
25
 Espesor crítico de aislación
Calor transferido por conducción:
𝑻∞ 𝑻𝒊 − 𝑻 𝒍𝒏 𝒓𝟎 𝒓𝒊
𝑸 𝑸 = 𝑻𝒊 − 𝑻 = 𝑸
𝒍𝒏 𝒓𝟎 𝒓𝒊 𝟐𝝅𝒌𝑳
𝑻 𝟐𝝅𝒌𝑳
𝒓𝟎 Calor transferido por convección:
𝑻𝒊
𝒓𝒊 𝑸
𝑸 = 𝒉 𝟐𝝅𝒓𝟎 𝑳 𝑻 − 𝑻∞ 𝑻 − 𝑻∞ =
𝟐𝝅𝒓𝟎 𝑳 𝒉
aislante

26
 Espesor crítico de aislación
Sumando las 𝟐𝝅𝑳 𝑻𝒊 − 𝑻∞
𝑸=
𝑻∞ expresiones: 𝒍𝒏 𝒓𝟎 𝒓𝒊 𝟏
+
𝑸 𝒌 𝒉𝒓𝟎

𝑻 El espesor óptimo de aislante se determina

𝒓𝟎 mediante el valor de r0 que minimiza Q (o bien
𝑻𝒊
maximiza resistencia total). Para ello:
𝒓𝒊
𝒅𝑸
aislante =𝟎
𝒅𝒓𝟎

27
Espesor crítico de aislación
𝟏 𝟏
−𝟐𝝅𝑳 𝑻𝒊 − 𝑻∞ −
𝒅𝑸 𝒌𝒓𝟎 𝒉𝒓𝟐𝟎
Por tanto: = 𝟐
=𝟎
𝒅𝒓𝟎 𝒍𝒏 𝒓𝟎 𝒓𝒊 𝟏
+
𝒌 𝒉𝒓𝟎

𝟏 𝟏 𝒌
Finalmente: − 𝟐=𝟎 𝒓𝟎 =
𝒌𝒓𝟎 𝒉𝒓𝟎 𝒉

28