Laboratorio de Sistemas Térmicos
Métodos numéricos para la conducción de calor en estado estacionario -
Práctica I
Docente: Omar Ardila – omar2198186@[Link]
August 20, 2025
Nota: se presenta la solución manual al problema planteado.
Problema 3
Considere una placa grande de uranio de espesor L = 4 cm y conductividad térmica k = 28 W/m◦ C
en la cual se genera calor de manera uniforme a una velocidad constante de e = 7 × 105 W/m3 . El
lado izquierdo de la placa es sometido a un flujo de calor de q0 = 700W/m2 , mientras que el derecho
está sujeto a convección hacia un medio ambiente a T∞ = 5°C, con un coeficiente de transferencia
de calor de h = 110W/m2 °C, como se muestra en la figura. Si considera un total de cinco nodos
igualmente espaciados, dos en las fronteras y tres internos, estime el perfil de temperaturas en la
placa mediante el procedimiento de diferencias finitas.
Placa de
Uranio (U238) Ambiente
Flujo de calor
convectivo
k = 28 W/m °C
h = 110 W/m2 °C
e = 7×105 W/m3 T∞ = 5°C
q0 = 700 W/m2
L = 0.04 m
Consideraciones
• Transferencia de calor en estado estacionario
• Transferencia de calor unidimensional
• Conductividad térmica constante
• Transferencia de calor por radiación es despreciable
• La frontera izquierda está sometida a flujo de calor (q0 )
Análisis
Se conoce el número de nodos, n = 5, y el espesor de pared, L = 0.04m, por lo tanto, se calcula el
espaciamiento nodal como
L 0.04 m
∆x = = = 0.01 m
n−1 5−1
Se enumeran los nodos desde el 1 hasta el 5 como se indica en la Figura 1.
Figure 1: Malla unidimensional y condiciones de contorno
Placa de
Uranio (U238) Ambiente
Flujo de calor
convectivo
k = 28 W/m °C
h = 110 W/m2 °C
e = 7×105 W/m3 T∞ = 5°C
q0 = 700 W/m2
1 2 3 4 5
∆x = 0.01 m
L = 0.04 m
La condición de frontera izquierda implica que existe un flujo de calor constante hacia la derecha
de q0 = 700 W/m2 , y su formulación en diferencias finitas puede obtenerse directamente a partir
de la Ecuación 1, para n = 1
Formulación para frontera izquierda con flujo de calor
T2 − T1
q0 A + kA + e(A∆x/2) = 0 (1)
∆x
Al reemplazar y eliminar el área (A) se obtiene
T2 − T1
q0 + k + e(∆x/2) = 0
∆x
Los nodos 2, 3 y 4 son internos, y su formulación en diferencias finitas se obtiene directamente a
partir de la Ecuación 2, para n = 2, 3, 4, respectivamente.
Formulación para nodos internos
Tn−1 − 2Tn + Tn+1 e
+ =0 (2)
(∆x)2 k
Reemplazando se obtienen las siguientes ecuaciones
T1 − 2T2 + T3 e
+ =0 ; n=2
(∆x)2 k
T2 − 2T3 + T4 e
+ =0 ; n=3
(∆x)2 k
T3 − 2T4 + T5 e
+ =0 ; n=4
(∆x)2 k
Hasta ahora se disponen de 4 ecuaciones con 5 incógnitas, que son T1 , T2 , T3 , T4 , T5 . Se requiere
de una quinta ecuación, que puede obtenerse con la condición de contorno derecha (convección),
reemplazando n = 5 en la Ecuación 3.
Formulación para convección en frontera derecha
Tn−1 − Tn
hA(T∞ − Tn ) + kA + e(A∆x/2) = 0 (3)
∆x
Al reemplazar y eliminar el área (A) obtiene
T4 − T5
h(T∞ − T5 ) + k + e(∆x/2) = 0
∆x
De este modo ya se tiene un sistema de 5 ecuaciones con 5 incógnitas, el cual puede resolverse
con ayuda de diferentes métodos y herramientas computacionales, como Matlab, Python, EES,
etc. Al resolver este sistema, se obtienen los resultados mostrados en la Tabla 1, y el perfil de
temperaturas mostrado en la Figura 2.
Table 1: Solución por diferencias finitas
Nodo 1 2 3 4 5
x [m] 0.00 0.01 0.02 0.03 0.04
T [°C] 286.89 285.39 281.38 274.39 265.87
Figure 2: Perfil de temperatura en la placa utilizando 5 nodos
Posición nodal [m]
