Download - 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Transcript
Page 1: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

ECUACION DE LA CONDUCCION DE CALOR

Q = Calor

Q=Velocidad de transferenciade calor

QA

=Flujo decalor

Conducción de calos en coordenadas cartesianas es tridimensional y depende del tiempo

T=T ( x , y , z ,t )

La transferencia de calor tiene dirección y magnitud, es cantidad vectorial.

Coordenadas cartesianas (x, y, z)

Coordenadas cilíndricas (r, θ, z)

Coordenadas esféricas (r, φ, θ)

Conducción de calor en estado estable: no hay cambio con el tiempo en cualquier punto dentro del medio.

Conducción de calor en estado transitorio: variación con el tiempo.

CONDUCCION DE CALOR MULTIDIMENSIONAL

Ley de Fourier en la conducción de calor unidimensional.

Qcond=−KA dTdx

El signo negativo garantiza que la transferencia de calor en la dirección positiva de x sea una cantidad positiva.

Distribución de temperaturas es tridimensional el vector de flujo de calor en un punto P sobre la superficie debe ser perpendicular a ella y debe apuntar en la dirección de la temperatura decreciente.

n es la normal a la superficie isotérmica.

Page 2: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Qn=−KA ∂T∂ x

Qn=Q x i+Q y j+Q z k

i , j , k vectoresunitarios

Q x ,Q y ,Q zmagnitudes de velocidadesde transferenciade calor en la direccion x , y , z

Q x=−A xK∂T∂ x,Q y=−A yK

∂T∂ y

,QZ=−A z K∂T∂z

Ax, Ay, Az áreas normales a las direcciones

Se considera materiales de naturaleza isotrópica

GENERACION DE CALOR

Temperatura de una resistencia de alambre se eleva cuando pasa una corriente eléctrica.

Los circuitos electrónicos eliminan calor.

Plantas nucleares de generación eléctrica.

Reacciones químicas exotérmicas, fuentes de calor.

Reacciones químicas endotérmicas, sumidero de calor.

Absorción de la energía solar.

Rayos gamma penetra en el cuerpo.

Page 3: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Generación de calor es un fenómeno volumétrico es decir ocurre en el medio, g(Wm3 )Velocidad de generación de calor G

G=∫V

g dV

Si es uniforme la generación, G= gV

ECUACION UNIDIMENSIONAL DE LA CONDUCCION DE CALOR

Ejemplos pared plana grande, vidrio de una ventana de una sola hoja, la placa metálica de la base de una plancha, un tubo para vapor de agua de fierro fundido, un elemento cilíndrico de combustible nuclear, una resistencia eléctrica de alambre, la pared de un recipiente esférico, una bola metálica que está siendo templada, en todos ellos la conducción de calor es dominante en una dirección y despreciable en las demás.

ECUACION DE LA CONDUCCION DE CALOR EN UNA PARED PLANA GRANDE

Para un elemento delgado de espesor ∆x,

Ρ = densidad

C = calor especifico

A = área de la pared perpendicular a la direccion de transferencia de calor

∆t = pequeño intervalo de tiempo

Page 4: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

[ Velocidaddeconducciondecalor en x

]−[ Velocidadde

conducciondecalor en x+∆ x

]+[ Velocidad degeneraciondecalor enelinterior delelemento

]=[Velocidad decambiodelcontenido deenergiadelelemento

] Q x−Qx+∆x+Gelemento=

∆ Eelemento∆ t

∆ Eelemento=Et+∆t−E t=mc (T t+∆ t−T t )=ρcA∆ x (T t+∆t−T t )

Gelemento= gV elemento= g A∆ x

Q x−Qx+∆x+ g A ∆ x=ρcA∆ x (T t+∆t−T t )

∆ t

Dividiendo entre A∆x

−1A {Q x+∆ x−Q x

∆x }+ g=ρc (T t+∆t−T t∆ t )

Para ∆x→0 ,∆t→0 por definición de derivada y ley de Fourier

lim∆ x→0 (Q x+∆x−Q x

∆ x )❑

=∂Q∂ x

=∂(−KA ∂T∂ x )

∂ x

Page 5: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

1A

∂(KA ∂T∂ x )∂ x

+ g=ρc ∂T∂ t

Como A es constante

∂(K ∂T∂x )∂ x

+ g= ρc ∂T∂ t

La conductividad térmica del material es variable K(T), depende de la temperatura pero se puede suponer que puede permanecer constante en algún valor promedio.

Para conductividad térmica constante

∂2T∂x2 + g

K= 1α∂T∂ t

∝= Kρc

α= Difusividad térmica del material representa la velocidad con la que se propaga el calor a través del mismo.

Sin generación, gK

=0

Estado estable, 1α∂T∂ t

=0

1.- ESTADO ESTABLE, ∂T∂ t

=0

∂2T∂x2 + g

K=0

2.- REGIMEN TRANSITORIO SIN GENERACION DE CALOR, g=0

∂2T∂x2 =1

α∂T∂ t

Page 6: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

3.- ESTADO ESTABLE, SIN GENERACION DE CALOR, ∂T∂ t

=0 y g=0

∂2T∂x2 =0

ECUACION DE LA CONDUCCION DE CALOR EN UN CILINDRO LARGO

Para un espesor ∆r,

Ρ = densidad

C = calor especifico

A = área del cilindro normal a la direccion de transferencia de calor, A = 2πrL

∆t = pequeño intervalo de tiempo

[ Velocidaddeconducciondecalor enr

]−[ Velocidadde

conduccionde calor enr+∆r

]+[ Velocidad degeneraciondecalor enelinteriordelcilindro

]=[Velocidad decambio delcontenidodeenergiadelelemento

]

Page 7: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Qr−Qr+∆r+Gelemento=∆Eelemento∆ t

∆ Eelemento=Et+∆t−E t=mc (T t+∆ t−T t )=ρcA∆ r (T t+∆t−T t )

Gelemento= gV elemento= g A∆ r

Q x−Qx+∆x+ g A ∆ r=ρcA ∆r (T t+∆t−T t )

∆ t

−1A {Qr+∆r−Qr

∆r }+ g=ρc (T t+∆ t−T t∆ t )

Para ∆r→0 ,∆t→0 por definición de derivada y ley de Fourier

lim∆r →0 (Q r+∆r−Q r

∆r )❑

=∂ Q∂ r

=∂(−KA ∂T∂r )

∂ r

1A

∂(KA ∂T∂r )∂r

+ g=ρc ∂T∂ t

A = 2πrL

Para conductividad térmica variable

1r

∂(rK ∂T∂r )∂r

+ g= ρc ∂T∂ t

Para conductividad térmica constante

Page 8: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

1r

∂(r ∂T∂r )∂r

+ gK

= 1∝∂T∂ t

1.- ESTADO ESTABLE, ∂T∂ t

=0

1r

∂(r ∂T∂r )∂r

+ gK

=0

2.- ESTADO TRANSITORIO SIN GENERACION DE CALOR, g=0

1r

∂(r ∂T∂r )∂r

= 1∝∂T∂ t

3.- ESTADO ESTABLE, SIN GENERACION DE CALOR, ∂T∂ t

=0 y g=0

∂(r ∂T∂r )∂ r

=0

ECUACION DE LA CONDUCCION DE CALOR EN UNA ESFERA

Para un casquete esférico de espesor ∆r,

ρ = densidad

C = calor especifico

A = área del casquete normal a la dirección de transferencia de calor, A = 4πr2

∆t = pequeño intervalo de tiempo

Page 9: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Para conductividad térmica variable

1

r2

∂(r2K∂T∂r )

∂r+ g=ρc ∂T

∂ t

Para conductividad térmica constante

1

r2

∂(r2 ∂T∂r )

∂r+ gK

= 1∝∂T∂ t

1.- ESTADO ESTABLE, ∂T∂ t

=0

1

r2

∂(r2 ∂T∂r )

∂r+ gK

=0

2.- REGIMEN TRANSITORIO SIN GENERACION, g=0

1

r2

∂(r2 ∂T∂r )

∂r= 1∝∂T∂ t

3.- ESTADO ESTABLE, SIN GENERACION DE CALOR, ∂T∂ t

=0 y g=0

Page 10: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

∂(r2 ∂T∂r )∂r

=0

ECUACION UNIDIMENSIONAL COMBINADA DE LA CONDUCCION DE CALOR

1

rn

∂(r nK ∂T∂r )∂r

+g=ρc ∂T∂ t

Pared plana n=0 ∂(K ∂T∂x )∂ x

+ g= ρc ∂T∂ t

Para cilindro n=1 1r

∂(rK ∂T∂r )∂r

+ g= ρc ∂T∂ t

Para esfera n=2 1

r2

∂(r2K∂T∂r )

∂r+ g=ρc ∂T

∂ t

ECUACION GENERAL DE CONDUCCION DE CALOR EN COORDENADAS RECTANGULARES

Page 11: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

[ Velocidaddeconduccionde calor enx , y , z

]−[Velocidad

deconducciondecalor en

x+∆ x , y+∆ y ,z+∆ z

]+[ Velocidad degeneraciondecalor enelinterior delsistema

]=[Velocidad decambiodelcontenido deenergia delelemento

]Q x+Q y+Q z−Q x+∆x−Q y+∆ y−Qz+∆ z+Gelemento=

∆ Eelemento∆ t

V = ∆x∆y∆z

∆ Eelemento=Et+∆t−E t=mc (T t+∆ t−T t )=ρc ∆x ∆ y ∆ z (T t+∆ t−T t )

Gelemento= gV elemento= g ∆ x ∆ y ∆ z

Q x+Q y+Q z−Q x+∆x−Q y+∆ y−Qz+∆ z+ g∆ x ∆ y ∆ z=ρc∆ x ∆ y ∆ z∆T∆ t

Dividir entre ∆x∆y∆z

−1∆ y ∆ z {Q x+∆ x−Q x

∆x }− 1∆ x ∆z {Q y+∆ y−Q y

∆ y }− 1∆ x ∆ y {Q z+∆ z−Q z

∆ z }+ g=ρc (T t+∆t−T t∆ t )

Ax= ∆y*∆z, Ay= ∆x*∆z, Az= ∆x*∆y

∆t, ∆x, ∆y, ∆z → 0

Page 12: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

− lim∆ x→0

1

∆ y ∆ z (Q x+∆ x−Q x

∆ x )−lim∆ y→0

1

∆ x ∆ z ( Q y+∆ y−Q y

∆ y )−lim∆ z→0

1

∆ x ∆ y ( Q z+∆ z−Q z

∆ z )+ g=ρc lim∆t →0 (T t+∆t−T t∆t )

lim∆x→0

1

∆ y ∆ z ( Q x+∆ x−Q x

∆x )= 1∆ y ∆ z

∂Q∂ x

= 1∆ y ∆ z

∂(−K ∆ y ∆z ∂T∂x )∂ x

=−∂(K ∂T∂ x )

∂x

lim∆ y→0

1

∆ x ∆ z ( Q y+∆ y−Q y

∆ y )=−∂(K ∂T∂ y )

∂ y

lim∆z→0

1

∆ x ∆ y ( Q z+∆z−Qz

∆ z )=−∂(K ∂T∂ z )

∂ z

Para K constante:

ECUACION DE FOURIER-BIOT

∂2T∂x2 + ∂

2T∂ y2 +

∂2T∂ z2 + g

K= 1∝∂T∂ t

ECUACION DE POISSON (ESTADO ESTABLE)

∂2T∂x2 + ∂

2T∂ y2 +

∂2T∂ z2 + g

K=0

ECUACION DE DIFUSION (REGIMEN TRANSITORIO, SIN GENERACION DE CALOR)

∂2T∂x2 + ∂

2T∂ y2 +

∂2T∂ z2 = 1

∝∂T∂ t

Page 13: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

ECUACION DE LAPLACE (ESTADO ESTABLE Y SIN GENERACION DE CALOR)

∂2T∂x2 + ∂

2T∂ y2 +

∂2T∂ z2 =0

ECUACION GENERAL DE CONDUCCION DE CALOR EN COORDENADAS CILINDRICAS

X = z

Y = rsenφ

X = r cosφ

X = z

Y = rsenφ

X = r cosφ

1r

∂(rK ∂T∂r )∂r

+

˙

1

r2

∂(rK ∂T∂∅ )∂∅

+∂(K ∂T∂ z )∂ z

+ g=ρc∂T∂ t

Page 14: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

ECUACION GENERAL DE CONDUCCION DE CALOR EN COORDENADAS ESFERICAS

1

r2

∂(r2K∂T∂r )

∂r+

˙

1

r2 sen2θ

∂(K ∂T∂φ )∂φ

+1

r2 sin θ

∂(Ksenθ ∂T∂θ )∂θ

+g=ρc∂T∂ t

CONDICIONES DE FRONTERA E INICIALES

1.- Condición de frontera de temperatura especifica

T(0,t) = T1 en x = 0

T(L,t) = T2 en x = L

Page 15: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

2.- Condición de frontera de flujo especifico de calor

Cuando se puede determinar la velocidad de transferencia de calor q, el flujo de calor sobre la superficie.

q=−K ∂T∂x

Flujo de calor en la dirección positiva de x (W/m2)

2a.- Caso especial: Frontera aislada

El aislamiento reduce la transferencia de calor pero no lo elimina en su totalidad, la transferencia de calor se puede tomar como cero

q=K∂T (0 , t )∂ x

=0

La función de la temperatura debe ser perpendicular a la superficie aislada.

Page 16: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

2b.- Caso especial: Simetría térmica

Pendiente cero

q=K∂T ( L2 , t)∂ x

=0

3.- Condiciones de convección de frontera

Balance de energía superficial

(Conduccion decalor enla superficie enunadireccion seleccionada)=(Conveccion decalor enla superficieen la

mismadireccion )

−K∂ ( 0 ,t )∂ x

=h1 (T ∞1−T (0 ,t ) )

Page 17: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

−K∂ (L , t )∂x

=h2 (T (L, t )−T ∞2 )

4.- Condición de radiación de frontera

En el vacío no hay convección, radiación sería el único mecanismo

−K∂T (0 , t )

∂ x=ε1σ (T alred , 14 −T (0 ,t )

4 )

−K∂T (L, t )

∂x=ε2σ (T (L,t )

4 −T alred ,24 )

σ = 5.67*10-8 W/m2-K4

T alred,1, Talred,2, K, ºR

5.- Condiciones de frontera en la interfase

Page 18: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Cuerpos formados por capas de materiales diferentes

a)Los dos cuerpos en contacto tienen la misma temperatura en el área de contacto

b)Una interfase no puede almacenar energía, flujo de calor sobre ambos lados de la interfase debe ser el mismo.

T A (x0 ,t )=T B (x0 ,t )

−K A

∂T A ( x0 , t )

∂ x=−K B

∂T B( x 0 ,t )

∂x

6.- Condiciones de frontera generalizadas

Una superficie puede comprender convección, radiación y flujo superficial de calor simultáneamente

( Transferencia decalor haciala superficie , en todos losmodos)=(Transferencia decalor desdelasuperficie en todos losmodos )

Page 19: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

GENERACION DE CALOR EN UN SOLIDO

La temperatura de un medio se eleva durante la generación de calor esto hasta que alcanzan las condiciones de operación estables y la velocidad de generación de calor es igual a la velocidad de la transferencia de calor a los alrededores.

La temperatura máxima ocurre en el plano medio de una pared plana, en la línea central de un cilindro hueco y en el punto medio en una esfera.

(Velocidad de la transferenciade calor desde el solido )=( Velocidadde lageneraciondeenerg iadentrodel solido )

Page 20: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Q= gV

Q=h A s (T s−T ∞ )

T S=T ∞+gVh A s

Q=Egen

Para las diversas geometrías:

Pared plana 2L, As = Apared, V =2L Apared

Cilindro ro, AS=2 πr oL ,V=π r02 L

Esfera r0, AS=4 πr 02 ,V= 4

3π r0

3

Temperaturas de superficie para las geometrías:

T s , pared plana=T ∞+gLh

T s , cilindro=T ∞+g r 0

2h

T s , esfera=T ∞+g r0

3h

Calor generado a través de un casco cilíndrico es igual al calor generado dentro de el

−K A rdTdr

= gV r

Page 21: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Ar=2πrL ,V r=πr2L

−K (2πrL ) dTdr

= g π r2L

dT=− g r2K

dr

∫T 0

T S

dT=−∫r=0

r=r0

g r2K

dr→T S−T 0=− g4K

(r 02−02 )

Para cilindro

T 0−T s=g r0

2

4 K=∆T max

T0 = Tcentro

T centro=T S+∆T max

Para placa plana

∆T max , pared=g L2

2K

Para esfera solida

∆T max , pared=g r 2

6K

Page 22: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

CONDUCTIVIDAD TERMICA VARIABLE K(T)

Cuando se conoce la variación de la conductividad térmica con la temperatura K(T), se puede determinar el valor promedio de la conductividad térmica en el rango de temperaturas entre T1 y T2 a partir de:

K prom=∫T1

T2

K (T )dT

T 2−T1

Si K(T) = K → Kprom = K

Si K(T) = Ko*( 1+βT)

K prom=∫T1

T2

K0 (1+βT )dT

T 2−T 1

=K0[1+β {T2+T 1

2 }]

Page 23: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

PROBLEMA

La distribución de temperaturas a través de una pared de 1m de espesor en cierto instante está

dada como T x=a+bx+c x2 donde T está dada en grados Celsius y x en metros, mientras

a=900°C, b= -300°C/m y c = -50°C/m2

Una generación de calor uniforme g = 1000W/m3, está presente en la pared de área 10m2 que tiene las propiedades ρ = 1600 Kg/m3, K = 40 W/m-K, Cp = 4 Kj/Kg-K

a)Determine la rapidez de transferencia de calor que entra a la pared x = 0 y sale de la pared.

b)Determine la rapidez de cambio de almacenamiento de energía en la pared.

c)Determine la rapidez con respecto al tiempo del cambio de temperatura en x = 0, x=0.25 y x=0.50m.

T x=900−300 x−50 x2→dTdx

=−300−100 x

a)Qentra=Q x=0=−KA dTdx

=−40∗10 (−300−100∗0 )=120000W=120KW

Qsale=Q x=L=−KA dTdx

=−40∗10 (−300−100∗1 )=160000W=160KW

Page 24: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

b) Balance de energía

Eentra+ Egenera−Esale=Ealmacenamiento

E genera = g*A*L = 1000*10*1 = 10000W = 10 KW

E almacenamiento = 120+10-160 = -30 KW

c)De ecuación general de conducción

∂2T∂x2 + g

K=ρCp

K∂T∂ t

∂T∂ t

= KρCp

∂2T∂ x2 + g

ρC p

∂2T

∂x2=∂( ∂ y∂x )∂x

=∂ (−300−100 x )

∂ x=−100

C

m2

∂T∂ t

= 401600∗4000

∗(−100 )+ 10001600∗4000

=−4.69∗10−4 Cs

Hallar la ecuación de Transferencia de calor por conducción entre dos caras paralelas, sin generación de calor y estado estacionario.

Para la dirección x

∂2T∂x2 =0

Como hay una sola variable, integrando: dTdx

=C1

Volviendo a integrar: T=C1 x+C2

Hallar C1 y C2 con las condiciones de contorno X = 0 → T = T1 y X = L → T = T2

T 1=C1∗0+C2→T1=C2

T 2=C1∗L+T1→C1=T2−T 1

L

Page 25: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Por lo tanto:

T x=(T 2−T 1

L ) x+T 1

dTdx

=T 2−T 1

L

Q=−KA dTdx

=−KA(T 2−T 1

L )

Q=−T 2−T 1

LKA

=T 1−T2

LKA

Hallar la ecuación de Transferencia de calor por conducción en una pared cilíndrica, sin generación de calor y estado estacionario.

Para la dirección radial

1r

∂(r ∂T∂r )∂r

=0→∂(r ∂T∂r )∂ r

=0

Integrando r∂T∂r

=C1→dTdr

=C1

r

Integrando nuevamente T=C1 lnr+C2

Para r = r1 → T = T1 y además r = r2 → T = T2

T 1=C1 ln r1+C2

T 2=C1 ln r2+C2 , ahora restamos

Page 26: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

T 1−T2=C1 (ln r1−ln r2 )=C1 lnr1

r2

→C1=T 1−T 2

lnr1

r2

Despejando C2

C2=T1−C1 ln r1=T 1−[T1−T2

lnr1

r2] ln r1

T=[T1−T2

lnr1

r2] ln r+¿T1−[T 1−T 2

lnr1

r2] lnr 1¿

Por lo tanto la ecuación de distribución de temperaturas en un cilindro es:

T r=T 1+[ T1−T 2

lnr1

r 2] ln

rr1

Calculo de la transferencia de calor:

Q=−K A r=r1

dTdrparar=r1

dTdr

=C1

r=

T1−T 2

lnr 1

r 2

r1

Q=−K2 π r1L

T1−T2

lnr1

r2

r1

=−2πKLT 1−T2

lnr1

r2

=T1−T2

lnr2

r1

2πKL

Hallar la ecuación de Transferencia de calor por conducción en una pared esférica, sin generación de calor y estado estacionario.

Page 27: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

1

r2

∂(r2 ∂T∂r )

∂r=0

Integrando r2 ∂T∂r

=C1→∂T∂r

=C1

r 2

Integrando nuevamente T=C1

r+C2

Para r = r1 → T = T1 y además r = r2 → T = T2

T 1=C1

r1

+C2

T 2=C1

r2

+C2 y restando ambas

T 1−T2=C1( 1r1

− 1r2 )→C1=

T1−T2

( 1r1

− 1r2

)

T 1=

T 1−T 2

( 1r1

− 1r2

)r1

+C2→C2=T1+T 1−T2

r1( 1r2

−1r1

)

T=

T1−T 2

( 1r1

− 1r2

)r

+T 1+T1−T 2

r1( 1r 2

−1r1

)=T 1−(T 1−T 2)

[1− r1

r ][1− r1

r2]

Calculo de la transferencia de calor:

Page 28: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Q=−KA dTdr

=−K 4 π r2 dTdr→Q

dr

r2=−4 πKdT

Integrando:

Q∫r1

r2

dr

r2=−∫

T 1

T 2

4 πKdT

−Q [ 1r2

− 1r1 ]=−4 πK (T2−T1 )→Q=

T 1−T 2

( 1r1

− 1r 2

) 14 πK

Q=T 1−T2

(r 2−r1 )4 πK r1 r2

PROBLEMA

Una pared homogénea de área “A” y espesor “L” tiene temperaturas en la superficie izquierda y derecha de 0°C y 40°C, respectivamente. Determinar las temperaturas en el centro de la pared. ¿Cuánto material debe agregarse y en cuál de los lados de la pared para que la temperatura en el centro de esta ascienda 4°C? ¿Cuánto material debe quitarse y en cuál de los lados para para que la temperatura en el centro de la pared baje 4°C? Exprese sus respuestas en términos de A y L

d2Td x2 =0

dTdx

=C1

T=C1 x+C2(1)

Page 29: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Para x = 0 → T = 0°C

en (1) 0=C1∗0+C2→C2=0

Para X = L → T = 40°C

en (1) 40=C1L+0→C1=40Lx

en (1) T=40L

Para x = L/2 en la mitad T=

40L

∗L

2→T=20℃

a) Que material debe agregarse para que el centro ascienda 4°C? (24°C)

Se toma el (a)

−K (40−0 )L+x

=−K ( 40−24 )

L2

L+x=

4016

∗L

2→x=

L4

b) Que material debe quitarse para que el centro descienda 4°C? (16°C)

Page 30: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Se toma el caso (a)

−K (40−0 )L−x

=−K ( 40−16 )

L2

L−x=

4024

∗L

2=

56L

x=L−56L→ x=L

6

CONDUCCION DE CALOR EN ESTADO ESTACIONARIO

Page 31: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 32: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 33: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 34: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 35: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 36: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 37: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor
Page 38: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

PAREDES PLANAS DE CAPAS MULTIPLES

Pared que consta de dos capas

Page 39: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Q=T ∞1−T ∞2

Rtotal

Rtotal=Rconv1+Rp1+Rp2+Rconv 2

Rtotal=1h1 A

+L1

K1 A+L2

K2 A+ 1h2 A

T1, T2 y T3 temperaturas de las superficies

Si queremos calcular la temperatura de la interfase 2, conociendo Q

Q=T ∞1−T 2

1h1 A

+L1

K1 A

→T 2=T ∞1+Q∗( 1h1 A

+L1

K 1 A )

RESISTENCIA TERMICA POR CONTACTO

Las superficies vistas al microscopio muestran picos y valles, al juntarse dos superficies hay espacios con aire el cual actúa como aislamiento por tener el aire conductividad baja. Por lo tanto la interfase ofrece resistencia a la transferencia de calor y esta resistencia por unidad de área se denomina resistencia térmica por contacto Rc.

Page 40: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

De la ley de newton de enfriamiento: Q=hc A ∆T interfase

Conductancia térmica por contacto directo hc=

QA

∆T interfase

La resistencia térmica por contacto es la inversa de la conductancia

Rc=1hc

=∆T interfaseQA

(m2−CW )

Esta resistencia térmica depende de la aspereza de la superficie, propiedades de los materiales y de la presión en la interfase y del tipo de fluido atrapado en esta.

Se puede minimizar esta resistencia aplicando grasa térmica o también insertando hojas metálicas suaves.

REDES GENERALIZADAS DE RESISTENCIAS TERMICAS

Para paredes en paralelo o combinadas serie – paralelo se puede obtener soluciones aproximadas suponiendo conducción unidimensional de calor

Page 41: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Dos placas paralelas aisladas exteriormente.

Q=Q 1+Q2=T 1−T 2

R1

+T1−T 2

R2

=(T 1−T2 )( 1R1

+ 1R2 )=

T1−T2

11R1

+ 1R2

Rtotal=1

1R1

+1R2

→1Rtotal

= 1R1

+ 1R2

Rtotal=R1R2

R1+R2

R1=L1

K1 A1

, R2=L2

K2 A2

Page 42: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Dos placas paralelas seguida de una pared (paralelas –serie)

Q=T 1−T3

Rtotal

Rtotal=R12+R3=R1R2

R1+R2

+R3

R12=R1R2

R1+R2

, R3=L3

K3 A3

CILINDROS CON CAPAS MULTIPLES

Flujo radial de calor a través de cilindros de diferentes conductividades térmicas.

Para condiciones de estado estable, el flujo de calor a través de cada sección será la misma.

Un cilindro compuesto de dos capas de longitud L, con convección en ambos lados.

Page 43: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Q=T ∞1−T ∞2

1h1 A1

+

ln( r2

r1)

2πL K1

+

ln( r3

r2)

2πL K2

+ 1h3 A3

Donde A1 = 2πr1L y A3 = 2πr3L

U = Coeficiente total de transferencia de calor para este sistema puede basarse sobre cualquier área, su valor numérico dependerá del área seleccionada

Como el diámetro mayor es el más fácil de medir, se escoge A3= 2πr3L, como el área base y la rapidez de transferencia de calor:

Q=U 0 A0 (T ∞1−T ∞2 )

Page 44: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

r3

2π r3 Lr1h1

=r 3

A0 r1h1

r3 ln( r 2

r1)

2π r3K1L=

r3 ln( r2

r1)

A3K1

r3 ln( r3

r 2)

2π r3K2L=

r3 ln( r3

r2)

A3K2

12π r3 Lh3

= 1A3h3

A0 = A3

U 0=1

r3

r1h1

+

r3 ln( r2

r1)

K 1

+

r 3 ln( r3

r2)

K2

+ 1h3

RADIO CRITICO DE AISLAMIENTO

Q=T ∞1−T ∞2

1h12 πr1L

+

ln(r 2

r 1)

2 πLK1

+ln( r❑r2

)2 πLK 2

+ 1h32 πrL

Para valor fijo de r1 y r2, la rapidez de flujo es una función de r y será un máximo para un valor de r

Page 45: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

dQdr

=0→d (∑ Ri )dr

=0 seraunMINIMO

∑ R i=1

h12 πr1L+

ln( r2

r1)

2 πLK1

+

ln( rr2)

2 πLK2

+ 1h32 πrL

∑ R i=1

h12 πr1L+

ln( r2

r1)

2 πLK1

+ ln r2 πLK2

−ln r2

2πL K2

+ 1h3 2πrL

d∑ R idr

= 1r 2π LK 2

− 1r2h32 πL

=0

rc=K2

h3

paraunCILINDRO

rc=2K2

h3

paraunaESFERA

Page 46: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

CALCULO DEL ESPESOR ECONOMICO DEL AISLANTE

Page 47: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

qa¿= Calor perdido con aislante por unidad de área

qb¿= Calor perdido sin aislante por unidad de área

θ = Tiempo (hr/año)

C2 = Costo de combustible ($/W-hr)

L2 = Espesor económico del aislante

C = costo de instalación ($/m3)

C1 = Depreciación anual

An = Ahorro neto anual por unidad de área ($/m2-año)

An = Cantidad de calor ahorrado al poner aislante por unidad de área menos depreciación

An=(qb} - {q} rsub {a} rsup { )θC2−L2CC1 (1 )

qa} = {{T} rsub {i} - {T} rsub {o}} over {{1} over {{h} rsub {i}} + {{L} rsub {1}} over {{K} rsub {1}} + {{L} rsub {2}} over {{K} rsub {2}} + {1} over {{h} rsub {o}}} left (2 right ¿

qb} = {{T} rsub {i} - {T} rsub {o}} over {{1} over {{h} rsub {i}} + {{L} rsub {1}} over {{K} rsub {1}} + {1} over {{h} rsub {o}}} left (3 right ¿

Page 48: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

(2) y (3) en (1)

An=¿

Derivando con respecto al aislante L2

d And L2

=(T i−T o )θC2

( 1hi

+L1

K1

+L2

K2

+ 1ho )

2

K2

−CC1=0

(T i−T o )θC2

( 1hi

+L1

K1

+L2

K2

+1ho )

2

K2

=CC1→1hi

+L1

K 1

+L2

K 2

+ 1ho

=√ (T i−T o )θC2

C C1 K2

L2={√ (T i−T o )θC2

CC1K2

− 1hi

+L1

K1

+ 1ho }K2

L2=√ K 2 (T i−T o )θC2

CC1

−( 1h i

+L1

K1

+ 1ho )K2

Ejercicio

Un cable eléctrico de 1.27cm de diámetro, va ser aislado con hule (K=0.1557 W/m-K); el cable estará instalado en el aire (h∞=8.5 W/m-K) a 21.1°C. Investigar el efecto que el espesor del aislamiento produce sobre la

Page 49: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

disipación de calor. Suponiendo una temperatura de 65°C en la superficie del cable.

r1=0.00635m

Radio critico

rc=Ka

h∞=0.1557

8.5=0.018m

Calor disipado sin aislante

Q=h∞2π r1L (T 1−T ∞ )→QL

=8.5∗π∗0.00635 (65−21.1 )=14.88Wm

Calor disipado con aislante

Q=T 1−T ∞

ln( rr1)

2π Kh L+

1h∞2πrL

→QL

=T 1−T ∞

ln( rr1)

2π Kh

+1

h∞2πr

Calor cedido con radio critico

QL

= 65−21.1

ln( 0.0180.00635 )

2 π∗0.1557+ 1

8.5∗2π 0018

=20.85Wm

Page 50: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

r (m) q/L (W/m)0.00635 14.880.008 17.0370.010 18.7870.0183 20.850.02 20.810.04 18.680.06 16.830.08 15.540.095 18.810.0935 14.880.10 14.6080.20 12.120.50 9.701.00 8.47

Ejercicio

Un horno de fierro con paredes blancas tiene un espesor de 0.025m y opera a tiempo completo θ= 8760hr/año a una temperatura interna T = 300°C, la

Page 51: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

temperatura ambiente es 26°C, el coeficiente combinado interior es hp-gas = 284 W/m2-K y la conductividad del fierro 62 W/m-K.

El horno será aislado externamente con un aislante de conductividad 0.087W/m-K, cuyo costo de instalación es de 1100$/m3, la depreciación anual (mantenimiento, impuestos) es estimado en 15%. el costo de combustible para la operación del horno es 0.8*10-5 $/W-hr, si el coeficiente pelicular externo es 9W/m2-K. Calcular

a) El espesor económico del aislanteb) El ahorro neto anual por unidad de área.

L2=√ K 2 (T i−T o )θC2

CC1

−( 1h i

+L1

K1

+ 1ho )K2

L2=√ 0.087 (300−26 )8760∗0.8∗10−5

CC1

−( 1284

+ 0.02562

+ 19 )0.087

L2=0.0906m

Page 52: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Ahorro neto

An=(qb} - {q} rsub {a} rsup { )θC2−L2CC1

qa} = {300-26} over {{1} over {284} + {0.025} over {62} + {0.0906} over {0.087} + {1} over {9}} =236.94 {W} over {{m} ^ {2}¿

qb} = {300-26} over {{1} over {284} + {0.025} over {62} + {1} over {9}} =2381.8 {W} over {{m} ^ {2}} ¿

An=(2381.8−236.94 ) 8760∗0.8∗10−5−0.096∗1100∗0.15

An=135.37$

m2−año

CONDUCCION DE CALOR A TRAVES DE TRAYECTORIAS CURVAS

Para ciertas condiciones el flujo de calor no es lineal, pero si unidimensional, por lo que se puede adaptar Fourier.

Q=−KA dTdS

(1 )

dTdS

gradiente de temperaturas en la trayectoria curvilínea

dTdr

=dTdZ

=0

Page 53: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

dTdS

=dTdθdθdS

(2 )

dS=rdθ→ dθdS

=1r

(3 )

(3) en (2)

dTdS

=dTdθ

1r

(4 )

(4) en (1)

d Q=−KdA dTdθ

1r

dA=Zdr

Page 54: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

d Q=−KZdr dTdθ

1r=−KZ∫

r1

r2

drr [ dTdθ ]

Q=−KZ ln( r2

r1) dTdθ ] (5 )

calculo dT/dθ usando coordenadas cilíndricas

1r2 ( ∂2T

∂θ2 )=0 integrando→dTdθ

=C1 (6 )→dT=C1dθ

volviendo a integrar

T=C1θ+C2 (7 )

condiciones de contorno

θ= 0 → T = T1 →T1 = C1*0 + C2 → T1 = C2 (8)

θ = θ2 → T = T2 →T2 = C1 θ2 + T1 →C1=T 2−T 1

θ2

(9 )

( 8 ) y ( 9 ) en ( 7 )

T=(T 2−T 1

θ2)θ+T1 (10 )

( 9 ) en ( 6 )

dTdθ

=T 2−T 1

θ2

(11)

( 11 ) en ( 5 )

Page 55: 2 2014 Curso Transferencia de Calor Ecuacion de Conduccion de Calor

Q=−KZ ln( r2

r1)(T 2−T 1

θ2)=¿KZ ln (r 1

r 2)(T 2−T 1

θ2)¿