CONDUCCION EN ESTADO

NO ESTACIONARIO
 Ley de Fourier
y=Y
y En este capitulo, se considera la
t<0
x variación de la temperatura con el
y=0
T0 tiempo así como con la posición
en sistemas Multidimensionales.

t=0 A diferencia de los procesos de

conducción de calor en estado
T0 T1
estable, en los de tipo transitorio
Q hay un aumento o disminución
en la energía interna del sistema
T (t , y )
mientras ocurre el proceso.
t>0

T0 T1

Q
En diferentes procesos de la
transferencia de calor, la temperatura
del sistema depende del tiempo.
Es el caso del calentamiento y
enfriamiento del techo de una casa
expuesta a la radiación solar; durante
el proceso de templado de vidrio ó de
una pieza de acero; en el proceso de
cocción de productos alimenticios,
etc…..
En todos estos casos la temperatura no
solo esta condicionada por la
distancia, sino también por el tiempo.
ANÁLISIS DE SISTEMAS CONCENTRADOS

En primer lugar, se consideran los sistemas

concentrados, en los que la temperatura varia con el
tiempo pero permanece uniforme a través del sistema
en todo momento.
La cantidad de transferencia de calor alcanza
su limite superior cuando el cuerpo llega a la
temperatura Tα de los alrededores.
Por lo tanto, la transferencia de calor máxima
entre el cuerpo y los alrededores es,

Q = mc p (Ti − Tα )
Consideremos un cuerpo como se muestra en la figura.
El cuerpo está inicialmente a una temperatura Ti.
En el instante t = 0, el cuerpo será colocado en un medio
a una temperatura Tα
Un balance de energía del solido

Ti > Tα
Ein − Eout = Eacum

− Eout = Eacum
Ti > Tα

dT
− As h(T(t ) − Tα ) = ρVc p
dt
 dT
− As h(T(t ) − Tα ) = ρVc p
dt
As h t T dT
−
ρVc p ∫0 dt = ∫Ti (T(t ) − Tα )

T − Tα  hAs 
= exp − ( )t 
Ti − Tα  ρc pV 
La temperatura de un cuerpo concentrado de
forma arbitraria, que se expone a convección
en el instante t =0, en el cual T =Ti, en un
medio a temperatura T∞, con un coeficiente
de transferencia de calor h, se expresa
como,
hAs
 b=( )
T − Tα hAs  ρc pV
= exp − ( )t 
Ti − Tα  ρc pV 

T − Tα
= exp[− bt ]
Ti − Tα
T − Tα
= exp[− bt ]
Ti − Tα
Criterios para el Análisis de Sistemas Concentrados
El primer paso en el establecimiento de un criterio para la
aplicabilidad del análisis de sistemas concentrados es definir la
longitud característica como,

Numero de Biot.

Rconduccion
Bi =
Rconveccion
En otras palabras:
Se aplica el método de sistema concentrado
cuando Bi≤0.1
 EJEMPLO:
 Unos balines de acero de 12 mm de diámetro se
templan mediante calentamiento a 1150ºK, luego se
enfrían lentamente hasta alcanzar los 400ºC en aire
ambiental que se encuentra a 325ºK, con h=20W/m2
ºK. Si las propiedades del acero son k=40W/mºK,
ρ=7800 Kg/m3 y cp=600J/KgºK.
 Estimar el tiempo que se requiere para el
enfriamiento.
Y cuando Bi>0.1 ….?

 Para el caso de una placa infinita:

∂ 2T ∂T
T = T (t , x ) k 2
= ρc p
∂x ∂t

T1 k ∂ 2T ∂T
=
Condiciones frontera: ρc p ∂x 2
∂t
t≤0 T=To -b≤x≤+b
To t=0 T=T1 x=-b
x=+b ∂T ∂ 2T
t>0 T=T(x,t) -b≤x≤+b =α
∂t ∂x 2

x
b
 2
∂T ∂ T
=α 2
∂t ∂y
 Solución Analítica
1 αt
T1 − T ∞
( −1) n −( n + ) 2 π 2 2 πy
= 2∑ e 2 b
cos(n + 1 2)
T1 − T0 n = 0 ( n + 1 2)π b

 Solución por números adimensionales (gráficas)

T − T0 y αt
Y = X = m=
T1 − T0 b b2
Fenómenos de transporte : Bird
 Considerando flujo convectivo en la pared:
 Grafica de Gurney Lurie

T1 − T
Y =
T1 − T0

αt
X =
x12

k
m=
hx1

x
n=
x1
 Una barra rectangular de mantequilla con 46 mm de
espesor y temperatura de 277.6 ºK (4.4 ºC) se extrae de
la nevera y se coloca en un medio ambiente a 297 ºK
(24ºC). (Puede considerarse que los lados de la
mantequilla están aislados) , por tanto, el área expuesta
al medio ambiente son las superficies planas superior e
inferior de la mantequilla. Calcular:
 A) La temperatura de la mantequilla a 15 mm por debajo de
la superficie y en el centro de la misma, después de una
exposición de 1.5 h.
 B) Si coeficiente convectivo tiene un valor de 8.52
W/m2ºK., idem. inciso A) .
 Las propiedades de la mantequilla son k = 0.197 W/m ºK, cp = 2.30
kJ/kgºK y ρ = 998 kg/m3.
Sólidos semi-infinitos
 Se entiende por un cuerpo idealizado
que tiene una sola superficie plana y
se extiende hacia el infinito en todas
las direcciones
 Este cuerpo idealizado se usa para
indicar que el cambio de temperatura
en la parte del cuerpo en la que se
interesa (región cercana a la
superficie) se debe a las condiciones
térmicas en una sola superficie.
Sólidos semi-infinitos
Por ejemplo una pared gruesa se
puede estimar como un medio
semiinfinito si lo único que
interesa es en la variación de la
temperatura en la región cercana a
una de las superficies. En este caso,
la temperatura en la región central
de la pared permanece inalterable.
Sólidos semi-infinitos
Durante periodos cortos, la mayor
parte de los cuerpos pueden
modelarse como sólidos
semiinfinitos, ya que el calor no tiene
el tiempo suficiente para penetrar a la
profundidad del cuerpo y por esta
razón el espesor de un cuerpo no
entra en el análisis de la transferencia
de calor.
Sólidos semi-infinitos
∂T ∂ 2T
=α 2
∂t ∂x
Condiciones frontera:
t≤0 T=To 0≤x≤∞

t=0 T=T1 x=0

t>0 T=T(x,t) 0<x≤ ∞

T − T0 x T − T0 x
= 1 − erf = erfc
T1 − T0 4αt T1 − T0 4αt
T − T0 x k (T1 − T0 )
= erfc q=
T1 − T0 4αt παt
T − T0 x  hx h 2αt   x h αt 
= erfc − exp  + 2 
erfc  + 
T1 − T0 4αt k k   4αt k 
 En invierno el piso se cubre con una capa de nieve a
-10ºC durante un periodo continuo de 3 meses. Se
puede suponer una temperatura inicial uniforme del
suelo de 15ºC y las propiedades promedio del suelo
son k = 0,4 W/mºC y a = 0,15x10-6 m2/s. Determine
la profundidad mínima de entierro para los tubos
que transportan agua, necesaria para impedir
que el agua transportada se congele.
Diferencias finitas: ecuaciones parabólicas
 Las ecuaciones parabólicas se emplean para caracterizar
problemas dependientes del tiempo y el espacio

ECUACIÓN DE CONDUCCIÓN DE CALOR

Caliente Frio
2
∂ T ∂T
α 2 =
∂x ∂t
Diferencias finitas: ecuaciones parabólicas

 Las EDP parabólicas pueden ser resueltas sustituyendo las derivadas

parciales por las diferencias divididas finitas
 Sin embargo, ahora hay que considerar cambios en el tiempo así como
en el espacio
 Mientras las ecuaciones elípticas están acotadas en todas las
dimensiones, las parabólicas están temporalmente abiertas en los
extremos
 Existen dos aproximaciones fundamentales para la solución de EDP
parabólicas:
 Esquema explícito
 Esquema implícito
Métodos explícitos

 La ecuación de conducción de calor requiere aproximaciones para la

segunda derivada en el espacio y para la primera derivada en el tiempo

∂ 2T Ti +l 1 − 2Ti l + Ti −l 1 ∂T Ti l +1 − Ti l
= =
∂x 2 ∆x 2 ∂t ∆t
 Sustituyendo

∂ 2T ∂T Ti +l 1 − 2Ti l + Ti −l 1 Ti l +1 − Ti l
α 2 = α =
∂x ∂t (∆x ) 2
∆t

α∆t
(
Ti l +1 = Ti l + λ Ti +l 1 − 2Ti l + Ti −l 1 ) λ=
(∆x )2
 α∆t
Ti l +1 l
(
= Ti + λ T − 2Ti + T l l l
) λ=
i +1 i −1
(∆x )
2

α∆t 1
= (T l
+ Ti −l 1 )
(∆x )
l +1
2
2 Ti = i +1
2

Ti −l 1
Ti l +1

Ti l Ti +l 1
Convergencia y estabilidad de los métodos
explícitos
 Convergencia: significa que conforme Δx y Δt tienden a cero, los
resultados de la técnica por diferencias finitas se aproximan a la
solución verdadera
 Estabilidad: significa que los errores en cualquier etapa del cálculo no
son amplificados, sino que son atenuados conforme el cálculo avanza
 Se puede demostrar que el método explícito es convergente y estable si
λ < 1/2, o

1 (∆x )
2
∆t ≤
2 k
 Si λ ≤ 1/2  la solución oscila
 Si λ ≤ 1/4  la solución no oscila
 Si λ ≤ 1/6  los errores por truncamiento se minimizan
 Una barra rectangular de mantequilla con 46 mm de
espesor y temperatura de 277.6 ºK (4.4 ºC) se
extrae de la nevera y se coloca en un medio
ambiente a 297 ºK (24ºC). (Puede considerarse que
los lados de la mantequilla están aislados) , por
tanto, el área expuesta al medio ambiente son las
superficies planas superior e inferior de la
mantequilla. Calcular:
 A) Las temperaturas en el interior de la mantequilla
luego de 1.0 h.