CAPITULO 3

CONDUCCIÓN NO PERMANENTE

Diferencias entre conducción permanente y no permanente

Conducción permanente Conducción no permanente

Las temperaturas y el calor son constantes, con el Las temperaturas y el calor no son constantes,
transcurso del tiempo. con el transcurso del tiempo.
Utiliza la ecuación de Laplace, cuya solución es: Utiliza la ecuación de Fourier, cuya solución el:
𝑇 = 𝑓(𝑥, 𝑦, 𝑧) 𝑇 = 𝑓(𝑥, 𝑦, 𝑧, 𝜃)
Se produce cuando existe un equilibrio térmico. Se produce cuando existe una ganancia o pérdida
de calor, donde con el transcurso del tiempo se
llega al equilibrio térmico (calor entregado igual a
calor perdido) o temperatura constante.
Aplicaciones en ingeniería: Proceso del templado,
elaboración de tablas para el temple, fabricación
de vidrios, cerámica vulcanizada, secado de una
tabla de madera, etc.

3.1 Introducción

Si un cuerpo sólido se somete a un cambio repentino en el medio ambiente, debe pasar un cierto
tiempo antes que se establezca en él una condición de temperatura en equilibrio. Nos referimos a la
condición de equilibrio como el estado estacionario o permanente y calculamos la distribución de
temperatura y la transferencia de calor por los métodos descritos en el Capítulo 2.

En el proceso de calentamiento o enfriamiento que tiene lugar en el periodo antes que se establezca
el equilibrio, el análisis debe tomar en cuenta el cambio de la energía interna del cuerpo con el
tiempo, y por tanto se da lugar a la conducción no permanente. La conducción no permanente es
aplicable en tratamientos térmicos como ser: temple, fabricación de vidrio, cerámica, vulcanizado,
entre otros.

3.2 Ecuación de Fourier

𝜕2𝑇 𝜕2𝑇 𝜕2𝑇 1 𝜕𝑇

+ + =
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2 𝛼 𝜕𝜃

Solución: 𝑇 = 𝑓(𝑥, 𝑦, 𝑧, 𝜃)

3.3 Ecuación de Fourier para flujo unidireccional

𝜕2𝑇 1 𝜕𝑇 𝜕2𝑇 1 𝜕𝑇 𝜕2𝑇 1 𝜕𝑇

2
= 2
= 2
=
𝜕𝑥 𝛼 𝜕𝜃 𝜕𝑦 𝛼 𝜕𝜃 𝜕𝑧 𝛼 𝜕𝜃
Aplicando condiciones de contorno, obtenemos la solución particular de estas ecuaciones.

Existen métodos para resolver problemas de conducción no permanente, basado en soluciones

particulares
3.4 Aplicaciones

Los problemas que se estudian en este caso son de dos tipos:

• Cuerpos que evolucionan hacia un equilibrio térmico (calentamientos o enfriamientos).

Ejemplo: el tratamiento térmico de materiales (en particular, el proceso de templado)

• Cuerpos que están sometidos a variaciones periódicas de Temperatura.

Ejemplo: problemas térmicos de regeneradores, máquinas de vapor, motores a explosión, etc.

DISTRIBUCION DE TEMPERATURAS:

Para resolver problemas de conducción no permanente, es necesario conocer cómo se distribuye la

temperatura en un cuerpo. Existen dos tipos de distribución de temperaturas:

Considerando unos planos concéntricos en un objeto, se tiene

1)Temperatura distribuida uniformemente

Todos los puntos del cuerpo tienen la misma temperatura

2)Temperatura estratificada

El número de Biot nos ayuda a determinar si la temperatura está distribuida uniformemente o

estratificada.
3.5 Numero de Biot

El número de Biot es un número adimensional utilizado para cálculos de transferencia de calor en

estado transitorio. Este número relaciona la transferencia de calor por conducción dentro de un
cuerpo y la transferencia de calor por convección en la superficie de dicho cuerpo.
𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑎𝑙 𝑓𝑙𝑢𝑗𝑜 𝑖𝑛𝑡𝑒𝑟𝑛𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑅𝑖
𝐵𝑖 = =
𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑎𝑙 𝑓𝑙𝑢𝑗𝑜 𝑒𝑥𝑡𝑒𝑟𝑛𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑅𝑜
DEDUCCION DEL NUMERO DE BIOT (para un cilindro hueco de pared delgada)

Corriente de un fluido

𝑅𝑖 = 𝐶𝑜𝑛𝑑𝑢𝑐𝑐𝑖𝑜𝑛 𝑎𝑡𝑟𝑎𝑣𝑒𝑠 𝑑𝑒𝑙 𝑒𝑠𝑝𝑒𝑠𝑜𝑟 𝑑𝑒 𝑝𝑎𝑟𝑒𝑑

𝑅𝑜 = 𝐶𝑜𝑛𝑣𝑒𝑐𝑐𝑖𝑜𝑛 𝑎𝑡𝑟𝑎𝑣𝑒𝑠 𝑑𝑒 𝑙𝑎 𝑝𝑒𝑙𝑖𝑐𝑢𝑙𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜, 𝑒𝑛 𝑐𝑜𝑛𝑡𝑎𝑐𝑡𝑜 𝑐𝑜𝑛 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑑𝑒𝑙 𝑐𝑙𝑖𝑛𝑑𝑟𝑜
CONDUCCION:
∆𝑇 𝑇𝑖 − 𝑇𝑜 ∆𝑟
𝑞𝑐𝑜𝑛𝑑 = = 𝑅𝑖 =
𝑅𝑖 𝑅𝑖 𝑘 𝐴𝑙𝑜𝑔

CONVECCION:
𝑇𝑜 − 𝑇∞ 𝑇𝑜 − 𝑇∞ 1
𝑞𝑐𝑜𝑛𝑣 = 𝐴𝑜 ℎ𝑜 (𝑇𝑜 − 𝑇∞ ) = = 𝑅𝑜 =
1 𝑅𝑜 𝐴𝑜 ℎ𝑜
𝐴𝑜 ℎ𝑜
𝑇∞ : 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜
𝐴𝑜 : 𝐴𝑟𝑒𝑎 𝑒𝑥𝑡𝑒𝑟𝑛𝑎
Reemplazando
∆𝑟
𝑘 𝐴𝑙𝑜𝑔 𝐴𝑜 ℎ𝑜 ∆𝑟
𝐵𝑖 = =
1 𝑘 𝐴𝑙𝑜𝑔
𝐴𝑜 ℎ𝑜
𝐴𝑜 ≈ 𝐴𝑙𝑜𝑔
ℎ𝑜 ∆𝑟
𝐵𝑖 = 𝑁𝑢𝑚𝑒𝑟𝑜 𝑑𝑒 𝐵𝑖𝑜𝑡
𝑘
NUMERO DE BIOT PARA CUALQUIER GEOMETRIA
ℎ𝐿
𝐵𝑖 =
𝑘
ℎ ∶ 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛.
𝑘 ∶ 𝐶𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑣𝑖𝑑𝑎𝑑 𝑡é𝑟𝑚𝑖𝑐𝑎 𝑑𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜 𝑠ó𝑙𝑖𝑑𝑜.
𝐿 = 𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑐𝑎𝑟𝑎𝑐𝑡𝑒𝑟í𝑠𝑡𝑖𝑐𝑎.
𝑉
𝐿=
𝐴𝑠
𝑉 ∶ 𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜.

𝐴𝑠 ∶ Á𝑟𝑒𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑎𝑙 𝑑𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜.

3.6 Métodos de solución

Dependiendo del valor calculado del número de Biot, podemos tener dos tipos de problemas, que
son los siguientes:

𝐵𝑖 ≤ 0.1 → 𝐶𝑢𝑒𝑟𝑝𝑜 𝑐𝑜𝑚𝑜 𝑏𝑙𝑜𝑞𝑢𝑒(𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑖𝑠𝑡𝑟𝑖𝑏𝑢𝑖𝑑𝑎 𝑢𝑛𝑖𝑓𝑜𝑟𝑚𝑒𝑚𝑒𝑛𝑡𝑒)

𝐵𝑖 > 0.1 → 𝐶𝑢𝑒𝑟𝑝𝑜 𝑐𝑜𝑛 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑒𝑠𝑡𝑟𝑎𝑡𝑖𝑓𝑖𝑐𝑎𝑑𝑎.
Caso Bi ≤ 0.1

3.6.1 MÉTODO ANÁLISIS DE BLOQUE

Medio
fluido

Como el cuerpo pierde calor, se produce una pérdida del calor sensible o una pérdida de energía
interna.

Desde todos los elementos que se encuentran en la superficie, sale energía térmica del cuerpo; esto
esta expresado mediante una sola flecha para mayor simplicidad.

Aplicando la primera ley de la termodinámica

Calor que sale del cuerpo Disminución de la energía interna del cuerpo
=
en el tiempo dθ en el tiempo dθ
Como la temperatura del cuerpo es uniforme, es decir, independiente de la posición dentro del
cuerpo, solo depende del tiempo. Por tanto, la ecuación puede escribirse de la siguiente manera:

ℎ𝐴𝑠 (𝑇 − 𝑇∞ ) 𝑑𝜃 = −𝜌𝐶𝑃 𝑉 𝑑𝑇
Aplicando condiciones iniciales:𝑇(0) = 𝑇𝑖
𝑇 𝜃
𝑑𝑇 ℎ𝐴𝑠
∫ =∫ − 𝑑𝜃
𝑇𝑖 𝑇 − 𝑇∞ 0 𝜌𝐶𝑃 𝑉

Asumimos ℎ , 𝜌 , 𝐶𝑃 como constantes

𝑇 − 𝑇∞ ℎ𝐴𝑠
𝐼𝑛 ( )=− 𝜃
𝑇𝑖 − 𝑇∞ 𝜌𝐶𝑃 𝑉
𝜌𝐶𝑃 𝑉 𝑇 − 𝑇∞
𝜃=− 𝐼𝑛 [ ] 𝜃 = 𝑓(𝑇)
ℎ𝐴𝑠 𝑇𝑖 − 𝑇∞

ℎ𝐴𝑠
[−( ) 𝜃]
𝑇 = 𝑇∞ + (𝑇𝑖 − 𝑇∞ )𝑒 𝜌𝐶𝑃 𝑉 𝑇 = 𝑓(𝜃)
La caída exponencial de temperatura es similar a la caída de voltaje en la descarga de un
condensador eléctrico, que es:
𝐸
= 𝑒 −𝑡⁄𝑅𝐶
𝐸𝑖
Para completar la analogía, se define la constante térmica de tiempo como:
1
(𝑅𝐶)𝑡𝑒𝑟𝑚𝑖𝑐𝑜 = ( ) (𝜌𝐶𝑃 𝑉) = (𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑡𝑒𝑟𝑚𝑖𝑐𝑎)(𝐶𝑎𝑝𝑎𝑐𝑖𝑡𝑎𝑛𝑐𝑖𝑎 𝑡𝑒𝑟𝑚𝑖𝑐𝑎)
ℎ𝐴𝑠
𝑇 − 𝑇∞ ℎ𝐴𝑠
[−( ) 𝜃]
= 𝑒 𝜌𝐶𝑃 𝑉 = 𝑒 −𝐵𝑖∗𝐹𝑜
𝑇𝑖 − 𝑇∞
Demostración:
ℎ𝐿 𝛼
𝐵𝑖 = 𝐹𝑜 = 𝜃
𝑘 𝐿2
ℎ𝐿 𝛼
𝐵𝑖 ∗ 𝐹𝑜 = ∗ 𝜃
𝑘 𝐿2
ℎ 𝛼
𝐵𝑖 ∗ 𝐹𝑜 = ∗ 𝜃
𝑘 𝐿
ℎ 𝑘 1
𝐵𝑖 ∗ 𝐹𝑜 = ∗ ∗ 𝜃
𝑘 𝜌 𝐶𝑝 𝑉
𝐴𝑠
ℎ𝐴𝑠
𝐵𝑖 ∗ 𝐹𝑜 = 𝜃
𝜌𝐶𝑃 𝑉
𝑇 − 𝑇∞ −𝐿2 𝑇 − 𝑇∞
= 𝑒 −𝐵𝑖∗𝐹𝑜 𝜃= 𝐼𝑛 ( )
𝑇𝑖 − 𝑇∞ 𝛼 ∗ 𝐵𝑖 𝑇𝑖 − 𝑇∞
Caso Bi > 0.1

• Método analítico
• Método analítico - grafico
• Método grafico(incluido con el método de la tabla numerica)

3.6.2 MÉTODO ANALÍTICO

Basado en la solución matemática de la ecuación de Fourier, para flujo unidimensional.

Este método tiene las siguientes restricciones:

1. Solo puede ser aplicado a placas planas de gran superficie, es decir:

𝐴𝑟𝑒𝑎
=∞
𝐸𝑠𝑝𝑒𝑠𝑜𝑟

2. No debe existir resistencia superficial a la transferencia de calor, la resistencia a la

convección debe ser igual a cero.
𝑅𝑐𝑜𝑛𝑣 = 0 → ℎ=∞
Teniendo la siguiente figura:

Donde:

∆𝑥 = 2𝑟𝑚 𝑟𝑚 = 𝐸𝑠𝑝𝑒𝑠𝑜𝑟 𝑚𝑒𝑑𝑖𝑜

𝑇𝑏 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑖𝑛𝑖𝑐𝑖𝑎𝑙 𝑑𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜
𝑇 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑒𝑛 𝑢𝑛 𝑝𝑙𝑎𝑛𝑜 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑑𝑜
𝑇∞ > 𝑇𝑏
𝑇∞ 𝑇∞
Condiciones de contorno:

𝑇 = 𝑇∞ → 𝑥 = 0 𝑥 = 2𝑟𝑚
Flujo de calor por ambas caras

𝑇 = 𝑇𝑏 → 𝜃 = 0
𝑇 = 𝑇∞ → 𝜃 = ∞
Ubicar el plano medio es de gran importancia debido a que el flujo de calor ocurre por ambas caras
de la placa.

Como el flujo de calor es unidireccional:

𝜕2𝑇 1 𝜕𝑇
=
𝜕𝑥 2 𝛼 𝜕𝜃
La solución de la ecuación diferencial es:
𝑇∞ − 𝑇 4 𝜋𝑥 1 −9𝛼 X 3𝜋𝑥 1 −25𝛼 X 5𝜋𝑥
𝑌= = [𝑒 −𝛼1 X sin + 𝑒 1 sin + 𝑒 1 sin +⋯]
𝑇∞ − 𝑇𝑏 𝜋 2𝑟𝑚 3 2𝑟𝑚 5 2𝑟𝑚
 𝜋 2 𝛼
𝛼1 = ( ) X= 𝜃
2 𝑟𝑚 2
𝑥 = 𝑝𝑜𝑠𝑖𝑐𝑖𝑜𝑛 𝑑𝑒𝑙 𝑝𝑙𝑎𝑛𝑜
Por tanto: 𝑇 = 𝑓(𝑥, 𝜃)

La solución de la ecuación diferencial, es una serie de infinitos términos se define como rápidamente
convergente. Es decir, que los primeros términos de la serie son los más significativos y los demás
términos pueden despreciarse, siendo el error mínimo.
𝑇∞ − 𝑇 4 𝜋𝑥 1 −9𝛼 X 3𝜋𝑥
𝑌= = [𝑒 −𝛼1 X sin + 𝑒 1 sin ]
𝑇∞ − 𝑇𝑏 𝜋 2𝑟𝑚 3 2𝑟𝑚
PROBLEMAS QUE SE PRESENTAN:

1) 𝑇 = 𝑓(𝑥, 𝜃) Datos: 𝑥, 𝜃
Se resuelve sin iteraciones

2) 𝑥 = 𝑓(𝑇, 𝜃) Datos: 𝑇, 𝜃
Se resuelve con iteraciones, debido a que x se encuentra dentro de la función sin.

Encontrar 𝑥𝑟𝑒𝑓 con el primer término de la ecuación:

𝑇∞ − 𝑇 4 𝜋𝑥
= [𝑒 −𝛼1 X sin ] 𝑥 = 𝑥𝑟𝑒𝑓
𝑇∞ − 𝑇𝑏 𝜋 2𝑟𝑚
Se hace variar un poco el valor 𝑥𝑟𝑒𝑓
Se calcula 𝑌𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 tomado los dos primeros valores de la ecuación.
4 𝜋𝑥 1 −9𝛼 X 3𝜋𝑥
𝑌𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = [𝑒 −𝛼1 X sin + 𝑒 1 sin ]
𝜋 2𝑟𝑚 3 2𝑟𝑚

Para verificar el resultado comparamos, 𝑌 con 𝑌𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 tomando en cuenta el error

admisible de ±5%

3) 𝜃 = 𝑓(𝑥, 𝑇) Datos: 𝑥, 𝑇
Se resuelve con iteraciones, debido a que 𝜃 es difícil de despejar.
Seguir el mismo procedimiento del problema 2

CALOR GANADO O PERDIDO POR EL CUERPO:

𝑄 8 1 1
= 1 − 2 (𝑒 −𝛼1 X + 𝑒 −9𝛼1 X + 𝑒 −25𝛼1 X + ⋯ )
2𝑟𝑚 𝐴𝜌𝐶𝑃 (𝑇∞ − 𝑇𝑏 ) 𝜋 9 25
𝑄 = 𝑓(𝜃)
3.6.3 MÉTODO ANALÍTICO-GRAFICO

Se basa en las soluciones analíticas de cuerpos geométricos, tales como placas, laminas, esferas,
cilindros. Dichas soluciones analíticas tienen sus expresiones gráficas.

La resolución analítica se facilita si

expresamos las variables en forma
adimensional.

Cambio relativo de temperatura:

𝑇∞ − 𝑇
𝑌=
𝑇∞ − 𝑇𝑏
Tiempo relativo:
𝛼
X= 𝜃
𝑟𝑚 2
Resistencia relativa:
𝑘
𝑚=
𝑟𝑚 ℎ
Posición relativa:
𝑟
𝑛=
𝑟𝑚
Posición relativa para placas planas:
𝑥
𝑛=
𝑥1

Posición relativa para esferas:

Para el centro de la esfera: 𝑟 = 0

Para la superficie externa de la esfera: 𝑟 = 𝑟𝑚

Posición relativa para cilindros:

Para el eje del cilindro: 𝑟 = 0

Para la superficie externa del cilindro: 𝑟 = 𝑟𝑚

CASOS EN LOS QUE ESTE MÉTODO NO ES APLICABLE

1. La línea del grafico intercepta con la línea trazada fuera de los límites de la gráfica.
Solución: Se puede extender el límite superior de la coordenada X manualmente.
2. La línea del grafico intercepta con la línea trazada en la zona de valores negativos de Y.
No existe solución, ya que el rango de valores para Y está comprendido desde cero a valores
positivos.
3. No existe la familia de líneas m para el valor calculado.
Solución: Si el valor de m se encuentra entre dos valores de m existentes en la gráfica, se
puede realizar una interpolación gráfica. Si no existen estos valores, el problema no puede
resolverse a través de este método.
4. El valor de n calculado no coincide o no se aproxima al valor de la gráfica.
Solución: Se puede realizar una interpolación gráfica.

INTERPOLACION GRAFICA

Ejemplo:

Si 𝑚 = 0.8 tomamos la media entre 1 y 0.5:

0.5 + 1
𝑚= = 0.75 ≅ 0.8
2
Entonces trazamos una línea media entre los el valor de 𝑛 que queremos calcular para m=0.5 y m=1

3.6.4 MÉTODO GRAFICO

Dependiendo de las condiciones que se tomen. Se considera un método aproximado o exacto.

Su uso se recomienda cuando los métodos analíticos y analítico-grafico no resuelven un problema

determinado de conducción no permanente. El método gráfico que también se llama método de
Schmidt, se basa en el llamado Principio de Desinverné.

DEDUCCION DE LA EXPRESION DEL PRINCIPIO DE DUSINBERRE

Consideramos una sección de un cuerpo, la misma que se ha enmallado; es decir, se han trazado
líneas verticales y horizontales constituyendo cuadrados de dimensiones ∆x = ∆y, a manera de una
malla. Se considera un flujo de calor bidireccional.

Análisis utilizando elementos finitos

Elegimos el nodo n para establecer un balance de energía:

Calor que se transmite del nodo 1 al nodo n + Calor que se transmite del nodo 2 al nodo n + Calor
que se transmite del nodo 3 al nodo n + Calor que se transmite del nodo 4 al nodo n = Calor que
queda en el nodo n

𝑞1→𝑛 + 𝑞2→𝑛 + 𝑞3→𝑛 + 𝑞4→𝑛 = 𝑞𝑛

Donde, por Fourier:
(𝑇1 − 𝑇𝑛 )
𝑞1→𝑛 = 𝐴𝑥 𝑘𝑥 𝐴𝑥 = ∆𝑦𝐿
∆𝑥
(𝑇2 − 𝑇𝑛 )
𝑞2→𝑛 = 𝐴𝑦 𝑘𝑦 𝐴𝑦 = ∆𝑥𝐿
∆𝑦
(𝑇3 − 𝑇𝑛 )
𝑞3→𝑛 = 𝐴𝑥 𝑘𝑥 𝐴𝑦 = ∆𝑦𝐿
∆𝑥
(𝑇4 − 𝑇𝑛 )
𝑞4→𝑛 = 𝐴𝑦 𝑘𝑦 𝐴𝑦 = ∆𝑥𝐿
∆𝑦
La energía que queda en el nodo n se transforma en calor sensible, y a su vez la temperatura en el
nodo n cambia en un incremento de tiempo definido.
(𝑇𝑛′ − 𝑇𝑛 )
𝑞𝑛 = 𝑚𝐶𝑝
∆𝜃
𝑇𝑛′ = La nueva temperatura del nodo n en un incremento de tiempo ∆𝜃
Conociendo:

𝐴𝑥 = 𝐴𝑦 = 𝐴𝑧 = 𝐴

𝑘𝑥 = 𝑘𝑦 = 𝑘𝑧 = 𝑘

∆𝑥 = ∆𝑦 = ∆𝑧
Reemplazando:
(𝑇1 − 𝑇𝑛 ) (𝑇2 − 𝑇𝑛 ) (𝑇3 − 𝑇𝑛 ) (𝑇4 − 𝑇𝑛 ) (𝑇𝑛′ − 𝑇𝑛 )
𝐴𝑥 𝑘𝑥 + 𝐴𝑦 𝑘𝑦 + 𝐴𝑥 𝑘𝑥 + 𝐴𝑦 𝑘𝑦 = 𝑚𝐶𝑝
∆𝑥 ∆𝑦 ∆𝑥 ∆𝑦 ∆𝜃
(𝑇1 − 𝑇𝑛 ) (𝑇2 − 𝑇𝑛 ) (𝑇3 − 𝑇𝑛 ) (𝑇4 − 𝑇𝑛 ) (𝑇𝑛′ − 𝑇𝑛 )
𝐴𝑘 +𝐴𝑘 +𝐴𝑘 +𝐴𝑘 = 𝑚𝐶𝑝
∆𝑥 ∆𝑥 ∆𝑥 ∆𝑥 ∆𝜃
∆𝑥 (𝑇𝑛′ − 𝑇𝑛 )
𝑇1 + 𝑇2 + 𝑇3 +𝑇4 − 4𝑇𝑛 = 𝑚𝐶𝑝
𝐴𝑘 ∆𝜃
Como:
𝑚 𝑘
𝜌= ∆𝑉 = 𝐴∆𝑥 𝛼=
∆𝑉 𝜌𝐶𝑝
∆𝑥 𝑚 𝐶𝑝 ∆𝑥(𝜌∆𝑉)𝐶𝑝 ∆𝑥 2 ∆𝑥 2 1
= = 𝜌𝐶𝑝 = ( )
𝐴 𝑘 ∆𝜃 ∆𝑉 ∆𝜃𝑘 ∆𝜃 𝛼
( ∆𝑥 ) 𝑘∆𝜃

∆𝒙𝟐
𝑴= → 𝒎𝒐𝒅𝒖𝒍𝒐
𝜶∆𝜽
Reemplazando, tenemos:

𝑇1 + 𝑇2 + 𝑇3 +𝑇4 − 4𝑇𝑛 = 𝑴(𝑇𝑛′ − 𝑇𝑛 )

𝑇1 + 𝑇2 + 𝑇3 +𝑇4 − 4𝑇𝑛 + 𝑴 𝑇𝑛 = 𝑴 𝑇𝑛′
𝑇1 + 𝑇2 + 𝑇3 +𝑇4 + (𝑀−4) 𝑇𝑛 = 𝑀 𝑇𝑛′
𝑇1 + 𝑇2 + 𝑇3 +𝑇4 + (𝑀−4) 𝑇𝑛
𝑇𝑛′ =
𝑀
A esta expresión se conoce como Principio de Dusinberre, en base a esta expresión se ha
desarrollado el Método Gráfico o Método de Schmidt.

MODULO

El valor del módulo M, adopta diferentes valores, según el tipo de problema:

𝑴 = 𝟐 𝒑𝒂𝒓𝒂 𝑪𝒐𝒏𝒅𝒖𝒄𝒄𝒊ó𝒏 𝒖𝒏𝒊𝒅𝒊𝒎𝒆𝒏𝒔𝒊𝒐𝒏𝒂𝒍

𝑴 = 𝟒 𝒑𝒂𝒓𝒂 𝑪𝒐𝒏𝒅𝒖𝒄𝒄𝒊ó𝒏 𝒃𝒊𝒅𝒊𝒎𝒆𝒏𝒔𝒊𝒐𝒏𝒂𝒍

𝑴 = 𝟔 𝒑𝒂𝒓𝒂 𝑪𝒐𝒏𝒅𝒖𝒄𝒄𝒊ó𝒏 𝒕𝒓𝒊𝒅𝒊𝒎𝒆𝒏𝒔𝒊𝒐𝒏𝒂𝒍

Por ejemplo, la expresión obtenida para conducción bidimensional, M = 4. Reemplazando en la

expresión tenemos:

DEDUCCION DE LA EXPRESION DEL PRINCIPIO DE DUSINBERRE (FLUJO UNDIRECCIONAL)

𝑇0 + 𝑇2 + (𝑀−2) 𝑇1
𝑇1′ =
𝑀
𝑀=2
𝑇0 + 𝑇2
𝑇1′ =
2

Como el análisis realizado es para un flujo Bidireccional(M=4):

𝑇1 + 𝑇2 + 𝑇3 +𝑇4
𝑇𝑛′ =
4
Esta expresión se interpreta: La nueva temperatura 𝑇𝑛 ′ que adquiere el nodo n transcurrido un
incremento de tiempo ∆𝜽, es la media aritmética de las temperaturas de los nodos adyacentes al
nodo n.

TIPOS DE PROBLEMAS:

CASO a) Incógnita: temperatura

Datos: Posición del plano, tiempo transcurrido, propiedades del cuerpo, dimensiones, espesor de
división ∆x o ∆r, resistencia superficial

CASO b) Incógnita: posición del plano

Datos: Temperatura del plano, tiempo transcurrido, propiedades del cuerpo, dimensiones, espesor
de división ∆x o ∆r, resistencia superficial

CASO c) Incógnita: tiempo

Datos: Posición del plano, temperatura, propiedades del cuerpo, dimensiones, espesor de división ∆x
o ∆r, resistencia superficial

PROCEDIMIENTO PARA UN FLUJO UNIDIRECCIONAL(M=2)

1) Se divide el cuerpo en partes iguales (número par), cada parte tiene el mismo espesor ∆x,
identificando los planos con números
∆𝒙𝟐
2) Se calcula el incremento de tiempo ∆𝜃, con la expresión ∆𝜃 = 𝑴𝜶

3) Se calcula el número de incrementos de tiempo 𝑁∆θ , de acuerdo a las condiciones del problema
𝜃
𝑁∆θ = , donde ɵ es el tiempo del proceso (calentamiento o enfriamiento)
∆𝜃

4) Se procede a efectuar los trazos, uno por cada incremento de tiempo hasta completar los 𝑁∆θ , en
el CASO a), o hasta llegar en la gráfica a la temperatura establecida como dato, en los CASOS b) y c)

NUMERO DE INCREMENTOS DE TIEMPO

Calentamiento:

A cada grafica le corresponde un valor de tiempo 𝜃

Como el tiempo es una variable continua, al dividir el tiempo

del proceso, en incrementos de tiempo iguales obtenemos
una división de tiempos (significa obtener un numero de
incrementos de tiempo), en la cual a cada incremento de
tiempo le corresponde una gráfica.

GRAFICA 1: 𝜃1 = ∆𝜃 𝑁∆θ = 1

GRAFICA 2: 𝜃2 = 𝜃1 + ∆𝜃 𝑁∆θ = 2
ENFRIAMIENTO:

TRAZO DE LINEAS

TRAZO 1: La primera línea debe cortar por lo menos un plano.

TRAZO 2: El trazo se sobrepone hasta 𝑡1′ , y se traza una línea hasta la el punto más cercano.
TRAZO 3: Repetimos el procedimiento para el trazo 1.

TRAZO 4:

𝑇4′ = 𝑇3 𝑒𝑙 𝑒𝑠𝑝𝑒𝑠𝑜𝑟 ∆𝑥 𝑒𝑠 𝑔𝑟𝑎𝑛𝑑𝑒

Para obtener valores precisos, debemos hacer que ∆𝑥 sea lo más pequeño posible.
METODO GRAFICO PARA UN CILINDRO

METODO GRAFICO PARA UNA ESFERA

CASOS PARTICULARES

a) h = ∞ y las caras externas expuestas a diferentes temperaturas.

Calentamiento

b) h << ∞ ,existe resistencia superficial y las caras externas expuestas a diferentes

temperaturas.
Enfriamiento.

Para ℎ = ∞ , se puede utilizar el método gráfico y el método de la tabla numérica.

Para ℎ ≪ ∞ , solo se puede utilizar el método gráfico. Debido a que este método no toma en cuenta
los espesores ficticios o reales.
3.7 Conducción no permanente en cuerpos semi infinitos.

Un cuerpo semi infinito es aquel cuya dimensión principal no es necesario conocer para resolver un
problema de conducción no permanente, debido a que solo nos interesa una cierta longitud del
cuerpo y no así toda la longitud del cuerpo.

El suelo de la tierra, es un ejemplo de un cuerpo semi infinito, donde la temperatura aumenta con la
profundidad.

Algunos problemas de conducción en estado inestable pueden tratarse aproximadamente

considerando que el cuerpo tiene una temperatura inicial uniforme cuando la temperatura en la
superficie cambia de manera brusca y luego se mantiene en un valor constante distinto a la inicial.
Ejemplos de este tipo de cuerpos son el suelo y una pared de gran espesor; la profundidad del suelo
no se considera, tampoco el espesor respectivamente

Donde

q = calor que gana o pierde el cuerpo semiinfinito

y = posición del plano de interés

Tb = temperatura inicial del cuerpo semiinfinito (varía con el tiempo)

Ts = temperatura de la superficie del cuerpo en contacto con el medio ambiente

𝑇∞ = temperatura del medio ambiente

h = coeficiente de convección (interface superficie-fluido)

En este caso, el calor se transfiere solo en una dirección, por lo que la ecuación de Fourier es:

𝜕2𝑇 1 𝜕𝑇
2
=
𝜕𝑥 𝛼 𝜕𝜃
SOLUCION T=f(y,𝜃)

Para resolver esta ecuación se plantean las siguientes condiciones de contorno:

𝑻 (𝒚 = 𝟎, 𝜽) = 𝑻𝒔 𝒑𝒂𝒓𝒂 𝜽 > 𝟎

𝑻(𝒚 = ∞,𝜽) = 𝑻 𝒑𝒂𝒓𝒂 𝜽 > 𝟎

Condición para el tiempo: 𝑻 (𝒚, 𝜽 = 𝟎) = 𝑻𝒃

CASOS:

a) 𝒉 = ∞ (no existe resistencia superficial)

𝑇𝑠 = 𝑇∞
Solución:
𝑇 − 𝑇𝑠 𝑇 − 𝑇∞ 𝑦
= = 𝑓𝑒𝑟 ( ) 𝑇 = 𝑓(𝑦, 𝜃)
𝑇𝑏 − 𝑇𝑠 𝑇𝑏 − 𝑇∞ √4𝛼𝜃
También se escribe:
𝑦
=𝜉 (𝑠𝑒 𝑙𝑒𝑒 𝑧𝑒𝑡𝑎)
√4𝛼𝜃
Donde la función error de Gauss se define como:
2 𝑢 −𝑛2
𝑓𝑒𝑟(𝜇) = ∫ 𝑒 𝑑𝑛 (∗)
√𝜋 0
La función error, se puede hallar de 3 maneras:
- Analíticamente, a partir de la ecuación (∗)
- De la gráfica fer vs. ξ

Para valores de ξ >2,5 se asume valores de fer=1

- De tablas que existen en los libros de Transferencia de Calor (Ejemplo Ocon y Tojo, Tabla 2-
3, pág. 143)
 b) 𝒉 << ∞ (resistencia superficial importante)
𝑇𝑠 ≠ 𝑇∞

𝑇 − 𝑇𝑏 𝑦 ℎ𝑦 ℎ2 𝛼𝜃 𝑦 ℎ√𝛼𝜃
+
= 1 − 𝑓𝑒𝑟 ( ) − [𝑒 𝑘 𝑘2 ] [1 − 𝑓𝑒𝑟 ( + )] 𝑇 = 𝑓(𝑦, 𝜃)
𝑇∞ − 𝑇𝑏 √4𝛼𝜃 √4𝛼𝜃 𝑘
Planteamiento de problemas
1) Incógnita= T, datos: y, 𝜃, Tb, 𝑇∞, h, 𝛼, 𝑘
2) Incógnita= 𝜃, datos: y, T, Tb, 𝑇∞, h, 𝛼, 𝑘
3) Incógnita= y, datos: 𝜃, T, Tb, 𝑇∞, h, 𝛼, 𝑘

Para ambos casos las soluciones tienen expresiones adimensionales.

Para el planteamiento de problemas 1 y 2 utilizar el primer valor calculado como referencia, para
ambos casos.