CONDUCION EN REGIMEN PERMANENTE TRIDIMENSIONAL

La transferencia de calor tiene dirección y magnitud. La razón de la transferencia de calor

por conducción en una dirección específica es proporcional al gradiente de temperatura, el
cual es la razón del cambio de la temperatura con respecto a la distancia, en esa dirección.
En general, la conducción de calor en un medio es tridimensional y depende del tiempo, y
la temperatura en un medio varía con la posición y con el tiempo; es decir, T= T (x, y, z, t).
En algunos problemas de transferencia de calor en ingeniería se requiere la determinación
de la distribución de temperatura (la variación de la temperatura) de uno a otro lado del
medio para calcular algunas cantidades de interés, como la razón local de transferencia de
calor, la expansión térmica y el esfuerzo térmico, en algunos lugares críticos en momentos
específicos. La especificación de la temperatura en un punto en un medio requiere en
primer lugar la determinación de la ubicación de ese punto. Esto se puede hacer al elegir un
sistema adecuado de coordenadas, como las rectangulares, cilíndricas o esféricas,
dependiendo de la configuración geométrica que intervenga, y un punto conveniente de
referencia (el origen). La ubicación de un punto se especifica cómo (x, y, z), en
coordenadas rectangulares, como (r, , z), en coordenadas cilíndricas, y como (r,  ), en
coordenadas esféricas, en donde las distancias x, y, z y r, y los ángulos  y  son como se
muestran en la figura. Entonces, la temperatura en un punto (x, y, z) en el instante t, en
coordenadas rectangulares, se expresa como T (x, y, z, t). El mejor sistema de coordenadas
para una configuración geométrica dada es la que describe mejor las superficies en dicha
configuración. Por ejemplo, un paralelepípedo se describe de la mejor manera en
coordenadas rectangulares, ya que cada una de las superficies se puede describir por un
valor constante de las coordenadas x, y o z. Un cilindro es lo más apropiado para las
coordenadas cilíndricas, ya que su superficie lateral se puede describir por un valor
constante del radio. De modo análogo, toda la superficie exterior de un cuerpo esférico se
puede describir del mejor modo por un valor constante del radio en coordenadas esféricas.
Para un cuerpo con forma arbitraria, lo normal es usar coordenadas rectangulares, ya que es
más fácil tratar con distancias que con ángulos. La notación que acaba de describirse
también se usa para identificar las variables que intervienen en un problema de
transferencia de calor. Por ejemplo, la notación T (x, y, z, t) implica que la temperatura
varía con las variables espaciales x, y y z, así como con el tiempo. Por otra parte, la
notación T (x, y, z) indica que la temperatura varía en las tres direcciones y que no se tiene
variación con el tiempo.
Se dice que la conducción en un medio es estacionaria (algunos autores emplean el término
estable o permanente) cuando la temperatura no varía con el tiempo, y no estacionaria o
transitoria, cuando lo hace. Se dice que la conducción de calor en un medio es
unidimensional cuando la transferencia de calor por conducción es significativa sólo en una
dimensión y despreciable en las otras dos direcciones primarias, bidimensional cuando la
conducción en la tercera dimensión es despreciable y tridimensional cuando la conducción
en todas las dimensiones es significativa. En este caso, nos centraremos en tridimensional.
Los problemas de transferencia de calor también se clasifican como unidimensionales,
bidimensionales o tridimensionales, dependiendo de las magnitudes relativas de las razones
de transferencia en las diferentes direcciones y del nivel de exactitud deseado. En el caso
más general la transferencia de calor a través de un medio es tridimensional. Es decir, la
temperatura varía a lo largo de las tres direcciones primarias dentro del medio durante el
proceso de transferencia de calor. En este caso general, la distribución de temperatura de
uno a otro lado del medio en un momento específico, así como la razón de la transferencia
de calor en cualquier ubicación se pueden describir por un conjunto de tres coordenadas,
tales como x, y y z, en el sistema de coordenadas rectangulares (o cartesianas), la r,  y z,
en el sistema de coordenadas cilíndricas, y la r, y , en el sistema de coordenadas esféricas
(o polares). En este caso, la distribución de temperatura se expresa como T (x, y, z, t), T (r,
, z, t) y T (r,  , t), en los respectivos sistemas de coordenadas.
La razón de la transferencia de calor a través de un medio en una dirección específica (por
ejemplo, en la dirección x) es proporcional a la diferencia de temperatura entre uno y otro
lado del medio y al área perpendicular a la dirección de la transferencia de calor, pero es
inversamente proporcional a la distancia en esa dirección. Esto se expresó en forma
diferencial por la ley de Fourier de la conducción del calor en forma unidimensional, como

Donde:
Qx = Tasa de flujo de calor a través del área A en la
dirección positiva.

k = Conductividad térmica  w/ m, k , BT u/ h pie.R 

A = Área de sección transversal de la transferencia de calor

dT/dx = Gradiente de temperatura

Con el fin de obtener una relación general para la ley de Fourier de la conducción del calor,
considere un medio en el cual la distribución de temperatura es tridimensional. El vector de
flujo de calor en un punto P sobre esta superficie debe ser perpendicular a ella y debe
apuntar en la dirección de la temperatura decreciente. Si n es la normal a la superficie
isotérmica en el punto P, la razón de la conducción de calor en ese punto se puede expresar
por la ley de Fourier como:

En coordenadas rectangulares el vector de conducción del calor

se puede expresar en términos de sus componentes como

en donde i →, j →y k →son los vectores unitarios, y Qx, Qy y Qz son las magnitudes de

las razones de transferencia de calor en las direcciones x, y y z, las cuales una vez más se
pueden determinar a partir de la ley de Fourier como:

Aquí, Ax, Ay y Az son las áreas de conducción del calor normales a las direcciones x, y y
z, respectivamente
ECUACION DE CONDUCCION
Una solución de esta ecuación, sujeta a condiciones iniciales y de frontera dados, produce
la distribución de temperatura en un sistema sólido. Una vez que se conoce la distribución
de temperatura, la tasa de transferencia de calor en el modo de conducción se puede evaluar
aplicando la ley de Fourier
La ecuación de conducción es una expresión matemática de la conservación de energía en
una sustancia sólida. Para deducir esta ecuación se efectúa un equilibrio de energía en un
volumen elemental de material en el que el calor se transfiere sólo por conducción. 

Coordenadas Rectangulares
Para ilustrar el enfoque analítico, primero se deducirá la ecuación de conducción para un sistema
coordenado rectangular unidimensional.

Considere un elemento delgado de espesor Δx en una pared plana grande, como se muestra
en la figura

Suponga que la densidad de la pared es r, el calor específico es C y el área de la pared

perpendicular a la dirección de transferencia de calor es A. Un balance de energía sobre
este elemento delgado, durante un pequeño intervalo de tiempo Δt, se puede expresar como

O bien:
Pero el cambio en el contenido de energía interna del elemento y la velocidad de
generación de calor dentro del elemento se pueden expresar como

Al sustituir, se tiene:

Al dividir entre A Δx

Al tomar el límite cuando Δx→0 y Δt→0 se obtiene

por la definición de derivada y a partir de la ley de Fourier de la conducción del calor,

Dado que el área A es constante para una pared plana, la ecuación unidimensional de
conducción de calor en régimen transitorio en una pared de ese tipo queda

Esta ecuación se aplica sólo al flujo de calor unidimensional debido a que se dedujo con base en la
suposición de que la distribución de temperatura es unidimensional. Si ahora se remueve esta
restricción y se supone que la temperatura es una función de las tres coordenadas, así como del
tiempo, es decir, T= T (x, y, z, t), entonces la forma tridimensional de la ecuación de conducción
se convierte en:
Si K es constante:

El operador Laplaciano en coordenadas cartesianas:

Donde α es la difusividad térmica, que es un grupo de propiedades del material definido

como

Si la temperatura de un material no es una función de tiempo, el sistema está en régimen

permanente y no almacena energía. La forma en régimen permanente de una ecuación de
conducción tridimensional en coordenadas rectangulares es
Si el sistema está en régimen permanente y no se genera calor internamente, la ecuación de
conducción se simplifica aún más a:

En coordenadas cilíndricas

Considere ahora un elemento delgado con forma de casco

cilíndrico, de espesor Δr, en un cilindro largo, como se muestra
en la figura 2-15. Suponga que la densidad del cilindro es r, el
calor específico es C y la longitud es L. El área del cilindro,
normal a la dirección de transferencia de calor en cualquier
lugar, es A = 2πrL, en donde r es el valor del radio en ese lugar.
Note que el área A de la transferencia de calor depende de r en
este caso y, por tanto, varía con el lugar. Un balance de energía
sobre este elemento delgado con forma de casco cilíndrico,
durante un pequeño intervalo de tiempo Δt, se puede expresar
como:

O bien,
El cambio en el contenido de energía del elemento y la velocidad de generación de calor
dentro del mismo se pueden expresar como

Al sustituir, se obtiene:

donde A = 2πrL. El lector puede sentirse tentado a expresar el área localizada a la mitad del
elemento, usando el radio promedio como A =2π (r + Δr/2) L. Pero nada hay que se pueda
ganar a partir de esta complicación, ya que, más adelante en el análisis, se tomará el límite
cuando Δr→0 y, por tanto, se cancelará el término Δr/2. Ahora se divide la ecuación
anterior entre AΔr y da

Si se toma el límite cuando Δr→0 y Δt→0, se obtiene

por la definición de derivada y a partir de la ley de Fourier de la conducción del calor,

Puesto que el área de transferencia de calor en este caso es A=2 πrL, la ecuación
unidimensional de conducción de calor en régimen transitorio en un cilindro queda:

 Para un problema tridimensional transitorio general

 Considerar el pequeño elemento cilíndrico de control

Δr, Δz, rΔ, de   densidad y cp  calor específico

Balance de Energía sobre

este elemento durante un pequeño intervalo de tiempo t  se puede expresar como

Reemplazando valores

Siendo el volumen del elemento V  rrz. El contenido de energía en dicho elemento y

la velocidad de generación de calor dentro del mismo se pueden expresar como:

Operando en la ecuación (4) y dividiendo entre r.z.r, se tiene:

Dado que el área de transferencia de calor del elemento para la conducción de ese calor en
las direcciones r,, z son

Ar  r.z , A  r.z, Az  r.r

Tomamos el límite cuando r, z, r y t tiende a cero se obtiene por definición de
derivada y de la Ley de Fourier de la conducción de calor.

Reemplazando

Si la conductividad térmica(k), es constante:

Si el sistema está en régimen permanente y no se genera calor internamente, la ecuación de

conducción se simplifica aún más a

=0
En coordenadas esféricas

Considere ahora una esfera con densidad r, calor

específico C y radio exterior R. El área de la esfera
normal a la dirección de transferencia de calor, en
cualquier lugar, es A=4 πr 2, en donde r es el valor del
radio en ese lugar. Note que, en este caso, el área de
transferencia de calor A, depende de r y, por tanto, varía
con la ubicación. Al considerar un elemento con forma
de casco esférico delgado de espesor Δr y repetir el
procedimiento descrito con anterioridad para el cilindro,
usando A=4 πr 2 en lugar de A=2 π rL, la ecuación
unidimensional de conducción de calor en régimen transitorio para una esfera:

la cual, en el caso de conductividad térmica constante, se reduce a

Para coordenadas esféricas, la temperatura

es una función de las tres coordenadas espaciales:
r = radial   = Polar, cenital o colatitud

V  r 2 sen = azimutal o longitud

A  rr; A  rsenr; Cp : Calor Específico ;   Densidad

 Balance de energía

Reemplazando:

 El contenido de energía en dicho elemento y la velocidad de generación de calor

dentro del mismo se pueden expresar como:

V= Volumen del elemento = V  r senr

2

Reemplazando se tiene:

 Dividiendo entre V, se tiene:

 Tomando límites y reemplazando la ecuación de FOURIER:

 Ordenando, se obtiene la ecuación diferencial de la conducción de calor en

coordenadas esféricas:
Si el sistema está en régimen permanente y no se genera calor internamente, la ecuación de
conducción se simplifica aún más a

=0

METODOS DE SOLUCION DE ´PROBLEMAS

CONDUCCIÓN TRIDIMENSIONAL (MÉTODO NUMÉRICO)
Considerar una cámara cúbica de paredes gruesas y aristas externas 2a, 2b, 2c, la octava
parte se representa en la figura

Nodo interno en un sistema

1) Considerar que la cámara cúbica es la de un horno, conteniendo en sus paredes

resistencias eléctricas calefactoras.

2) Se requiere tres coordinas x, y, z para expresar la distribución de temperaturas en sus

paredes, para tal fin se ha dividido en celdillas cúbicas de dimensiones Δx =Δy =Δz de tal
forma que sus nodos extremos coinciden con sus límites.

3) La aplicación de Balance de energía para el nodo interno (i, j, k) ver la figura

Se tiene Q0  0

Nodo Interno (i, j, k) de una

cámara cúbica (con generación


interna de calor q0  W / m
3
,
en estado estacionario).

4) Sean las celdillas identificadas como puntos nodales

Nodo 0  i, j, k Nodo 4  i, j 1, k 

Nodo 1  i 1, j, k Nodo 5  i, j, k 1

Nodo 2  i 1, j, k  Nodo 6  i, j, k 1

Nodo 3  i, j 1, k 

5) Reemplazando las cantidades de calor en el balance de energía

Simplificando

Si en la cámara no hay generación interna de calor

PROBLEMAS
 Hallar la ecuación nodal para un nodo en un sistema tridimensional ubicado en una
arista exterior, en el cual existe transferencia por convección y se produce
generación interna de calor.

Nodo situado en una

arista convectiva exterior
de un sistema para
evaluar la ecuación nodal
de temperaturas para los
nodos: 0, con generación
interna de calor, en
estado estacionario

SOL.
 Balance de energía para el sistema mostrado en la figura

Reemplazando las cantidades de calor y simplificando se tiene.

Ahora, si: Δx =Δy =Δz

. Si multiplicamos todo por 4/k Δx, obtenemos:

 Obtener la ecuación nodal para un nodo situado en la cara interior de un horno que
se encuentra expuesta a los gases calientes cuya T= Tf de coeficiente h, de
conductividad térmica k, con generación de energía interna. En donde las
amplitudes para este horno en las direcciones x, y, z son Δx, Δy, Δz; iguales.
SOL.

Nodo situado en la cara interna de un

horno hueco, sistema tridimensional

 Balance de energía

 Reemplazando la ecuación de Fourier en cada una de las cantidades de calor por

conducción, el calor convectivo y generado y simplificando se tiene,
REFERENCIAS BIBLIOGRAFICAS

Bohn M., Kreith F. & Manglik M. (2012). Principios de Transferencia de Calor. (7va edic.).
Editorial: Cengage Learning. México.

Cengel Y. (2007). Transferencia de Calor y Masa. (3era edic.). Editorial: McGraw-Hill

Interamericana. México.

Henriquez O. (2009). Transferencia de Calor. Perú: Lima.