apuntes transferencia de calor_2010

Notas de Transferencia de Calor, 2010, Dr. Juan Carlos Prince Avelino

1

TRANSFERENCIA DE CALOR

OBJETIVO: conocer y comprender los diversos mecanismos de transferencia de calor y aplicar sus leyes en la solución de problemas de ingeniería. UNIDAD I FUNDAMENTOS DE TRANSFERENCIA DE CALOR

1.1 Conceptos Básicos, definiciones, objetivo, mecanismos de transferencia de calor. 1.2 Conducción, Ley de Fourier. 1.3 Convección, Ley de Newton. 1.4 Radiación, Ley de Stefan-Boltzman.

UNIDAD II CONDUCCIÓN UNIDIMENSIONAL ESTACIONARIA

2.1 Superficies planas simples y compuestas. 2.2 Superficies cilíndricas simples y compuestas. 2.3 Superficies extendidas ( Aletas ).

UNIDAD III CONDUCCIÓN BIDIMENSIONAL ESTACIONARIA

3.1 Ecuación de conducción bidimensional, condiciones de frontera. 3.2 Solución numérica, método de diferencias finitas.

UNIDAD IV CONDUCCIÓN TRANSITORIA

4.1 Análisis Global, temperatura en función del tiempo. 4.2 Análisis gráfico de temperatura en función del tiempo en paredes planas,

cilindros y esferas. UNIDAD V CONVECCION FORZADA

5.1 Fundamentos de convección forzada, números adimensionales. 5.2 Convección forzada en flujo laminar. 5.3 Convección forzada en flujo turbulento.

UNIDAD VI CONVECCIÓN NATURAL

6.1 Fundamentos de convección natural, números adimensionales. 6.2 Convección natural en flujo laminar. 6.3 Convección forzada en flujo turbulento.

UNIDAD VII INTERCAMBIADORES DE CALOR

7.1 Generalidades de intercambiadores de calor. 7.2 Análisis de intercambiadores de calor por el método LMTD. 7.3 Análisis de intercambiadores de calor por el método eficiencia-NTU.

UNIDAD VIII RADIACIÓN

8.1 Fundamentos de radiación 8.2 Factores de forma e intercambio de calor radiante

Bibliografía: 1.- Introduction to Thermal Systems Engineering: Heat Transfer. Moran-Shapiro-Munson-DeWitt.Holman. Edit. John Wiley.


2

UNIDAD I FUNDAMENTOS DE TRANSFERENCIA DE CALOR DEFINICIONES Y CONCEPTOS BÁSICOS TERMODINAMICA: Estudia La relación entre calor y otras formas de energía. TRANSFERENCIA DE CALOR: Analiza la razón del calor transmitido en un sistema,

siempre que exista un gradiente de temperatura. En la primera Ley de la Temperatura hay dos tipos de energía de transferencia: el calor (Q) y el trabajo (W) y se relaciona mediante la siguiente expresión: Q-W=∆E La primera ley de la termodinámica es conocida también como ley de la conservación de la energía. TEMPERATURA: Es la propiedad del sistema que nos determina si este sistema esta en

equilibrio térmico con otro sistema. En otras palabras, los sistemas A y B en equilibrio térmico si sus temperaturas son idénticas, esto es:

TA = TB

LEY CERO DE LA TERMODINAMICA: Si los sistemas D y C están cada uno de manera

separada en equilibrio con un tercer sistema A, entonces B y C están en equilibrio.

TA = TB TB = TC, por lo tanto TC = TA

OBJETIVO DE LA TRANSFERENCIA DE CALOR: Es describir la manera en que la diferencia de temperaturas, digamos TA y TB gobierna la magnitud de la razón de transferencia de calor entre el sistema A y sus alrededores B, es decir, Q = Q [TA, TB, tiempo, propiedades termofísicas, geometría, flujo].


3

AREAS DE APLICACIÓN DE TRANSFERENCIA DE CALOR

a) Aislamiento térmico. b) Aumento de transferencia de calor. c) Control de temperatura.

MODOS DE TRANSFERENCIA DE CALOR

a) Conducción. b) Convección c) Radiación.

a) CONDUCCION Aplicando la 1ª Ley de la termodinámica al elemento diferencial

[ ]dtdEWxQxQx =−∆+−

Considerando W = 0

dtdExQxQx =∆+−

E = energía interna E = mU U = energía interna especifica m = masa del sistema m = Vρ V = A∆xρ


4

E = ρA∆xU Y a su vez: dU = C d T E= m C ∆T

C = calor especifico

][dtdTxC

dtdE

Α∆= ρ

LEY DE FOURIER EXPERIMENTAL

Qx = -KA dxdT

K = conductividad térmica A = área transversal

dxdT = gradiente de temperatura

Qx = -KA ⎢⎣

⎡⎥⎦⎤

∆−xTT 12

xdx

dQxQxxQx ∆+=∆+

xdxdTKA

dxd

dxdTKAxQx ∆⎢

⎣

⎡⎥⎦⎤−+−=∆+


5

dtdExQxQx =∆+−

dxdTxCAx

dxTdKA

dxdTKA

dxdTKA ∆=

⎥⎥⎦

⎤

⎢⎢⎣

⎡∆⎜⎜

⎝

⎛⎟⎟⎠

⎞−+−−− ρ2

2

tTxCAx

xTKA

∂∂

∆=∆∂∂ ρ2

2

tTC

xTK

∂∂

=∂∂ ρ2

2

Tratamos de obtener una solución de la forma T(x,t) para derivar xT∂∂

dxdTKAQx −=

tTC

xTK

∂∂

=∂∂ ρ2

2

Si divido entre K:

tT

KC

xT p

∂∂

=∂∂ ρ

2

2

α1

=K

pCKρ

α = [Difusividad térmica]

Ecuación de conducción con generación interna de calor

tTq

xT

∂∂

=+∂∂

α1

2

2

K, conductividad térmica; KmW° ó Cm

W°

Ley de Fourier:


6

dxdTKAq −=

dTKAdxqT

T∫∫ −=2

1

1

0

) )(( 12

2

0

2

1

TTKALqKATqxT

T−−=−⇒−= θ

⎢⎣

⎡⎥⎦⎤−

−=L

TTKAq 12

NOTAS:

Las temperaturas pueden ser en ºC ó ºK al final nos dará el mismo resultado. Cuando se multiplica por un diferencial de temperaturas las unidades de

temperatura no importan.


7

CONVECCCION Ley de Enfriamiento de Newton (?) La transferencia de calor por convección o simplemente convección es el que se lleva acabo por el flujo de un fluido. El fluido actúa como un acarreador de energía que es tomada (o entregada a ) una pared sólida. Tw= Temperatura de pared Tα = temperatura de fluido El flujo de calor por convección obedece la Ley de Newton:

)( ATTwhq ∞−=

Donde h: Coeficiente de transferencia de calor por convección )( Kmw o2

La h, es el coeficiente de mayor importancia de la convección. Existen dos clasificaciones para la convección:

a) Libre o Natural b) Forzada

Convección Natural: Cuando un fluido es llevado a través de una pared por efectos de

flotación.


8

Convección Forzada: Cuando un fluido es obligado a pasar a través de las paredes. La convección puede ser monofásica, bifásica ó puede darse con un cambio de fase

(ebullición y condensación). RADIACCION En este caso hablaríamos de la energía radiante emitida por un cuerpo, debido a su temperatura trasmitida en el espacio. Según Planck la radiación es trasmitida en forma de fotones discretos. La radiación puede ser transmitida en el vació mientras que la conducción y convección requieren de un medio material. El máximo flujo ( )2m

w en que el calor puede ser radiado de una superficie es: 4" sTq σ=

Donde: TS= Temperatura absoluta en °K de la superficie σ = Constante de stephan Boltman )( 4

2867.5 Km

wE −

La razón neta de intercambio de calor entre una superficie y sus alrededores por unidad de área es:

)( 4sup

4" TTwTq −= ε Donde: ε = Emisividad 0 1≤≤ ε


9

Para un cuerpo negro 1=ε • Cuando no me dan el valor ε , considero 1=ε • Siempre la temperatura debe estar dada en °K

NOTA: 273+°C=°K


10

1.-La superficie exterior de un muro de concreto de 0.2 m de espesor se mantiene a

una temperatura de -5 °C, mientras que la interior se mantiene a 20°C. la conductividad térmica del concreto es 1.2 w/m°K. Determinar la perdida de calor a través del muro de 10m de largo y 3m de alto.

Solución: Q= 4500 W


11

2.- Se transfiere calor a razón de 0.1 KW a través de un aislante de fibra de vidrio (ρ=100 kg/m2) de 5cm de espesor y 2m2 de área. Si la superficie caliente esta a 70 °C, determine la temperatura de la superficie fría.

Solución: T2= 0.55 °C


12

3.- Cuando la superficie de una placa plana de 0.1x0.5m se mantiene a 30°C, la tasa de transferencia de calor del aire caliente a 100°C que fluye a su lado es de 125W. cual es el coeficiente de transferencia de calor por convección entre la placa y el aire.

Solución: h = 35 W/m² °C


13

4.-Durante un invierno, la superficie de un rió desarrolla una capa de hielo de espesor desconocido l. se conoce la temperatura del agua del rió de es de 4°C, la temperatura del aire atmosférico es de -30°C, y la temperatura debajo de la capa de hielo de 0°C. La conductividad térmica del hielo es K=2.25 w/m°K, los coeficientes de transferencia de calor sobre el lado del agua y del aire de la capa de hielo son h1=500w/m2°K y h2=100 w/m2°K respectivamente. Calcule la temperatura sobre la superficie superior de la capa de hielo y el espesor L del hielo. Solución: L= 0.01125 m T2= -10 °C


14

5.- Un calentador radiante de franjas metálicas tiene 6mm de ancho, una longitud total de 3m. la emisividad de la superficie de la franja es 0.85. A que temperatura debe estar la franja para disipar 1600W de calor a un cuarto que esta a 25°C. Solución: Tw = 893 °C


15

6.- Un cilindro de 5 cm de diámetro se calienta hasta una temperatura de 200°C (473°K), mientras una corriente de aire a 30°C (303°K) y con una velocidad de 50 m/s le sopla transversalmente. Si la emisividad de la superficie es 0.7. Calcular la perdida total de calor por unidad de longitud si las paredes de la habitación en la que esta colocado el cilindro esta a 10°C. Solución: Q/L = 5078 W/m

7.- Una pared de 5cm de asbesto poco compacto esta colocada entre dos placas a 100 y 200 °C, calcular el calor transferido a través de la capa de asbesto (poco compacto)


16

Solución: q/L = 226 W/m²


17

8.- Un aislante tiene una conductividad térmica de 1 W/m°C ¿Cuál será el espesor necesario para que exista una caída de temperatura de 500°C para un flujo de calor de 400W/m2? Solución: ∆x = 0.04 m


18

9.- Las temperaturas de las caras de una pared plana de 15cm de espesor son 370°C y 93°C. La pared esta construida con un vidrio especial que tiene las siguientes propiedades: K=0.78W/m°C, ρ=2.700 kg/m3, Cp=0.84 KJ/kg°C, ¿Cuál es el flujo de calor a través de la pared en condiciones estacionarias? Solución: q/L = 1440 W/m²


19

10.- una de las caras de una pared plana se mantiene a 100°C mientras que la otra se expone al ambiente que esta a 10°C, siendo h=10W/m2°C el coeficiente de convección. La pared tiene una conductividad térmica K=1.6 W/m°C y un espesor de 40 cm. Calculese el flujo de calor a través de la pared.

q/L = 257 W/m²


20

11.- Un oleoducto de 50 cm de diámetro transporta en el ártico petróleo 30°C y esta expuesto a una temperatura ambiente de -20°C. Un aislante especial de polvo de 5 cm de espesor y de conductividad térmica 7mW/m°C cubre la superficie del oleoducto, el coeficiente de transferencia de calor es 12 W/m2°C. Estímese la perdida de energía del oleoducto por unidad de longitud. Q/L = 10.77 W/m


21

UNIDAD II CONDUCCION UNIDIMENSIONAL ESTACIONARIA PAREDES PLANAS SIMPLES O COMPUESTAS

T0 TL

KALRt =

El fenómeno se rige por al siguiente Ecuación:

tT

xT

∂∂

=∂∂

α1

2

2

Por ser un fenómeno estacionario

02

2

=dtdT


22

Por lo tanto:

02

2

=dx

Td 1CdxdT

=⇒ 21 CXCTL

o

T

T

L

o

+

0=⎢⎣

⎡⎥⎦⎤

dxdT

dxd 21 CdxCdT L

oTo

L += ∫∫ ( ) 21 CoLCTT oL +−=−

0=⎢⎣

⎡⎥⎦⎤

dxdTd Condiciones de frontera 21

21

CxCTCLCTT oL

+=+=−

1CdxdTd =⎢⎣

⎡⎥⎦⎤∫

L

o

TTLxTTox

=→==→=

Aplicando oTTox =→=

( ) 21 CoCTo +=

2CTo = Aplicando LTTLx =→=

oTXCT += 1

( ) oL TLCT += 1

oL TLCT += 1

LTT

C oL −=1

21 CXCT +=

ooL TX

LTT

T +−

=

LTT

CdxdT oL −== 1


23

⎥⎦

⎤⎢⎣

⎡ −=Κ−=

L

oL

TTT

q"

[ ]oL TTL

qq −=ΚΑ

−=Α= "

[ ]Lo TTL

−ΚΑ

=

*Al término KAL se le conoce como resistencia térmica.

KALRt =

RtTT

q Lo −=


24

PAREDES COMPUESTAS Pared Simple Pared Compuesta

RtTT

q Lo −= Rt

TTq FC −=

FC RRRRRRt ++++= 321


25

⎥⎦

⎤⎢⎣

⎡+

Κ+

Κ+

Κ+

Α

−=

Fc

FC

hLLL

h

TTq

111

3

3

2

2

1

1

Fc hKL

KL

KL

h

U11

1

3

3

2

2

1

1 ++++=

U = Coeficiente global de transferencia de calor.

( )FC TTUQ −Α=


26

1.- La caja aislada mostrada en la figura esta diseñada para mantener el aire atrapando en ella a una alta temperatura de 50°C. La temperatura afuera es de 10°C. Para mantener constante la temperatura del aire, La transferencia de calor que sale por el aislante se recupera mediante un calentador de resistencia eléctrica colocado en el centro de la caja. Calcular la potencia eléctrica disipada por el calentador; las dimensiones del espacio interno (aire) son x = 1m, y = 0.4m y z = 0.3m. La pared aislada consiste de una placa de 10 cm. de espesor de fibra de vidrio que se encuentra entre dos placas de madera de abeto cada una de 1 cm. de espesor. Los coeficientes de transferencia de calor sobre las superfi-cies interna y externa de la pared son Κo25 m

w y Κo215 mw respectivamente.

U = 0.304 W/m2°K Am = 2.125 m2

Q = 26 W


27

PLACAS CILINDRICAS Simples y Compuestas a) Cilindro Simple Q, se transfiere del centro de la tubería hacia la pared externa del mismo. Ley de Fourier

AqQ "= Tuberías: Ae: área exterior Ai: área interior

iiee AqAqQ "" == q”e=(2πreL)=q”i(2πriL) q”ere=q” iri

e

i

i

e

rr

qq

=""

Para un cilindro simple, las temperaturas se distribuyen de la siguiente manera:


28

)ln(

)ln()(

i

e

ieii

rr

rr

TTTT −−=

Calculamos de la siguiente manera:

( )ei

i

eTT

rrKLQ −=

)(ln

2π

Nota: La temperatura y el calor están en función del radio.


29

CILINDRO COMPUESTO

t

FC

RTTQ −

=

ee

e

i

ii hLKr

r

LKr

r

LKr

r

hRt

Α+

⎟⎠⎞⎜

⎝⎛

+⎟⎠⎞⎜

⎝⎛

+⎟⎠⎞⎜

⎝⎛

+Α

=1

2

ln

2

ln

2

ln1

3

2

2

1

2

1

1

πππ

Lrii π2=Α Lree π2=Α

Nota: Los problemas se resuelven hasta donde me lo permitan los datos. Ejemplo: Datos:

Twext t

Fwext

RTTQ −

=

TF como solo tengo he

he ee

t hR

Α=

1 ( )Fwextee

ee

Fwext

TThh

TTQ−Α

=

Α

−=

11

Convección

( ) ( )Fwexteew TThQTThAq −Α=⇔−= ∞ Convección


30

2.- Una tubería de 2 in cedula 40, tiene una conductividad e 47 w/m°C. el fluido dentro del tubo tiene un coeficiente de convección de 170 w/m2°C, la superficie exterior de la tubería es cubierta con un aislante de fibra de vidrio de 12.5 mm de espesor y K=0.04 w/m°C. el coeficiente de convección para la superficie exterior del aislante es 12 w/m2°C. la temperatura del fluido dentro de la tubería es 160 °C y la temperatura ambiente es 21°C. Calcular:

a) la perdida de calor por metro de longitud b) la temperatura en cada una de las interfaces de la tubería y el aislante.

⎥⎥⎥⎥

⎦

⎤

⎢⎢⎢⎢

⎣

⎡

+⎟⎠⎞⎜

⎝⎛

+⎟⎠⎞⎜

⎝⎛

+

−=

ee

e

i

ii

ei

rhKr

r

Kr

r

rhL

TTQ

ππππ 21

2

ln

2

ln

211

2

1

1

1

ee

e

i

ii hLKr

r

LKr

r

hRt

Α+

⎟⎠⎞⎜

⎝⎛

+⎟⎠⎞⎜

⎝⎛

+Α

=1

2

ln

2

ln1

2

1

1

1

ππ

Q = 80.6 W/m Rt = 1.73 W/m2 °C T1 = 157 °C T2 = 157 °C T3 = 46 °C T4 = 21 °C


31


32

SUPERFICIES EXTENDIDAS ALETAS Al desarrollar el balanceo de energía se obtiene la siguiente ecuación diferencial de la aleta.

( ) 02

2

=∞−−Α TThPdx

TdK C

Resolviéndose para:

1) Aleta larga 2) Aleta corta aislada en el extremo 3) Aleta corta con convección en el extremo

ALETA LARGA

∞−= TTSi x)(θ ( ) mxbx eTTTT −

∞∞ −+=)(

∞−= TTbbθ Si x = 0

eKAhPm =

( )

bx

bx

mbx

TTTTTT

eTTTT

=−+=−+=

∞∞

−∞∞

)1(

)0(


33

( )( ) ( )mxTTTT

mx

bx

b

−−=−−=

∞∞ expexp

)(

θθ

Nota: x debe ser en el punto donde se desea conocer la temperatura ( ) 2

1ρθ hKAq Cbb =

ALETA CORTA AISLADA EN EL EXTREMO

( )[ ]( )

( ) ( )mLtghhPKAq

mLxLm

ebb

b

21

coshcosh

θ

θθ

=

−=

ALETA CORTA CON CONVECCION EN EL EXTREMO

( )[ ] ( )[ ]

( ) [ ]

( ) ( )Cebb

b

mLtghhPKAq

mLsenhmKhmL

xLmsenhmKhxLm

21

cosh

cosh

θ

θθ

=

⎟⎠⎞

⎜⎝⎛+

−⎟⎠⎞

⎜⎝⎛+−

=

Donde: PALL C

C +=


34

EFICIENCIA DE ALETA

b

b

Ahq

calordeciatransferenMaxrealcalordeciatransferen

θη

exp

==

Donde: Aexp= área expuesta

bb Ahq θη exp=


35

3.- El contenedor cilíndrico de una motocicleta es construido de duraluminio y tiene una altura de 0.15 mts. Y un diámetro externo de 50 mm. bajo condiciones típicas de operación la superficie exterior del cilindro esta a una temperatura de 500 ºK y expuesta al aire ambiente que se encuentra a 300ºK, con un coeficiente de convección de Km

wº50 2

.Aletas anulares de perfil rectangular son agregadas para incrementar la transferencia de calor a los alrededores. Se agregan 5 de estas aletas con un espesor de 6 mm. y una longitud de 20 mm., igualmente espaciadas ¿Cuál es el incremento de transferencia de calor debido a la adición de las aletas?

Qaletas = 413 W Qbase = 189 W Qcilindro = 236 W Am = 0.00014 m2

Aexp = 0.043 m2


36


37

4.- Una tubería de vapor de radio externo de 4 cm esta cubierta con una placa de asbesto de 1 cm de espesor, que a su vez esta cubierto con fibra de vidrio con 3 cm de espesor. La superficie de la tubería de vapor esta a 330°C y la superficie exterior de la fibra de vidrio esta a 30°C. CALCULAR:

a) La temperatura de interfase entre el asbesto y la fibra de vidrio b) La transferencia de calor por metro de longitud de la tubería.

Kasbesto = 0.149 W/m °C Kvidrio = 0.038 W/m °C Tinterfase = 298 °C q/L = 136 W/m


38

UNIDAD III CONDUCCION BIDIMENSIONAL ESTABLE

02

2

2

2

=∂∂

+∂∂

yT

xT

Condiciones de frontera

0=x bTT = 0=y ∞= TT Lx = ∞= TT Hy = ∞= TT

Cambio de variable

∞−= TT yxyx ),(),(θ

TTTT

22 ∂=∂

∂=∂−= ∞

θ

θθ

Sustituyendo (3) en (1)

02

2

2

2

=∂∂

+∂∂

yT

xT


39

0=x ∞−= TTbθ bθθ =

Lx = 0=θ 0=y 0=θ Hy = 0=θ

Separación de variables

yxdx

'=∂θ 'xy

y=

∂∂θ

yxx

''2

2

=∂∂ θ ''2

2

xyy

=∂∂ θ

0""=

+y

xyyx

0

"

=+

xy

xyx

0""=+

yy

xx

2"" λ=−=

yy

xx

2" λ=

xx

xx 2" λ=

xx 2" λ=

0" 2 =− xx λ ……….(1)

02 =+′′ yy λ ………(2)


40

[ ]

λαλα

λα

λα

α

α

α

αα

α

α

α

±==

=−

=−

=

=

=

22

22

22

2

00

"'

x

xx

x

x

x

eee

exex

ex

[ ]

λαλα

λα

λα

α

α

α

αα

α

α

α

i

eee

eyey

ey

y

yy

y

y

y

±=−=

=+

=+

=

=

=

22

22

22

2

00

"'

( )

( )yiyi

y

xxx

eCeCy

eCeCxλλ

λλ

−

−

+=

+=

43

21

( ) ( ) ( )

( ) [ ][ ]yiyixxyx

yxyx

eCeCeCeC

yxλλλλθ

θ−− ++=

=

4321,

,

Identidades de Euler

( ) ( )[ ] ( )[ ][ ]( )[ ] ( )[ ][ ]

( ) ( ) ( )[ ] ( ) ( )[ ]yBCosySenxACoshxSenhK

yiSenhyCoshCyiSenhyCoshC

xSenhxCoshCxSenhxCoshC

yx

yx

λλλλθ

λλλλ

λλλλθ

++=

−++

+++=

,

43

21,

Donde:

3

4

1

2

31

CCB

CCA

CCK

=

=

=

Aplicar condiciones de frontera


41

x = L θ = 0

[ ][ ]

LCoshLSenhA

LACoshLSenhyCoshBySenhLACoshLSenhK

λλ

λλλλλλ

−=

=+++=

00

Sustituimos

( ) [ ]

( )

( )[ ]( )

( )( )[ ]( ) [ ]yCosBySen

LCoshLxSenhK

xLCoshLxSenh

LCoshLCoshLSenhLLCoshSenh

aCoshbSenhbCoshaSenhbaSenh

LCoshLLCoshSenhLSenh

yCosBySenLCoshLCoshLSenhLSenhK

yx

yx

λλλ

λθ

λ

λλλλλ

λλλλ

λλλλλλθ

+−

=

−

−

±=±

−

+⎥⎦

⎤⎢⎣

⎡+⎟⎟

⎠

⎞⎜⎜⎝

⎛−+=

,

,

y = 0 θ = 0

( )[ ]( ) [ ])0()0(0 λλλ

λ CosBSenLCosh

LxSenhK +−

=

B=0

( )[ ]( ) ySen

LCoshLxSenhKyx λ

λλθ −

=),(


42

∏=

mH2

y = H θ = 0

( )[ ]( )

( )

( )y

HnSen

LHnCosh

LxHnSenh

K

Hn

nHHSen

HSenLCosh

LxSenhK

nnyx

ππ

π

θ

πλ

πλλ

λλ

λ

−=

=

==

−=

∑∞

=1,

0

0

x = 0 θ = θb Sustituyendo y para una Kn en particular

( )

∫ ∫ ⎟⎠⎞

⎜⎝⎛−=⎟

⎠⎞

⎜⎝⎛

−=

−=

H

y

H

nyb

nb

nb

ydH

mSenyHnSenL

HntghKyd

HmSen

YHnLSen

HntghK

yHnSen

LHnCosh

LHnSenh

K

0 0

0

ππππθ

ππθ

ππ

π

θ

[ ]

[ ]1

000

−−=

−−=−=∫

ππ

ππ

ππ

π

CosmmH

CosCosmmHy

HmCos

mHyd

HmSen

H

y

H

Con valores de: 0 para n par

πmH2 para n impar


43

y

H

ydH

mySenHnSen ππ

∫0

Con valores de: O para nm ≠

2H para m = n

22 H

HLntghK

mH

nbπ

πθ −= m = n = impar

HLntghn

K bn ππ

θ4−=

( )

( )y

HnSen

LHnCosh

LxHnSenh

HLnntgh

Kn

nb

yxπ

π

π

ππθ

θ−

−= ∑∞

=

141

,


44

Ecuación de condición

tT

yT

xT

∂∂

=∂∂

+∂∂

α1

2

2

2

2

02

2

2

2

=∂∂

+∂∂

yT

xT

Para x son valores en y Para y son valores en x Condiciones en la frontera

0=x bθθ = 0=y 0=θ Lx = 0=θ Hy = 0=θ

Separación de variables

( ) ( ) ( )yxyx yx=,θ

yxx

'=∂∂θ 'xy

y=

∂∂θ

yxx

''2

2

=∂∂ θ ''2

2

xyy

=∂∂ θ

Sustituimos en Ec.1

0"" =+ xyyx


45

Dividido entre “y”

0"" =+y

xyx

Dividido entre “x”

0""=+

yy

xx

2"" λ=−=yy

xx

000

2

2

2

≤

≥

=

λ

λ

λ

2" λ=

xx

xx 2" λ=

0" 2 =− xx λ

2" λ=−

yy

yy 2" λ=−

0"2 =+ yyλ

Resolvemos 0" 2 =− xx λ

x

x

x

exex

ex

α

α

α

α

α2"

'=

=

=

[ ] 00

22

22

=−

=−

λα

λαα

αα

x

xx

eee


46

λαλα

λα

±==

=−2

22 0

( )xx

x eCeCx λλ −+= 21

Resolvemos 0" 2 =+ yy λ

[ ]

λαλα

λα

λα

α

α

α

αα

α

α

α

i

eee

eyey

ey

y

yy

y

y

y

±=−=

=+

=+

=

=

=

22

22

22

2

00

"'

( )

yiyiy eCeCy λλ −+= 43

Decimos entonces:

( ) ( ) ( )

( ) [ ][ ]yCosCySenCeCeC

yxxx

yx

yxyx

λλθ

θλλ

4321,

,

++=

=−

Condiciones de Frontera

∞→= Lx 0=θ

[ ][ ]yCosCySenCeCeC LL λλλλ43210 ++= −

1=λ 01 =C

( ) [ ]yCosCySenCeC xyx λλθ λ

432, +=∴ − Y=0 θ=0


47

[ ]( ) ( )[ ]

( ) [ ]( )

( ) ySenKe

KCC

ySeneCC

ySenCeC

CCCeC

CosCSenCeC

xyx

xyx

xyx

x

x

λθ

λθ

λθ

λ

λ

λ

λ

λ

−

−

−

−

−

=

=

=

=

=+=

+=

,

32

32,

32,

4

432

432

0100

000

Y=H θ=0

( ) yHnSeneK

Hn

nHHSenHSenKe

n

XHn

nyx

x

πθ

πλ

πλλλ

π

λ

∑∞

=

−

−

=

=

=∴==

1,

00

X=0 θ=0

yHnSenK

nnb

πθ ∑∞

=

=1

Para una n en particular

] [ ]10

00

00

−−=−=

⎥⎦⎤

⎢⎣⎡−==

⎥⎦⎤

⎢⎣⎡=⎥

⎦

⎤⎢⎣

⎡

∫∫

∫∫

ππ

ππ

ππ

πππ

πππθ

CosmmHy

HmCos

mH

yH

mCosmHyd

HmSen

Hmyd

HmSen

ydH

mSenyH

mSenKydH

mSen

H

y

H

y

H

y

H

ny

H

b

Con valores de: 0 para m par

πmH2 para m impar


48

Integral de Fourier

⎥⎦⎤

⎢⎣⎡∫ y

H

ydH

mSenyH

mSen ππ

0

Con valores de: 0 para nm ≠

2H para m = n

( )πθ

πθπ

θ

πθ

nHnH

HnH

K

HKm

H

bbb

n

nb

422

2

22

2

===

−=

( ) n

yHnSene

XHn

n

byx

π

πθ

θ

π−

∞

=∑−=

3,2,1,

4


49

METODO DE DIFERENCIAS FINITAS

T1 = 63.6 °C T2 = 72.1 °C T3 = 55.09 °C T4 = 62.47 °C T5 = 31.81 °C T6 = 36.07 °C


50

Ejemplo 2


51

UNIDAD IV CONDUCCION TRANSITORIA ANALISIS POR BLOQUES Este es el caso más sencillo de conducción transitoria ya que se considera que el gradiente de temperatura dentro de un cuerpo es despreciable y por lo tanto la única variable independiente es el tiempo (t). Aplicando la primera Ley de la Termodinámica Razón del flujo de = Razón de incremento de energía Calor hacia el sólido interna del sólido

( )[ ]dtdTVCTThA t ρ=−∞

( ) ( )[ ]dtdTTT

VChA

t =−− ∞ρ1

( )[ ]dtdTTT

VChA

t =+−− ∞ρ

Si VChAmρ

=


52

( )

( ) ∞−=−

TTdT

mdtt

t

( )

( )

( )

∫ ∫ ∞−=−

t

o

T

Tt

tt

i TTdT

dtm

( )[ ] ( )∞∞ −−−=− TTTTmt it lnln

( )

∞−

∞−=−

TTTT

mti

tln

( )

∞−

∞−=−

TTTT

ei

tmt

( ) ( ) mt

it eTTTT −∞−+∞= Temperatura en función del tiempo. Nota:

a)

yxe

yxz

z =

=−

−

ln

b) yxyx lnlnln =−

c) 1º=e Definimos la longitud característica como:

AVLs =

Donde: V, volumen de la pieza A, área de la pieza

( )( )[ ] ( )t

i

T

Ttto TTmt

+∞−=− ln


53

Definimos al Número de Biot como:

KhLB s

i = Donde: K: Conductividad térmica de la pieza

Si 1.0≤iB entonces aplicaremos la teoría anterior. CONDUCCION TRANSITORIA UNIDIMENSIONAL

( ) 01, 2

2

=−→dtdT

dxTdtxT

α

Solución ( )txT , En la figura se observa una placa de espesor 2L que inicialmente se encuentra a la temperatura iT y que es inmensa en un medio que esta a temperatura ∞T .


54

Si definimos: ( ) ( ) ∞−= TTtx tx,,θ

el problema se puede expresar como:

dtd

xdxd θθ 1

2

2

=

Esta expresión es la ecuación de conducción unidimensional considerando al tiempo como variable, con condición inicial.

iθθ = en 0=t Donde ∞−= TTiiθ Y condiciones de frontera

0=dxdθ En x = 0

θθ hdxdK =− En x = 1

Utilizado el método de separación de variables y colocándoles las condiciones de frontera iniciales, se llega a una solución ( )tx,θ Para la cual se emplean series de Fourier para las siguientes 3 casos:

a) Placa plana b) Cilindro largo c) Esfera

Graficado mediante tres parámetros

1) Número de Biot

Khró

KhLB o

i =


55

2) número de Fourier

22o

o rTó

LTF αα

=

3) Parámetro Geométrico

orró

Lx

Tomar L de tabla 2.3 Utilizar tablas 2.3, 2.38 y 2.39

Nota: Esfera 2

3

34

DA

rV

π

π

=

=


56

PROBLEMAS UNIDAD IV

1.- La temperatura de una corriente de gas es medida con un termocople, el acoplamiento puede ser aproximado como una esfera de diámetro 1mm, con una conductividad de 25 W/m°C entre el acoplamiento, la densidad es igual a 8400 kg/m3 y un Cp= 400 J/kg°C. con un coeficiente de transferencia de calor de 560 W/m2°C entre el acoplamiento y el gas ¿Cuánto tiempo le toma al termocople registrar el 99% de la diferencia de temperaturas aplicada?

Bi = 0.0037 t = 4.6 s


57

2.- Una placa de hierro de 5cm de espesor, esta inicialmente a 225°C de repente la placa se expone a un ambiente que se encuentra a 25°C y con un coeficiente de transferencia de calor de 500 W/m2°C. calcular la temperatura en el centro de la placa y a 1cm de profundidad a 2 minutos de comenzar el enfriamiento. Determinar la energía removida de la placa por metro cuadrado durante los 2 minutos.

F0 = 3.9 T0 = 125 °C Q/A = 15349 Kj/ m2

1/Bi = 5.16


58

3.- Una barra de acero dulce de 6cm de diámetro a 38°C de repente se sumerge en un liquido a 93°C con un coeficiente de transferencia de calor de 110 W/m2°K. Determinar el tiempo requerido para que la barra se caliente hasta 88°C.

K = 43 W/ m °C ρ = 7800 kg/ m3

Cp = 473 J/Kg °C t = 121 s


59

4.- Una esfera de cobre de 5cm de diámetro esta inicialmente a una temperatura uniforme de 250°C. La esfera se expone de forma rápida a un ambiente a 30°C y con un coeficiente de transferencia de calor de 28 W/m2°C. Utilizando el método de análisis de la capacidad global, calcular el tiempo necesario para que la esfera alcance una temperatura de 90°C.

t = 1327 s m = 0.000979 1/s K = 374 W/m °C


60

UNIDAD V CONVECCION FORZADA Conceptos básicos de Convección. La convección es el mecanismo de transferencia de calor debido al movimiento de un fluido. El problema fundamental en la transferencia de calor por convección consiste en determinar la relación entre el flujo de calor a través de la pared de un sólido (q”) y la diferencia de temperaturas ( )∞−TTw que existe entre la pared y el fluido. El coeficiente que relaciona las dos cantidades anteriores han sido ya establecidos como: q”=h ( )∞−TTw Donde h: Es el coeficiente de transferencia de calor por convección. La convección puede ser libre o forzada; el flujo es forzado cuando otro mecanismo lo empuja a pasar por un cuerpo sólido. En convección libre o natural el movimiento del fluido ocurre por si mismo sin asistencia de un mecanismo externo, se debe a la flotación relativa en diferentes regiones del fluido. Los flujos pueden ser internos o externos, laminares o turbulentos. Los fluidos pueden estar en una fase, bifásico o en cambio de fase (condensación y ebullición). PARAMETROS ADIMENCIONALES

Reynolds: v

uLRe =

Prandtl: K

CP p

r

µ=

Nusselt: KhLNu =

Donde: u, velocidad del fluido L, longitud característica Cp, capacidad calorífica K, conductividad térmica ν, viscosidad cinemática µ, viscosidad dinámica del fluido NOTA: Si hablamos de tubería L = Ø Si hablamos de una placa L = longitud de la placa


61

FLUJO LAMINAR EN DUCTOS Factor de fricción y caída de presión.

eDRf 64=

2

2u

DfLP

h

ρ=∆

PDh

Α=

4

Α= um ρ

vuDR h

eD =

Donde: hD = diámetro hidráulico A = área de la sección transversal P = perímetro ρ = densidad


62

DD

DD

D

D

D

PDh ===

⎥⎦⎤

⎢⎣⎡

=Α

=ππ

π

π

π

π

44

1

44

44

4 2

22

DDh = K

hDKhLNu h==

( )( )

aD

aaa

aaaD

h

h

=

===44

44 2

( )( )

32

64

132

4

32

4 2

2

aaa

a

a

a

aa

Dh ===⎥⎦⎤

⎢⎣⎡

=

TRANSFERENCIA DE CALOR EN FLUJO LAMINAR Cuando el flujo esta completamente desarrollado y es laminar, el valor del coeficiente de transferencia de calor h, se puede obtener a través de una tabla (tabla 6.1). Existen dos casos:

a) Tubería a temperatura constante b) Tubería con flujo de calor constante

La evaluación de las propiedades físicas del fluido se hacen a la temperatura promedio:


63

( )oib TTT +=21

Thq ∆=" (Definición básica de convección) ( )∞−=∆ TTT w

Nota: si me dan longitud y me piden To, propongo una temperatura. DIFERENCIA DE TEMPERATURAS MEDIA LOGARITMICA

2

1

21

lnTT

TTT

LTMDT

LM

LM

∆∆∆−∆

=∆

=∆

FLUJO TURBULENTO A TRAVES DE UN DUCTO

a) Factor de fricción y caída de presión

2

2u

DfLp ρ=∆

- diagrama de moody (tuberías lisas y rugosas) f 4

1316.0

−= eRf 42ERe <

- Lisas 51

184.0−

= eRf 64 12 ERE e <<


64

a) Fórmula analítica de Colburn

328

r

tP

fs =

Donde: St = es el número de Stanton

uChs

pt ρ=

b) Fórmula pata tubería lisa.

3

18.0023.0 reu PRN = Considerando 64 12 ERE e <<

c) Fórmula de Dittus - Boelter

nreu PRN 8.0023.0=

Considerando:

524.12500

1207.0

ER

P

e

r

≤≤

≤≤

n = 0.4, si el fluido se esta calentando, bw TT > n = 0.3, si el fluido se esta enfriando, bw TT <

d) Fórmula de Seider – Tate

14.0

318.0027.0 ⎥

⎦

⎤⎢⎣

⎡=

wreu PRN

µµ

Preferentemente para fluidos muy viscosos. Donde: wµ , es la viscosidad del fluido a la temperatura de la pared (Tw) µ , viscosidad dinámica a Tb En todos los casos:

KhdNu = , de donde se despeja la h


65

PROBLEMAS UNIDAD V

1.- Por tubo de 60mm de diámetro entra agua a una temperatura de 20°C y un flujo masico de 0.01 Kg/s. considere que por la pared del tubo se transfieren 2000 W/m2 Calcular:

a) la longitud del tubo requerida para obtener una temperatura de salida de 80°C b) la temperatura superficial en la salida del tubo

RESPUESTAS: Re = 377 Nu = 3.657 h = 47 w/m² ºC L = 6.64 m Twsal = 131ºC


66

2.- Fluye agua a una razón de 0.5 kg/seg a través de una tubería de 10cm de longitud con un diámetro interior de 2 cm. La tubería esta siendo calentada con flujo uniforme de calor en la pared de 5E10 W/m2. Evalué las propiedades del agua a 20°C. Calcular:

a) la caída de presión a lo largo del tubo b) el coeficiente retransferencia de calor basado en la teoría de Colburn c) el coeficiente de transferencia de calor con la teoría de Dittus-Beelter d) la diferencia de temperaturas de la pared y la temperatura local media del agua e) el incremento de la temperatura experimentada por la temperatura media del agua

desde la entrada hasta la salida.

Solución: Re = 32 330 a) ∆P = 14 520 Pa b) h = 5 212 w/m² ºC c) h = 6052 w/m² ºC d) ( Tw –Tb)m = 8.2 ºC e) ( Ti – To ) = 15 ºC


67


68

3.- fluye agua a una razón de 5kg/seg a través de una tubería de 5 cm de diámetro y 10 m de longitud estando la pared a 80°C. Si el agua entra a 20°C, calcular la temperatura de salida del agua.

Solución: Re = 186 912 h = 7890 w/ m² ºC Segunda interacción; Nu = 596 To = 46 ºC


69


70


71

4.- Considere el flujo de agua a una razón de 0.015 kg/seg a través de un ducto cuadrado de 2 cm de lado cuyas paredes se mantienen a una temperatura unidor me de100°C. Considere que es un flujo hidrodinámicamente y térmicamente desarrollado. Determine la longitud del ducto requerido para calentar agua desde 30°C hasta 70°C.

Solución: Re = 1 334 Nu = 3 h = 96 w/ m² ºC L = 7 mts.


72


73

5.- Considere el calentamiento de aire atmosférico fluyendo a una velocidad de 0.5 m/seg dentro de un tubo de pared delgada de 2.5 cm de diámetro en una región térmicamente desarrollada. El calentamiento puede ser hecho ya sea por condensamiento de vapor sobre la superficie inferior del tubo, de manera que se mantiene la temperatura constante o por calentamiento por resistencia eléctrica de manera que se mantiene una superficie con flujo de calor uniforme. Calcular el coeficiente de transferencia de calor para ambas condiciones de calentamiento, considerando propiedades del aire a 350°K en promedio.

Solución: Re = 602 Nu = 3.657 a) h = 4.4 w/ m² ºC Nu = 4.35 b) h = 5.24 |w/ m² ºC


74


75

UNIDAD VI CONVECCION NATURAL Convección: fluido con una pared sólida, existe transferencia de calor entre ellos ( )Fc TT − .

- Genera una h de valor pequeño. - El aire es el elemento que presenta mayor convección natural. - Debido a la diferencia de densidades el aire y una pared presentan

convección natural. En algunas situaciones de convección el movimiento del fluido es establecido sin forzar la velocidad. Consideremos una placa vertical caliente colocada en un fluido en reposo que esta a temperatura uniforme menor que la placa. La transferencia de calor se realizara primero por conducción pura y un gradiente de temperatura se establecerá en el fluido. Esta variación de temperatura dentro del fluido generara un gradiente de densidad provocando un movimiento convectivo como resultado de las fuerzas de flotación y es denominando convección natural o convección libre. El movimiento en convección natural generalmente es menor que en convección forzada y por lo tanto la transferencia de calor por convección natural es menor que en la forzada.


76

En el esquema anterior se muestra el desarrollo del campo de velocidades enfrente de la placa vertical caliente debido a fuerzas de flotación. El fluido calentado enfrente de la placa, sube entrando fluido de la región exterior en reposo. Una capa limite de velocidades es desarrollada con una velocidad pico que cae en algún lugar de la capa limite. La velocidad es cero, tanto en la superficie de la placa como en el borde de la capa limite. PARAMETROS ADIMENSIONALES DE CONVECCION LIBRE Número de Grashof (Gr)

( )2

3

vTTLg

G wr

∞−=

β

L, es una longitud característica g, es la gravedad β , es el coeficiente volumétrico de expansión térmica. La ecuación anterior representa el radio de las fuerzas de flotación a las fuerzas viscosas actuando sobre el fluido.

ρρρ

β ⎥⎦⎤

⎢⎣⎡−=dTd1

asvisFuerzasFuerzasRe cos,

=

En convección libre la transición de flujo laminar a turbulento es gobernada por el valor critico del número de Grashof. Algunas veces el parámetro adimensional denominado número de Rayleig esta definido como:

[ ]α

βv

TTLgPGR w

rra∞−

==3

α , difusividad térmica Esta ecuación es usada en lugar de Grashof para correlacionar la transferencia de calor en convección natural, esto es:

( )rau PRfN ,=


77

Mientras más viscoso sea un fluido mas pequeño es Reynolds.

EVALUACION DE PROPIEDADES

[ ]⇒−= ∞TTT wf 21 Para encontrar propiedades del fluido

T1 (Para gases donde la temperatura es dad en ºK)

β Ver apéndice para líquidos. Transición de Laminar a Turbulento

Con 33 11 EPE r ≤≤− y ( )

2

3

vTTyg

G wry

∞−=

β

Para el aire:

72.0=rP , la ecuación anterior se reduce a:

( ) 41

515.068.0 ayuy RN +=

Nota: Si 91EPR

Gr

ayry ≈= y para [ ]33 11 EPE r ≤≤− el régimen es laminar.

Si ,1 9EGry > el régimen es turbulento


78

PARED ISOTERMICA VERTICAL (Placa vertical)

a) Laminar

( ) 41

41

21

492.0986.0671.0 ay

rr

ru R

PP

PN⎥⎥⎦

⎤

⎢⎢⎣

⎡

++=

Donde:

( )α

βv

TTygR w

ay∞−

=3

Donde: Y = altura de la placa


79

b) Turbulento La siguiente ecuación cubre todos los rangos (lamina y turbulento)

( )

2

278

169

61

492.01

387.0825.0

⎥⎥⎥⎥⎥⎥⎥

⎦

⎤

⎢⎢⎢⎢⎢⎢⎢

⎣

⎡

⎥⎥

⎦

⎤

⎢⎢

⎣

⎡⎟⎟⎠

⎞⎜⎜⎝

⎛+

+=

r

ayuy

P

RN

Para el aire: 72.0=rP se reduce a :

( )[ ]261

325.0825.0 ayuy RN +=

Nota: KhLNu = ;

LKNh u=

Si la placa es vertical L = y Si la placa es horizontal

PL Α=

PARED VERTICAL CON FLUJO DE CALOR UNIFORME Cuando wq" es constante, la temperatura en la pared crece monotonicamente en la

dirección “y”. La diferencia de temperaturas ( )∞−TTw se incrementa en función 51

y . Una correlación valida para todos los números de Rayleig y Prandtl es la siguiente:

( )

2

278

169

61

437.01

387.0825.0

⎥⎥⎥⎥⎥⎥⎥

⎦

⎤

⎢⎢⎢⎢⎢⎢⎢

⎣

⎡

⎥⎥

⎦

⎤

⎢⎢

⎣

⎡⎟⎟⎠

⎞⎜⎜⎝

⎛+

+=

r

ayuy

P

RN

En esta expresión, i ayR es basada en la diferencia de temperaturas promedio en “y”, denominada ∞−TTw .


80

PLACAS HORIZONTALES ISOTERMICAS

( )( )1073/1

7441

1115.0

1154.0

ERERN

ERERN

aLaLuL

aLaLuL

<<=

<<= A

( )10541

1127.0 ERERN aLaLuL <<= B

( )α

βv

LTTgR waL

3∞−

=

Siendo P

L Α=

Donde: A, área de la superficie plana P, perímetro del ár


81

CILINDRO HORIZONTAL

( )

2

278

169

61

559.01

387.06.0

⎥⎥⎥⎥⎥⎥⎥

⎦

⎤

⎢⎢⎢⎢⎢⎢⎢

⎣

⎡

⎥⎥

⎦

⎤

⎢⎢

⎣

⎡⎟⎟⎠

⎞⎜⎜⎝

⎛+

+=

r

auD

P

RN 125 11 ERE aD <<−

ESFERA

( )9

416

9

41

469.01

589.02

⎥⎥

⎦

⎤

⎢⎢

⎣

⎡⎟⎟⎠

⎞⎜⎜⎝

⎛+

+=

r

aDuD

P

RN 7.0≥rP y "1ERaD ≤

En ambos casos:

( )v

TTDgR w

aD αβ ∞−

=3

Área de la Esfera: A= πD²


82

PROBLEMAS UNIDAD VI

1.- Un calentador de inversión para agua consiste en una placa vertical delgada de forma rectangular de 8cm de altura y 15 cm de ancho, la placa es calentada eléctricamente y mantenida a 55°C, mientras que la temperatura promedio del agua que la rodea es de 15°C. Calcular: la razón de transferencia de calor eliminado por el calentador dentro del agua.

Formula q” = hA∆T

Respuestas Con interpolación: Ray = 5.93 E h = 735 w/m² ºC q”= 705 W

8


83


84

2.- Flujo de aire a través de un ducto de calentamiento rectangular largo de 0.75 de ancho y 0.3 m de alto, mantiene la superficie del ducto a 45°C. Si el ducto no estuviera aislado y esta expuesto al aire que esta a 15°C ¿Cuál es la perdida de calor del ducto por metro de longitud? Formula q” = hA∆T

Respuestas Con interpolación: Placas Verticales

Ray = 7.25 E 8 Placas Horizontales

Ray = 1.42 E 8 q”/L= 247 w/m


85

UNIDAD VII

INTERCAMBIADORES DE CALOR: Se engloba a todos aquellos dispositivos utilizados para transferir energía de un medio a otro. Los radiadores de calefacción, los calentadores de agua, las baterías frigoríficas... son algunos ejemplos de cambiadores de calor.

CLASIFICACIÓN: La clasificación más general que puede realizarse de los cambiadores de calor, se efectúa atendiendo al grado de contacto entre los fluidos. Según este criterio, los cambiadores de calor se dividen en dos grandes grupos:

Intercambiadores de contacto directo.

Intercambiadores de contacto indirecto.

Estos últimos pueden a su vez dividirse en alternativos y de superficie.

Los intercambiadores de contacto directo, también conocidos como cambiadores de mezcla, son aquellos dispositivos en los que los fluidos sufren una mezcla física completa, realizándose, como consecuencia, la transferencia energética entre ellos.

Pertenecen a este grupo, entre otros tipos de cambiadores, las denominadas torres de refrigeración o torres húmedas, así como los enfriadores de gases.

En cuanto a los intercambiadores alternativos, ambos fluidos recorren un mismo espacio de forma alternada, sin coincidencia entre ellos, de forma tal que la mezcla física de ambos fluidos puede considerarse despreciable. El elemento fundamental de este subgrupo de cambiadores es la superficie que alternativamente recibe y cede la energía térmica.

Por otra parte, se denominan intercambiadores de superficie a aquellos equipos o dispositivos en los que la transferencia térmica se realiza a través de una superficie, plana o cilíndrica, que separa físicamente las corrientes de ambos fluidos, no existiendo por tanto ninguna posibilidad de contacto directo o contaminación entre dichos fluidos, salvo en el caso de rotura de la antedicha superficie de separación.

CLASIFICACIÓN DE LOS CAMBIADORES DE CALOR DE SUPERFICIE: La clasificación más usual de este grupo de cambiadores, se realiza en base a la dirección relativa de los flujos de ambos fluidos, pudiéndose hablar entonces de cambiadores de flujos paralelos y de cambiadores de flujos cruzados, según sus direcciones sean paralelas en el espacio o formen cualquier ángulo en él.

Los cambiadores de flujos paralelos, son generalmente utilizados en el intercambio térmico líquido-líquido, mientras que los de flujos cruzados se utilizan generalmente en el intercambio líquido-gas.

Como se decía anteriormente se denomina cambiadores de calor de flujos paralelos a aquellos cambiadores en los que circulan ambos fluidos con direcciones paralelas en el espacio, si además de tener ambos flujos la misma dirección, tienen el mismo sentido, reciben el nombre de en “equicorriente”, denominándose en “contracorriente” a aquellos en los que los flujos tienen sentidos contrarios.


86

Dentro del subgrupo de cambiadores de flujos paralelos, se emplean entre otros, los denominados “de placas”, “de tubo”, también llamados “de doble tubo”, de “inmersión”, “multitubulares” o de “carcasa y tubos”.

En los cambiadores multitubulares o de carcasa y tubos, es normal combinar la clasificación anterior con otra, basada en el número de veces que cada partícula de los fluidos recorre el cambiador, recibiendo el nombre de paso cada recorrido, así, un cambiador en equicorriente o contracorriente, sería un cambiador de un paso por carcasa y un paso por tubos.

INTERCAMBIADORES DE HAZ TUBULAR Y CARCASA: Este tipo de aparato es el más extendido. En principio, el haz alojado en el interior de una carcasa se ha mantenido en la construcción de condensadores y rehervidores igualmente. El aparato está constituido por un haz de tubos montados sobre dos placas tubulares que llevan un número determinado de placas deflectoras. Por cada extremo se fijan las cajas de distribución que aseguran la circulación del fluido por el interior del haz, en varias fases. El haz está alojado en una carcasa provisto de una tobera de entrada y otra de salida para el segundo fluido que circula por el exterior de los tubos, siguiendo el camino impuesto por las placas deflectoras. Todos los elementos que entran en la construcción de los intercambiadores, han sido objeto de una normalización publicada por T.E.M.A (Estándar of Tubular Exchanger Manufacturers Association), que especifica las características mecánicas y térmicas correspondientes a las diversas condiciones de funcionamiento.

Carcasa: El material más usado para la construcción de las carcasas es el acero al carbono. Para diámetros inferiores a 24”, en la carcasa se emplea un tubo de acero L.P.S (Schedule 30 hasta 12” y 1 cm. de espesor entre 12” y 24”), si la presión de servicio es inferior a 20 Kg./cm2.

Para más de 24” la carcasa se realiza con planchas de acero enrolladas y soldadas. Por cada extremo se sueldan las bridas que llevarán las tapas y las cajas de distribución. Las toberas de entrada y salida se sueldan, o no, con una placa de refuerzo según la presión de servicio. Por último la carcasa se podrá equipar con anillos para poder levantarla y llevará, además, la placa de identidad del aparato.

Las perforaciones de los huecos en estas placas están normalizadas, efectuándose según una disposición, ya sea de paso cuadrado o paso triangular. Teniendo en cuenta la orientación del haz en relación a la dirección general del fluido que circula en la carcasa, se obtienen las cuatro disposiciones de las figuras siguientes.


87

El paso triangular permite colocar alrededor de un 10% de tubos más que en el paso cuadrado sobre una placa tubular de diámetro dado pero, en contrapartida, la disposición de los tubos

hace imposible la limpieza exterior, introduciendo rascadores a través del haz. Para estos aparatos, es necesario recurrir a la limpieza química y reservar su empleo a productos limpios.

El haz de tubos lleva deflectoras transversales que tienen por finalidad alargar el camino del fluido que circular por la carcasa y mejorar así, la transmisión por el exterior de los tubos. Estas deflectoras están constituidas, generalmente, por un disco que tiene un diámetro ligeramente inferior al de la carcasa y que posee un segmento libre igual al 25% del diámetro interior de Dc de la carcasa. El espaciado B, entre deflectoras, que condiciona directamente la velocidad del fluido, está comprendida entre Dc/5 y Dc.

Además estas reflectoras aseguran la rigidez del haz y son solidarias de la placa tubular fija por medio de unos tirantes.

En ciertos casos se emplean deflectoras longitudinales constituidas por una simple chapa inserta en el medio del haz. Esta disposición obliga a efectuar al fluido un ir y venir en la carcasa y se tiene en este caso, un aparato que se denomina “dos pasos del lado carcasa”.


88


89

ANÁLISIS DE INTERCAMBIADORES DE CALOR Estudiaremos dos métodos diferentes de análisis para los intercambiadores de calor.

a) Método de diferencia de temperatura media logarítmica (LMTD) b) Método de número de unidades de transferencia (NTU).

Método LMTD Considerando un intercambiador de calor de doble tubería, los fluidos pueden circular tanto en corrientes paralelas como a contracorrientes y se calcularía la transferencia de calor en este dispositivo mediante: mTUq Α∆=" Si consideramos un intercambiador de calor distinto del de doble tubería, la transferencia de calor se calculara utilizando un factor de corrección que se aplica a la LMTD para un dispositivo de doble tubería a contracorriente con las mismas temperaturas fría y caliente para el fluido, entonces la ecuación anterior se expresa de la siguiente manera:

FTUq mΑ∆=" Donde:

ei hKL

h

U11

1

++=

Determinando los coeficientes de convección ( )ei hh , obtendremos una U. Donde: U, coeficiente global de transferencia de calor A, superficie de transferencia de calor consistente con la definición de U

mT∆ , diferencia media de temperaturas apropiada a través del cambiador de calor F, factor de corrección para diversos tipos de intercambiadores de calor. (Ver tablas 10.8 a 10.11).

Nota: Cuando interviene un cambio de fase, como el caso de la condensación o la ebullición (evaporación), el fluido permanece normalmente a una temperatura prácticamente constante y las relaciones se simplifican, entonces: F = 1, para ebullición o condensación El factor de corrección depende también del tipo de intercambiador de calor.

a) Flujo paralelo b) Contra flujo c) Flujo cruzado


90

La diferencia de temperatura media logarítmica se expresaría de la siguiente manera:

2

1

21

lnTT

TTTm

∆∆∆−∆

=∆

Balance de Energía

( )TmCq p∆= Fluido frío = ( )TmCp∆ Fluido caliente


91

ANALISIS DE INTERCAMBIADORES DE CALOR POR EL METODO EFICIENCIA (Є)-NTU El método LMTD para el análisis de intercambiadores de calor, es útil cuando las temperaturas de entrada y salida son conocidas o se pueden obtener con facilidad. Cuando tenemos que evaluar las temperaturas de entrada o salida de un intercambiador determinado, el análisis se debe realizar mediante un procedimiento iterativo, debido a la función logarítmica que aparece en la LMTD. En estos casos podemos realizar un análisis más fácil basado en la eficiencia del intercambiador de calor durante la transferencia de una entrada determinada de calor. Este método ofrece también ventajas para el análisis de problemas donde se tiene que comparar varios tipos de intercambiadores de calor, para elegir el más adecuado para un objetivo de transferencia de calor particular. Definiremos la eficiencia como:

maxposiblecalor de ncia transfereMáxima

realcalor de ciaTransferenq

q==ε

Donde podemos calcular a q como: ( ) ( )

CalientepFriap TmCTmCq ∆∆= Mientras maxq para el intercambiador es conocido. Debemos primero reconocer que este máximo valor podría ser alcanzado si uno de los fluidos desarrolla un cambio de temperatura igual a la máxima diferencia de temperaturas presente en el intercambiador de calor; que es la diferencia de las temperaturas de entrada entre los fluidos frío y caliente. El fluido que podría alcanzar esta máxima diferencia de temperatura es el que tiene un mínimo valor de ( pmC ) ya que el balance de energía requiere que la energía recibida por un fluido sea igual al del otro fluido que se esta suministrando. ( ) [ ]entradaFentradaCp TTmCq −=

minmax Las expresiones para el cálculo de la eficiencia en función del flujo son:

a) Flujo Paralelo b) Contra Flujo

( )[ ]C

CNTU+

+−−=

11exp1ε ( )[ ]

( )[ ]CNTUCCNTU+−−+−−

=1exp1

1exp1ε


92

Para Ambas:

minCUANTU =

max

min

CCC =

( )( )

maxmax

minmin

p

p

mCC

mCC

=

=

Donde: NTU,= Number of Thermal Units (número de unidades de transferencia). A= área de transferencia de calor. U= coeficiente global de transferencia de calor.


93

PROBLEMAS UNIDAD VII

1.- Un intercambiador de calor de tubo y coraza debe ser diseñado para calentar 2.5 Kg./seg. de agua desde 15°C hasta 85°C. El calentamiento se realiza pasando aceite de motor caliente que entra a 160°C. El aceite tiene un coeficiente de convección promedio de 400 W/m2°K sobre la parte exterior de los tubos.10 tubos pasan el agua a través de la coraza, cada tubo es de pared delgada con un diámetro de 25 mm. y hace 8 pasos a través de la coraza, si el aceite deja el intercambiador de calor a 100°C ¿Cuál es la razón de flujo de aceite? ¿Qué longitud deben tener los tubos para alcanzar el calentamiento deseado?

Solución: segkgmaceite /5=

22600Re = 116=Nu

CmWhi º/3000 2=9.0=F mL 36=


94


95

2.- Gases calientes de escape entran a 300°C a un intercambiador de calor de flujo cruzado y salen a 100°C. Estos gases son usados para calentar agua que entran a una razón de 1Kg. /seg. desde 35°C hasta 125°C. El calor especifico de los gases de escape es aproximadamente 1KJ/kg°K, el coeficiente global de transferencia de calor basado en la superficie del lado del gas es igual a 100 W/m2°K. Determinar el área requerida usando el método NTU.

Solución: 75.0=ε 45.0=C

46.2=NTU25.46 mA =


96


97

3.- Se utiliza aceite (Cp = 2.1 KJ/Kg. ºK) para calentar agua en un intercambiador de calor de tubos y casco con un solo recorrido por el casco y dos por los tubos. El coeficiente total de transferencia de calor es de 525 W/m2 ºK. Las velocidades de flujo másico son de 7 Kg./s para el aceite y de 10 Kg./s para el agua. El aceite y el agua entran en el intercambiador de calor a 240ºC y 20ºC, respectivamente. Se tiene que diseñar el intercambiador de calor de modo que el agua salga de el a una temperatura mínima de 80ºC. Calcule el área para la transferencia de calor necesaria para lograr esta temperatura.

Solución: To,aceite = 69.46° C A = 24.70 m=


98

4.- En un intercambiador de calor de tubos y casco se introduce agua a 35 ºC para calentarla utilizando aceite hasta 75º C, que entra a 110 ºC y sale a 75 ºC. El intercambiador esta dispuesto a contra flujo, donde el agua tiene un recorrido por el casco y el aceite tiene dos recorridos por los tubos. Si la razón del flujo del agua es de 68 Kg. por minuto y el coeficiente total de transferencia de calor se estima en 320 W/m2 ºK de acuerdo con la tabla 8.1 calcule el área de intercambio necesaria.

Solución: maceite=2.3 Kg/seg q = 189448 W F= 0.85 A=1.8 m2

apuntes transferencia de calor_2010

Documents