Sección 3: Ecuación general de conducción de calor

En la sección anterior vimos cómo se puede calcular temperaturas en las superficies o los flujos
de calor que se transfieren a lo largo de un sistema unidireccional, sin embargo en muchas
aplicaciones estamos interesados en saber que pasa dentro de los volúmenes de control; con la
ecuación general de conducción podemos obtener el perfil o distribución de temperaturas a lo
largo de todo el sólido, además se tienen en cuenta las posibles fuentes o generadores de calor
dentro del sistema.

3.1 Distribución de temperatura para un sólido tridimensional

Conocer la distribución de temperaturas es útil cuando se quiere comprobar la integridad

estructural de un material y determinar los esfuerzos térmicos, las expansiones y deflexiones a
las que se puede ver sometido, especialmente para especificaciones de diseño o de seguridad.
En sistemas que requieren un aislante también es útil para determinar espesores óptimos.

Balance por conservación de energía

𝐸̇𝑒𝑛𝑡 𝐸̇𝑠𝑎𝑙𝑒
𝐸̇𝑔 𝐸̇𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑑𝑜

𝑒𝑛𝑒𝑟𝑔í𝑎 𝐽
𝐸̇ [=] = = 𝑊, Velocidad de energía térmica
𝑡𝑖𝑒𝑚𝑝𝑜 𝑠

𝐸𝑒𝑛𝑡 + 𝐸𝑔 − 𝐸𝑠𝑎𝑙𝑒 = 𝐸𝑎𝑐𝑢𝑚 , En un intervalo de tiempo ∆𝑡

𝐸[=] 𝐽 → Cantidades de los cambios de energía

Metodología:

1. Definir volumen de control finito o diferencial-infinitesimal.

2. Identificar base de tiempo.
3. Identificar procesos de energía apropiados que entran o salen del volumen de control.
4. Escribir ecuaciones de conservación de energía y sustituir expresión de flujo apropiadas.
Se espera al final de estos cuatro pasos obtener las expresiones que permitan:
1. Calcular el flujo de calor a través del volumen de control
2. Obtener la distribución de temperatura respecto a variables dimensionales, por ejemplo
respecto a 𝑥 si es unidireccional.

En el siguiente esquema se representa un sistema sencillo con transferencia de calor

unidireccional en dirección 𝑥:

∆𝑇
𝑇1 𝑞𝑥 𝛼 𝐴
∆𝑥
𝑇(𝑥) ∆𝑇
𝑞𝑥 = −𝑘𝐴
∆𝑥

𝑥 𝑇2

𝑞𝑥

Para un ∆𝑡 → 0
𝑑𝑇
𝑞𝑥 = −𝑘𝐴
𝑑𝑥

El calor debe transferirse en dirección al decrecimiento de la temperatura, sin embargo cuando

no tenemos un sistema simple como el anterior, es necesario expresar la dirección del calor por
conducción en términos de todas las coordenadas posibles y se debe reescribir la ecuación en
forma vectorial:
𝜕𝑇 𝜕𝑇 𝜕𝑇
𝑞 = −𝑘∇𝑇 = −𝑘 �𝑖 +𝑗 +𝑘 �
𝜕𝑥 𝜕𝑦 𝜕𝑧

Donde,
∇: Operador nabla, tridimensional.
𝑇(𝑥, 𝑦, 𝑧): Campo escalar de temperaturas.

Es de anotar que el flujo de calor siempre es perpendicular a las superficies isotérmicas cuando
el medio es isotrópico, es decir, que sus características físicas y químicas no dependen de la
dirección.

Aunque partiendo de la ecuación en forma vectorial se puede obtener la expresión de calor casi
para cualquier geometría, existen expresiones básicas para las geometrías más comunes:
- Rectangulares o cartesianas
- Cilíndricas
- Esféricas
3.1.1 Coordenadas cartesianas

En una geometría tipo rectangular, el calor por conducción está compuesto por tres
componentes, cada una en una dirección dada (𝑥, 𝑦 y 𝑧).

𝑞′′ = 𝑖𝑞′′𝑥 + 𝑗𝑞′′𝑦 + 𝑘𝑞′′𝑧

𝜕𝑇 𝜕𝑇 𝜕𝑇
𝑞′′𝑥 = −𝑘 , 𝑞′′𝑦 = −𝑘 , 𝑞′′𝑧 = −𝑘
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝑞𝑧+𝑑𝑧
𝑞𝑦+𝑑𝑦

𝐸̇𝑔
̇
𝐸𝑎𝑐𝑢𝑚
𝑑𝑧
𝑞𝑥 𝑞𝑥+𝑑𝑥

𝑧 𝑦
𝑑𝑦
𝑥
𝑑𝑥
𝑞𝑦
𝑞𝑧

Si realizamos un balance de calor en el volumen de control infinitesimal (𝑑𝑥 𝑑𝑦 𝑑𝑧)

representado por la gráfica anterior y para cada una de sus caras, podemos obtener expresiones
de los componentes de calor en cada dirección, por ejemplo, el componente 𝑥 de la rapidez de
transferencia de calor en 𝑥 + 𝑑𝑥 es igual al valor de este componente en 𝑥 más la cantidad por
la que cambia 𝑑𝑥, haciendo el mismo análisis en cada cara:

𝜕𝑞𝑥
𝑞𝑥+𝑑𝑥 = 𝑞𝑥 + 𝑑𝑥
𝜕𝑥

𝜕𝑞𝑦
𝑞𝑦+𝑑𝑦 = 𝑞𝑦 + 𝑑𝑦
𝜕𝑦
 𝜕𝑞𝑧
𝑞𝑧+𝑑𝑧 = 𝑞𝑧 + 𝑑𝑧
𝜕𝑧

Además se debe realizar un balance de energía global:

𝐸̇𝑒𝑛𝑡 + 𝐸̇𝑔 − 𝐸̇𝑠𝑎𝑙𝑒 = 𝐸̇𝑎𝑐𝑢𝑚

Donde 𝐸̇𝑔 es la energía generada dentro del volumen de control, la existencia de esta energía
depende de factores o condiciones del sistema, que puede involucrar reacciones exotérmicas o
endotérmicas, cambio de fases, entre otros. Y se representa como el calor generado 𝑞̇
multiplicado por el volumen de control:

𝐸̇𝑔 = 𝑞̇ (𝑑𝑥 𝑑𝑦 𝑑𝑧)

Para que esta relación sea coherente el calor generado debe tener unidades de volumen
𝑞̇ [=]𝑊/𝑚3 .

Si no hay cambios de fase en el volumen de control y asumiendo que la capacidad calorífica del
material y la densidad del mismo no cambian significativamente con la temperatura, la energía
acumulada 𝐸̇𝑎𝑐𝑢𝑚 puede calcularse en términos de la variación de la temperatura del sistema
respecto al tiempo:
𝜕𝑇
𝐸̇𝑎𝑐𝑢𝑚 = 𝜌𝐶𝑝 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑡
𝑘𝑔 𝐽
𝜌[=] , 𝐶𝑝 [=]
𝑚3 𝑘𝑔.𝐾

Reemplazando estas expresiones en el balance de energía:

𝜕𝑇
𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 + 𝑞̇ 𝑑𝑥 𝑑𝑦 𝑑𝑧 − 𝑞𝑥+𝑑𝑥 − 𝑞𝑦+𝑑𝑦 − 𝑞𝑧+𝑑𝑧 = 𝜌𝐶𝑝 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑡

Sustituyendo 𝑞𝑥+𝑑𝑥 , 𝑞𝑦+𝑑𝑦 , 𝑞𝑧+𝑑𝑧

𝜕𝑞𝑥 𝜕𝑞𝑦 𝜕𝑞𝑧 𝜕𝑇

𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 + 𝑞̇ 𝑑𝑥 𝑑𝑦 𝑑𝑧 − 𝑞𝑥 − 𝑑𝑥 − 𝑞𝑦 − 𝑑𝑦 − 𝑞𝑧 − 𝑑𝑧 = 𝜌𝐶𝑝 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡
𝜕𝑞𝑥 𝜕𝑞𝑦 𝜕𝑞𝑧 𝜕𝑇
𝑞̇ 𝑑𝑥 𝑑𝑦 𝑑𝑧 − 𝑑𝑥 − 𝑑𝑦 − 𝑑𝑧 = 𝜌𝐶𝑝 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡

Evaluando la rapidez de calor en términos de la Ley de Fourier, se tiene:

 𝜕𝑇 𝜕𝑇 𝜕𝑇
𝑞𝑥 = −𝑘 𝑑𝑦 𝑑𝑧 , 𝑞𝑦 = −𝑘 𝑑𝑥 𝑑𝑧 , 𝑞𝑧 = −𝑘 𝑑𝑥 𝑑𝑦
𝜕𝑥 𝜕𝑦 𝜕𝑧

Sustituyendo y dividiendo por 𝑑𝑥 𝑑𝑦 𝑑𝑧, obtenemos:

𝜕 𝜕𝑇 𝜕 𝜕𝑇 𝜕 𝜕𝑇 𝜕𝑇
𝑞̇ + �𝑘 � + �𝑘 � + �𝑘 � = 𝜌𝐶𝑝
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑡

Esta última ecuación es la forma general de la ecuación de difusión de calor por conducción en
coordenadas cartesianas. Al resolver la ecuación diferencial se obtiene la distribución de
temperaturas 𝑇(𝑥, 𝑦, 𝑧) como función del tiempo. La ecuación se puede reescribir así:

𝜕 2 𝑇 𝜕 2 𝑇 𝜕 2 𝑇 𝑞̇ 1 𝜕𝑇
+ + + =
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2 𝑘 𝛼 𝜕𝑡

La constante 𝛼 se define entonces como la difusividad térmica, y es propia de cada material,

para este caso en particular se asume que las propiedades del material no varían mucho con la
temperatura. La difusividad térmica se define entonces como:

𝑘 𝑚2
𝛼= [=]
𝜌𝐶𝑝 𝑠

Siguiendo con la expresión para la distribución de temperatura, en estado estable no hay

acumulación de energía, por lo tanto se puede asumir que:

𝜕 2 𝑇 𝜕 2 𝑇 𝜕 2 𝑇 𝑞̇
+ + + =0
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2 𝑘

Si la ecuación anterior además la sometemos a condiciones tales que la transferencia de calor es

unidimensional y no hay generación de energía, el calor deberá ser constante en la dirección de
transferencia:

𝑑 𝜕𝑇
�𝑘 � = 0
𝑑𝑥 𝜕𝑥

𝑑𝑞𝑥′′ 𝑞𝑥
=0 𝑞𝑥′′ =
𝑑𝑥 𝐴
3.1.2 Coordenadas cilíndricas

Al igual que en coordenadas rectangulares, el calor por conducción en un cilindro debe

expresarse en términos de coordenadas cilíndricas:

𝜕𝑇 1 𝜕𝑇 𝜕𝑇
𝑞 ′′ = −𝑘∇𝑇 == −𝑘 �𝑖 +𝑗 +𝑘 �
𝜕𝑟 𝑟 𝜕∅ 𝜕𝑧

𝜕𝑇 1 𝜕𝑇 𝜕𝑇
𝑞′′𝑟 = −𝑘 , 𝑞′′∅ = −𝑘 , 𝑞′′𝑧 = −𝑘
𝜕𝑟 𝑟 𝜕∅ 𝜕𝑧

Tomando un volumen de control tal y como se define en la gráfica siguiente, y realizando los
respectivos cálculos, se llega a la ecuación general de transferencia de calor por conducción en
un sistema cilíndrico:

𝑟
𝑇(𝑟, ∅, 𝑧)

1 𝜕 𝜕𝑇 1 𝜕 𝜕𝑇 𝜕 𝜕𝑇 𝜕𝑇
�𝑘𝑟 � + 2 �𝑘 � + �𝑘 � + 𝑞̇ = 𝜌𝐶𝑝
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕∅ 𝜕∅ 𝜕𝑧 𝜕𝑧 𝜕𝑡

3.1.3 Coordenadas esféricas

En una esfera el flujo de calor esta dado en términos de coordenadas según la siguiente
ecuación:
 𝜕𝑇 1 𝜕𝑇 1 𝜕𝑇
𝑞 ′′ = −𝑘∇𝑇 == −𝑘 �𝑖 +𝑗 +𝑘 �
𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝑠𝑒𝑛𝜃 𝜕∅

𝜕𝑇 1 𝜕𝑇 𝑘 𝜕𝑇
𝑞′′𝑟 = −𝑘 , 𝑞′′𝜃 = −𝑘 , 𝑞′′∅ = −
𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝑠𝑒𝑛𝜃 𝜕∅

Análogamente se toma un volumen de control, se realizan los balances correspondientes y se

llega a la ecuación general de conducción para un sistema esférico:

𝜃 𝑇(𝑟, 𝜃, ∅)
𝑟
𝑦

1 𝜕 𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕 𝜕𝑇 𝜕𝑇
2
�𝑘𝑟 2 � + 2 2
�𝑘 � + 2 2
�𝑘 𝑠𝑒𝑛 𝜃 � + 𝑞̇ = 𝜌𝐶𝑝
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑒𝑛 𝜃 𝜕∅ 𝜕∅ 𝑟 𝑠𝑒𝑛 𝜃 𝜕𝜃 𝜕𝜃 𝜕𝑡

Ejemplo 3.1:

Derive una expresión para 𝑘(𝑥) con 𝐴(𝑥) = (1 − 𝑥)[=]𝑚2 , 𝑇(𝑥) = 300(1 − 2𝑥 − 𝑥 3 ),
𝑞 = 6000 𝑊

𝑞̇ 𝑥
Ejemplo 3.2:

La distribución de temperatura en estado estable en una pared unidimensional con 𝑘 =

1 𝑊/𝑚. 𝐾 y un espesor de 0.3 m, esa dada por la expresión:

𝑇(°𝐶) = 𝑎 + 𝑏𝑥 + 𝑐𝑥 2
𝑎 = 200 °𝐶
°𝐶
𝑏 = −200
𝑚
°𝐶 𝐾
𝑐 = 30 2 = 303 2
𝑚 𝑚
𝑥[=]𝑚

Hallar:
a) La rapidez de generación de calor 𝑞̇ en la pared. Velocidad de 𝑞 hacia adentro y hacia afuera
del área unitaria.
b) Si la superficie fría se expone a un fluido a 100 °C cual es ℎ?