Seccin2
Seccin2
Conduccinrgimen
Conduccinrgimen
transitorio
transitorio
1. Introduccin
2. Transferencia de calor por
conduccin
2.1 Introduccin
2.2 Ecuacin general
2.3 Condiciones de frontera
2.4 Conduccin en estado estable
unidimensional sin y con
generacin
2. Continuacin
2.5 Paredes compuestas
2.6 Superficies extendidas (aletas)
2.7 Conduccin en estado estable
bidimensional sin generacin
Mtodo analtico.
2.8 Conduccin en estado transitorio
unidimensional. Mtodo analtico.
2.9 Conduccin en estado transitorio
bidimensional con y sin generacin
2.10 Mtodos numricos
Contenido
1. 1. Entender y resolver problemas relacionados con problemas de
conduccin de calor en rgimen transitorio
2. Aprender a resolver ecuaciones en derivadas parciales por
medio de mtodos numricos. Aplicadas a la resolucin de
fenmenos de conduccin (difusin) en rgimen permanente y
transitorio.
OBJETIVOS
1. [Link]
2. [Link]
3. [Link]
4. wikipedia...
MANUALES DE PROGRAMACIN
En el curso se ha ido incrementando gradualmente la
complejidad de las situaciones estudiadas. Pero no se han
estudiado casos en los que la temperatura vare con el tiempo
(transitorio)
Los problemas de transferencia de calor dependientes del
tiempo son frecuentes. Surgen debido a la variacin de las
condiciones de frontera
Al inicio se considerar que la temperatura dentro del slido
permanece constante T(t), luego la temperatura podr variar
con la posicin y con el tiempo T(x,t) en 1D o T(x,y,t) en 2D
CONDUCCIN TRANSITORIA
Los sistemas concentrados son aquellos en los que la
temperatura interior no vara con la posicin: T(t)==> todo el
sistema se encuentra a la misma temperatura
Bola de cobre vs trozo de carne
Ocurre cuando la resistencia al interior del cuerpo es muy
pequea comparada con la resistencia a la transferencia de
calor entre el sistema y los alrededores.
La transferencia de calor se determina por medio de un
balance de calor global
(Transf. calor desde el cuerpo)=(Increm. energaenel tiempo)
hA(TT
)=m
c
dT
dt
Sist. concentrados/resistencia interna despreciable Sist. concentrados/resistencia interna despreciable
Reemplazando y haciendo un cambio de variable
O tambin
El inverso de la cte en el exp
se conoce como cte de tiempo
trmica t
m=pV
0=TT
-
hA0=pV c
d 0
dt
-
pV c
hA
d 0
0
=dt
pV c
hA
0
i
0
d 0
0
=
0
t
dt - t =
pV c
hA
ln(0/ 0
i
)
0
0
i
=e
t
hA
pV c
Sistemas concentrados Sistemas concentrados
Conociendo la temperatura se puede evaluar el intercambio de
calor con los alrededores y por tanto tambin el cambio en la
energa interna del sistema
La cantidad total de calor transferida entre el sistema y el medio
desde t=0 hasta t=t ser
El calor mximo transferido se alcanzar cuando la temperatura del
cuerpo sea la de los alrededores
Q(t )=hA(T (t )T
) W|
Q(t )=m
c (T(t )T
I
) kJ |
Q
max
=m
c(T
T
I
) kJ |
Sistemas concentrados Sistemas concentrados
Dada la sencillez del mtodo resulta importante establecer
cuando es adecuado aplicarlo
Existe un parmetro que cuantifica cuan exacto es el mtodo.
Para obtenerlo se considera la conduccin en una placa que es
enfriada por conveccin.
kA
(T
s , 1
T
s , 2
)
L
=h A(T
s , 2
T
)
(T
s , 1
T
s , 2
)
(T
s , 2
T
)
=
(L/ kA)
(1/ hA)
=
R
cond
R
conv
(T
s , 1
T
s , 2
)
(T
s , 2
T
)
=
hL
k
Bi
Validez del mtodo Validez del mtodo
Bi es el nmero de Biot. Por conveniencia
se define L=V/As (fcil geom complejas)
Puede ser interpretado como la relacin
entre las resistencias. Si Bi<<1 la resistencia
de conduccin es mucho menor que la de
conveccin. El mtodo puede ser aplicado.
Tambin se puede entender como la relacin entre la diferencia de
temperaturas en el slido, respecto a la de los alrededores
Si Bi<<1 la variacin de la temperatura dentro del cuerpo es
pequea
R
cond
R
conv
=
hL
k
Bi
Validez del mtodo Validez del mtodo
En general se considera adecuado el uso del mtodo de resistencia
despreciable o modelo de sistema concentrado si
Bi<0.1
T(t) en funcin del Bi
Placas: Bi<0.1
Cilindro: Bi<0.05
Esferas: Bi<0.03
Validez del mtodo Validez del mtodo
Reemplazando L=V/As en la ecuacin de se obtiene
Fo es el nmero de Fourier. El cual es un tiempo adimensional
que caracteriza problemas transitorios
hAt
V c
=
ht
c L
=
hL
k
k
c
t
L
2
=
h L
k
t
L
2
hAt
V c
=BiFo , Fo=
t
L
2
Validez del mtodo Validez del mtodo
Bolas de acero al carbn (=7833kg/m, k=54W/mC,
Cp=0.465kJ/kgC y =1.474x10 m/s) de 8 mm de dimetro se
recuecen calentndolas primero hasta 900C en un horno y a
continuacin, dejndolas enfriar con lentitud hasta 100C en aire
ambiente a 35C. Si h=75W/m C, determine cuanto tardar el proceso
de recocido. Si se deben recocer 2500 bolas por hora, determine la
razn total de transferencia de calor de las bolas al aire ambiente.
Ejercicio Ejercicio
Considere una plancha de 1000 W cuya base est hecha de aleacin de
aluminio 2024-T6 (=2770 Kg/m, Cp=875 J/Kg* C, =7.3*10 m/s). La
placa base tiene un rea superficial de 0,03 m. En un principio la placa esta
en equilibrio trmico con el ambiente a 22C. Si el coeficiente de transferencia
de calor por conveccin en la superficie de la placa base es de 12 W/m C y
se supone que el 85% del calor generado en los alambres de resistencia se
transfiere a la placa, determine el tiempo que pasar para que la temperatura
de sta llegue a 140 C. Es realista suponer que la temperatura de la placa
es uniforme en todo momento?
Ejercicio Ejercicio
Pared plana T(x,t)
Considerar: sin generacin, 1D, prop. fsicas ctes, espesor
2L, inicialmente a Ti expuesta al ambiente a T con un h
Adimensionalizando con X=x/L y
2
T
x
+
e
gen
k
=
1
o
T
t
2L
Q
CONV
2
T
x
=
1
o
T
t
, T ( x , 0)=T
i
T (0, t )
x
=0 , k
T ( L, t )
x
=hT ( L, t )T
|
T=f (x , L, t , k , o, h)
0=(TT
)/ (T
i
T
)
Conduccin unidimensional transitoria Conduccin unidimensional transitoria
Reemplazando:
Introduciendo el nmero de Fourier y Bi
0
X
=
0
( x / L)
=
L
T
i
T
T
x
,
2
0
X
2
=
L
2
T
i
T
2
T
x
2
0
t
=
1
T
i
T
T
t
2
0
X
2
=
L
2
o
0
t
,
0(1, t )
X
=
hL
k
0(1, t )
2
0
X
2
=
0
Fo
,
0(1, Fo)
X
=Bi 0(1, Fo)
0(0, Fo)
X
=0 , 0( X , 0)=1
Conduccin unidimensional placa plana Conduccin unidimensional placa plana
Conduccin unidimensional placa plana Conduccin unidimensional placa plana
0( X , Fo)-Temperatura adimensional
X=
x
L
-Distancia adimensional
Bi=
hL
k
-Coeficiente adimensional de transferencia
Fo=
ot
L
2
-Tiempo adimensional
0=f ( X , Bi , Fo)
Conduccin unidimensional placa plana Conduccin unidimensional placa plana
Solucin exacta
aplicando el mtodo de separacin de variables con
imponiendo las condiciones de contorno se obtiene:
0( X , Fo)=F( X)G(Fo)
0=
n=1
C
n
e
\
n
2
Fo
cos(\
n
X)
C
n
=
4sin \
n
2\
n
+sin(2\
n
)
\
n
tan\
n
=Bi
Conduccin unidimensional placa plana Conduccin unidimensional placa plana
Solucin aproximada
si Fo>0.2 la serie se puede aproximar usando slo el primer
termino (error < 2%)
C1 y 1 se obtienen de tablas en funcin del Bi.
En el centro de la pared
Las temperaturas en cualquier posicin estn relacionadas con la
del centro as:
La variacin de la temperatura en cualquier punto respecto al
tiempo es la misma que la del centro de la pared
0=C
1
e
\
1
2
Fo
cos(\
1
X)
0
0
=C
1
e
\
1
2
Fo
0=0
0
cos(\
1
X)
Conduccin unidimensional transitorio Conduccin unidimensional transitorio
Solucin aproximada
para el caso de cilindros y esferas se puede realizar el mismo
procedimiento.
Cilindro:
Esfera:
0
cil
=C
1
e
\
1
2
Fo
Jo(\
1
r / r
0
) , 0
0, cil
=C
1
e
\
1
2
Fo
0
cil
0
0, cil
=Jo(\
1
r / r
0
)
0
esf
=C
1
e
\
1
2
Fo
sin(\
1
r / r
0
)
\
1
r / r
0
|
, 0
0, esf
=C
1
e
\
1
2
Fo
0
esf
0
0, esf
=
sin(\
1
r / r
0
)
\
1
r / r
0
Conduccin unidimensional transitorio Conduccin unidimensional transitorio
Solucin aproximada flujo de calor
Placa plana:
Cilindro:
Esfera:
Q
cil
Q
max , cil
=120
0, cil
J
1
(\
1
)
\
1
Q
esf
Q
max , esf
=130
0, esf
sin(\
1
)\
1
cos(\
1
)
\
1
3
Q
max
=pV c
p
(T
T
I
) kJ |
Q
placa
Q
max , placa
=10
0, placa
sin(\
1
)
\
1
Conduccin unidimensional Conduccin unidimensional
transitorio transitorio
A1=C1
Funciones de Bessel de cero y
primer orden
Conduccin unidimensional transitorio Conduccin unidimensional transitorio
Solucin aproximada de Heisler Solucin aproximada de Heisler
Solucin aproximada
para Fo>0.2 se desarrollaron grficas que permiten aproximar la
temperatura y el flujo de calor en placas, cilindros y esferas en
funcin del Bi y de Fo.
Son grficas difciles de leer para Fo<1
Grficas de Heisler Grficas de Heisler
To: temperatura en el plano central de la placa=T(x=0,t)
x/L
[ ] 1 , 0
L
x
Grficas de Heisler Grficas de Heisler
) (
) (
0
0
7 .
0
=
T T C V Q
Q
Q
Q
t Q
i p
FIG
Heisler cilindros Heisler cilindros
r/r
0
[ ] 1 , 0
0
r
r
To: temperatura en el eje
del cilindro=T(r=0,t)
3 . FIG
i
o
4 . FIG
o i
T T
T T
T T
T T
T T
T ) t , r ( T
Heisler esferas Heisler esferas
5 . FIG
i
o
6 . FIG
o i
T T
T T
T T
T T
T T
T ) t , r ( T
r/r
0
[ ] 1 , 0
0
r
r
To: temperatura en el
centro de la
esfera=T(r=0,t)
Solido Semiinfinito Solido Semiinfinito
Es un cuerpo idealizado que tiene una sola superficie plana y se
extiende haca el infinito en todas las otras direcciones
Se usa para indicar que el cambio de temperatura en la regin cercana
a la superficie se debe slo a las condiciones en la superficie plana
A una cierta profundidad las condiciones en los alrededores no afectan
la temperatura
Se considera transferencia de calor unidimensional perpendicular a la
superficie plana
Solido Semiinfinito Solido Semiinfinito
Solido Semiinfinito Solido Semiinfinito
Aproximacin vlida para intervalos de tiempo cortos
Cambios sbitos de las condiciones trmicas en la superficie
Haciendo un cambio de variables
2
T
x
=
1
o
T
t
, T (0, t )=T
s
T ( x-, t )=T
i
, T ( x , 0)=T
i
p=
x
.
4 ot
-
2
T
p
2
=2p
T
p
T (p=0)=T
s
, T (p -)=T
i
Solido Semiinfinito Solido Semiinfinito
Soluciones
1. Temperatura superficial conocida y constante
2. Flujo de calor superficial constante
3. Conveccin superficial
Donde erf es la funcin de error y erfc=1-erf funcin error complementaria
T ( x , t )T
i
T
s
T
i
=erfc
(
x
2
.
ot
)
,
q
s
(t )=
k(T
s
T
i
)
.
ot
T ( x, t )T
i
=
q
s
k
.
4 ot
exp
(
x
2
4 ot
)
x erfc
(
x
2
.
ot
)
|
T ( x , t )T
i
T
T
i
=erfc
(
x
2
.
ot
)
exp
(
hx
k
+
h
2
ot
k
2
)
erfc
(
x
2
.
ot
+
h
.
ot
k
)
Solido Semiinfinito Solido Semiinfinito
Funcin de error
Funcin de error complementaria
La variable u es una variable ficticia
erf (w)=
2
.
0
w
e
u
2
du
erfc(w)=1erf (w)
Solido Semiinfinito Solido Semiinfinito
Variacin de la temperatura con la posicin y el tiempo
Resumen Resumen
Ejercicio Ejercicio
Una tabla gruesa de madera (k=0.17W/mC y =1.28 x
10 m/s) que esta inicialmente a una temperatura uniforme de
25C se expone a gases calientes a 550C durante un periodo
de 5min. El coeficiente de transferencia de calor entro los gases
y la tabla es de 35W/m C. Si la temperatura de ignicin de la
tabla es de 450C, determine si se encender.
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
Algunas configuraciones pueden ser estudiadas usando los
procedimientos explicados anteriormente y superponindolos
==>solucin producto (multiplicacin de casos 1D)
Deben estar rodeados por el mismo fluido a una T y h constantes
Un cilindro corto (2D), por ejemplo, puede ser considerado como una
interseccin entre un cilindro y una placa plana infinitos que se cruzan
perpendicularmente
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
La temperatura en el cilindro corto sera:
Las soluciones 1D se denotan
(
T (r , x , t )T
T
i
T
)
cil. corto
=
(
T ( x , t )T
T
i
T
)
placa
(
T (r , t )T
T
i
T
)
cil
0
placa
( x , t )=
(
T ( x , t )T
T
i
T
)
placa
0
cil
(r , t )=
(
T (r , t )T
T
i
T
)
cil
0
semiinf
( x , t )=
(
T (x , t )T
T
i
T
)
sol. semiinfinito
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
La transferencia de calor en el cilindro corto sera:
Para un cuerpo 3D es:
(
Q
Q
max
)
total 2D
=
(
Q
Q
max
)
1
+
(
Q
Q
max
)
2
1
(
Q
Q
max
)
1
|
(
Q
Q
max
)
total 2D
=
(
Q
Q
max
)
1
+
(
Q
Q
max
)
2
1
(
Q
Q
max
)
1
|
+
(
Q
Q
max
)
3
1
(
Q
Q
max
)
1
|
1
(
Q
Q
max
)
2
|
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
La temperatura una barra rectangular slida es:
0
barra
( x, y , t )=
(
T ( x , y , t )T
T
i
T
)
barra
=0
placa
( x , t )0
placa
( y , t )
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
La temperatura un prisma slido de medidas 2Hx2Lx2W:
prisma
(x , y , z , t )=
placa
(x , t )
placa
( y , t )
placa
( z , t )
h
L
h
L
h
H
h
H
h
w
h
w
T
2
H
2L
x
y
y
H
0
W
0
z
0 L
x
z
2
W
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
Conduccin en sistemas multidimensionales Conduccin en sistemas multidimensionales
Ejercicio Ejercicio
Un bloque rectangular de hielo (k=2.22W/mC y =0.124x10 m/s)
de 5cm de alto y base cuadrada de 4 cm X 4cm inicialmente a -20C se
coloca sobre una mesa en una habitacin a 18C. El coeficiente de
transferencia de calor sobre las superficies expuestas del bloque de
hielo es de 12 W/m C. Si se descarta toda transferencia de calor de la
base hacia la mesa, determine cunto tiempo transcurrir antes que el
hielo se empiece a fundir. En dnde se formarn las primeras gotas de
lquido?
Seccin2
Seccin2
Mtodosnumricos
Mtodosnumricos
Conduccin transitoria Conduccin transitoria
La disponibilidad de computadores rpidos y software poderosos y
fciles de usar han tenido un alto impacto en la prctica ingenieril
Hoy da los ingenieros tienen que conocer el fenmeno fsico, el
modelo matemtico que lo rige, aplicar correctamente las
herramientas computacionales y ser capaces de interpretar los
resultados
Es un error usar un cdigo sin entender cmo se esta resolviendo el
problema ==> el cdigo no debe ser una caja negra
El know-how adquirido al resolver problemas sencillos es
indispensable
Introduccin Introduccin
Introduccin: ciclo produccin asistido
CAD: computer-assisted design
CAE: computer-assisted engineering
CSM: computational solid mechanics
CFD: computational fluid dynamics
CAA: computational aero-acoustic
CEM: computational electro-magnetics
CAM: computer assisted manufacturing
Introduccin: ciclo produccin asistido Introduccin: ciclo produccin asistido
Introduccin: CFD ventajas
Introduccin: CFD ventajas Introduccin: CFD ventajas
Introduccin: CFD BOEING Introduccin: CFD BOEING
Introduccin: estructura simulacin CFD
Introduccin: estructura simulacin CFD Introduccin: estructura simulacin CFD
Los mtodos numricos convierten las ecuaciones diferenciales en
ecuaciones algebraicas
El dominio es dividido en un conjunto de puntos, volmenes, etc.
La solucin analtica provee una funcin continua, mientras la solucin
numrica es discreta
Introduccin Introduccin
Introduccin Introduccin
Los mtodos numricos son bsicamente de tres tipos: elementos
finitos, diferencias finitas y volmenes finitos.
Elementos finitos generalmente usados en resistencia de materiales y
diseo de mquinas.
En el mtodo de las diferencias finitas las derivadas son reemplazadas
por diferencias.
En volmenes finitos se realiza una integracin considerando que el
valor en el punto es representativo del volumen de control que lo rodea
En diferencias finitas se usa la definicin de la derivada
d f ( x)
d x
=
f ( x+A x)f ( x)
A x
Introduccin Introduccin
