Conducción de calor en

estado transitorio
Considere la conducción bidimensional
en estado estacionario es una sección
cuadrada con las temperaturas
mostradas en las superficies.
Establezca la ecuación de diferencias
finitas que corresponde al nodo 1

Nodo 1:
Análisis: las cuatro superficies transfieren calor por Δx
conducción 100

Balance de energía 0

Qizq + Qsup + Qder + Qinf + e1. Δ V = o

T1
k Δy (50-T1)/Δx +k Δx (100-T1) Δy + k Δy (T2-T1)/Δx 50 T2 Δy
+ k Δx (T3-T1) Δy = 0
Δy
Δ x = Δy y dividiendo entre k
50 – T1 +100 -T1 +T2 -T1 + T3 – T1 = 0 Δx
– 4T1 +T2 + T3 = -150 T3
 Objetivos
• Al finalizar este unidad el estudiante debe ser capaz de:
– Visualizar el cambio de temperatura en un pared plana con
conducción unidimensional en estado transitorio
– Distinguir para la conducción unidimensional en estado transitorio
en placas planas, cilindros y esferas el uso de sistema
concentrado, soluciones generales y aproximadas, y sólido
seminifinito
– Resolver problemas prácticos de transferencia de calor por
conducción unidimensional en estado transciente para placas
planas, cilindros y esferas
– Aplicar el método de la superposición para resolver problemas
prácticos de transferencia de calor por conducción para sistemas
multidimensionales en estado transitorio
 ¿Cuáles son los perfiles de temperatura?

Experimento mental transferencia de

t= 0 calor unidimensional por
conducción en una pared expuesta
Ti Ti
a convección en ambas paredes que
Inicialmente
T = Ti se encuentra a una temperatura
h uniforme
T∞ T∞
h Ti t0
Finalmente t=∞
T = T∞ t1
t2
t= 0
O
L x t3

T∞
espesor = 2L
tf
O L
Ti >T∞ x
 . . t o < t1
Q Q
Sistemas concentrados T o < T1

To To To to t1 T1 T1 T1
T1 T1 T1
To To To
Se comportan como un T1 T2 T1
To To To
bulto, su temperatura
interior permanece To T1
espacialmente uniforme
en todo momento T
( C)
0
durante el proceso de
transferencias de calor Tf

método de resistencia
interna despreciable o T1
sistema concentrados

T0
Temperatura es función
del tiempo y no de la
posición
T(t) t0 t1 t
para cualquier posición (min)
en el sólido
 . T∞ t=t
t = t0 = 0 h
Q
T = T∞
T = T∞ m
V Ts = T(t)
T s = Ti AS
T∞ > T i ρ
Cp

Un intervalo de tiempo dt  Un incremento dT en la temperatura del

cuerpo

Transferencia de calor hacia Incremento de la energía

el cuerpo por convección = interna en el cuerpo durante
durante un intervalo dt un intervalo dt

hAs(T∞–T) d t = mdu = mcp dT (1)

Sí,
m = ρV constante (2)

dT = d(T- T∞) (3)

Sustituyendo (2) y (3) en (1) y reordenando
d(T- T∞) hAs
= dt
(T∞-T) ρVcp
d(T- T∞) -hAs (4)
= dt
(T - T∞) ρVcp
Integrando (4)

t = 0, T = Ti (5)
t = t, T = T(t)
(6)
T(t)- T∞ hAs
ln = - t (7)
Ti - T∞ ρVcp constante de
tiempo [=] s-1
hAs
b= - (8)
ρVcp
Sustituyendo (8) en (7) y aplicando exponencial a (7)

T(t) - T∞ =
θ= e –bt (9)
Ti - T ∞
 (Ti - T∞)
T(t) = T∞ + (10)
e bt
 T1 t1

T2 t2

h
h T3 t3
Tf
Tf

Tf t1

x
 Razón de transferencia para la T(t)

.
Q(t) = hAs[T(t)–T∞) (10)

Cantidad total Transferencia de

de energía calor máxima o
intercambiada para límite superior
un valor de Δt

Q = mcp[T(t)–Ti] (11) Q = mcp[T∞–T(t)] (12)

 Criterios para la aplicabilidad
de sistemas concentrados

Longitud Aplicable
característica Número de
del cuerpo Biot
Cuerpos pequeños, V

Bi =hLc/k Conductividad alta, k

Lc = V/As (13) (14)
Bi = (Lc/k)/(1/h) Medio mal conductor, h
conducción Formas geométricas resistencia a resistencia a la
unidimensional simples la conducción convección en la
interna superficie externa

Pared Cilindro Esfera Exacto Aproximado Aplicable

(espesor (r0) (r0)
2L)

Lc = L Lc = r0/2 Lc = r0/3 Bi = 0 Bi > 0 Bi ≤ 0.1

(error
del 5%)
 Supuestos:
• Pared plana infinita
• Transferencia de calor
unidimensional en dirección x
• Todas las superficies inicialmente a
las mismas condiciones térmicas a la
temperatura uniforme Ti

T∞ T∞ • Cuerpo no genera calor

• Espesor de la pared 2L
h h
• En el tiempo t = 0 se baja la
0 L temperatura a T = T∞
• Pared en estado transitorio

La ecuación de balance de energía es:

(18)
Condiciones frontera:
En el centro
de la pared (19)

En la
superficie (20)

Condición inicial:
(21)

(22)

(23)
Distancia
adimensional
X= 0 x=0

X= 1x=L

Temperatura
adimensional (24)

θ(x,t) = 0  T(x, t) = T∞
θ(x,t) = 1  T(x, t) = Ti
 ∂ T/ ∂ x = ∂[(T- T∞)/ ∂x

∂[(T- T∞)/ (Ti -T∞) ](Ti –T∞)

(25)
∂ T/ ∂x =
∂(x/L) L
(26) (Ti –T∞) ∂θ
∂ T/ ∂x =
L ∂(X)
(27)
L
∂θ/∂(X) = ∂T/ ∂x
(Ti –T∞)

Sustituyendo (26) y (27) en (18) y (25) en (20)

(28)

(29)
 Se define un coeficiente adimensional de transferencia de calor o número
de Biot:
Bi = hL/k
(30)
Y un tiempo adimensional o número de Fourier):
τ = αt/L2= Fo (31)

En consecuencia la Ecuación de conducción unidimensional estado transitorio

adimensional toma la forma siguiente;

∂2θ/∂X2 = ∂θ /∂τ (32)

las condiciones frontera adimensionales;

∂ θ(0, τ)/∂X = 0 (33)

∂ θ(1, τ)/∂X = -Biθ(1, τ) (34)

y la condición inicial :

θ(X,0) = 1 (35)

θ = f (X, Bi, Fo) (36)

Se aplica el método de separación de variables para resolver la
ecuación diferencial y se basa en el desarrollo de una función
arbitraria en términos de series de Fourier. Se asume que la variable
dependiente es el producto de varias funciones cada una de ellas
función de una sola variable independiente; con lo cual se transforma
el sistema de ecuaciones diferenciales parciales a diferenciales
ordinarias

Para la pared plana (ver restricciones adicionales en el texto)

entonces:

θ(X, τ) = F(X)G(τ )
(37)
Sustituyendo (37) en (32)

1/F(d2F/dX2) = 1/G(dG/dτ) = -λ2 = cte

(38)
Dando origen a

(39)
 (40)

Cuyas soluciones son :

F = C1 cos (λX) +C2 sen (λX) (41)

G = C3exp(-λ 2τ) (42)

(43)

Aplicando condiciones frontera (33) y (34)

(44)
 (45)

función periódica con un

número infinito de raíces

ecuación
característica

raíces o valores
característicos o
eigenvalores

Entonces hay un infinito número de soluciones de la forma

Por lo tanto, la solución es la combinación lineal de todas ellas:

(46)
An se obtiene a partir de la condición inicial:

(47)
 Se observa
que el aporte
del cuarto
término de la
serie es cero
Términos de la
suma decrecen
rápidamente
conforme
aumenta n

Basta con evaluar algunos de los

primeros términos de la serie infinita (al
menos uno) si el tiempo adimensional,
τ, es grande, para determinar la
temperatura adimensional
Aproximada
a un término

Conserva primer termino

de la serie

Error < 2 %

Soluciones
aproximadas

Tabular

Gráfica
 (47)

(48)

(49)

Recuerde

A1 y λ1 son funciones de Bi
En la tabla 4-2 (Cengel, 2011: 2039) encuentran los
valores A1 y λ1 en función de Bi

Jo es la función de Bessel de primera especie y de

orden cero

En la tabla 4-3 (Cengel, 2011: 2039) encuentran los

valores Jo
Solución aproximada para el centro de las figuras geométricas

(50)

(51)

(52)

La temperatura adimensionales están relacionadas con la

temperatura en el centro de la forma siguiente:

(53)

(54)

(55)
La energía transferida desde la cámara anterior del ojo a
través de la cornea varía considerablemenete dependiendo
de que se usa o no un lente de contacto. Trate el ojo como un
sistema esférico y asuma que se encuentra en estado
estable. El coeficiente de convección no cambia con o sin el
lente de contacto. La cornea y el lente cubre un tercio del
area superficial esférica.

cámara Construya el circuito

anterior
térmico
Rcon,1 = 1/ (4/3πr12)

Resf,1 = (r2 - r1)// (4/3πkr1r2)

Resf,2 = (r3 - r2)// (4/3πkr2r3)

Rcon,2 = 1/ (4/3πr32)
 D = 5 cm

Huevo
T∞ = 95 oC

h = 1200 W/m2.oC

Ti = 5 oC Tf = 70 oC
• Supuestos:
– Conducción en estado transitorio, no depende del tiempo
– Conducción unidimensional en dirección del radio, r,
debido a su simetría térmica con respecto al centro de la
esfera
– La forma del huevo puede ser aproximado a una esfera
con radio igual a 2.5 cm
– Propiedades del huevo y la constante de transferencia de
calor por convección son constantes ( k, α y h cte)
– Dado que el huevo puede considerarse formado
mayoritariamente por agua sus pripiedades seran las del
agua
Análisis: primero se determinará si se puede aplicar la solución de
.sólido concentrado, parea ello se debe cumplir que Bi < 0.1
:Bi viene dado por
Bi = h.Lc /k
Sí h = 1200 W/m2.oC , k = 0.627 W/m.oC y Lc = r/3 = (0.05/2)/3 = 0.0083 m el
:número de Biot sera
Bi = (1200 W/m2.oC x 0.0083 m)/ 0.627 W/m.oC = 15.9 < 0.1 no cumple con el
modelo

En segundo lugar se revisará si cumple con el modelo de aproximación a

.un término, y que es valida para τ > 0.2, donde τ = αt/L2
La ecuación de transferencia de calor unidimensional y transitoria para
:el centro de una esfera es

En este situación debe calcularse nuevamente Bi con la diferencia que se utiliza

r en lugar de Lc, donde r = 0.05/2 m = 0.025 m

Bi = (1200 W/m2.oC x 0.025 m)/ 0.627 W/m.oC = 47.8

Para este valor de Biot se obtienen λ1 y A1 para una esfera en la tabla
:4.2 de Cengel
Bi = 40, A1 = 1.9942 y λ1 = 3.0632
Bi = 50, A1 = 1.9962 y λ1 = 3.0788
Para Bi = 47.8
A1 = (47.8- 40/(50 - 40) x (1.9962 - 1.9942) + 1.9942 = 1.9958

λ1 = (47.8- 40/(50 - 40) x (3.0788 – 3.0632) + 3.0632 = 3.0754

θ0, esf = (70- 95)/(5-95) = 0.277 = 1.9958 exp( - 3.0754 2τ)

ln(0.1392) = -9.4580τ
Τ = 0.209, en consecuencia cumple con la restricción del modelo de
:aproximación a un termino, y podemos pasar a calcular el tiempo, t
t = Τr2/α = (0.209x 0.0252m2)/(0.151 x 10-6 m2/s) = 865 s = 14,4 min

Por lo tanto, el huevo necesitará estar sumergido alrededor de 15 min en

.agua a 95 oC si inicialmente se encontraba a 5 oC
 T∞ = 500 oC
h =120 W/m2.0C
Latón
7 min
Ti = 20 oC
y Ts

x
4 cm
• Supuestos:
– Conducción en estado transitorio, no depende
del tiempo
– Conducción unidimensional en dirección x,
debido a su espesor es pequeño en relación a la
profundidad y altura de la placa y además hay
simetría térmica con respecto al plano central
– Propiedades y la constante de transferencia de
calor por convección son constantes
• A partir de la tabla A-3 de Cengel se
obtuvieron las siguientes propiedades del
latón a 300 K:
– k = 110 W/m2.oC
– ρ = 8530 kg/m3
– Cp= 380 J/kg/oC
– α = 33.9 x 10-6 m2/s
:Análisis
Primero se determinará si se puede aplicar la solución de sólido
.concentrado, para ello se debe cumplir que Bi < 0.1
:Bi viene dado por
Bi = h.Lc /k
Sí h = 120 W/m2.oC , k = 110 W/m.oC y Lc = 2L/2 = (0.04/2) = 0.02 m el número
:de Biot será
Bi = (120 W/m2.oC x 0.02 m)/ 110 W/m.oC = 0.0218 < 0.1 cumple con el modelo

Por lo tanto la temperatura se puede determinar por la siguiente ecuación:

Donde,

b = hAs/ρCp V
As = 2 A (transferencia por ambas caras)
V = espesor x área = 2L . A
b = h(2A)/(ρCp2LA) = h/ρCpL =120 W/m2.oC/ ( 8350kg/m3)(380J/kg.oC)(0.02m)
b = 0.00185 s -1
T (t) = 279 0C

En segundo lugar se puede utilizara el modelo de aproximación a un

,término sí τ > 0.2
τ = α.t/L2 = [(33.9 x 10-6 m2/s).(7min x 60 s/min)]/(0.022m2) = 35.6
.Por lo tanto se puede utilizar este modelo
Si se utilizan las graficas de Heisler debemos realizar el siguiente
: proceso
a) Se obtiene la temperatura adimensional en el centro de la placa
para el tiempo de 7 min
b) Se obtiene la relación entre la temperatura adimensional en la
superficie y la temperatura adimensional en el centro de la placa
para el tiempo de 7 min
c) Encuentre el valor de T en la superficie usando las definiciones de
temperatura adimensional.
Para 1/Bi = 1/0.0218 = 45.8 y τ = 35.6 de la figura 4-16.a se obtiene

Para 1/Bi = 45.8 y x/L = 1/1 =1 de la figura 4-16.b se obtiene

Consecuentemente,
To - T∞ = 0.46 (Ti - T∞)
To = 0.46 (Ti - T∞) +T∞
To = 0.46 Ti +0.54T∞
y
T - T∞ = 0.99(To - T∞)
T = 0.99(To - T∞) + T∞
T =0.99 To - 0.01T∞
Sustituyendo To
T = 0.99 [0.46 Ti + 0.54T∞] - 0.01T∞
T = 0.455 Ti + 0.545T∞
T = 0.455 (20) + 0.545 (500) = 281.6 oC
Se observa que la diferencia es entre los dos modelos es de 2.6 oC

También podemos utilizar las ecuaciones

Es necesario obtener A1 y λ1 de la tabla 4-2 y evaluar la ecuación, dado que

ambos dependen de Bi. Para Bi = 0.0218 hay que interpolar

λ 1 = 0.1461
A1 = 1.0036
Si x/L = 1/1 = 1
θpared = 1.0036x exp(- 0.14612 x 35.6) x cos (1.00362 x 1)
θpared = 0.47

T - T∞ = 0.47(Ti - T∞)
T = 0.47 (Ti - T∞) +T∞
T = 0.47 Ti +0.53T∞
T = 0.47(20) +0.53(500)
T = 274 oC

Se observa que la diferencia es entre los dos modelos es de 5 oC