Conducción en Estado Transitorio

Escuela Ingeniería Mecánica
Sección 2
Conducción régimen
transitorio
Contenido
1. Introducción 2. Continuación
2. Transferencia de calor por 2.5 Paredes compuestas
conducción 2.6 Superficies extendidas (aletas)
2.1 Introducción 2.7 Conducción en estado estable
2.2 Ecuación general bidimensional sin generación
Método analítico.
2.3 Condiciones de frontera
2.8 Conducción en estado transitorio
2.4 Conducción en estado estable
unidimensional. Método analítico.
unidimensional sin y con
generación 2.9 Conducción en estado transitorio
bidimensional con y sin generación
2.10 Métodos numéricos
OBJETIVOS
1. Entender y resolver problemas relacionados con problemas de

conducción de calor en régimen transitorio
2. Aprender a resolver ecuaciones en derivadas parciales por

medio de métodos numéricos. Aplicadas a la resolución de
fenómenos de conducción (difusión) en régimen permanente y
transitorio.
MANUALES DE PROGRAMACIÓN
1. http://www.programatium.com/c.htm
2. http://www.cplusplus.com/doc/tutorial/
3.
http://www.cprogramming.com/tutorial.html
4. wikipedia...
CONDUCCIÓN TRANSITORIA
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 varíe con el tiempo
(transitorio)
Los problemas de transferencia de calor dependientes del
tiempo son frecuentes. Surgen debido a la variación de las
condiciones de frontera
Al inicio se considerará que la temperatura dentro del sólido
permanece constante T(t), luego la temperatura podrá variar
con la posición y con el tiempo T(x,t) en 1D o T(x,y,t) en 2D
Sist. concentrados/resistencia interna despreciable
Los sistemas concentrados son aquellos en los que la

temperatura interior no varía con la posición: 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
pequeña 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. energía en el tiempo)
dT
−hA (T −T ∞ )=m c̄
dt
Sistemas concentrados
Reemplazando m= V y haciendo un cambio de variable

=T −T ∞ 
d V c d
−hA = V c  − =dt
dt hA 
 t
− V c d  −V c
∫
hA   0
=∫ dt  t =
hA
ln / i 
i
O también hA
 −t
V c
=e
i
El inverso de la cte en el exp

se conoce como cte de tiempo
térmica τt
Sistemas concentrados
Conociendo la temperatura se puede evaluar el intercambio de calor

con los alrededores y por tanto también el cambio en la energía
interna del sistema
Q̇ t =hA T t −T ∞  [W ]
La cantidad total de calor transferida entre el sistema y el medio
desde t=0 hasta t=t será
Q(t)=m ̄
c (T (t)−T I ) [kJ ]
El calor máximo transferido se alcanzará cuando la temperatura del

cuerpo sea la de los alrededores
Qmax =m ̄c (T ∞ −T I ) [kJ ]
Validez del método
Dada la sencillez del método resulta importante establecer

cuando es adecuado aplicarlo
Existe un parámetro que cuantifica cuan exacto es el método.
Para obtenerlo se considera la conducción en una placa que es
enfriada por convección.
(T s , 1−T s , 2 )
kA =h A (T s , 2−T ∞ )
L
(T s ,1−T s ,2 ) (L/ kA) R cond
= =
(T s ,2 −T ∞ ) (1/hA ) R conv
(T s ,1−T s ,2 ) hL
= ≡ Bi
(T s , 2−T ∞ ) k
Validez del método
R cond hL
= ≡Bi
R conv k
Bi es el número de Biot. Por conveniencia
se define L=V/As (fácil geom complejas)
Puede ser interpretado como la relación

entre las resistencias. Si Bi<<1 la resistencia
de conducción es mucho menor que la de
convección. El método puede ser aplicado.
También se puede entender como la relación entre la diferencia de

temperaturas en el sólido, respecto a la de los alrededores
Si Bi<<1 la variación de la temperatura dentro del cuerpo es pequeña

Validez del método
En general se considera adecuado el uso del método de resistencia

despreciable o modelo de sistema concentrado si
Bi<0.1
Lc: distancia de la máxima
diferencia de temperatura
T(t) en función del Bi

Validez del método
Reemplazando L=V/As en la ecuación de θ se obtiene
hA t ht hL k t h L α t
= = =
ρV c ρ c L k ρ c L 2
k L2
hA t αt
=Bi⋅Fo , Fo= 2
ρV c L
Fo es el número de Fourier. El cual es un tiempo adimensional

que caracteriza problemas transitorios
Ejercicio
Bolas de acero al carbón (ρ=7833kg/m³, k=54W/mºC,

Cp=0.465kJ/kgºC y α=1.474x10⁻⁶ m²/s) de 8 mm de diámetro se
recuecen calentándolas primero hasta 900ºC en un horno y a
continuación, dejándolas enfriar con lentitud hasta 100ºC en aire
ambiente a 35ºC. Si h=75W/m² ºC, determine cuanto tardará el proceso
de recocido. Si se deben recocer 2500 bolas por hora, determine la
razón total de transferencia de calor de las bolas al aire ambiente.
Ejercicio
Considere una plancha de 1000 W cuya base está hecha de aleación 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² y un espesor de 5mm. En un
principio la placa esta en equilibrio térmico con el ambiente a 22ºC. Si el
coeficiente de transferencia de calor por convección 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?
Conducción unidimensional transitoria
Pared plana T(x,t) ∂2 T ė gen 1 ∂ T

 =
∂ x² k  ∂t
Considerar: sin generación, 1D, prop. físicas ctes, espesor
2L, inicialmente a Ti expuesta al ambiente a Tꝏ con un h QCONV
∂2 T 1 ∂ T
= , T  x , 0=T i
∂ x²  ∂ t
∂ T 0, t  ∂ T  L , t
=0 , −k =h[T  L , t −T ∞ ]
∂x ∂x 2L
T =f x , L , t , k ,  , h
=T −T ∞ /T i −T ∞ 
Adimensionalizando con X=x/L y
Conducción unidimensional placa plana
∂ ∂ L ∂T ∂2  L 2 ∂2 T
= = , =
∂ X ∂x / L T i −T ∞ ∂ x ∂ X T i −T ∞ ∂ x 2
2
∂ 1 ∂T
=
∂ t T i−T ∞ ∂ t
Reemplazando:
∂ 2  L2 ∂  ∂ 1, t  hL
= , =  1,t 
∂X 2  ∂t ∂X k
Introduciendo el número de Fourier y Bi

∂2  ∂ ∂  1, Fo
= , =−Bi 1, Fo
∂ X ∂ Fo ∂X
2
∂ 0, Fo
=0 ,   X , 0=1
∂X
  X , Fo  Temperatura adimensional
x
X =  Distancia adimensional
L
hL
Bi= Coeficiente adimensional de transferencia
k
t
Fo= 2  Tiempo adimensional
L
=f  X , Bi , Fo
Solución exacta
aplicando el método de separación de variables con
 X , Fo = F  X ⋅G Fo
imponiendo las condiciones de contorno se obtiene:
∞
− 2n Fo
=∑ C n e cos n X 
n=1
4sin  n
C n=
2  n sin 2  n 
 n tan  n =Bi
Solución aproximada
si Fo>0.2 la serie se puede aproximar usando sólo el primer
termino (error < 2%)
2
− 1 Fo
=C 1 e cos 1 X 
C1 y λ1 se obtienen de tablas en función del Bi.
2
−1 Fo
En el centro de la pared  0 =C 1 e
Las temperaturas en cualquier posición están relacionadas con la
del centro así:
=0 cos 1 X 
La variación de la temperatura en cualquier punto respecto al

tiempo es la misma que la del centro de la pared
Conducción unidimensional transitorio
para el caso de cilindros y esferas se puede realizar el mismo
procedimiento.
Cilindro: 2
− 1 Fo
2
−1 Fo
cil =C 1 e Jo 1 r /r 0  , 0, cil =C 1 e
cil
=Jo 1 r /r 0 
0, cil
Esfera:
2
− 1 Fo
esf =C 1 e
[sin 1 r /r 0 
 1 r /r 0 ] , 0, esf =C 1 e
2
− Fo
1
esf sin  1 r /r 0 
=
0, esf 1 r /r 0
Solución aproximada flujo de calor

Qmax = V c p T ∞ −T I  [kJ ]
Placa plana:
Q placa sin  1 
=1−0, placa
Q max , placa 1
Cilindro: Q cil J 1  1 
=1−2 0, cil
Q max , cil 1
Esfera: Q esf sin 1 − 1 cos 1 

=1−3 0, esf 3
Q max , esf 1
A1=C1
Funciones de Bessel de cero y

primer orden
Solución aproximada de Heisler
para Fo>0.2 se desarrollaron gráficas que permiten aproximar la
temperatura y el flujo de calor en placas, cilindros y esferas en
función del Bi y de Fo.
Son gráficas difíciles de leer para Fo<1
Gráficas de Heisler
x/L
x
  0,1
L
To: temperatura en el plano central de la placa=T(x=0,t)

Gráficas de Heisler
 Q
Q(t )     Q0
 Q0  FIG .7
Q0   V  C p  (Ti  T )
Heisler cilindros
r/r0
r
  0,1
r0
T (r, t )  T  T  T   To  T 
     
Ti  T  To  T  FIG.4  Ti  T  FIG.3
To: temperatura en el eje

del cilindro=T(r=0,t)
Heisler esferas
r/r0
r
  0,1
r0
T( r, t )  T  T  T   To  T 
     
Ti  T T
 o  T  FIG.6
T
 i  T  FIG.5
To: temperatura en el
centro de la
esfera=T(r=0,t)
Solido Semiinfinito
Es un cuerpo idealizado que tiene una sola superficie plana y se

extiende hacía el infinito en todas las otras direcciones
Se usa para indicar que el cambio de temperatura en la región cercana
a la superficie se debe sólo 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
Aproximación válida para intervalos de tiempo cortos

Cambios súbitos de las condiciones térmicas en la superficie
∂2 T 1 ∂ T
= , T 0, t =T s
∂ x²  ∂ t
T  x  ∞ , t=T i , T  x , 0=T i
Haciendo un cambio de variables
2
x ∂ T ∂T
=  2
=−2 
 4  t ∂ ∂
T  =0=T s , T   ∞=T i
Solido Semiinfinito
Soluciones
1. Temperatura superficial conocida y constante
T  x ,t −T i k T s−T i 
T s−T i
=erfc
x
2  t   , q̇ s t =
  t
2. Flujo de calor superficial constante
T  x , t −T i =
q̇ s
k [ 4 t

exp − x2
4 t 
− x erfc
x
2t  ]
3. Convección superficial
T  x ,t −T i
T ∞ −T i
=erfc
x
2  t
−exp
k  
hx h2  t
 2 erfc
k
x
2t

ht
k   
Donde erf es la función de error y erfc=1-erf función error complementaria
Solido Semiinfinito
Función de error
w
2 2
erf w = ∫
 0
e−u
du
Función de error complementaria
erfc w =1−erf w
La variable u es una variable ficticia

Solido Semiinfinito
Variación de la temperatura con la posición y el tiempo

Resumen
Ejercicio
Una tabla gruesa de madera (k=0.17W/mºC y α=1.28 x

10⁻⁷m²/s) que esta inicialmente a una temperatura uniforme de
25ºC se expone a gases calientes a 550ºC 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 ignición de la
tabla es de 450ºC, determine si se encenderá.
Conducción en sistemas multidimensionales
Algunas configuraciones pueden ser estudiadas usando los

procedimientos explicados anteriormente y superponiéndolos
==>solución producto (multiplicación 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
intersección entre un cilindro y una placa plana infinitos que se cruzan
perpendicularmente
La temperatura en el cilindro corto sera:
 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
Las soluciones 1D se denotan
 placa  x ,t =

T  x ,t −T ∞
T i −T ∞  placa
cil r ,t =

T r , t −T ∞
T i −T ∞ cil
 semiinf  x ,t =

T x , t −T ∞
T i −T ∞  sol. semiinfinito
La transferencia de calor en el cilindro corto sera:
Q
( )
Q max total 2D
( ) ( )[ ( )]
=
Q
+
Q
Q max 1 Q max 2
1−
Q
Qmax 1
Para un cuerpo 3D es:
( ) ( ) ( )[ ( )]
Q
Q max total 2D
=
Q
+
Q
Q max 1 Q max 2
1−
Q
Q max 1
+
( Q ) [ 1−( Q ) ][ 1−( Q ) ]
Q
max 3
Q Q
max 1 max 2
La temperatura una barra rectangular sólida es:

barra  x , y ,t =
T  x , y ,t −T ∞
T i −T ∞ 
barra
= placa  x ,t  placa  y ,t 
La temperatura un prisma sólido de medidas 2Hx2Lx2W:
θ prisma (x , y , z ,t )=θ placa (x , t)θ placa ( y ,t )θ placa ( z ,t )
hw
T y
y hw
H hH
x 0
2H
z hH
2W
hL hL
2L
0
x z W
0 L


Ejercicio
Un bloque rectangular de hielo (k=2.22W/mºC y α=0.124x10⁻⁷ m²/s)

de 5cm de alto y base cuadrada de 4 cm X 4cm inicialmente a -20ºC se
coloca sobre una mesa en una habitación a 18ºC. 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 cuánto tiempo transcurrirá antes que el
hielo se empiece a fundir. ¿En dónde se formarán las primeras gotas de
líquido?
Ejercicio
Considere un bloque cúbico cuyos lados tienen 5cm y un bloque

cilíndrico cuya altura y diámetro también son de 5cm. Los dos bloques
se encuentran al principio a 20°C y están hechos de granito
(k=2.5W/m°C, alpha=1.15x10-6 m2/s). Ahora los dos bloques se
exponen en un horno a gases calientes a 500°C sobre todas sus
superficies, con un coeficiente de transferencia de calor de 40W/m2
°C. Determine la temperatura en el centro de la superficie superior de
una de las configuraciones al cabo de 10min.
Conducción transitoria
Sección 2
Métodos numéricos
Introducción
La disponibilidad de computadores rápidos y software poderosos y

fáciles de usar han tenido un alto impacto en la práctica ingenieril
Hoy día los ingenieros tienen que conocer el fenómeno físico, el
modelo matemático que lo rige, aplicar correctamente las
herramientas computacionales y ser capaces de interpretar los
resultados
Es un error usar un código sin entender cómo se esta resolviendo el
problema ==> el código no debe ser una caja negra
El know-how adquirido al resolver problemas sencillos es
indispensable
Introducción: ciclo producción 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

Introducción: CFD ventajas

Introducción: CFD BOEING

Introducción: estructura simulación CFD

Introducción
Los métodos numéricos convierten las ecuaciones diferenciales en

ecuaciones algebraicas
El dominio es dividido en un conjunto de puntos, volúmenes, etc.
La solución analítica provee una función continua, mientras la solución
numérica es discreta
Introducción
Los métodos numéricos son básicamente de tres tipos: elementos

finitos, diferencias finitas y volúmenes finitos.
Elementos finitos generalmente usados en resistencia de materiales y
diseño de máquinas.
En el método de las diferencias finitas las derivadas son reemplazadas
por diferencias.
En volúmenes finitos se realiza una integración considerando que el
valor en el punto es representativo del volumen de control que lo rodea
Introducción
En diferencias finitas se usa la definición de la derivada

d f  x  f  x x −f  x
=
dx x

Conducción en Estado Transitorio

Cargado por

Copyright:

Formatos disponibles

Conducción en Estado Transitorio

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Conducción en Estado Transitorio

Cargado por

Copyright:

Formatos disponibles

Escuela Ingeniería Mecánica

1. Entender y resolver problemas relacionados con problemas de

2. Aprender a resolver ecuaciones en derivadas parciales por

En el curso se ha ido incrementando gradualmente la

Sist. concentrados/resistencia interna despreciable

Los sistemas concentrados son aquellos en los que la

Reemplazando m= V y haciendo un cambio de variable

El inverso de la cte en el exp

Conociendo la temperatura se puede evaluar el intercambio de calor

El calor máximo transferido se alcanzará cuando la temperatura del

Validez del método

Dada la sencillez del método resulta importante establecer

Validez del método

Puede ser interpretado como la relación

También se puede entender como la relación entre la diferencia de

Si Bi<<1 la variación de la temperatura dentro del cuerpo es pequeña

Validez del método

En general se considera adecuado el uso del método de resistencia

T(t) en función del Bi

Validez del método

Reemplazando L=V/As en la ecuación de θ se obtiene

Fo es el número de Fourier. El cual es un tiempo adimensional

Bolas de acero al carbón (ρ=7833kg/m³, k=54W/mºC,

Considere una plancha de 1000 W cuya base está hecha de aleación de

Conducción unidimensional transitoria

Pared plana T(x,t) ∂2 T ė gen 1 ∂ T

Conducción unidimensional placa plana

Introduciendo el número de Fourier y Bi

Conducción unidimensional placa plana

Conducción unidimensional placa plana

Conducción unidimensional placa plana

La variación de la temperatura en cualquier punto respecto al

Conducción unidimensional transitorio

Conducción unidimensional transitorio

Solución aproximada flujo de calor

Esfera: Q esf sin 1 − 1 cos 1 

Conducción unidimensional transitorio

Conducción unidimensional transitorio

Funciones de Bessel de cero y

Solución aproximada de Heisler

To: temperatura en el plano central de la placa=T(x=0,t)

To: temperatura en el eje

Es un cuerpo idealizado que tiene una sola superficie plana y se

Aproximación válida para intervalos de tiempo cortos

Función de error complementaria

erfc w =1−erf w

La variable u es una variable ficticia

Variación de la temperatura con la posición y el tiempo

Una tabla gruesa de madera (k=0.17W/mºC y α=1.28 x

Conducción en sistemas multidimensionales

Algunas configuraciones pueden ser estudiadas usando los

Conducción en sistemas multidimensionales

La temperatura en el cilindro corto sera:

Las soluciones 1D se denotan

Conducción en sistemas multidimensionales

La transferencia de calor en el cilindro corto sera:

Para un cuerpo 3D es:

Conducción en sistemas multidimensionales