Transferencia de Calor
Prof. Dr. Hernán Olguín Astudillo
e-mail:
[email protected] Oficina C-208B
Horario de Atención: Martes: 13:00-13:45
Condiciones de Frontera y
Condiciones Iniciales
Condiciones de Frontera e Iniciales
¿Para qué las necesitamos?
Ecuación de Laplace:
𝑑2𝑇
𝛻 2𝑇 = 0 2
=0
𝑑𝑥
Primera integración:
𝑑 𝑑𝑇 𝑑𝑇
=0 = 𝐶1
𝑑𝑥 𝑑𝑥 𝑑𝑥
Segunda integración integración:
𝑇 = 𝐶1 𝑥 + 𝐶2 Es una familia de soluciones
C.F. Temperatura Específica
Caso unidimensional:
𝑇 0, 𝑡 = 𝑇1 𝑇 𝐿, 𝑡 = 𝑇2
Condiciones especificadas en un punto
Caso bidimensional:
𝑇 0, 𝑦, 𝑡 = 𝑇1 𝑇 𝐿, 𝑦, 𝑡 = 𝑇2
Condiciones especificadas en una linea
Caso tridimensional:
𝑇 0, 𝑦, 𝑧, 𝑡 = 𝑇1 𝑇 𝐿, 𝑦, 𝑧, 𝑡 = 𝑇2
Condiciones especificadas en una superficie
C.F. Temperatura Específica- Ejemplo
Una pared plana (k=20 W/(mK), L=1 m) tiene una temperatura de 100 K en
una de sus caras y 50 K en la otra. Calcular el flujo de calor, utilizando la
ecuación de calor en forma diferencial unidimensional y suponiendo estado
estacionario
𝜕𝑇
𝛻 ∙ 𝑘𝛻𝑇 + 𝑔 = 𝜌𝐶 𝛻 ∙ 𝑘𝛻𝑇 = 0 𝛻 2𝑇 = 0
𝜕𝑡
𝜕2𝑇 𝜕𝑇
2
=0 = 𝐶1 𝑇 = 𝐶1 𝑥 + 𝐶2
𝜕𝑥 𝜕𝑥
BC:
𝑇 0 = 100 = 𝐶1 ∗ 0 + 𝐶2 100 = 𝐶2
𝑇 𝐿 = 50 = 𝐶1 ∗ 𝐿 + 100 −50
𝐶1 =
𝐿
C.F. Temperatura Específica- Ejemplo
• Perfil de temperatura:
50
𝑇(𝑥) = − 𝑥 + 100
𝐿
• Se puede evaluar T en distintas ubicaciones:
50
𝑇(𝐿) = − 𝐿 + 100 = 50
𝐿
50 𝐿
𝑇(𝐿/2) = − + 100 = 75
𝐿 2
• Flujo de calor:
𝑊
𝑞 = 20 ∗ 50 = 1000 2
𝑚
¿Que rol juegan L y k?
C.F. Flujo Específico
Caso general:
𝜕𝑇 𝑊
𝑞𝑖 = −𝑘
𝜕𝑥𝑖 𝑚2
Para un caso unidimensional:
𝜕𝑇(0, 𝑡)
−𝑘 = 50
𝜕𝑥
𝜕𝑇 𝐿, 𝑡
−𝑘 = −50
𝜕𝑥
C.F. Flujo Específico
Frontera aislada:
𝜕𝑇(0, 𝑡) 𝜕𝑇(0, 𝑡)
−𝑘 =0 =0
𝜕𝑥 𝜕𝑥
Simetría térmica:
𝐿
𝜕𝑇 ,𝑡
2 =0
𝜕𝑥
C.F. Flujo Específico- Ejemplo
Una pared plana (k=20 W/(mK), L=1 m) tiene una fuente de calor que
entrega 1000 W/m2 en una de sus caras y una temperatura de 50 K en el
lado opuesto. Calcular la temperatura en la cara donde la fuente de calor es
aplicada. Utilizar la ecuación de calor en forma diferencial unidimensional y
suponer estado estacionario.
𝜕𝑇
𝛻 ∙ 𝑘𝛻𝑇 + 𝑔 = 𝜌𝐶 𝛻 ∙ 𝑘𝛻𝑇 = 0 𝛻 2𝑇 = 0
𝜕𝑡
BC
𝜕2𝑇 𝜕𝑇 𝜕𝑇
2
=0 = 𝐶1 −𝑘 = −𝑘𝐶1 = 1000
𝜕𝑥 𝜕𝑥 𝜕𝑥
1000
𝐶1 = − = −50 𝑇 = −50𝑥 + 𝐶2
20
C.F. Flujo Específico- Ejemplo
𝑇 = −50𝑥 + 𝐶2
BC: 𝑇(𝐿) = −50 + 𝐶2 = 50 𝐶2 = 100
Reemplazando en la expresión original:
𝑇(𝑥) = −50𝑥 + 100
Evaluando en el origen:
𝑇(0) = −50 ∗ 0 + 100 = 100 𝐾
C.F. Convección/Radiación
Condición de Convección:
𝜕𝑇 0, 𝑡 𝜕𝑇 𝐿, 𝑡
−𝑘 = ℎ[𝑇∞ − 𝑇 0, 𝑡 ] −𝑘 = ℎ[𝑇 𝐿, 𝑡 − 𝑇∞ ]
𝜕𝑥 𝜕𝑥
Condición de Radiación:
𝜕𝑇 0, 𝑡 𝜕𝑇 𝐿, 𝑡
−𝑘 = 𝜀𝜎 𝑇𝑠4 − 𝑇∞4 −𝑘 = 𝜀𝜎 𝑇∞4 − 𝑇𝑠4
𝜕𝑥 𝜕𝑥
C.F. Convección- Ejemplo
Una pared plana (k=20 W/(mK), L=1 m) tiene una temperatura de 350 K en
una de sus caras y una condición de borde convectiva en la otra. Si la
temperatura exterior es de 300 K y h=100 W/(m2K), calcular la temperatura
en la cara exterior de la pared. Utilizar la ecuación de calor en forma
diferencial unidimensional y suponer estado estacionario.
𝜕𝑇
𝛻 ∙ 𝑘𝛻𝑇 + 𝑔 = 𝜌𝐶 𝛻 ∙ 𝑘𝛻𝑇 = 0 𝛻 2𝑇 = 0
𝜕𝑡
BC1
𝜕2𝑇 𝜕𝑇 𝜕𝑇
2
=0 = 𝐶1 −𝑘 = −𝑘𝐶1 = ℎ[𝑇 𝐿 − 𝑇∞ ]
𝜕𝑥 𝜕𝑥 𝜕𝑥
100 𝑇 𝐿 − 𝑇∞
𝐶1 = − = −5 [𝑇 𝐿 − 𝑇∞ ]
20
C.F. Convección- Ejemplo
Seguimos con la integración
𝜕𝑇
= 𝐶1 𝑇(𝑥) = 𝐶1 𝑥 + 𝐶2
𝜕𝑥
BC2:
𝑇 0 = 350 = 𝐶1 ∗ 0 + 𝐶2 𝐶2 = 350
Reemplazando en la BC1:
𝐶1 = −5 𝑇 𝐿 − 𝑇∞ 𝐶1 = −5 [𝐶1 𝐿 + 𝐶2 − 𝑇∞ ]
𝐶1 = −5[𝐶1 + 350 − 300] 𝐶1 = −5[𝐶1 + 50]
−250
𝐶1 = = 41, 6 𝑇(𝑥) = −42 𝑥 + 350
6
C.F. Convección- Ejemplo
Temperatura en la cara exterior de la pared:
𝑇(𝑥) = −42 𝑥 + 350 𝑇(𝐿) = −42 + 350 = 308 K
Si no me cree, la pueden verificar!
𝑞𝑐𝑜𝑛𝑣 = 𝑞𝑐𝑜𝑛𝑑
𝜕𝑇
−𝑘 = ℎ∆𝑇 20 ∗ 42 = 100 ∗ 8
𝜕𝑥
840 ≠ 800 ¿¿Donde está el error??
C.F. Convección- Ejemplo
Temperatura en la cara exterior de la pared usando valores „exactos“
𝑇(𝑥) = −41, 6𝑥 + 350 𝑇(𝐿) = −41, 6 + 350 = 308, 3 K
Calculamos todo de nuevo:
𝑞𝑐𝑜𝑛𝑣 = 𝑞𝑐𝑜𝑛𝑑
𝜕𝑇
−𝑘 = ℎ∆𝑇 20 ∗ 41, 6 = 100 ∗ 8, 3
𝜕𝑥
833, 3 = 833, 3 𝑇(𝐿) = 308, 3 K
C.F. Interfase
Condición 1:
𝑇𝐴 𝑥𝑐 , 𝑡 = 𝑇𝐵 (𝑥𝑐 , 𝑡)
Condición 2:
𝜕𝑇𝐴 𝑥𝑐 , 𝑡 𝜕𝑇𝐵 𝑥𝑐 , 𝑡
−𝑘𝐴 = −𝑘𝐵
𝜕𝑥 𝜕𝑥
Generalización:
𝑇𝑜𝑑𝑜 𝑙𝑜 𝑞𝑢𝑒 𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜 𝑔𝑎𝑛𝑎/𝑝𝑖𝑒𝑟𝑑𝑒 = 𝑇𝑜𝑑𝑜 𝑙𝑜 𝑞𝑢𝑒 𝑒𝑙 𝑒𝑛𝑡𝑜𝑟𝑛𝑜 𝑔𝑎𝑛𝑎/𝑝𝑖𝑒𝑟𝑑𝑒
